Parasitoid fig-wasp evolutionary diversification and variation in ecological opportunity.

Ecological processes are manifest in the evolution and form of phenotype diversity. The great abundance of parasitoid species has led to speculation whether rates of speciation and extinction are dependent on parasitoid diversity. If these factors are mutually exclusive, species diversity should fluctuate instead of remaining relatively constant over time. It is not known whether radiations constrained by coevolutionary interactions conform to density-dependent diversification processes. Here we test the prediction that parasitoid fig wasp diversification responds to changes in ecological opportunity and density-independent processes. A phylogenetic approach is used to estimate relative divergence times and infer diversification rate changes using gamma-statistics. Monte Carlo constant rates tests that accommodate incomplete sampling could not reject constant rates diversification. Parasitoid fig wasp diversification is consistent with a more complex explanation than density-dependent cladogenesis. The results suggest contemporary African parasitoid fig wasp diversity remains a legacy of an ancient ecological opportunity facilitated by fig tree diversification following the breakup of Pan-African forests and evolution of the savanna biome over the last 55 Ma and the more recent aridification of the African continent in the last 5 Ma. These results imply that amplified phenotypic differentiation of specialist insects coevolving with plants is coupled to evolutionarily infrequent changes in ecological opportunity.

Parallel diversification of Australian gall-thrips on Acacia.

The diversification of gall-inducing Australian Kladothrips (Insecta: Thysanoptera) on Acacia has produced a pair of sister-clades, each of which includes a suite of lineages that utilize virtually the same set of 15 closely related host plant species. This pattern of parallel insect-host plant radiation may be driven by cospeciation, host-shifting to the same set of host plants, or some combination of these processes. We used molecular-phylogenetic data on the two gall-thrips clades to analyze the degree of concordance between their phylogenies, which is indicative of parallel divergence. Analyses of phylogenetic concordance indicate statistically-significant similarity between the two clades. Their topologies also fit with a hypothesis of some degree of host-plant tracking. Based on phylogenetic and taxonomic information regarding the phylogeny of the Acacia host plants in each clade, one or more species has apparently shifted to more-divergent Acacia host-plant species, and in each case these shifts have resulted in notable divergence in aspects of the phenotype including morphology, life history and behaviour. Our analyses indicate that gall-thrips on Australian Acacia have undergone parallel diversification as a result of some combination of cospeciation, highly restricted host-plant shifting, or both processes, but that the evolution of novel phenotypic diversity in this group is a function of relatively few shifts to divergent host plants. This combination of ecologically restricted and divergent radiation may represent a microcosm for the macroevolution of host plant relationships and phenotypic diversity among other phytophagous insects.

Specific inhibitors of Plasmodium falciparum thioredoxin reductase as potential antimalarial agents.

Plasmodium falciparum thioredoxin reductase (PfTrxR: NADPH+Trx(S)2+H+<-->NADP++Trx(SH)2) is a high Mr flavin-dependent TrxR that reduces thioredoxin (Trx) via a CysXXXXCys pair located penultimately to the C-terminal Gly. In this respect, PfTrxR differs significantly from its human counterpart which bears a Cys-Sec redox pair at the same position. PfTrxR is essentially involved in antioxidant defense and redox regulation of the parasite and has been previously validated by knock-out studies as a potential drug target for malaria chemotherapy. Moreover, human TrxR is present in most cancer cells at levels tenfold higher than in normal cells. Here we report the discovery of a series of potent inhibitors of PfTrxR. The three most promising inhibitors, 3(IC50(PfTrxR)=2 microM and IC50(hTrxR)=50 microM), 7(IC50(PfTrxR)=2 microM and IC50(hTrxR)=140 microM), and 11(IC50(PfTrxR)=0.5 microM and IC50(hTrxR)=4 microM) were selective for the parasite enzyme. Detailed mechanistic characterization of the effects of these compounds on the PfTrxR-catalyzed reaction showed clear uncompetitive inhibition with respect to both substrate and cofactor. For the most specific PfTrxR inhibitor 7, an alkylation mechanism study based on a thiol conjugation model was performed. Furthermore, all three compounds were active in the lower micromolar range on the chloroquine-resistant P. falciparum strain K1 in vitro.

Inbreeding ancestors: the role of sibmating in the social evolution of gall thrips.

We used microsatellite data to estimate levels of inbreeding in four species of solitary gall thrips that are in the same clade as the six species with soldier castes. Three of the four species were highly inbred (Fis 0.54-0.68), and the other apparently mated randomly (Fis near zero). These estimates, combined with previous data from species with soldiers, suggest that inbreeding is a pervasive life-history feature of the gall-inducing thrips on Australian Acacia. Mapping of inbreeding estimates onto the phylogeny of the gall inducers showed that the ancestral lineage that gave rise to soldiers was apparently highly inbred, and therefore, inbreeding could have played a role in the origin of sociality within this group. Moreover, there was a trend from high levels of inbreeding at the origin of soldiers to low levels in the most derived species with soldiers, which exhibits the highest levels of reproductive division of labor and soldier altruism. These patterns are consistent with considerations from population genetics, which show that the likelihood of the origin of soldier altruism is higher in inbreeding populations but that, once soldiers have evolved, a reduction in inbreeding levels may facilitate the evolution of enhanced division of labor and reproductive skew.

Getting the adrenaline going: crystal structure of the adrenaline-synthesizing enzyme PNMT.

BACKGROUND: Adrenaline is localized to specific regions of the central nervous system (CNS), but its role therein is unclear because of a lack of suitable pharmacologic agents. Ideally, a chemical is required that crosses the blood-brain barrier, potently inhibits the adrenaline-synthesizing enzyme PNMT, and does not affect other catecholamine processes. Currently available PNMT inhibitors do not meet these criteria. We aim to produce potent, selective, and CNS-active PNMT inhibitors by structure-based design methods. The first step is the structure determination of PNMT. RESULTS: We have solved the crystal structure of human PNMT complexed with a cofactor product and a submicromolar inhibitor at a resolution of 2.4 A. The structure reveals a highly decorated methyltransferase fold, with an active site protected from solvent by an extensive cover formed from several discrete structural motifs. The structure of PNMT shows that the inhibitor interacts with the enzyme in a different mode from the (modeled) substrate noradrenaline. Specifically, the position and orientation of the amines is not equivalent. CONCLUSIONS: An unexpected finding is that the structure of PNMT provides independent evidence of both backward evolution and fold recruitment in the evolution of a complex enzyme from a simple fold. The proposed evolutionary pathway implies that adrenaline, the product of PNMT catalysis, is a relative newcomer in the catecholamine family. The PNMT structure reported here enables the design of potent and selective inhibitors with which to characterize the role of adrenaline in the CNS. Such chemical probes could potentially be useful as novel therapeutics.

An unusually low pK(a) for Cys282 in the active site of human muscle creatine kinase.

All phosphagen kinases contain a conserved cysteine residue which has been shown by crystallographic studies, on both creatine kinase and arginine kinase, to be located in the active site. There are conflicting reports as to whether this cysteine is essential for catalysis. In this study we have used site-directed mutagenesis to replace Cys282 of human muscle creatine kinase with serine and methionine. In addition, we have replaced Cys282, conserved across all creatine kinases, with alanine. No activity was found with the C282M mutant. The C282S mutant showed significant, albeit greatly reduced, activity in both the forward (creatine phosphorylation) and reverse (MgADP phosphorylation) reactions. The K(m) for creatine was increased approximately 10-fold, but the K(m) for phosphocreatine was relatively unaffected. The V and V/K pH-profiles for the wild-type enzyme were similar to those reported for rabbit muscle creatine kinase, the most widely studied creatine kinase isozyme. However, the V/K(creatine) profile for the C282S mutant was missing a pK of 5.4. This suggests that Cys282 exists as the thiolate anion, and is necessary for the optimal binding of creatine. The low pK of Cys282 was also determined spectrophotometrically and found to be 5.6 +/- 0.1. The S284A mutant was found to have reduced catalytic activity, as well as a 15-fold increase in K(m) for creatine. The pK(a) of Cys282 in this mutant was found to be 6.7 +/- 0.1, indicating that H-bonding to Ser284 is an important, but not the sole, factor contributing to the unusually low pK(a) of Cys282.

Phenylethanolamine N-methyltransferase kinetics: bovine versus recombinant human enzyme.

Epinephrine (Epi) acts as a neurotransmitter in the brain, but its function therein is not well understood. Phenylethanolamine N-methyltransferase (PNMT) catalyzes the final step in the biosynthesis of Epi and is thus a pharmacological target to investigate the function of Epi in the central nervous system. The kinetic differences between bovine adrenal PNMT and human brain PNMT for a number of substrates and inhibitors are examined and the results reported.

Mutagenesis of two acidic active site residues in human muscle creatine kinase: implications for the catalytic mechanism.

Creatine kinase (CK) catalyzes the reversible phosphorylation of the guanidine substrate, creatine, by MgATP. Although several X-ray crystal structures of various isoforms of creatine kinase have been published, the detailed catalytic mechanism remains unresolved. A crystal structure of the CK homologue, arginine kinase (AK), complexed with the transition-state analogue (arginine-nitrate-ADP), has revealed two carboxylate amino acid residues (Glu225 and Glu314) within 2.8 A of the proposed transphosphorylation site. These two residues are the putative catalytic groups that may promote nucleophilic attack by the guanidine amino group on the gamma-phosphate of ATP. From primary sequence alignments of arginine kinases and creatine kinases, we have identified two homologous creatine kinase acidic amino acid residues (Glu232 and Asp326), and these were targeted for examination of their potential roles in the CK mechanism. Using site-directed mutagenesis, we have made several substitutions at these two positions. The results indicate that of these two residues the Glu232 is the likely catalytic residue while Asp326 likely performs a role in properly aligning substrates for catalysis.

Spectroscopic detection of transient thiamin diphosphate-bound intermediates on benzoylformate decarboxylase.

Thiamin diphosphate (ThDP)-dependent enzymes catalyze a range of transformations, such as decarboxylation and ligation. We report a novel spectroscopic assay for detection of some of the ThDP-bound intermediates produced on benzoylformate decarboxylase. Benzoylformate decarboxylase was mixed with its alternate substrate p-nitrobenzoylformic acid on a rapid-scan stopped-flow instrument, resulting in formation of three absorbing species (lambda(max) in parentheses): I(1) (a transient at 620 nm), I(2) (a transient at 400 nm), and I(3) (a stable absorbance with lambda(max) > 730 nm). Analysis of the kinetics of the two transient species supports a model in which a noncovalent complex of the substrate and the enzyme is converted to the first covalent intermediate I(1); the absorbance corresponding to I(1) is probably a charge-transfer band arising from the interaction of the thiamin diphosphate-p-nitrobenzoylformic acid covalent adduct (2-p-nitromandelylThDP) and the enzyme. The rate of disappearance of I(1) parallels the rate of formation of I(2). Chemical models suggest the lambda(max) of I(2) (near 400 nm) to be appropriate to the enamine, a key intermediate in ThDP-dependent reactions resulting from the decarboxylation of the thiamin diphosphate-p-nitrobenzoylformic acid covalent adduct. Therefore, the rate of disappearance of I(1) and/or the appearance of I(2) directly measure the rate of decarboxylation. A relaxation kinetic treatment of the pre-steady-state kinetic data also revealed a hitherto unreported facet of the mechanism, alternating active-sites reactivity. Parallel studies of the His70Ala BFD active-site variant indicate that it cannot form the complex reported by the charge-transfer band (I(1)) at the level of the wild-type protein.

Studies on the conformational properties of CP-10(42-55), the hinge region of CP-10, using circular dichroism and RP-HPLC.

The conformational properties of CP-10(42-55), a peptide corresponding to the hinge region of CP-10, were investigated using circular dichroism spectroscopy and reverse-phase high-performance liquid chromatography (RP-HPLC). The circular dichroism studies indicated that CP-10(42-55) formed considerable secondary structure in the presence of hydrophobic solution environments including 50% acetonitrile, 50% trifluoroethanol and 200 mM sodium dodecyl sulfate, which comprised a mixture of alpha-helix and beta-sheet. The effect of temperature on the conformation of CP-10(42-55) was investigated between 5 and 40 degrees C, with very small changes in the spectra being observed. RP-HPLC was then used to investigate the effect of temperature on the conformation of CP-10(42-55) in the presence of a hydrophobic surface. Using a C18-adsorbent, CP-10(42-55) exhibited a conformational transition at 25 degrees C, which was associated with an increase in the chromatographic contact area and the binding affinity of the peptide for the stationary phase. In addition, near-planar bandbroadening behaviour indicated that conformational species interconverted with rapid rate constants compared with the chromatographic time scale. These results indicated that the conformational change at 25 degrees C in the RP-HPLC system most likely corresponds to the unfolding of an alpha-helical and/or beta-sheet structure to an extended coil structure. Therefore, the strong chemotactic properties of this peptide may be attributed to its ability to form considerable secondary structure in the presence of a hydrophobic environment.

A comparative study of human muscle and brain creatine kinases expressed in Escherichia coli.

We report the expression of the human muscle (CK-MM) and brain (CK-BB) creatine kinases in Escherichia coli. The proteins have been purified to apparent homogeneity and several of their physical and kinetic properties investigated. In the process, we have conclusively verified the correct DNA sequence of the genes encoding the respective isozymes, and determined the correct primary structure and mass of the gene products. Alignment of the primary sequences of these two enzymes shows 81% sequence identity with each other, and no obvious gross structural differences. However, Western blot analyses demonstrated the general lack of antigenic cross-reactivity between these isozymes. Preliminary kinetic analyses show the K(m) and k(cat) values for the creatine and MgATP substrates are similar to values reported for other isozymes from various tissues and organisms. The human muscle and brain CKs do not, however, exhibit the synergism of substrate binding that is observed, for example, in rabbit muscle creatine kinase.

Involvement of electrostatic interactions in the mechanism of peptide folding induced by sodium dodecyl sulfate binding.

Sodium dodecyl sulfate (SDS) has consistently been shown to induce secondary structure, particularly alpha-helices, in polypeptides, and is commonly used to model membrane and other hydrophobic environments. However, the precise mechanism by which SDS induces these conformational changes remains unclear. To examine the role of electrostatic interactions in this mechanism, we have designed two hydrophilic, charged amphipathic alpha-helical peptides, one basic (QAPAYKKAAKKLAES) and the other acidic (QAPAYEEAAEELAKS), and their structures were studied by CD and NMR. The design of the peptides is based on the sequence of the segment of residues 56-70 of human platelet factor 4 [PF4(56-70), QAPLYKKIIKKLLES]. Both peptides were unstructured in water, and in the presence of neutral, zwitterionic, or cationic detergents. However, in SDS at neutral pH, the basic peptide folded into an alpha-helix. By contrast, the pH needed to be lowered to 1.8 before alpha-helix formation was observed for the acidic peptide. Strong, attractive electrostatic interactions, between the anionic groups of SDS and the cationic groups of the lysines, appeared to be necessary to initiate the folding of the basic peptide. NMR analysis showed that the basic peptide was fully embedded in SDS-peptide micelles, and that its three-dimensional alpha-helical structure could be superimposed on that of the native structure of PF4(56-70). These results enabled us to propose a working model of the basic peptide-SDS complex, and a mechanism for SDS-induced alpha-helical folding. This study demonstrates that, while the folding of peptides is mostly driven by hydrophobic effects, electrostatic interactions play a significant role in the formation and the stabilization of SDS-induced structure.

Identification of initiation sites for T4 lysozyme folding using CD and NMR spectroscopy of peptide fragments.

Using CD and 2D (1)H NMR spectroscopy, we have identified potential initiation sites for the folding of T4 lysozyme by examining the conformational preferences of peptide fragments corresponding to regions of secondary structure. CD spectropolarimetry showed most peptides were unstructured in water, but adopted partial helical conformations in TFE and SDS solution. This was also consistent with the (1)H NMR data which showed that the peptides were predominantly disordered in water, although in some cases, nascent or small populations of partially folded conformations could be detected. NOE patterns, coupling constants, and deviations from random coil Halpha chemical shift values complemented the CD data and confirmed that many of the peptides were helical in TFE and SDS micelles. In particular, the peptide corresponding to helix E in the native enzyme formed a well-defined helix in both TFE and SDS, indicating that helix E potentially forms an initiation site for T4 lysozyme folding. The data for the other peptides indicated that helices D, F, G, and H are dependent on tertiary interactions for their folding and/or stability. Overall, the results from this study, and those of our earlier studies, are in agreement with modeling and HD-deuterium exchange experiments, and support an hierarchical model of folding for T4 lysozyme.

Structural analysis of the heparin-binding site of the NC1 domain of collagen XIV by CD and NMR.

Type XIV collagen, a fibril-associated collagen with interrupted triple helices (FACIT), interacts with the surrounding extracellular matrix and/or with cells via its binding to glycosaminoglycans (GAGs). To further characterize such interactions in the NC1 domain of chicken collagen XIV, we identified amino acids essential for heparin binding by affinity chromatography analysis after proteolytic digestion of the synthetic peptide NC1(84-116). The 3D structure of this peptide was then obtained using circular dichroism and NMR. The NC1(84-116) peptide appeared poorly structured in water, but the stabilization of its conformation by the interaction with hydrophobic surfaces or by using cosolvents (TFE, SDS) revealed a high propensity to adopt an alpha-helical folding. A 3D structure model of NC1(84-116), calculated from NMR data recorded in a TFE/water mixture, showed that the NC1-heparin binding site forms a amphipathic alpha-helix exhibiting a twisted basic groove. It is structurally similar to the consensus spatial alpha-helix model of heparin-binding [Margalit et al. (1993) J. Biol. Chem. 268, 19228-19231], except that the GAG binding domain of NC1 may be extended over 18 residues, that is, the NC1(94-111) segment. In addition, the formation of a hydrophobic groove upon helix formation suggests the contribution of additional sequences to ensure the stability of the GAG-binding domain. Overall the NC1(84-116) model exhibits a nativelike conformation which presents suitably oriented residues for the interaction with a specific GAG.

The crystal structure of benzoylformate decarboxylase at 1.6 A resolution: diversity of catalytic residues in thiamin diphosphate-dependent enzymes.

The crystal structure of the thiamin diphosphate (ThDP)-dependent enzyme benzoylformate decarboxylase (BFD), the third enzyme in the mandelate pathway of Pseudomonas putida, has been solved by multiple isomorphous replacement at 1.6 A resolution and refined to an R-factor of 15.0% (free R = 18.6%). The structure of BFD has been compared to that of other ThDP-dependent enzymes, including pyruvate decarboxylase. The overall architecture of BFD resembles that of the other family members, and cofactor- and metal-binding residues are well conserved. Surprisingly, there is no conservation of active-site residues not directly bound to the cofactor. The position of functional groups in the active site may be conserved, however. Three classes of metal ions have been identified in the BFD crystal structure: Ca2+ bound to the cofactor in each subunit, Mg2+ on a 2-fold axis of the tetramer, and Ca2+ at a crystal contact. The structure includes a non-proline cis-peptide bond and an unusually long and regular polyproline type II helix that mediates the main contact between tetramers in the crystal. The high-quality electron-density map allowed the correction of errors totaling more than 10% of the amino acid sequence, which had been predicted from the reported sequence of the mdlC gene. Analysis of the BFD structure suggests that requirements for activation of the cofactor, the nature of the reaction intermediates, and architectural considerations relating to the protein fold have been dominant forces in the evolution of ThDP-dependent enzymes.

Conformational analysis of LYS(11-36), a peptide derived from the beta-sheet region of T4 lysozyme, in TFE and SDS.

The solution conformation of a peptide LYS(11-36), which corresponds to the beta-sheet region in T4 lysozyme, has been examined in aqueous solution, TFE, and SDS micelles by CD and 1H NMR spectroscopy. Secondary structure predictions suggest some beta-sheet and turn character in aqueous solution but predict a helical conformation in a more hydrophobic environment. The predictions were supported by the CD and NMR studies which showed the peptide to be relatively unstructured in aqueous solution, although there was some evidence of a beta-turn conformer which was maintained in 200 mM SDS and, to a lesser extent, in 50% TFE. The peptide was significantly helical in the presence of either 50% TFE or 200 mM SDS. TFE and SDS titrations showed that the peptide could form helical, sheet, or extended structure depending on the TFE or SDS concentration. The studies indicate that peptide environment is the determining factor in secondary structure adopted by LYS(11-36).

Recombinant human phenylethanolamine N-methyltransferase: overproduction in Escherichia coli, purification, and characterization.

The gene encoding phenylethanolamine N-methyltransferase (PNMT) has been amplified from a human adrenal medulla cDNA library. Following ligation of the gene into a pET3a-derived expression vector and transformation into Escherichia coli BL21(DE3)pLysS, PNMT was expressed, yielding about 10% of the soluble protein. The enzyme was purified to homogeneity by ammonium sulfate fractionation followed by ion-exchange chromatography and gel filtration. The Km for phenylethanolamine and S-adenosyl-L-methionine were determined to be 130 and 16 microM, respectively. The enzyme could be inhibited by reagents expected to modify cysteine, arginine, tyrosine, and histidine residues, but not by methyl acetimidate, a reagent expected to modify lysine residues.

Pharmacokinetics of thiopental and pentobarbital enantiomers after intravenous administration of racemic thiopental.

We studied the pharmacokinetics of thiopental enantiomers in 14 healthy patients aged 37-73 yr receiving racemic thiopental by intravenous (IV) bolus or IV infusion. Plasma concentration of each enantiomer was measured by chiral high-performance liquid chromatography. After IV bolus, the total plasma clearance (CL) (295 +/- 132 mL/min) and volume of distribution at steady state (Vss) (139 +/- 38.5 L) of R-thiopental were significantly greater than those of S-thiopental (230 +/- 104 mL/min and 114 +/- 47.5 L, respectively). The plasma unbound fraction (fu) was determined by ultrafiltration of plasma from six healthy volunteers. The fu of R-thiopental (12.4% +/- 0.6%) was significantly greater than that of S-thiopental (10.0% +/- 1.0%). When the CL and Vss of the two enantiomers were corrected for the difference in mean fu, there were no significant differences between enantiomers for these variables. As the 20%-30% difference between the enantiomers in total CL and total Vss could be accounted for by stereoselectivity in fu, these differences are not likely to be clinically significant. During 105-180 min IV infusion of racemic thiopental to the other patients, there was no difference between enantiomers in mean plasma concentrations of total or unbound thiopental or total pentobarbital, a major metabolite of thiopental (P > 0.05). Therefore, it is appropriate to relate pharmacodynamic effects to racemic plasma concentrations of thiopental during IV infusion of racemic thiopental.

Determination of (R)-(+)- and (S)-(-)-isomers of thiopentone in plasma by chiral high-performance liquid chromatography.

A method for the determination of R-(+)- and (S)-(-)-isomers of thiopentone in plasma was developed. Following liquid-liquid extraction, the separation of enantiomers of thiopentone and the internal standard (racemic ketamine) was achieved by high-performance liquid chromatography on an alpha1-acid glycoprotein (AGP) column with ultraviolet detection at 280 nm. The mobile phase consisted of 20 mM KH2PO4 buffer-2-propanol-methanol (93.5:5.0:1.5) at pH 5.0. The flow-rate was 0.9 ml/min. The limit of quantification for each isomer was approximately 10 ng/ml. The assay is suitable for pharmacokinetic studies of (R)-(+)- and (S)-(-)-isomers of thiopentone, following usual bolus intravenous clinical doses of the racemic drug.

CD and NMR determination of the solution structure of a peptide corresponding to T4 lysozyme residues 38-51.

Solid phase methods have been used to synthesise a peptide corresponding to residues 38-51 of T4 lysozyme. The peptide, LYS(38-51), encompasses helix B in the crystal structure of T4 lysozyme. CD and 1H-NMR analysis showed that the peptide was unstructured in aqueous solution but adopted a helical conformation in the more hydrophobic environment provided by 50% TFE and SDS micelles. The solution structure derived from the NMR data was similar to that of the helix in the X-ray structure, although there was some fraying at the N-terminus.