Long-range structural effects in a second-site revertant of a mutant dihydrofolate reductase

Influence of solution ionic strength on the stabilities of M20 loop conformations in apo E. coli dihydrofolate reductase

Interactions among ions and their specific interactions with macromolecular solutes are known to play a central role in biomolecular stability. However, similar effects in the conformational stability of protein loops that play functional roles, such as binding ligands, proteins, and DNA/RNA molecules, remain relatively unexplored. A well-characterized enzyme that has such a functional loop is Escherichia coli dihydrofolate reductase (ecDHFR), whose so-called M20 loop has been observed in three ordered conformations in crystal structures. To explore how solution ionic strengths may affect the M20 loop conformation, we proposed a reaction coordinate that could quantitatively describe the loop conformation and used it to classify the loop conformations in representative ecDHFR x-ray structures crystallized in varying ionic strengths. The Protein Data Bank survey indicates that at ionic strengths (I) below the intracellular ion concentration-derived ionic strength in E. coli (I ≤ 0.237M), the ecDHFR M20 loop tends to adopt open/closed conformations, and rarely an occluded loop state, but when I is >0.237M, the loop tends to adopt closed/occluded conformations. Distance-dependent electrostatic potentials around the most mobile M20 loop region from molecular dynamics simulations of ecDHFR in equilibrated CaCl₂ solutions of varying ionic strengths show that high ionic strengths (I = 0.75/1.5M) can preferentially stabilize the loop in closed/occluded conformations. These results nicely correlate with conformations derived from ecDHFR structures crystallized in varying ionic strengths. Altogether, our results suggest caution in linking M20 loop conformations derived from crystal structures solved at ionic strengths beyond that tolerated by E. coli to the ecDHFR function.

Directed evolution of multiple genomic loci allows the prediction of antibiotic resistance

Antibiotic development is frequently plagued by the rapid emergence of drug resistance. However, assessing the risk of resistance development in the preclinical stage is difficult. Standard laboratory evolution approaches explore only a small fraction of the sequence space and fail to identify exceedingly rare resistance mutations and combinations thereof. Therefore, new rapid and exhaustive methods are needed to accurately assess the potential of resistance evolution and uncover the underlying mutational mechanisms. Here, we introduce directed evolution with random genomic mutations (DIvERGE), a method that allows an up to million-fold increase in mutation rate along the full lengths of multiple predefined loci in a range of bacterial species. In a single day, DIvERGE generated specific mutation combinations, yielding clinically significant resistance against trimethoprim and ciprofloxacin. Many of these mutations have remained previously undetected or provide resistance in a species-specific manner. These results indicate pathogen-specific resistance mechanisms and the necessity of future narrow-spectrum antibacterial treatments. In contrast to prior claims, we detected the rapid emergence of resistance against gepotidacin, a novel antibiotic currently in clinical trials. Based on these properties, DIvERGE could be applicable to identify less resistance-prone antibiotics at an early stage of drug development. Finally, we discuss potential future applications of DIvERGE in synthetic and evolutionary biology.

The cytochrome b Zn binding amino acid residue histidine 291 is essential for ubihydroquinone oxidation at the Qo site of bacterial cytochrome bc1

The ubiquinol:cytochrome (cyt) c oxidoreductase (or cyt bc₁) is an important membrane protein complex in photosynthetic and respiratory energy transduction. In bacteria such as Rhodobacter capsulatus it is constituted of three subunits: the iron-sulfur protein, cyt b and cyt c₁, which form two catalytic domains, the Q_o (hydroquinone (QH₂) oxidation) and Q_i (quinone (Q) reduction) sites. At the Q_o site, the pathways of bifurcated electron transfers emanating from QH₂ oxidation are known, but the associated proton release routes are not well defined. In energy transducing complexes, Zn²⁺ binding amino acid residues often correlate with proton uptake or release pathways. Earlier, using combined EXAFS and structural studies, we identified Zn coordinating residues of mitochondrial and bacterial cyt bc₁. In this work, using the genetically tractable bacterial cyt bc₁, we substituted each of the proposed Zn binding residues with non-protonatable side chains. Among these mutants, only the His291Leu substitution destroyed almost completely the Q_o site catalysis without perturbing significantly the redox properties of the cofactors or the assembly of the complex. In this mutant, which is unable to support photosynthetic growth, the bifurcated electron transfer reactions that result from QH₂ oxidation at the Q_o site, as well as the associated proton(s) release, were dramatically impaired. Based on these findings, on the putative role of His291 in liganding Zn, and on its solvent exposed and highly conserved position, we propose that His291 of cyt b is critical for proton release associated to QH₂ oxidation at the Q_o site of cyt bc₁.

The structure of Na⁺-translocating of NADH:ubiquinone oxidoreductase of Vibrio cholerae: implications on coupling between electron transfer and Na⁺ transport

The Na⁺-translocating NADH:ubiquinone oxidoreductase (Na⁺-NQR) of Vibrio cholerae is a respiratory complex that couples the exergonic oxidation of NADH to the transport of Na⁺ across the cytoplasmic membrane. It is composed of six different subunits, NqrA, NqrB, NqrC, NqrD, NqrE, and NqrF, which harbor FAD, FMN, riboflavin, quinone, and two FeS centers as redox co-factors. We recently determined the X-ray structure of the entire Na⁺-NQR complex at 3.5-Å resolution and complemented the analysis by high-resolution structures of NqrA, NqrC, and NqrF. The position of flavin and FeS co-factors both at the cytoplasmic and the periplasmic side revealed an electron transfer pathway from cytoplasmic subunit NqrF across the membrane to the periplasmic NqrC, and via NqrB back to the quinone reduction site on cytoplasmic NqrA. A so far unknown Fe site located in the midst of membrane-embedded subunits NqrD and NqrE shuttles the electrons over the membrane. Some distances observed between redox centers appear to be too large for effective electron transfer and require conformational changes that are most likely involved in Na⁺ transport. Based on the structure, we propose a mechanism where redox induced conformational changes critically couple electron transfer to Na⁺ translocation from the cytoplasm to the periplasm through a channel in subunit NqrB.

Systematic study of the genetic response of a variable virus to the introduction of deleterious mutations in a functional capsid region

We have targeted the intersubunit interfaces in the capsid of foot-and-mouth disease virus to investigate the genetic response of a variable virus when individual deleterious mutations are systematically introduced along a functionally defined region of its genome. We had previously found that the individual truncation (by mutation to alanine) of 28 of the 42 amino acid side chains per protomer involved in interactions between capsid pentameric subunits severely impaired infectivity. We have now used viral RNAs individually containing each of those 28 deleterious mutations (or a few others) to carry out a total of 96 transfections of susceptible cells, generally followed by passage(s) of the viral progeny in cell culture. The results revealed a very high frequency of fixation in the capsid of second-site, stereochemically diverse substitutions that compensated for the detrimental effect of primary substitutions at many different positions. Most second-site substitutions occurred at or near the capsid interpentamer interfaces and involved residues that are spatially very close to the originally substituted residue. However, others occurred far from the primary substitution, and even from the interpentamer interfaces. Remarkably, most second-site substitutions involved only a few capsid residues, which acted as "second-site hot spots." Substitutions at these hot spots compensated for the deleterious effects of many different replacements at diverse positions. The remarkable capacity of the virus to respond to the introduction of deleterious mutations in the capsid with the frequent fixation of diverse second-site mutations, and the existence of second-site hot spots, may have important implications for virus evolution.

Deterministic, compensatory mutational events in the capsid of foot-and-mouth disease virus in response to the introduction of mutations found in viruses from persistent infections

The evolution of foot-and-mouth disease virus (FMDV) (biological clone C-S8c1) in persistently infected cells led to the emergence of a variant (R100) that displayed increased virulence, reduced stability, and other modified phenotypic traits. Some mutations fixed in the R100 genome involved a cluster of highly conserved residues around the capsid pores that participate in interactions with each other and/or between capsid protomers. We have investigated phenotypic and genotypic changes that occurred when these replacements were introduced into the C-S8c1 capsid. The C3007V and M3014L mutations exerted no effect on plaque size or viral yield during lytic infections, or on virion stability, but led to a reduction in biological fitness; the D3009A mutation caused drastic reductions in plaque size and viability. Remarkably, competition of the C3007V mutant with the nonmutated virus invariably resulted in the fixation of the D3009A mutation in the C3007V capsid. In turn, the presence of the D3009A mutation invariably led to the fixation of the M3014L mutation. In both cases, two individually disadvantageous mutations led, together, to an increase in fitness, as the double mutants outcompeted the nonmutated genotype. The higher fitness of C3007V/D3009A was related to a faster multiplication rate. These observations provide evidence for a chain of linked, compensatory mutational events in a defined region of the FMDV capsid. Furthermore, they indicate that the clustering of unique amino acid replacements in viruses from persistent infections may also occur in cytolytic infections in response to changes caused by previous mutations without an involvement of the new mutations in the adaptation to a different environment.

Searching sequence space: two different approaches to dihydrofolate reductase catalysis

There are numerous examples of proteins that catalyze the same reaction while possessing different structures. This review focuses on two dihydrofolate reductases (DHFRs) that have disparate structures and discusses how the catalytic strategies of these two DHFRs are driven by their respective scaffolds. The two proteins are E. coli chromosomal DHFR (Ec DHFR) and a type II R-plasmid-encoded DHFR, typified by R67 DHFR. The former has been described as a very well evolved enzyme with an efficiency of 0.15, while the latter has been suggested to be a model for a "primitive" enzyme that has not yet been optimized by evolution. This comparison underlines what is important to catalysis in these two enzymes and concurrently highlights fundamental issues in enzyme catalysis.

Laboratory-evolved vanillyl-alcohol oxidase produces natural vanillin

The flavoenzyme vanillyl-alcohol oxidase was subjected to random mutagenesis to generate mutants with enhanced reactivity to creosol (2-methoxy-4-methylphenol). The vanillyl-alcohol oxidase-mediated conversion of creosol proceeds via a two-step process in which the initially formed vanillyl alcohol (4-hydroxy-3-methoxybenzyl alcohol) is oxidized to the widely used flavor compound vanillin (4-hydroxy-3-methoxybenzaldehyde). The first step of this reaction is extremely slow due to the formation of a covalent FAD N-5-creosol adduct. After a single round of error-prone PCR, seven mutants were generated with increased reactivity to creosol. The single-point mutants I238T, F454Y, E502G, and T505S showed an up to 40-fold increase in catalytic efficiency (kcat/Km) with creosol compared with the wild-type enzyme. This enhanced reactivity was due to a lower stability of the covalent flavin-substrate adduct, thereby promoting vanillin formation. The catalytic efficiencies of the mutants were also enhanced for other ortho-substituted 4-methylphenols, but not for p-cresol (4-methylphenol). The replaced amino acid residues are not located within a distance of direct interaction with the substrate, and the determined three-dimensional structures of the mutant enzymes are highly similar to that of the wild-type enzyme. These results clearly show the importance of remote residues, not readily predicted by rational design, for the substrate specificity of enzymes.

Correlated motion and the effect of distal mutations in dihydrofolate reductase

Dihydrofolate reductase (DHFR) catalyzes the reduction of dihydrofolate to tetrahydrofolate. The catalytic rate in this system has been found to be significantly affected by mutations far from the site of chemical activity in the enzyme [Rajagopalan, P. T. R, Lutz, S., and Benkovic, S. J. (2002) Biochemistry 41, 12618-12628]. On the basis of extensive computer simulations for wild-type DHFR from Escherichia coli and four mutants (G121S, G121V, M42F, and M42F/G121S), we show that key parameters for catalysis are changed. The parameters we study are relative populations of different conformations sampled and hydrogen bonds. We find that the mutations result in long-range structural perturbations, rationalizing the effects that the mutations have on the kinetics of the enzyme. Such perturbations also provide a rationalization for the reported nonadditivity effect for double mutations. We finally examine the role a structural perturbation will have on the hydride transfer step. On the basis of our new findings, we discuss the role of coupled motions between distant regions in the enzyme, which previously was reported by Radkiewicz and Brooks.

GABA(A) receptor M2-M3 loop secondary structure and changes in accessibility during channel gating

The gamma-aminobutyric acid type A (GABA(A)) receptor M2-M3 loop structure and its role in gating were investigated using the substituted cysteine accessibility method. Residues from alpha(1)Arg-273 to alpha(1)Ile-289 were mutated to cysteine, one at a time. MTSET(+) or MTSES(-) reacted with all mutants from alpha(1)R273C to alpha(1)Y281C, except alpha(1)P277C, in the absence and presence of GABA. The MTSET(+) closed-state reaction rate was >1000 liters/mol-s at alpha(1)N274C, alpha(1)S275C, alpha(1)K278C, and alpha(1)Y281C and was <300 liters/mol-s at alpha(1)R273C, alpha(1)L276C, alpha(1)V279C, alpha(1)A280C, and alpha(1)A284C. These two groups of residues lie on opposite sides of an alpha-helix. The fast reacting group lies on a continuation of the M2 segment channel-lining helix face. This suggests that the M2 segment alpha-helix extends about two helical turns beyond alpha(1)N274 (20'), aligned with the extracellular ring of charge. At alpha(1)S275C, alpha(1)V279C, alpha(1)A280C, and alpha(1)A284C the reaction rate was faster in the presence of GABA. The reagents had no functional effect on the mutants from alpha(1)A282C to alpha(1)I289C, except alpha(1)A284C. Access may be sterically hindered possibly by close interaction with the extracellular domain. We suggest that the M2 segment alpha-helix extends beyond the predicted extracellular end of the M2 segment and that gating induces a conformational change in and/or around the N-terminal half of the M2-M3 loop. Implications for coupling ligand-evoked conformational changes in the extracellular domain to channel gating in the membrane-spanning domain are discussed.

The binding site of aminergic G protein-coupled receptors: the transmembrane segments and second extracellular loop

In the current chapter, we review approaches to the identification of the residues forming the binding sites for agonists, antagonists, and allosteric modulators in the family of aminergic G protein-coupled receptors (GPCRs). We then review the structural bases for ligand binding and pharmacological specificity based on the application of these methods to muscarinic cholinergic, adrenergic, dopaminergic, serotonergic, and histaminergic receptors, using the high resolution rhodopsin structure as a template. Furthermore, we propose a critical role of the second extracellular loop in forming the binding site for small molecular weight aminergic ligands, much as this loop dives down into the binding-site crevice and contacts retinal in rhodopsin.

Identification of a region of strong discrimination in the pore of CFTR

The variety of methods used to identify the structural determinants of anion selectivity in the cystic fibrosis transmembrane conductance regulator Cl(-) channel has made it difficult to assemble the data into a coherent framework that describes the three-dimensional structure of the pore. Here, we compare the relative importance of sites previously studied and identify new sites that contribute strongly to anion selectivity. We studied Cl(-) and substitute anions in oocytes expressing wild-type cystic fibrosis transmembrane conductance regulator or 12-pore-domain mutants to determine relative permeability and relative conductance for 9 monovalent anions and 1 divalent anion. The data indicate that the region of strong discrimination resides between T338 and S341 in transmembrane 6, where mutations affected selectivity between Cl(-) and both large and small anions. Mutations further toward the extracellular end of the pore only strongly affected selectivity between Cl(-) and larger anions. Only mutations at S341 affected selectivity between monovalent and divalent anions. The data are consistent with a narrowing of the pore between the extracellular end and a constriction near the middle of the pore.

Evaluation of site-directed spin labeling for characterizing protein-ligand complexes using simulated restraints

Simulation studies have been performed to evaluate the utility of site-directed spin labeling for determining the structures of protein-ligand complexes, given a known protein structure. Two protein-ligand complexes were used as model systems for these studies: a 1.9-A-resolution x-ray structure of a dihydrofolate reductase mutant complexed with methotrexate, and a 1.5-A-resolution x-ray structure of the V-Src tyrosine kinase SH2 domain complexed with a five-residue phosphopeptide. Nitroxide spin labels were modeled at five dihydrofolate reductase residue positions and at four SH2 domain residue positions. For both systems, after energy minimization, conformational ensembles of the spin-labeled residues were generated by simulated annealing while holding the remainder of the protein-ligand complex fixed. Effective distances, simulating those that could be obtained from (1)H-NMR relaxation measurements, were calculated between ligand protons and the spin labels. These were converted to restraints with several different levels of precision. Restrained simulated annealing calculations were then performed with the aim of reproducing target ligand-binding modes. The effects of incorporating a few supplementary short-range (< or =5.0 A) distance restraints were also examined. For the dihydrofolate reductase-methotrexate complex, the ligand-binding mode was reproduced reasonably well using relatively tight spin-label restraints, but methotrexate was poorly localized using loose spin-label restraints. Short-range and spin-label restraints proved to be complementary. For the SH2 domain-phosphopeptide complex without the short-range restraints, the peptide did not localize to the correct depth in the binding groove; nevertheless, the orientation and internal conformation of the peptide was reproduced moderately well. Use of the spin-label restraints in conjunction with the short-range restraints resulted in relatively well defined structural ensembles. These results indicate that restraints derived from site-directed spin labeling can contribute significantly to defining the orientations and conformations of bound ligands. Accurate ligand localization appears to require either a few supplementary short-range distance restraints, or relatively tight spin-label restraints, with at least one spin label positioned so that some of the restraints draw the ligand into the binding pocket in the latter case.

Syntaxin 1A interacts with multiple exocytic proteins to regulate neurotransmitter release in vivo

Biochemical studies suggest that syntaxin 1A participates in multiple protein-protein interactions in the synaptic terminal, but the in vivo significance of these interactions is poorly understood. We used a targeted mutagenesis approach to eliminate specific syntaxin binding interactions and demonstrate that Drosophila syntaxin 1A plays multiple regulatory roles in neurotransmission in vivo. Syntaxin mutations that eliminate ROP/Munc-18 binding display increased neurotransmitter release, suggesting that ROP inhibits neurosecretion through its interaction with syntaxin. Syntaxin mutations that block Ca2+ channel binding also cause an increase in neurotransmitter release, suggesting that syntaxin normally functions in inhibiting Ca2+ channel opening. Additionally, we identify and characterize a syntaxin Ca2+ effector domain, which may spatially organize the Ca2+ channel, cysteine string protein, and synaptotagmin for effective excitation-secretion coupling in the presynaptic terminal.

Non-pseudogene-derived complex acid beta-glucosidase mutations causing mild type 1 and severe type 2 gaucher disease

Gaucher disease is an autosomal recessive inborn error of glycosphingolipid metabolism caused by the deficient activity of the lysosomal hydrolase, acid beta-glucosidase. Three phenotypically distinct subtypes result from different acid beta-glucosidase mutations encoding enzymes with absent or low activity. A severe neonatal type 2 variant who presented with collodion skin, ichthyosis, and a rapid neurodegenerative course had two novel acid beta-glucosidase alleles: a complex, maternally derived allele, E326K+L444P, and a paternally inherited nonsense mutation, E233X. Because the only other non-pseudogene-derived complex allele, D140H+E326K, also had the E326K lesion and was reported in a mild type 1 patient with a D140H+E326K/K157Q genotype, these complex alleles and their individual mutations were expressed and characterized. Because the E233X mutation expressed no activity and the K157Q allele had approximately 1% normal specific activity based on cross-reacting immunologic material (CRIM SA) in the baculovirus system, the residual activity in both patients was primarily from their complex alleles. In the type 1 patient, the D140H+E326K allele was neuroprotective, encoding an enzyme with a catalytic efficiency similar to that of the N370S enzyme. In contrast, the E326K+L444P allele did not have sufficient activity to protect against the neurologic manifestations and, in combination with the inactive E233X lesion, resulted in the severe neonatal type 2 variant. Thus, characterization of these novel genotypes with non-pseudogene-derived complex mutations provided the pathogenic basis for their diverse phenotypes.

Second-site mutations in cyclic AMP-sensitive revertants of a Ka mutant of S49 mouse lymphoma cells reduce the affinity of regulatory subunit of cyclic AMP-dependent protein kinase for catalytic subunit

Ka mutants of S49 mouse lymphoma cells are generally heterozygous for expression of wild-type and mutant regulatory (R) subunits of type I alpha cyclic AMP-(cAMP)-dependent protein kinase, where the mutant R subunit has a defect in cAMP-binding to one of two intrachain cAMP-binding sites. Several cAMP-sensitive revertants of such a Ka mutant were found previously to harbor second-site mutations in the mutant allele, and we have now identified three such mutations by sequence analysis of PCR-amplified cDNAs. The resulting amino acid changes were Ala98 to Thr, Gly179 to Arg, or Gly224 to Asp. The Ka mutation in these strains (Glu201 to Lys) eliminated cAMP-binding to the more aminoterminal cAMP-binding site (site A). None of the second-site mutations restored this activity in bacterially expressed recombinant R subunit. On the other hand, all three second-site mutations reduced the apparent affinity of the mutant R subunit for catalytic (C) subunit with the effects of the substitutions at Ala98 and Gly179 substantially greater than the effect of the substitution at Gly224. Patterns of phosphorylation and turnover of wild-type and mutant R subunits in intact revertant cells were consistent with reduced association of the doubly mutant subunits with C subunit, but the free mutant subunits apparently were more stable than free wild-type subunits. Differences in metabolic turnover of mutant and wild-type subunits did not correlate with the sensitivities of the isolated proteins to proteolytic cleavage.

Enzymes: chemistry and biochemistry

Basic principles underlying enzyme action are considered. Catalytic antibodies (abzymes), catalytic RNA (ribozymes), and non-biological counterparts of enzyme-catalyzed reactions are mentioned. Enzyme evolution is considered in terms of divergence, convergence, and lateral gene transfer.