Redesign of soluble fatty acid desaturases from plants for altered substrate specificity and double bond position

Acyl-acyl carrier protein (ACP) desaturases introduce double bonds at specific positions in fatty acids of defined chain lengths and are one of the major determinants of the monounsaturated fatty acid composition of vegetable oils. Mutagenesis studies were conducted to determine the structural basis for the substrate and double bond positional specificities displayed by acyl-ACP desaturases. By replacement of specific amino acid residues in a Δ6-palmitoyl (16:0)-ACP desaturase with their equivalents from a Δ9-stearoyl (18:0)-ACP desaturase, mutant enzymes were identified that have altered fatty acid chain-length specificities or that can insert double bonds into either the Δ6 or Δ9 positions of 16:0- and 18:0-ACP. Most notably, by replacement of five amino acids (A181T/A200F/S205N/L206T/G207A), the Δ6-16:0-ACP desaturase was converted into an enzyme that functions principally as a Δ9-18:0-ACP desaturase. Many of the determinants of fatty acid chain-length specificity in these mutants are found in residues that line the substrate binding channel as revealed by x-ray crystallography of the Δ9-18:0-ACP desaturase. The crystallographic model of the active site is also consistent with the diverged activities associated with naturally occurring variant acyl-ACP desaturases. In addition, on the basis of the active-site model, a Δ9-18:0-ACP desaturase was converted into an enzyme with substrate preference for 16:0-ACP by replacement of two residues (L118F/P179I). These results demonstrate the ability to rationally modify acyl-ACP desaturase activities through site-directed mutagenesis and represent a first step toward the design of acyl-ACP desaturases for the production of novel monounsaturated fatty acids in transgenic oilseed crops.

The change in hydrogen bond strength accompanying charge rearrangement: Implications for enzymatic catalysis

The equilibrium for formation of the intramolecular hydrogen bond (KHB) in a series of substituted salicylate monoanions was investigated as a function of ΔpKa, the difference between the pKa values of the hydrogen bond donor and acceptor, in both water and dimethyl sulfoxide. The dependence of log KHB upon ΔpKa is linear in both solvents, but is steeper in dimethyl sulfoxide (slope = 0.73) than in water (slope = 0.05). Thus, hydrogen bond strength can undergo substantially larger increases in nonaqueous media than aqueous solutions as the charge density on the donor or acceptor atom increases. These results support a general mechanism for enzymatic catalysis, in which hydrogen bonding to a substrate is strengthened as charge rearranges in going from the ground state to the transition state; the strengthening of the hydrogen bond would be greater in a nonaqueous enzymatic active site than in water, thus providing a rate enhancement for an enzymatic reaction relative to the solution reaction. We suggest that binding energy of an enzyme is used to fix the substrate in the low-dielectric active site, where the strengthening of the hydrogen bond in the course of a reaction is increased.

Unusual 1H NMR chemical shifts support (His) Cɛ1—H⋅⋅⋅O⩵C H-bond: Proposal for reaction-driven ring flip mechanism in serine protease catalysis

13C-selective NMR, combined with inhibitor perturbation experiments, shows that the Cɛ1—H proton of the catalytic histidine in resting α-lytic protease and subtilisin BPN′ resonates, when protonated, at 9.22 ppm and 9.18 ppm, respectively, which is outside the normal range for such protons and ≈0.6 to 0.8 ppm further downfield than previously reported. They also show that the previous α-lytic protease assignments [Markley, J. L., Neves, D. E., Westler, W. M., Ibanez, I. B., Porubcan, M. A. & Baillargeon, M. W. (1980) Front. Protein Chem. 10, 31–61] were to signals from inactive or denatured protein. Simulations of linewidth vs. pH demonstrate that the true signal is more difficult to detect than corresponding signals from inactive derivatives, owing to higher imidazole pKa values and larger chemical shift differences between protonated and neutral forms. A compilation and analysis of available NMR data indicates that the true Cɛ1—H signals from other serine proteases are similarly displaced downfield, with past assignments to more upfield signals probably in error. The downfield displacement of these proton resonances is shown to be consistent with an H-bond involving the histidine Cɛ1—H as donor, confirming the original hypothesis of Derewenda et al. [Derewenda, Z. S., Derewenda, U. & Kobos, P. M. (1994) J. Mol. Biol. 241, 83–93], which was based on an analysis of literature x-ray crystal structures of serine hydrolases. The invariability of this H-bond among enzymes containing Asp-His-Ser triads indicates functional importance. Here, we propose that it enables a reaction-driven imidazole ring flip mechanism, overcoming a major dilemma inherent in all previous mechanisms, namely how these enzymes catalyze both the formation and productive breakdown of tetrahedral intermediates.

Forced unfolding modulated by disulfide bonds in the Ig domains of a cell adhesion molecule

Cell adhesion molecules (CAMs) mediate cell attachment and stress transfer through extracellular domains. Here we forcibly unfold the Ig domains of a prototypical Ig superfamily CAM that contains intradomain disulfide bonds. The Ig domains of all such CAMs have conformations homologous to cadherin extracellular domains, titin Ig-type domains, and fibronectin type-III (FNIII) domains. Atomic force microscopy has been used to extend the five Ig domains of Mel-CAM (melanoma CAM)—a protein that is overexpressed in metastatic melanomas—under conditions where the disulfide bonds were either left intact or disrupted through reduction. Under physiological conditions where intradomain disulfide bonds are intact, partial unfolding was observed at forces far smaller than those reported previously for either titin's Ig-type domains or tenascin's FNIII domains. This partial unfolding under low force may be an important mechanism for imparting elasticity to cell–cell contacts, as well as a regulatory mechanism for adhesive interactions. Under reducing conditions, Mel-CAM's Ig domains were found to fully unfold through a partially folded state and at slightly higher forces. The results suggest that, in divergent evolution of all such domains, stabilization imparted by disulfide bonds relaxes requirements for strong, noncovalent, folded-state interactions.

Locking in alternate conformations of the integrin αLβ2 I domain with disulfide bonds reveals functional relationships among integrin domains

We used integrin αLβ2 heterodimers containing I domains locked open (active) or closed (inactive) with disulfide bonds to investigate regulatory interactions among domains in integrins. mAbs to the αL I domain and β2 I-like domain inhibit adhesion of wild-type αLβ2 to intercellular adhesion molecule-1. However, with αLβ2 containing a locked open I domain, mAbs to the I domain were subdivided into subsets (i) that did not inhibit, and thus appear to inhibit by favoring the closed conformation, and (ii) that did inhibit, and thus appear to bind to the ligand binding site. Furthermore, αLβ2 containing a locked open I domain was completely resistant to inhibition by mAbs to the β2 I-like domain, but became fully susceptible to inhibition after disulfide reduction with DTT. This finding suggests that the I-like domain indirectly contributes to ligand binding by regulating opening of the I domain in wild-type αLβ2. Conversely, locking the I domain closed partially restrained conformational change of the I-like domain by Mn2+, as measured with mAb m24, which we map here to the β2 I-like domain. By contrast, locking the I domain closed or open did not affect constitutive or Mn2+-induced exposure of the KIM127 epitope in the β2 stalk region. Furthermore, locked open I domains, in αLβ2 complexes or expressed in isolation on the cell surface, bound to intercellular adhesion molecule-1 equivalently in Mg2+ and Mn2+. These results suggest that Mn2+ activates αLβ2 by binding to a site other than the I domain, most likely the I-like domain of β2.

Activation of αvβ3-Vitronectin Binding Is a Multistage Process in which Increases in Bond Strength Are Dependent on Y747 and Y759 in the Cytoplasmic Domain of β3

Integrin receptors serve as mechanical links between the cell and its structural environment. Using αvβ3 integrin expressed in K562 cells as a model system, the process by which the mechanical connection between αvβ3 and vitronectin develops was analyzed by measuring the resistance of these bonds to mechanical separation. Three distinct stages of activation, as defined by increases in the αvβ3-vitronectin binding strength, were defined by mutational, biochemical, and biomechanical analyses. Activation to the low binding strength stage 1 occurs through interaction with the vitronectin ligand and leads to the phosphorylation of Y747 in the β3 subunit. Stage 2 is characterized by a 4-fold increase in binding strength and is dependent on stage1 and the phosphorylation of Y747. Stage 3 is characterized by a further 2.5-fold increase in binding strength and is dependent on stage 2 events and the availability of Y759 for interaction with cellular proteins. The Y747F mutant blocked the transition from stage 1 to stage 2, and the Y759F blocked the transition from stage 2 to stage 3. The data suggest a model for tension-induced activation of αvβ3 integrin.

Detecting the conformational change of transmembrane signaling in a bacterial chemoreceptor by measuring effects on disulfide cross-linking in vivo.

Transmembrane signaling by bacterial chemoreceptors is thought to involve relative movement among the four transmembrane helices of the homodimer. We assayed that movement by measuring effects of ligand occupancy on rates of oxidative cross-linking between cysteines introduced into neighboring helices of the transmembrane domain of chemoreceptor Trg from Escherichia coli. Measurements were done on chemoreceptors in their native environment, intact cells that were motile and chemotactically responsive. Receptor occupancy did not appear to cause drastic rearrangement of the four-helix structure since, among 67 cysteine pairs tested, the same 19 exhibited oxidative cross-linking in the presence or absence of saturating chemoattractant. However, occupancy did cause subtle changes that were detected as effects on rates of cross-linking. Among the seven disulfides appropriate for measurements of initial rates of formation, ligand occupancy had significant and different effects on all three cross-links that connected the two helices within a subunit but had minimal effects on the four that spanned the packing interface between subunits. This constitutes direct evidence that the conformational change of transmembrane signaling involves significant movement within a subunit and minimal movement between subunits, a pattern deduced from several previous studies and now documented directly. Among possible modes of movement between the two helices of a subunit, axial sliding of one helix relative to the other was the conformational change that best accounted for the observed effects on cross-linking.

The fifth epidermal growth factor-like domain of thrombomodulin does not have an epidermal growth factor-like disulfide bonding pattern.

The disulfide bonding pattern of the fourth and fifth epidermal growth factor (EGF)-like domains within the smallest active fragment of thrombomodulin have been determined. In previous work, this fragment was expressed and purified to homogeneity, and its cofactor activity, as measured by Kcat for thrombin activation of protein C, was the same as that for full-length thrombomodulin. CNBr cleavage at the single methionine in the connecting region between the domains and subsequent deglycosylation yielded the individual EGF-like domains. The disulfide bonds were mapped by partial reduction with tris(2-carboxyethyl)phosphine according to the method of Gray [Gray, W. R. (1993) Protein Sci. 2, 1732-1748], which provides unambiguous results. The disulfide bonding pattern of the fourth EGF-like domain was (1-3, 2-4, 5-6), which is the same as that found previously in EGF and in a synthetic version of the fourth EGF-like domain. Surprisingly, the disulfide bonding pattern of the fifth domain was (1-2, 3-4, 5-6), which is unlike that found in EGF or in any other EGF-like domain analyzed so far. This result is in line with an earlier observation that the (1-2, 3-4, 5-6) isomer bound to thrombin more tightly than the EGF-like (1-3, 2-4, 5-6) isomer. The observation that not all EGF-like domains have an EGF-like disulfide bonding pattern reveals an additional element of diversity in the structure of EGF-like domains.

NMR evidence for the participation of a low-barrier hydrogen bond in the mechanism of delta 5-3-ketosteroid isomerase.

Delta 5-3-Ketosteroid isomerase (EC 5.3.3.1) promotes an allylic rearrangement involving intramolecular proton transfer via a dienolic intermediate. This enzyme enhances the catalytic rate by a factor of 10(10). Two residues, Tyr-14, the general acid that polarizes the steroid 3-carbonyl group and facilitates enolization, and Asp-38 the general base that abstracts and transfers the 4 beta-proton to the 6 beta-position, contribute 10(4.7) and 10(5.6) to the rate increase, respectively. A major mechanistic enigma is the huge disparity between the pKa values of the catalytic groups and their targets. Upon binding of an analog of the dienolate intermediate to isomerase, proton NMR detects a highly deshielded resonance at 18.15 ppm in proximity to aromatic protons, and with a 3-fold preference for protium over deuterium (fractionation factor, phi = 0.34), consistent with formation of a short, strong (low-barrier) hydrogen bond to Tyr-14. The strength of this hydrogen bond is estimated to be at least 7.1 kcal/mol. This bond is relatively inaccessible to bulk solvent and is pH insensitive. Low-barrier hydrogen bonding of Tyr-14 to the intermediate, in conjunction with the previously demonstrated tunneling contribution to the proton transfer by Asp-38, provide a plausible and quantitative explanation for the high catalytic power of this isomerase.

Intermolecular disulfide bonds stabilize VirB7 homodimers and VirB7/VirB9 heterodimers during biogenesis of the Agrobacterium tumefaciens T-complex transport apparatus.

The Agrobacterium tumefaciens VirB7 lipoprotein contributes to the stabilization of VirB proteins during biogenesis of the putative T-complex transport apparatus. Here, we report that stabilization of VirB7 itself is correlated with its ability to form disulfide cross-linked homodimers via a reactive Cys-24 residue. Three types of beta-mercaptoethanol-dissociable complexes were visualized with VirB7 and/or a VirB7::PhoA41 fusion protein: (i) a 9-kDa complex corresponding in size to a VirB7 homodimer, (ii) a 54-kDa complex corresponding in size to a VirB7/VirB7::PhoA41 mixed dimer, and (iii) a 102-kDa complex corresponding to a VirB7::PhoA41 homodimer. A VirB7C24S mutant protein was immunologically undetectable, whereas the corresponding VirB7C24S::PhoA41 derivative accumulated to detectable levels but failed to form dissociable homodimers or mixed dimers with wild-type VirB7. We further report that VirB7-dependent stabilization of VirB9 is correlated with the ability of these two proteins to dimerize via formation of a disulfide bridge between reactive Cys-24 and Cys-262 residues, respectively. Two types of dissociable complexes were visualized: (i) a 36-kDa complex corresponding in size to a VirB7/VirB9 heterodimer and (ii) an 84-kDa complex corresponding in size to a VirB7/VirB9::PhoA293 heterodimer. A VirB9C262S mutant protein was immunologically undetectable, whereas the corresponding VirB9C262S::PhoA293 derivative accumulated to detectable levels but failed to form dissociable heterodimers with wild-type VirB7. Taken together, these results support a model in which the formation of disulfide cross-linked VirB7 dimers represent critical early steps in the biogenesis of the T-complex transport apparatus.

Bond orientational order, molecular motion, and free energy of high-density DNA mesophases.

By equilibrating condensed DNA arrays against reservoirs of known osmotic stress and examining them with several structural probes, it has been possible to achieve a detailed thermodynamic and structural characterization of the change between two distinct regions on the liquid-crystalline phase diagram: (i) a higher density hexagonally packed region with long-range bond orientational order in the plane perpendicular to the average molecular direction and (ii) a lower density cholesteric region with fluid-like positional order. X-ray scattering on highly ordered DNA arrays at high density and with the helical axis oriented parallel to the incoming beam showed a sixfold azimuthal modulation of the first-order diffraction peak that reflects the macroscopic bond-orientational order. Transition to the less-dense cholesteric phase through osmotically controlled swelling shows the loss of this bond orientational order, which had been expected from the change in optical birefringence patterns and which is consistent with a rapid onset of molecular positional disorder. This change in order was previously inferred from intermolecular force measurements and is now confirmed by 31P NMR. Controlled reversible swelling and compaction under osmotic stress, spanning a range of densities between approximately 120 mg/ml to approximately 600 mg/ml, allow measurement of the free-energy changes throughout each phase and at the phase transition, essential information for theories of liquid-crystalline states.

The in vitro ejection of zinc from human immunodeficiency virus (HIV) type 1 nucleocapsid protein by disulfide benzamides with cellular anti-HIV activity.

Several disulfide benzamides have been shown to possess wide-spectrum antiretroviral activity in cell culture at low micromolar to submicromolar concentrations, inhibiting human immunodeficiency virus (HIV) type 1 (HIV-1) clinical and drug-resistant strains along with HIV-2 and simian immunodeficiency virus [Rice, W. G., Supko, J. G., Malspeis, L., Buckheit, R. W., Jr., Clanton, D., Bu, M., Graham, L., Schaeffer, C. A., Turpin, J. A., Domagala, J., Gogliotti, R., Bader, J. P., Halliday, S. M., Coren, L., Sowder, R. C., II, Arthur, L. O. & Henderson, L. E. (1995) Science 270, 1194-1197]. Rice and coworkers have proposed that the compounds act by "attacking" the two zinc fingers of HIV nucleocapsid protein. Shown here is evidence that low micromolar concentrations of the anti-HIV disulfide benzamides eject zinc from HIV nucleocapsid protein (NCp7) in vitro, as monitored by the zinc-specific fluorescent probe N-(6-methoxy-8-quinoyl)-p-toluenesulfonamide (TSQ). Structurally similar disulfide benzamides that do not inhibit HIV-1 in culture do not eject zinc, nor do analogs of the antiviral compounds with the disulfide replaced with a methylene sulfide. The kinetics of NCp7 zinc ejection by disulfide benzamides were found to be nonsaturable and biexponential, with the rate of ejection from the C-terminal zinc finger 7-fold faster than that from the N-terminal. The antiviral compounds were found to inhibit the zinc-dependent binding of NCp7 to HIV psi RNA, as studied by gel-shift assays, and the data correlated well with the zinc ejection data. Anti-HIV disulfide benzamides specifically eject NCp7 zinc and abolish the protein's ability to bind psi RNA in vitro, providing evidence for a possible antiretroviral mechanism of action of these compounds. Congeners of this class are under advanced preclinical evaluation as a potential chemotherapy for acquired immunodeficiency syndrome.

Structure-activity relationships derived by machine learning: the use of atoms and their bond connectivities to predict mutagenicity by inductive logic programming.

We present a general approach to forming structure-activity relationships (SARs). This approach is based on representing chemical structure by atoms and their bond connectivities in combination with the inductive logic programming (ILP) algorithm PROGOL. Existing SAR methods describe chemical structure by using attributes which are general properties of an object. It is not possible to map chemical structure directly to attribute-based descriptions, as such descriptions have no internal organization. A more natural and general way to describe chemical structure is to use a relational description, where the internal construction of the description maps that of the object described. Our atom and bond connectivities representation is a relational description. ILP algorithms can form SARs with relational descriptions. We have tested the relational approach by investigating the SARs of 230 aromatic and heteroaromatic nitro compounds. These compounds had been split previously into two subsets, 188 compounds that were amenable to regression and 42 that were not. For the 188 compounds, a SAR was found that was as accurate as the best statistical or neural network-generated SARs. The PROGOL SAR has the advantages that it did not need the use of any indicator variables handcrafted by an expert, and the generated rules were easily comprehensible. For the 42 compounds, PROGOL formed a SAR that was significantly (P < 0.025) more accurate than linear regression, quadratic regression, and back-propagation. This SAR is based on an automatically generated structural alert for mutagenicity.

Evidence for DNA phosphate backbone alkylation and cleavage by pyrrolo[1,2-a]benzimidazoles: small molecules capable of causing base-pair-specific phosphodiester bond hydrolysis.

This report presents evidence that a reduced pyrrolo[1,2-a]benzimidazole (PBI) cleaves DNA as a result of phosphate alkylation followed by hydrolysis of the resulting phosphate triester. The base-pair specificity of the phosphate alkylation results from Hoogsteen-type hydrogen bonding of the reduced PBI in the major groove at only A.T and G.C base pairs. Alkylated phosphates were detected by 31P NMR and the cleavage products were detected by 1H NMR and HPLC. Evidence is also presented that a reduced PBI interacts with DNA in the major groove rather than in the minor groove or by intercalation.

Enhanced conformational sampling in Monte Carlo simulations of proteins: application to a constrained peptide.

A Monte Carlo simulation method for globular proteins, called extended-scaled-collective-variable (ESCV) Monte Carlo, is proposed. This method combines two Monte Carlo algorithms known as entropy-sampling and scaled-collective-variable algorithms. Entropy-sampling Monte Carlo is able to sample a large configurational space even in a disordered system that has a large number of potential barriers. In contrast, scaled-collective-variable Monte Carlo provides an efficient sampling for a system whose dynamics is highly cooperative. Because a globular protein is a disordered system whose dynamics is characterized by collective motions, a combination of these two algorithms could provide an optimal Monte Carlo simulation for a globular protein. As a test case, we have carried out an ESCV Monte Carlo simulation for a cell adhesive Arg-Gly-Asp-containing peptide, Lys-Arg-Cys-Arg-Gly-Asp-Cys-Met-Asp, and determined the conformational distribution at 300 K. The peptide contains a disulfide bridge between the two cysteine residues. This bond mimics the strong geometrical constraints that result from a protein's globular nature and give rise to highly cooperative dynamics. Computation results show that the ESCV Monte Carlo was not trapped at any local minimum and that the canonical distribution was correctly determined.