Protein engineering applications of industrially exploitable enzymes: Geobacillus stearothermophilus LDH and Candida methylica FDH.

Enzymes have become important tools in several industries due to their ability to produce chirally pure and complex molecules with interesting biological properties. The NAD(+)-dependent LDH (lactate dehydrogenase) [bsLDH [Geobacillus stearothermophilus (formerly Bacillus stearothermophilus) LDH] from G. stearothermophilus and the NAD(+)-dependent FDH (formate dehydrogenase) [cmFDH (Candida methylica FDH)] enzyme from C. methylica are particularly crucial enzymes in the pharmaceutical industry and are related to each other in terms of NADH use and regeneration. LDH catalyses the interconversion of pyruvate (oxo acid) and lactate (alpha-hydroxy acid) using the NADH/NAD(+) pair as a redox cofactor. Employing LDH to reduce other oxo acids can generate chirally pure alpha-hydroxy acids of use in the production of pharmaceuticals. One important use of FDH is to regenerate the relatively expensive NADH cofactor that is used by NAD(+)-dependent oxidoreductases such as LDH. Both LDH and FDH from organisms of interest were previously cloned and overproduced. Therefore they are available at a low cost. However, both of these enzymes show disadvantages in the large-scale production of chirally pure compounds. We have applied two routes of protein engineering studies to improve the properties of these two enzymes, namely DNA shuffling and site-directed mutagenesis. Altering the substrate specificity of bsLDH by DNA shuffling and changing the coenzyme specificity of cmFDH by site-directed mutagenesis are the most successful examples of our studies. The present paper will also include the details of these examples together with some other applications of protein engineering regarding these enzymes.

NGF receptor TrkAd5: therapeutic agent and drug design target.

Biochemical studies have shown that domain 5 of the TrkA (tropomyosin receptor kinase A) receptor is involved in the binding of NGF (nerve growth factor). Crystallographic studies have confirmed this, demonstrating that one homodimer of NGF binds to two TrkAd5 molecules. TrkAd5 has been made recombinantly in Escherichia coli, purified and shown to bind NGF with picomolar affinity. We have used the co-ordinates of the crystal structure of the NGF-TrkAd5 complex to screen approximately two million compounds in silico for the identification of small molecule agonists/antagonists. Selected hits were shown to be active in an in vitro ligand-binding assay; structure-activity relationships are now being investigated. In addition, TrkAd5 has been shown to be efficacious in preclinical models of inflammatory pain and asthma by the sequestration of excess levels of endogenous NGF, and therefore represents a novel therapeutic agent.

Amino acid pairing preferences in parallel beta-sheets in proteins.

Statistical approaches have been applied to examine amino acid pairing preferences within parallel beta-sheets. The main chain hydrogen bonding pattern in parallel beta-sheets means that, for each residue pair, only one of the residues is involved in main chain hydrogen bonding with the strand containing the partner residue. We call this the hydrogen bonded (HB) residue and the partner residue the non-hydrogen bonded (nHB) residue, and differentiate between the favorability of a pair and that of its reverse pair, e.g. Asn(HB)-Thr(nHB)versus Thr(HB)-Asn(nHB). Significantly (p < or = 0.000001) favoured pairings were rationalised using stereochemical arguments. For instance, Asn(HB)-Thr(nHB) and Arg(HB)-Thr(nHB) were favoured pairs, where the residues adopted favoured chi1 rotamer positions that allowed side-chain interactions to occur. In contrast, Thr(HB)-Asn(nHB) and Thr(HB)-Arg(nHB) were not significantly favoured, and could only form side-chain interactions if the residues involved adopted less favourable chi1 conformations. The favourability of hydrophobic pairs e.g. Ile(HB)-Ile(nHB), Val(HB)-Val(nHB) and Leu(HB)-Ile(nHB) was explained by the residues adopting their most preferred chi1 and chi2 conformations, which enabled them to form nested arrangements. Cysteine-cysteine pairs are significantly favoured, although these do not form intrasheet disulphide bridges. Interactions between positively and negatively charged residues were asymmetrically preferred: those with the negatively charged residue at the HB position were more favoured. This trend was accounted for by the presence of general electrostatic interactions, which, based on analysis of distances between charged atoms, were likely to be stronger when the negatively charged residue is the HB partner. The Arg(HB)-Asp(nHB) interaction was an exception to this trend and its favorability was rationalised by the formation of specific side-chain interactions. This research provides rules that could be applied to protein structure prediction, comparative modelling and protein engineering and design. The methods used to analyse the pairing preferences are automated and detailed results are available (http://www.rubic.rdg.ac.uk/betapairprefsparallel/).

Crystal structure of Plasmodium berghei lactate dehydrogenase indicates the unique structural differences of these enzymes are shared across the Plasmodium genus.

As Plasmodium rely extensively on homolactic fermentation for energy production, Plasmodium falciparum lactate dehydrogenase (PfLDH)--the key enzyme in this process--has previously been suggested as a novel target for antimalarials. This enzyme has distinctive kinetic and structural properties that distinguish it from its human homologues. In this study, we now describe the expression, kinetic characterisation and crystal structure determination of the LDH from Plasmodium berghei. This enzyme is seen to have a similar kinetic profile to its P. falciparum counterpart, exhibiting the characteristic lack of substrate inhibition that distinguishes plasmodial from human LDHs. The crystal structure of P. berghei lactate dehydrogenase (PbLDH) shows a very similar active site arrangement to the P. falciparum enzyme. In particular, an insertion of five amino acid residues in the active site loop creates an enlarged volume in the substrate binding site, and characteristic changes in the residues lining the NADH cofactor binding pocket result in displacement of the cofactor relative to its observed position in mammalian and all other LDH structures. These results imply the special features previously described for PfLDH may be shared across the Plasmodium genus, supporting the universal application of therapeutics targeting this enzyme.

The F-box-containing protein UFO and AGAMOUS participate in antagonistic pathways governing early petal development in Arabidopsis.

The UNUSUAL FLORAL ORGANS (UFO) gene is required for multiple processes in the developing Arabidopsis flower, including the proper patterning and identity of both petals and stamens. The gene encodes an F-box-containing protein, UFO, which interacts physically and genetically with the Skp1 homolog, ASK1. In this report, we describe four ufo alleles characterized by the absence of petals, which uncover another role for UFO in promoting second whorl development. This UFO-dependent pathway is required regardless of the second whorl organ to be formed, arguing that it affects a basic process acting in parallel with those establishing organ identity. However, the pathway is dispensable in the absence of AGAMOUS (AG), a known inhibitor of petal development. In situ hybridization results argue that AG is not transcribed in the petal region, suggesting that it acts non-cell-autonomously to inhibit second whorl development in ufo mutants. These results are combined into a genetic model explaining early second whorl initiation/proliferation, in which UFO functions to inhibit an AG-dependent activity.

Intercellular adhesion molecule-4 binds alpha(4)beta(1) and alpha(V)-family integrins through novel integrin-binding mechanisms.

The LW blood group glycoprotein, ICAM-4, is a member of the intercellular adhesion molecule (ICAM) family expressed in erythroid cells. To begin to address the function of this molecule, ligands for ICAM-4 on hemopoietic and nonhemopoietic cell lines were identified. Peptide inhibition studies suggest that adhesion of cell lines to an ICAM-4-Fc construct is mediated by an LDV-inhibitable integrin on hemopoietic cells and an RGD-inhibitable integrin on nonhemopoietic cells. Antibody inhibition studies identified the hemopoietic integrin as alpha(4)beta(1.) Antibody inhibition studies on alpha(4)beta(1)-negative, nonhemopoietic cell lines suggested that adhesion of these cells is mediated by alpha(V) integrins (notably alpha(V)beta(1) and alpha(V)beta(5)). The structure of ICAM-4 modeled on the crystal structure of ICAM-2 was used to identify surface-exposed amino acid residues for site-directed mutagenesis. Neither an unusual LETS nor an LDV motif in the first domain of ICAM-4 was critical for integrin binding. ICAM-4 is the first ICAM family member shown to be a ligand for integrins other than those of the beta(2) family, and the data suggest that ICAM-4 has a novel integrin-binding site(s). These findings suggest a role for ICAM-4 in normal erythropoiesis and may also be relevant to the adhesive interactions of sickle cells.

Protochlorophyllide oxidoreductase: a homology model examined by site-directed mutagenesis.

An homology model of protochlorophyllide reductase (POR) from Synechocystis sp. was constructed on a template from the tyrosine-dependent oxidoreductase family. The model showed characteristics appropriate to a globular, soluble protein and was used to generate a structure of the ternary complex of POR, nicotinamide adenine dinucleotide phosphate (NADPH), and protochlorophyllide. The POR ternary model was validated by mutagenesis experiments involving predicted coenzyme-binding residues and by chemical modification experiments. A core tryptophan residue was shown to be responsible for much of the protein's fluorescence. Both quenching of this residue by coenzyme and fluorescence resonance energy transfer (FRET) from the protein to the coenzyme allowed the binding constant of NADPH to be determined. Replacement of this residue by Tyr gave an active mutant with approximately halved fluorescence and a negligible FRET signal, consistent with the role of this residue in energy transfer to the NADPH at the active site and with the model. The mechanism of the enzyme is discussed in the context of the model and semiempirical molecular orbital calculations.

Planned neck dissection for advanced primary head and neck malignancy treated with organ preservation therapy: disease control and survival outcomes.

BACKGROUND: The role of planned neck dissection after organ preservation therapy with radiotherapy or chemotherapy/radiotherapy for advanced head and neck cancers presenting with clinically positive neck disease is still being elucidated. The aim of this study is to review the outcomes of such patients treated by organ preservation therapy at our institution. METHODS: A retrospective chart review of 33 patients who underwent planned neck dissections after organ preservation therapy for advanced primary head and neck malignancy. Endpoints measured were disease-free survival and local, regional, and distant control. SETTING: Tertiary metropolitan medical center. RESULTS: Two-year actuarial disease-free survival was 61%, and neck control was 92%, with only two failures in the neck. The use of neoadjuvant chemotherapy and total dose of radiotherapy did not correlate with neck control or disease-free survival. The presence of pathologically positive nodal disease at the time of neck dissection did not correlate with recurrent neck disease, but was a predictor of local recurrence (p = .0086). CONCLUSIONS: Our data suggest that for patients undergoing planned neck dissection after organ preservation therapy, neck control is obtained in almost all cases. The presence of pathologically positive nodal disease at the time of surgery may have implications for the incidence of local recurrence.

Structural basis for altered activity of M- and H-isozyme forms of human lactate dehydrogenase.

Lactate dehydrogenase (LDH) interconverts pyruvate and lactate with concomitant interconversion of NADH and NAD(+). Although crystal structures of a variety of LDH have previously been described, a notable absence has been any of the three known human forms of this glycolytic enzyme. We have now determined the crystal structures of two isoforms of human LDH-the M form, predominantly found in muscle; and the H form, found mainly in cardiac muscle. Both structures have been crystallized as ternary complexes in the presence of the NADH cofactor and oxamate, a substrate-like inhibitor. Although each of these isoforms has different kinetic properties, the domain structure, subunit association, and active-site regions are indistinguishable between the two structures. The pK(a) that governs the K(M) for pyruvate for the two isozymes is found to differ by about 0.94 pH units, consistent with variation in pK(a) of the active-site histidine. The close similarity of these crystal structures suggests the distinctive activity of these enzyme isoforms is likely to result directly from variation of charged surface residues peripheral to the active site, a hypothesis supported by electrostatic calculations based on each structure. Proteins 2001;43:175-185.

Ab initio protein structure prediction using physicochemical potentials and a simplified off-lattice model.

This study describes a computational method for ab inito protein structure prediction. Protein conformation has been modeled by using six optimized backbone torsion angles and fixed side chains approximating rotationally averaged real side chains. The approximations aim to keep complexity of the structure description to a minimum without seriously compromising the accuracy of the structural representation. An evolutionary Monte Carlo algorithm has been developed to search through this restricted conformational space to locate low-energy protein structures. A simple physicochemical force field has been developed to assess the energies of different conformations within this structural description. The corresponding residue interaction energies are based on hydrophobic, hydrophilic, steric, and hydrogen-bonding potentials. The search procedure has been used to locate native energy minima from primary sequence alone. The 3-D structures of polypeptides up to 38 residues with both beta and alpha secondary structural elements have been accurately predicted. The search procedure has been found to be highly efficient and follows an energetically and structurally plausible pathway to locate native populations. The simple force field described in the study has been compared with a more complex all-atom model and been found to be similarly effective in predicting the structures of proposed independent folding units. Proteins 2001;43:186-202.

Is there a conserved interaction between cardiolipin and the type II bacterial reaction center?

In a recent publication, the structural details of an interaction between the Rhodobacter sphaeroides reaction center and the anionic phospholipid diphosphatidyl glycerol (cardiolipin) were described (K. E. McAuley, P. K. Fyfe, J. P. Ridge, N. W. Isaacs, R. J. Cogdell, and M. R. Jones, 1999, Proc. Natl. Acad. Sci. U.S.A. 96:14706-14711). This was the first crystallographic description of an interaction between this biologically important lipid and an integral membrane protein and was also the first piece of evidence that the reaction center has a specific interaction with cardiolipin. We have examined the extent to which the residues that interact with the cardiolipin are conserved in other species of photosynthetic bacteria with this type of reaction center and discuss the possibility that this cardiolipin binding site is a conserved feature of these reaction centers. We look at how sequence variations that would affect the shape of the cardiolipin binding site might affect the protein-cardiolipin interaction, by modeling the binding of cardiolipin to the reaction center from Rhodopseudomonas viridis.

Ultrasound-guided fine-needle aspiration and thyroid disease.

BACKGROUND: Fine-needle aspiration represents a critical diagnostic test in determining proper management of thyroid disease and the use of ultrasound-guided fine-needle aspiration (USGFNA) has increased over the years. METHODS: A retrospective chart review of patients undergoing USGFNA. Two hundred fifteen patients underwent 234 procedures with 362 nodules aspirated within a 2 (1/2)-year period. RESULTS: The mean ages of women and men were 51.9 and 57.8, respectively. The average size of nodules was 2.1 cm. A difficult to assess gland or nodule was the most common indication for USGFNA (33%). The sensitivity was 88.2%, specificity was 80.0%, the PPV was 65.2%, the negative predictive value was 94.1%, and the accuracy was 82.5%. The cancer yield, inadequacy, and complication rates were 44%, 10.5%, and 8.5%, respectively. CONCLUSIONS: USGFNA aspiration is a safe and effective diagnostic modality in the management of thyroid disease, especially for nodules that are difficult to palpate.

The asparagine-stabilized beta-turn of apamin: contribution to structural stability from dynamics simulation and amide hydrogen exchange analysis.

Molecular dynamics simulations of bee venom apamin, and an analogue having an Asn to Ala substitution at residue 2 (apamin-N2A), were analyzed to explore the contribution of hydrogen bonds involving Asn2 to local (beta-turn residues N2, C3, K4, A5) and global stability. The wild-type peptide retained a stable conformation during 2.4 ns of simulation at 67 degrees C, with high beta-turn stability characterized by backbone-side chain hydrogen bonds involving beta-turn residues K4 and A5, with the N2 side chain amide carbonyl. The loss of stabilizing interactions involving the N2 side chain resulted in the loss of the beta-turn conformation in the apamin N2A simulations (27 or 67 degrees C). This loss of beta-turn stability propagates throughout the peptide structure, with destabilization of the C-terminal helix connected to the N-terminal region by two disulfide bonds. Backbone stability in a synthetic peptide analogue (apamin-N2A) was characterized by NMR and amide hydrogen exchange measurements. Consistent with the simulations, loss of hydrogen bonds involving the N2 side chain resulted in destabilization of both the N-terminal beta-turn and the C-terminal helix. Amide exchange protection factors in the C-terminal helix were reduced by 9-11-fold in apamin N2A as compared with apamin, corresponding to free energy (deltaDeltaG(uf)) of around 1.5 kcal M(-1) at 20 degrees C. This is equivalent to the contribution of hydrogen bond interactions involving the N2 side chain to the stability of the beta-turn. Together with additional measures of exchange protection factors, the three main contributions to backbone stability in apamin that account for virtually the full thermodynamic stability of the peptide have been quantitated.

New sequence data enable modelling of the fungal alternative oxidase and explain an absence of regulation by pyruvate.

Respiratory rates involving the alternative oxidase (AO) were studied in mitochondria from Tapesia acuformis. There was no evidence for regulation by pyruvate, in contrast with plant AO. The site of interaction of pyruvate with the plant AO is a conserved cysteine. The primary sequence was obtained for AO from Magnaporthe grisea and compared with four published sequences for fungal AO. In all cases this cysteine was absent. Sequence data were obtained for the C-terminal domain of a further five fungal AOs. In this region the fungal sequences were all consistent with a four-helix, di-iron binding structure as in the ferritin-fold family. A molecular model of this domain was deduced from the structure of Delta-9 desaturase. This is in general agreement with that developed for plant AOs, despite very low sequence identity between the two kingdoms. Further modelling indicated an appropriate active site for binding of ubiquinol, required in the AO redox reaction.

Controlled induction of GUS marked clonal sectors in Arabidopsis.

Stably transformed Arabidopsis lines in which GUS marked cell clones are readily produced in response to heat-shock have been established and characterized. Control of GUS activation is achieved by heat-shock-induced FLP recombinase activity which "switches on" expression of a GUS marker gene previously held transcriptionally silent. To obtain efficient GUS sectoring, single insert Arabidopsis lines carrying FLP recombinase under the control of a heat-shock-inducible promoter and an FLP-activatable GUS construct were generated. Analysis of GUS sectoring in lines hemizygous and homozygous for both inserts was conducted after various regimes of heat-shock were given at various developmental stages. It is shown that GUS sectoring events can be efficiently induced in most vegetative, aerial and sexual structures in Arabidopsis. Furthermore, the frequency of sectoring events, sector size and, to some extent, the tissues in which sectors are generated can be readily controlled by choice of the conditions and timing of heat-shock used.

Progressive stabilization of intermediate and transition states in protein folding reactions by introducing surface hydrophobic residues.

It can be argued from the principle of solvent exclusion that the introduction of hydrophobic residues onto the surface of a protein will not destabilize the folded state because the nonpolar side chain will be at least as exposed in the unfolded state as it is when the protein chain is folded. A comparison of the folding pathway of wild type and 11 site-directed mutants of CD2.d1 shows this to be true. In fact, owing to partial burial of nonpolar groups as folding proceeds, we find that the rapidly formed intermediate state and, to a greater extent, the transition state are generally stabilized by hydrophobic surface mutations. This effect is slightly moderated in the folded state presumably by the perturbation of van der Waals' contacts and/or local electrostatic interactions that have a greater influence in this fully compact structure. The fact that in all but one case we find that stabilization of the rapidly collapsed intermediate is accompanied by a faster acquisition of the folded state refutes the argument that I states are generally "off pathway" conformations or ensembles that lead to the inhibition of otherwise more rapid folding trajectories.

Molecular modelling and site-directed mutagenesis of the inositol 1,3,4,5-tetrakisphosphate-binding pleckstrin homology domain from the Ras GTPase-activating protein GAP1IP4BP.

GAP1(IP4BP) is a Ras GTPase-activating protein (GAP) that in vitro is regulated by the cytosolic second messenger inositol 1,3,4,5-tetrakisphosphate [Ins(1,3,4,5)P(4)]. We have studied Ins(1,3,4,5)P(4) binding to GAP1(IP4BP), and shown that the inositol phosphate specificity and binding affinity are similar to Ins(1,3,4,5)P(4) binding to Bruton's tyrosine kinase (Btk), evidence which suggests a similar mechanism for Ins(1,3,4,5)P(4) binding. The crystal structure of the Btk pleckstrin homology (PH) domain in complex with Ins(1,3,4,5)P(4) has shown that the binding site is located in a partially buried pocket between the beta 1/beta 2- and beta 3/beta 4-loops. Many of the residues involved in the binding are conserved in GAP1(IP4BP). Therefore we generated a model of the PH domain of GAP1(IP4BP) in complex with Ins(1,3,4,5)P(4) based on the Btk-Ins(1,3,4,5)P(4) complex crystal structure. This model had the typical PH domain fold, with the proposed binding site modelling well on the Btk structure. The model has been verified by site-directed mutagenesis of various residues in and around the proposed binding site. These mutations have markedly reduced affinity for Ins(1,3,4,5)P(4), indicating a specific and tight fit for the substrate. The model can also be used to explain the specificity of inositol phosphate binding.

Tall cell variant: an aggressive form of papillary thyroid carcinoma.

Twenty-four cases of the tall cell variant (TCV), a subset of papillary thyroid carcinoma, were identified in a group of 624 patients with thyroid cancer. All pathology specimens were reviewed, and each patient's carcinoma was categorized according to characteristics on presentation, local recurrence, distant metastases, follow-up, and tumor-related mortality. The TCV group was compared with a historical control group (Mazzaferri and Jhiang: 1355 patients). The TCV group had a statistically higher percentage of stage 3 and 4 carcinoma, extrathyroidal invasion, and tumor size less than 1.5 cm than the control group. There was no statistical relationship between age greater than 50 years and stage in the TCV group. No relationship could be found between TCV histology and recurrence or mortality. These findings, combined with those of studies that link stage on presentation to poor outcomes, have led to our conclusion that TCV is an aggressive malignancy warranting appropriate treatment and close follow-up.