Specificity of Loxosceles α clade phospholipase D enzymes for choline-containing lipids: Role of a conserved aromatic cage.

Spider venom GDPD-like phospholipases D (SicTox) have been identified to be one of the major toxins in recluse spider venom. They are divided into two major clades: the α clade and the ß clade. Most α clade toxins present high activity against lipids with choline head groups such as sphingomyelin, while activities in ß clade toxins vary and include preference for substrates containing ethanolamine headgroups (Sicarius terrosus, St_ßIB1). A structural comparison of available structures of phospholipases D (PLDs) reveals a conserved aromatic cage in the α clade. To test the potential influence of the aromatic cage on membrane-lipid specificity we performed molecular dynamics (MD) simulations of the binding of several PLDs onto lipid bilayers containing choline headgroups; two SicTox from the α clade, Loxosceles intermedia αIA1 (Li_αIA) and Loxosceles laeta αIII1 (Ll_αIII1), and one from the ß clade, St_ßIB1. The simulation results reveal that the aromatic cage captures a choline-headgroup and suggest that the cage plays a major role in lipid specificity. We also simulated an engineered St_ßIB1, where we introduced the aromatic cage, and this led to binding with choline-containing lipids. Moreover, a multiple sequence alignment revealed the conservation of the aromatic cage among the α clade PLDs. Here, we confirmed that the i-face of α and ß clade PLDs is involved in their binding to choline and ethanolamine-containing bilayers, respectively. Furthermore, our results suggest a major role in choline lipid recognition of the aromatic cage of the α clade PLDs. The MD simulation results are supported by in vitro liposome binding assay experiments.

Extreme divergence between one-to-one orthologs: the structure of N15 Cro bound to operator DNA and its relationship to the λ Cro complex.

The gene cro promotes lytic growth of phages through binding of Cro protein dimers to regulatory DNA sites. Most Cro proteins are one-to-one orthologs, yet their sequence, structure and binding site sequences are quite divergent across lambdoid phages. We report the cocrystal structure of bacteriophage N15 Cro with a symmetric consensus site. We contrast this complex with an orthologous structure from phage λ, which has a dissimilar binding site sequence and a Cro protein that is highly divergent in sequence, dimerization interface and protein fold. The N15 Cro complex has less DNA bending and smaller DNA-induced changes in protein structure. N15 Cro makes fewer direct contacts and hydrogen bonds to bases, relying mostly on water-mediated and Van der Waals contacts to recognize the sequence. The recognition helices of N15 Cro and λ Cro make mostly nonhomologous and nonanalogous contacts. Interface alignment scores show that half-site binding geometries of N15 Cro and λ Cro are less similar to each other than to distantly related CI repressors. Despite this divergence, the Cro family shows several code-like protein-DNA sequence covariations. In some cases, orthologous genes can achieve a similar biological function using very different specific molecular interactions.

Increased sequence hydrophobicity reduces conformational specificity: A mutational case study of the Arc repressor protein.

The amino-acid sequences of soluble, globular proteins must have hydrophobic residues to form a stable core, but excess sequence hydrophobicity can lead to loss of native state conformational specificity and aggregation. Previous studies of polar-to-hydrophobic mutations in the ß-sheet of the Arc repressor dimer showed that a single substitution at position 11 (N11L) leads to population of an alternate dimeric fold in which the ß-sheet is replaced by helix. Two additional hydrophobic mutations at positions 9 and 13 (Q9V and R13V) lead to population of a differently folded octamer along with both dimeric folds. Here we conduct a comprehensive study of the sequence determinants of this progressive loss of fold specificity. We find that the alternate dimer-fold specifically results from the N11L substitution and is not promoted by other hydrophobic substitutions in the ß-sheet. We also find that three highly hydrophobic substitutions at positions 9, 11, and 13 are necessary and sufficient for oligomer formation, but the oligomer size depends on the identity of the hydrophobic residue in question. The hydrophobic substitutions increase thermal stability, illustrating how increased hydrophobicity can increase folding stability even as it degrades conformational specificity. The oligomeric variants are predicted to be aggregation-prone but may be hindered from doing so by proline residues that flank the ß-sheet region. Loss of conformational specificity due to increased hydrophobicity can manifest itself at any level of structure, depending upon the specific mutations and the context in which they occur.

Evolutionary dynamics of origin and loss in the deep history of phospholipase D toxin genes.

BACKGROUND: Venom-expressed sphingomyelinase D/phospholipase D (SMase D/PLD) enzymes evolved from the ubiquitous glycerophosphoryl diester phosphodiesterases (GDPD). Expression of GDPD-like SMaseD/PLD toxins in both arachnids and bacteria has inspired consideration of the relative contributions of lateral gene transfer and convergent recruitment in the evolutionary history of this lineage. Previous work recognized two distinct lineages, a SicTox-like (ST-like) clade including the arachnid toxins, and an Actinobacterial-toxin like (AT-like) clade including the bacterial toxins and numerous fungal homologs. RESULTS: Here we expand taxon sampling by homology detection to discover new GDPD-like SMase D/PLD homologs. The ST-like clade now includes homologs in a wider variety of arthropods along with a sister group in Cnidaria; the AT-like clade now includes additional fungal phyla and proteobacterial homologs; and we report a third clade expressed in diverse aquatic metazoan taxa, a few single-celled eukaryotes, and a few aquatic proteobacteria. GDPD-like SMaseD/PLDs have an ancient presence in chelicerates within the ST-like family and ctenophores within the Aquatic family. A rooted phylogenetic tree shows that the three clades derived from a basal paraphyletic group of proteobacterial GDPD-like SMase D/PLDs, some of which are on mobile genetic elements. GDPD-like SMase D/PLDs share a signature C-terminal motif and a shortened ßα1 loop, features that distinguish them from GDPDs. The three major clades also have active site loop signatures that distinguish them from GDPDs and from each other. Analysis of molecular phylogenies with respect to organismal relationships reveals a dynamic evolutionary history including both lateral gene transfer and gene duplication/loss. CONCLUSIONS: The GDPD-like SMaseD/PLD enzymes derive from a single ancient ancestor, likely proteobacterial, and radiated into diverse organismal lineages at least in part through lateral gene transfer.

Variable Substrate Preference among Phospholipase D Toxins from Sicariid Spiders.

Venoms of the sicariid spiders contain phospholipase D enzyme toxins that can cause severe dermonecrosis and even death in humans. These enzymes convert sphingolipid and lysolipid substrates to cyclic phosphates by activating a hydroxyl nucleophile present in both classes of lipid. The most medically relevant substrates are thought to be sphingomyelin and/or lysophosphatidylcholine. To better understand the substrate preference of these toxins, we used (31)P NMR to compare the activity of three related but phylogenetically diverse sicariid toxins against a diverse panel of sphingolipid and lysolipid substrates. Two of the three showed significantly faster turnover of sphingolipids over lysolipids, and all three showed a strong preference for positively charged (choline and/or ethanolamine) over neutral (glycerol and serine) headgroups. Strikingly, however, the enzymes vary widely in their preference for choline, the headgroup of both sphingomyelin and lysophosphatidylcholine, versus ethanolamine. An enzyme from Sicarius terrosus showed a strong preference for ethanolamine over choline, whereas two paralogous enzymes from Loxosceles arizonica either preferred choline or showed no significant preference. Intrigued by the novel substrate preference of the Sicarius enzyme, we solved its crystal structure at 2.1 Å resolution. The evolution of variable substrate specificity may help explain the reduced dermonecrotic potential of some natural toxin variants, because mammalian sphingolipids use primarily choline as a positively charged headgroup; it may also be relevant for sicariid predatory behavior, because ethanolamine-containing sphingolipids are common in insect prey.

A role for indels in the evolution of Cro protein folds.

Insertions and deletions in protein sequences, or indels, can disrupt structure and may result in changes in protein folds during evolution or in association with alternative splicing. Pfl 6 and Xfaso 1 are two proteins in the Cro family that share a common ancestor but have different folds. Sequence alignments of the two proteins show two gaps, one at the N terminus, where the sequence of Xfaso 1 is two residues shorter, and one near the center of the sequence, where the sequence of Pfl 6 is five residues shorter. To test the potential importance of indels in Cro protein evolution, we generated hybrid variants of Pfl 6 and Xfaso 1 with indels in one or both regions, chosen according to several plausible sequence alignments. All but one deletion variant completely unfolded both proteins, showing that a longer N-terminal sequence was critical for Pfl 6 folding and a longer central region sequence was critical for Xfaso 1 folding. By contrast, Xfaso 1 tolerated a longer N-terminal sequence with little destabilization, and Pfl 6 tolerated central region insertions, albeit with substantial effects on thermal stability and some perturbation of the surrounding structure. None of the mutations appeared to convert one stable fold into the other. On the basis of this two-protein comparison, short insertion and deletion mutations probably played a role in evolutionary fold change in the Cro family, but were also not the only factors.

Protein salvage and repurposing in evolution: Phospholipase D toxins are stabilized by a remodeled scrap of a membrane association domain.

The glycerophosphodiester phosphodiesterase (GDPD)-like SMaseD/PLD domain family, which includes phospholipase D (PLD) toxins in recluse spiders and actinobacteria, evolved anciently in bacteria from the GDPD. The PLD enzymes retained the core (ß/α)8 barrel fold of GDPD, while gaining a signature C-terminal expansion motif and losing a small insertion domain. Using sequence alignments and phylogenetic analysis, we infer that the C-terminal motif derives from a segment of an ancient bacterial PLAT domain. Formally, part of a protein containing a PLAT domain repeat underwent fusion to the C terminus of a GDPD barrel, leading to attachment of a segment of a PLAT domain, followed by a second complete PLAT domain. The complete domain was retained only in some basal homologs, but the PLAT segment was conserved and repurposed as the expansion motif. The PLAT segment corresponds to strands ß7-ß8 of a ß-sandwich, while the expansion motif as represented in spider PLD toxins has been remodeled as an α-helix, a ß-strand, and an ordered loop. The GDPD-PLAT fusion led to two acquisitions in founding the GDPD-like SMaseD/PLD family: (1) a PLAT domain that presumably supported early lipase activity by mediating membrane association, and (2) an expansion motif that putatively stabilized the catalytic domain, possibly compensating for, or permitting, loss of the insertion domain. Of wider significance, messy domain shuffling events can leave behind scraps of domains that can be salvaged, remodeled, and repurposed.

Measuring 1HN temperature coefficients in invisible protein states by relaxation dispersion NMR spectroscopy.

A method based on the Carr-Purcell-Meiboom-Gill relaxation dispersion experiment is presented for measuring the temperature coefficients of amide proton chemical shifts of low populated 'invisible' protein states that exchange with a 'visible' ground state on the millisecond time-scale. The utility of the approach is demonstrated with an application to an I58D mutant of the Pfl6 Cro protein that undergoes exchange between the native, folded state and a cold denatured, unfolded conformational ensemble that is populated at a level of 6% at 2.5°C. A wide distribution of amide temperature coefficients is measured for the unfolded state. The distribution is centered about -5.6 ppb/K, consistent with an absence of intra-molecular hydrogen bonds, on average. However, the large range of values (standard deviation of 2.1 ppb/K) strongly supports the notion that the unfolded state of the protein is not a true random coil polypeptide chain.

Determinants of gas-phase disassembly behavior in homodimeric protein complexes with related yet divergent structures.

The overall structure of a protein-protein complex reflects an intricate arrangement of noncovalent interactions. Whereas intramolecular interactions confer secondary and tertiary structure to individual subunits, intermolecular interactions lead to quaternary structure--the ordered aggregation of separate polypeptide chains into multisubunit assemblies. The specific ensemble of noncovalent contacts dictates the stability of subunit folds, enforces protein-protein binding specificity, and determines multimer stability. Consequently, noncovalent architecture is likely to play a role in the gas-phase dissociation of these assemblies during tandem mass spectrometry (MS/MS). To further advance the applicability of MS/MS to analytical problems in structural biology, a better understanding of the interplay between the structures and fragmentation behaviors of noncovalent protein complexes is essential. The present work constitutes a systematic study of model protein homodimers (bacteriophage N15 Cro, bacteriophage λ Cro, and bacteriophage P22 Arc) with related but divergent structures, both in terms of subunit folds and protein-protein interfaces. Because each of these dimers has a well-characterized structure (solution and/or crystal structure), specific noncovalent features could be correlated with gas-phase disassembly patterns as studied by collision-induced dissociation, surface-induced dissociation, and ion mobility. Of the several respects in which the dimers differed in structure, the presence or absence of intermolecular electrostatic contacts exerted the most significant influence on the gas-phase dissociation behavior. This is attributed to the well-known enhancement of ionic interactions in the absence of bulk solvent. Because salt bridges are general contributors to both intermolecular and intramolecular stability in protein complexes, these observations are broadly applicable to aid in the interpretation or prediction of dissociation spectra for noncovalent protein assemblies.

Transitive homology-guided structural studies lead to discovery of Cro proteins with 40% sequence identity but different folds.

Proteins that share common ancestry may differ in structure and function because of divergent evolution of their amino acid sequences. For a typical diverse protein superfamily, the properties of a few scattered members are known from experiment. A satisfying picture of functional and structural evolution in relation to sequence changes, however, may require characterization of a larger, well chosen subset. Here, we employ a "stepping-stone" method, based on transitive homology, to target sequences intermediate between two related proteins with known divergent properties. We apply the approach to the question of how new protein folds can evolve from preexisting folds and, in particular, to an evolutionary change in secondary structure and oligomeric state in the Cro family of bacteriophage transcription factors, initially identified by sequence-structure comparison of distant homologs from phages P22 and lambda. We report crystal structures of two Cro proteins, Xfaso 1 and Pfl 6, with sequences intermediate between those of P22 and lambda. The domains show 40% sequence identity but differ by switching of alpha-helix to beta-sheet in a C-terminal region spanning approximately 25 residues. Sedimentation analysis also suggests a correlation between helix-to-sheet conversion and strengthened dimerization.