A plasmid-encoded dihydrofolate reductase from trimethoprim-resistant bacteria has a novel D2-symmetric active site

From a binding module to essential catalytic activity: how nature stumbled on a good thing

Enzymes are complex macromolecules capable of catalyzing a wide variety of chemical reactions with high efficiency. Nonetheless, biological catalysis can be rudimentary. Here, we describe an enzyme that is built from a simple protein fold. This short protein sequence - almost a peptide - belongs to the ancient SH3 family of binding modules. Surprisingly, this binding module catalyzes the specific reduction of dihydrofolate using NADPH as a reducing cofactor, making this a dihydrofolate reductase. Too small to provide all the required binding and catalytic machinery on its own, it homotetramerizes, thus creating a large, central active site environment. Remarkably, none of the active site residues is essential to the catalytic function. Instead, backbone interactions juxtapose the reducing cofactor proximal to the target imine of the folate substrate, and a specific motion of the substrate promotes formation of the transition state. In this feature article, we describe the features that make this small protein a functional enzyme capable of catalyzing a metabolically essential reaction, highlighting the characteristics that make it a model for the evolution of primitive enzymes from binding modules.

Epistasis between promoter activity and coding mutations shapes gene evolvability

The evolution of protein-coding genes proceeds as mutations act on two main dimensions: regulation of transcription level and the coding sequence. The extent and impact of the connection between these two dimensions are largely unknown because they have generally been studied independently. By measuring the fitness effects of all possible mutations on a protein complex at various levels of promoter activity, we show that promoter activity at the optimal level for the wild-type protein masks the effects of both deleterious and beneficial coding mutations. Mutations that are deleterious at low activity but masked at optimal activity are slightly destabilizing for individual subunits and binding interfaces. Coding mutations that increase protein abundance are beneficial at low expression but could potentially incur a cost at high promoter activity. We thereby demonstrate that promoter activity in interaction with protein properties can dictate which coding mutations are beneficial, neutral, or deleterious.

Discovery of Highly Trimethoprim-Resistant DfrB Dihydrofolate Reductases in Diverse Environmental Settings Suggests an Evolutionary Advantage Unrelated to Antibiotic Resistance

Type B dihydrofolate reductases (DfrB) are intrinsically highly resistant to the widely used antibiotic trimethoprim, posing a threat to global public health. The ten known DfrB family members have been strongly associated with genetic material related to the application of antibiotics. Several dfrB genes were associated with multidrug resistance contexts and mobile genetic elements, integrated both in chromosomes and plasmids. However, little is known regarding their presence in other environments. Here, we investigated the presence of dfrB beyond the traditional areas of enquiry by conducting metagenomic database searches from environmental settings where antibiotics are not prevalent. Thirty putative DfrB homologues that share 62 to 95% identity with characterized DfrB were identified. Expression of ten representative homologues verified trimethoprim resistance in all and dihydrofolate reductase activity in most. Contrary to samples associated with the use of antibiotics, the newly identified dfrB were rarely associated with mobile genetic elements or antibiotic resistance genes. Instead, association with metabolic enzymes was observed, suggesting an evolutionary advantage unrelated to antibiotic resistance. Our results are consistent with the hypothesis that multiple dfrB exist in diverse environments from which dfrB were mobilized into the clinically relevant resistome. Our observations reinforce the need to closely monitor their progression.

Computational Development of Inhibitors of Plasmid-Borne Bacterial Dihydrofolate Reductase

Resistance to trimethoprim and other antibiotics targeting dihydrofolate reductase may arise in bacteria harboring an atypical, plasmid-encoded, homotetrameric dihydrofolate reductase, called R67 DHFR. Although developing inhibitors to this enzyme may be expected to be promising drugs to fight trimethoprim-resistant strains, there is a paucity of reports describing the development of such molecules. In this manuscript, we describe the design of promising lead compounds to target R67 DHFR. Density-functional calculations were first used to identify the modifications of the pterin core that yielded derivatives likely to bind the enzyme and not susceptible to being acted upon by it. These unreactive molecules were then docked to the active site, and the stability of the docking poses of the best candidates was analyzed through triplicate molecular dynamics simulations, and compared to the binding stability of the enzyme–substrate complex. Molecule 32 ([6-(methoxymethyl)-4-oxo-3,7-dihydro-4H-pyrano[2,3-d]pyrimidin-2-yl]methyl-guanidinium) was shown by this methodology to afford extremely stable binding towards R67 DHFR and to prevent simultaneous binding to the substrate. Additional docking and molecular dynamics simulations further showed that this candidate also binds strongly to the canonical prokaryotic dihydrofolate reductase and to human DHFR, and is therefore likely to be useful to the development of chemotherapeutic agents and of dual-acting antibiotics that target the two types of bacterial dihydrofolate reductase.

Identification and validation of potent inhibitor of Escherichia coli DHFR from MMV pathogen box

The present study is conducted to find the solution of rising antimicrobial resistance (AMR) in Escherichia coli which is a pathogen responsible for fatal systemic infections in human and animals. The enzyme dihydrofolate reductase (DHFR) is found in all organisms. In this study DHFR of E. coli (ec-DHFR) and human DHFR (h-DHFR) is targeted by novel chemical entities (NCE) from the Pathogen box of Medicines for Malaria Venture, Switzerland (MMV) using molecular modelling. The in-silico studies were further validated by in-vitro assays. The virtual screening of 400 MMV compounds was conducted using PyRx standard tool followed by manual docking of selected compounds by Autodock vina and Ligplot program. The in-silico studies showed good binding energy and strong hydrogen bond in docking of MMV675968 with ec-DHFR and no hydrogen bond with h-DHFR. This was further validated by the Molecular dynamic studies that revealed high binding free energy in ec-DHFR and in-vitro assays that produced good synergy in combination study of MMV675968 with last line (meropenem) and last resort (colistin) antibiotics. The extensive MD simulation and energetic analysis thus concludes that MMV675968 targets ec-DHFR. The combination studies were conducted with MMV675968 and FDA approved drugs against a panel of multidrug resistant Escherichia coli isolates. The synergistic results obtained in combination studies concluded that in-vitro data is consistent with in-silico data and that MMV675968 is a potential lead for future process of antimicrobial drug development against the multidrug resistance E. coli.Communicated by Ramaswamy H. Sarma.

The Bacterial Genomic Context of Highly Trimethoprim-Resistant DfrB Dihydrofolate Reductases Highlights an Emerging Threat to Public Health

Type B dihydrofolate reductase (dfrb) genes were identified following the introduction of trimethoprim in the 1960s. Although they intrinsically confer resistance to trimethoprim (TMP) that is orders of magnitude greater than through other mechanisms, the distribution and prevalence of these short (237 bp) genes is unknown. Indeed, this knowledge has been hampered by systematic biases in search methodologies. Here, we investigate the genomic context of dfrbs to gain information on their current distribution in bacterial genomes. Upon searching publicly available databases, we identified 61 sequences containing dfrbs within an analyzable genomic context. The majority (70%) of those sequences also harbor virulence genes and 97% of the dfrbs are found near a mobile genetic element, representing a potential risk for antibiotic resistance genes. We further identified and confirmed the TMP-resistant phenotype of two new members of the family, dfrb10 and dfrb11. Dfrbs are found both in Betaproteobacteria and Gammaproteobacteria, a majority (59%) being in Pseudomonas aeruginosa. Previously labelled as strictly plasmid-borne, we found 69% of dfrbs in the chromosome of pathogenic bacteria. Our results demonstrate that the intrinsically TMP-resistant dfrbs are a potential emerging threat to public health and justify closer surveillance of these genes.

Differentiation of the binding of two ligands to a tetrameric protein with a single symmetric active site by 19 F NMR

R67 dihydrofolate reductase (R67 DHFR) is a plasmid‐encoded enzyme that confers resistance to the antibacterial drug trimethoprim. R67 DHFR is a tetramer with a single active site that is unusual as both cofactor and substrate are recognized by symmetry‐related residues. Such promiscuity has limited our previous efforts to differentiate ligand binding by NMR. To address this problem, we incorporated fluorine at positions 4, 5, 6, or 7 of the indole rings of tryptophans 38 and 45 and characterized the spectra to determine which probe was optimal for studying ligand binding. Two resonances were observed for all apo proteins. Unexpectedly, the W45 resonance appeared broad, and truncation of the disordered N‐termini resulted in the appearance of one sharp W45 resonance. These results are consistent with interaction of the N‐terminus with W45. Binding of the cofactor broadened W38 for all fluorine probes, whereas substrate, dihydrofolate, binding resulted in the appearance of three new resonances for 4‐ and 5‐fluoroindole labeled protein and severe line broadening for 6‐ and 7‐fluoroindole R67 DHFR. W45 became slightly broader upon ligand binding. With only two peaks in the ¹⁹F NMR spectra, our data were able to differentiate cofactor and substrate binding to the single, symmetric active site of R67 DHFR and yield binding affinities.

Crowders Steal Dihydrofolate Reductase Ligands through Quinary Interactions

Dihydrofolate reductase (DHFR) reduces dihydrofolate (DHF) to tetrahydrofolate using NADPH as a cofactor. Due to its role in one carbon metabolism, chromosomal DHFR is the target of the antibacterial drug, trimethoprim. Resistance to trimethoprim has resulted in a type II DHFR that is not structurally related to the chromosomal enzyme target. Because of its metabolic significance, understanding DHFR kinetics and ligand binding behavior in more cell-like conditions, where the total macromolecule concentration can be as great as 300 mg/mL, is important. The progress-curve kinetics and ligand binding properties of the drug target (chromosomal E. coli DHFR) and the drug resistant (R67 DHFR) enzymes were studied in the presence of macromolecular cosolutes. There were varied effects on NADPH oxidation and binding to the two DHFRs, with some cosolutes increasing affinity and others weakening binding. However, DHF binding and reduction in both DHFRs decreased in the presence of all cosolutes. The decreased binding of ligands is mostly attributed to weak associations with the macromolecules, as opposed to crowder effects on the DHFRs. Computer simulations found weak, transient interactions for both ligands with several proteins. The net charge of protein cosolutes correlated with effects on NADP⁺ binding, with near neutral and positively charged proteins having more detrimental effects on binding. For DHF binding, effects correlated more with the size of binding pockets on the protein crowders. These nonspecific interactions between DHFR ligands and proteins predict that the in vivo efficiency of DHFRs may be much lower than expected from their in vitro rates.

In Vivo Titration of Folate Pathway Enzymes

How enzymes behave in cells is likely different from how they behave in the test tube. Previous in vitro studies find that osmolytes interact weakly with folate. Removal of the osmolyte from the solvation shell of folate is more difficult than removal of water, which weakens binding of folate to its enzyme partners. To examine if this phenomenon occurs in vivo, osmotic stress titrations were performed with Escherichia coli Two strategies were employed: resistance to an antibacterial drug and complementation of a knockout strain by the appropriate gene cloned into a plasmid that allows tight control of expression levels as well as labeling by a degradation tag. The abilities of the knockout and complemented strains to grow under osmotic stress were compared. Typically, the knockout strain could grow to high osmolalities on supplemented medium, while the complemented strain stopped growing at lower osmolalities on minimal medium. This pattern was observed for an R67 dihydrofolate reductase clone rescuing a ΔfolA strain, for a methylenetetrahydrofolate reductase clone rescuing a ΔmetF strain, and for a serine hydroxymethyltransferase clone rescuing a ΔglyA strain. Additionally, an R67 dihydrofolate reductase clone allowed E. coli DH5α to grow in the presence of trimethoprim until an osmolality of ∼0.81 is reached, while cells in a control titration lacking antibiotic could grow to 1.90 osmol.IMPORTANCEE. coli can survive in drought and flooding conditions and can tolerate large changes in osmolality. However, the cell processes that limit bacterial growth under high osmotic stress conditions are not known. In this study, the dose of four different enzymes in E. coli was decreased by using deletion strains complemented by the gene carried in a tunable plasmid. Under conditions of limiting enzyme concentration (lower than that achieved by chromosomal gene expression), cell growth can be blocked by osmotic stress conditions that are normally tolerated. These observations indicate that E. coli has evolved to deal with variations in its osmotic environment and that normal protein levels are sufficient to buffer the cell from environmental changes. Additional factors involved in the osmotic pressure response may include altered protein concentration/activity levels, weak solute interactions with ligands which can make it more difficult for proteins to bind their substrates/inhibitors/cofactors in vivo, and/or viscosity effects.

Small Angle Neutron Scattering Studies of R67 Dihydrofolate Reductase, a Tetrameric Protein with Intrinsically Disordered N-Termini

R67 dihydrofolate reductase (DHFR) is a homotetramer with a single active site pore and no sequence or structural homology with chromosomal DHFRs. The R67 enzyme provides resistance to trimethoprim, an active site-directed inhibitor of Escherichia coli DHFR. Sixteen to twenty N-terminal amino acids are intrinsically disordered in the R67 dimer crystal structure. Chymotrypsin cleavage of 16 N-terminal residues results in an active enzyme with a decreased stability. The space sampled by the disordered N-termini of R67 DHFR was investigated using small angle neutron scattering. From a combined analysis using molecular dynamics and the program SASSIE ( http://www.smallangles.net/sassie/SASSIE_HOME.html ), the apoenzyme displays a radius of gyration (Rg) of 21.46 ± 0.50 Å. Addition of glycine betaine, an osmolyte, does not result in folding of the termini as the Rg increases slightly to 22.78 ± 0.87 Å. SASSIE fits of the latter SANS data indicate that the disordered N-termini sample larger regions of space and remain disordered, suggesting they might function as entropic bristles. Pressure perturbation calorimetry also indicated that the volume of R67 DHFR increases upon addition of 10% betaine and decreased at 20% betaine because of the dehydration of the protein. Studies of the hydration of full-length R67 DHFR in the presence of the osmolytes betaine and dimethyl sulfoxide find around 1250 water molecules hydrating the protein. Similar studies with truncated R67 DHFR yield around 400 water molecules hydrating the protein in the presence of betaine. The difference of ∼900 waters indicates the N-termini are well-hydrated.

X-ray crystal structure of Escherichia coli HspQ, a protein involved in the retardation of replication initiation

The heat shock protein HspQ (YccV) of Escherichia coli has been proposed to participate in the retardation of replication initiation in cells with the dnaA508 allele. In this study, we have determined the 2.5-Å-resolution X-ray structure of the trimer of HspQ, which is also the first structure of a member of the YccV superfamily. The acidic character of the HspQ trimer suggests an interaction surface with basic proteins. From these results, we discuss the cellular function of HspQ, including its relationship with the DnaA508 protein.

Aspects of Weak Interactions between Folate and Glycine Betaine

Folate, or vitamin B9, is an important compound in one-carbon metabolism. Previous studies have found weaker binding of dihydrofolate to dihydrofolate reductase in the presence of osmolytes. In other words, osmolytes are more difficult to remove from the dihydrofolate solvation shell than water; this shifts the equilibrium toward the free ligand and protein species. This study uses vapor-pressure osmometry to explore the interaction of folate with the model osmolyte, glycine betaine. This method yields a preferential interaction potential (μ23/RT value). This value is concentration-dependent as folate dimerizes. The μ23/RT value also tracks the deprotonation of folate's N3-O4 keto-enol group, yielding a pKa of 8.1. To determine which folate atoms interact most strongly with betaine, the interaction of heterocyclic aromatic compounds (as well as other small molecules) with betaine was monitored. Using an accessible surface area approach coupled with osmometry measurements, deconvolution of the μ23/RT values into α values for atom types was achieved. This allows prediction of μ23/RT values for larger molecules such as folate. Molecular dynamics simulations of folate show a variety of structures from extended to L-shaped. These conformers possess μ23/RT values from -0.18 to 0.09 m-1, where a negative value indicates a preference for solvation by betaine and a positive value indicates a preference for water. This range of values is consistent with values observed in osmometry and solubility experiments. As the average predicted folate μ23/RT value is near zero, this indicates folate interacts almost equally well with betaine and water. Specifically, the glutamate tail prefers to interact with water, while the aromatic rings prefer betaine. In general, the more protonated species in our small molecule survey interact better with betaine as they provide a source of hydrogens (betaine is not a hydrogen bond donor). Upon deprotonation of the small molecule, the preference swings toward water interaction because of its hydrogen bond donating capacities.

Tales of Dihydrofolate Binding to R67 Dihydrofolate Reductase

Homotetrameric R67 dihydrofolate reductase possesses 222 symmetry and a single active site pore. This situation results in a promiscuous binding site that accommodates either the substrate, dihydrofolate (DHF), or the cofactor, NADPH. NADPH interacts more directly with the protein as it is larger than the substrate. In contrast, the p-aminobenzoyl-glutamate tail of DHF, as monitored by nuclear magnetic resonance and crystallography, is disordered when bound. To explore whether smaller active site volumes (which should decrease the level of tail disorder by confinement effects) alter steady state rates, asymmetric mutations that decreased the half-pore volume by ∼35% were constructed. Only minor effects on k(cat) were observed. To continue exploring the role of tail disorder in catalysis, 1-ethyl-3-[3-(dimethylamino)propyl]carbodiimide-mediated cross-linking between R67 DHFR and folate was performed. A two-folate, one-tetramer complex results in the loss of enzyme activity where two symmetry-related K32 residues in the protein are cross-linked to the carboxylates of two bound folates. The tethered folate could be reduced, although with a ≤30-fold decreased rate, suggesting decreased dynamics and/or suboptimal positioning of the cross-linked folate for catalysis. Computer simulations that restrain the dihydrofolate tail near K32 indicate that cross-linking still allows movement of the p-aminobenzoyl ring, which allows the reaction to occur. Finally, a bis-ethylene-diamine-α,γ-amide folate adduct was synthesized; both negatively charged carboxylates in the glutamate tail were replaced with positively charged amines. The K(i) for this adduct was ∼9-fold higher than for folate. These various results indicate a balance between folate tail disorder, which helps the enzyme bind substrate while dynamics facilitates catalysis.

Identification and genomic analysis of temperate Pseudomonas bacteriophage PstS-1 from the Japan trench at a depth of 7000 m

Viruses play important roles in aquatic ecosystems, but deep-sea bacteriophages remain largely unexplored. A temperate bacteriophage (termed vB_PstS-1) was identified from the psychrotolerant gammaproteobacterium Pseudomonas stutzeri 1-1-1b, which was isolated from hadopelagic water (depth of 7000 m) of the Japan Trench in the Northwest Pacific Ocean. The genome size of PstS-1 was 48,666 bp; its genome displayed a 59.8% G + C content and a total of 79 coding sequences were identified in its genome. The PstS-1 phage belongs to the family Siphoviridae, but its genomic sequence and organization are distinct from those of any other well-known Siphoviridae phage. The mosaic genomic structure of PstS-1 suggests the occurrence of genetic exchange between distinct temperate phages in deep-sea Pseudomonas populations. The PstS-1 genome also harbors three distinct sequence regions corresponding to spacers within a single clustered regularly interspaced short palindromic repeat (CRISPR) locus in the rhizosphere-associated diazotrophic P. stutzeri A1501 genome. The extension of these spacers to the soil environment and the presence of many homologs of both the hadal deep-sea phage PstS-1 and terrestrial Pseudomonas phages suggest the early co-evolution of temperate phages and their host genus Pseudomonas prior to the divergence of their habitational and physiological adaptation.

Asymmetric mutations in the tetrameric R67 dihydrofolate reductase reveal high tolerance to active-site substitutions

Type II R67 dihydrofolate reductase (DHFR) is a bacterial plasmid-encoded enzyme that is intrinsically resistant to the widely-administered antibiotic trimethoprim. R67 DHFR is genetically and structurally unrelated to E. coli chromosomal DHFR and has an unusual architecture, in that four identical protomers form a single symmetrical active site tunnel that allows only one substrate binding/catalytic event at any given time. As a result, substitution of an active-site residue has as many as four distinct consequences on catalysis, constituting an atypical model of enzyme evolution. Although we previously demonstrated that no single residue of the native active site is indispensable for function, library selection here revealed a strong bias toward maintenance of two native protomers per mutated tetramer. A variety of such "half-native" tetramers were shown to procure native-like catalytic activity, with similar KM values but kcat values 5- to 33-fold lower, illustrating a high tolerance for active-site substitutions. The selected variants showed a reduced thermal stability (Tm ∼12°C lower), which appears to result from looser association of the protomers, but generally showed a marked increase in resilience to heat denaturation, recovering activity to a significantly greater extent than the variant with no active-site substitutions. Our results suggest that the presence of two native protomers in the R67 DHFR tetramer is sufficient to provide native-like catalytic rate and thus ensure cellular proliferation.

Variations in metabolic pathways create challenges for automated metabolic reconstructions: Examples from the tetrahydrofolate synthesis pathway

The availability of thousands of sequenced genomes has revealed the diversity of biochemical solutions to similar chemical problems. Even for molecules at the heart of metabolism, such as cofactors, the pathway enzymes first discovered in model organisms like Escherichia coli or Saccharomyces cerevisiae are often not universally conserved. Tetrahydrofolate (THF) (or its close relative tetrahydromethanopterin) is a universal and essential C1-carrier that most microbes and plants synthesize de novo. The THF biosynthesis pathway and enzymes are, however, not universal and alternate solutions are found for most steps, making this pathway a challenge to annotate automatically in many genomes. Comparing THF pathway reconstructions and functional annotations of a chosen set of folate synthesis genes in specific prokaryotes revealed the strengths and weaknesses of different microbial annotation platforms. This analysis revealed that most current platforms fail in metabolic reconstruction of variant pathways. However, all the pieces are in place to quickly correct these deficiencies if the different databases were built on each other's strengths.

Investigation of osmolyte effects on FolM: comparison with other dihydrofolate reductases

A weak association between osmolytes and dihydrofolate (DHF) decreases the affinity of the substrate for the Escherichia coli chromosomal and R67 plasmid dihydrofolate reductase (DHFR) enzymes. To test whether the osmolyte-DHF association also interferes with binding of DHF to FolM, an E. coli enzyme that possesses weak DHFR activity, ligand binding was monitored in the presence of osmolytes. The affinity of FolM for DHF, measured by kcat/Km(DHF), was decreased by the addition of an osmolyte. Additionally, binding of the antifolate drug, methotrexate, to FolM was weakened by the addition of an osmolyte. The changes in ligand binding with water activity were unique for each osmolyte, indicating preferential interaction between the osmolyte and folate and its derivatives; however, additional evidence provided support for further interactions between FolM and osmolytes. Binding of the reduced nicotinamide adenine dinucleotide phosphate (NADPH) cofactor to FolM was monitored by isothermal titration calorimetry as a control for protein-osmolyte association. In the presence of betaine (proposed to be the osmolyte most excluded from protein surfaces), the NADPH Kd decreased, consistent with dehydration effects. However, other osmolytes did not tighten binding to the cofactor. Rather, dimethyl sulfoxide (DMSO) had no effect on the NADPH Kd, while ethylene glycol and polyethylene glycol 400 weakened cofactor binding. Differential scanning calorimetry of FolM in the presence of osmolytes showed that both DMSO and ethylene glycol decreased the stability of FolM, while betaine increased the stability of the protein. These results suggest that some osmolytes can destabilize FolM by preferentially interacting with the protein. Further, these weak attractions can impede ligand binding. These various contributions have to be considered when interpreting osmotic pressure results.

Further studies on the role of water in R67 dihydrofolate reductase

Previous osmotic pressure studies of two nonhomologous dihydrofolate reductase (DHFR) enzymes found tighter binding of the nicotinamide adenine dinucleotide phosphate cofactor upon addition of neutral osmolytes. This result is consistent with water release accompanying binding. In contrast, osmotic stress studies found weaker binding of the dihydrofolate (DHF) substrate for both type I and type II DHFRs in the presence of osmolytes; this observation can be explained if dihydrofolate interacts with osmolytes and shifts the equilibrium from the enzyme-bound state toward the unbound substrate. Nuclear magnetic resonance experiments support this hypothesis, finding that osmolytes interact with dihydrofolate. To consider binding without added osmolytes, a high-pressure approach was used. In this study, the type II enzyme, R67 DHFR, was subjected to high hydrostatic pressure (HHP). Both enzyme activity and fluorescence measurements find the protein tolerates pressures up to 200 MPa. Binding of the cofactor to R67 DHFR weakens with increasing pressure, and a positive association volume of 11.4 ± 0.5 cm(3)/mol was measured. Additionally, an activation volume of 3.3 ± 0.5 cm(3)/mol describing k(cat)/K(m(DHF)) was determined from progress curve analysis. Results from these HHP experiments suggest water release accompanies binding of both the cofactor and DHF to R67 DHFR. In an additional set of experiments, isothermal titration calorimetry studies in H2O and D2O find that water reorganization dominates the enthalpy associated with binding of DHF to R67 DHFR·NADP(+), while no obvious effects occur for cofactor binding. The combined results indicate that water plays an active role in ligand binding to R67 DHFR.

Insights into the evolution of gene organization and multidrug resistance from Klebsiella pneumoniae plasmid pKF3-140

Plasmid-mediated transfer of drug-resistance genes among various bacterial species is considered one of the most important mechanisms for the spread of multidrug resistance. To gain insights into the evolution of gene organization and antimicrobial resistance in clinical bacterial samples, a complete plasmid genome of Klebsiella pneumoniae pKF3-140 is determined, which has a circular chromosome of 147,416bp in length. Among the 203 predicted genes, 142 have function assignment and about 50 appear to be involved in plasmid replication, maintenance, conjugative transfer, iron acquisition and transport, and drug resistance. Extensive comparative genomic analyses revealed that pKF3-140 exhibits a rather low sequence similarity and structural conservation with other reported K. pneumoniae plasmids. In contrast, the overall organization of pKF3-140 is highly similar to Escherichia coli plasmids p1ESCUM and pUTI89, which indicates the possibility that K. pneumoniae pKF3-140 may have a potential origin in E. coli. Meanwhile, interestingly, several drug resistant genes show high similarity to the plasmid pU302L in Salmonella enterica serovar Typhimurium U302 strain G8430 and the plasmid pK245 in K. pneumoniae. This mosaic pattern of sequence similarities suggests that pKF3-140 might have arisen from E. coli and acquired the resistance genes from a variety of enteric bacteria and underscores the importance of a further understanding of horizontal gene transfer among enteric bacteria.

Weak interactions between folate and osmolytes in solution

Previous osmotic stress studies on the role of solvent in two structurally unrelated dihydrofolate reductases (DHFRs) found weaker binding of dihydrofolate (DHF) to either enzyme in the presence of osmolytes. To explain these unusual results, weak interactions between DHF and osmolytes were proposed, with a competition between osmolyte and DHFR for DHF. High osmolyte concentrations will inhibit binding of the cognate pair. To evaluate this hypothesis, we devised a small molecule approach. Dimerization of folate, monitored by nuclear magnetic resonance, was weakened 2-3-fold upon addition of betaine or dimethyl sulfoxide (DMSO), supporting preferential interaction of either osmolyte with the monomer (as it possesses a larger surface area). Nuclear Overhauser effect (NOE) spectroscopy experiments found a positive NOE for the interaction of the C3'/C5' benzoyl ring protons with the C9 proton in buffer; however, a negative NOE was observed upon addition of betaine or DMSO. This change indicated a decreased tumbling rate, consistent with osmolyte interaction. Osmotic stress experiments also showed that betaine, DMSO, and sucrose preferentially interact with folate. Further, studies with the folate fragments, p-aminobenzoic acid and pterin 6-carboxylate, revealed interactions for both model compounds with betaine and sucrose. In contrast, DMSO was strongly excluded from the pterin ring but preferentially interacted with the p-aminobenzoyl moiety. These interactions are likely to be important in vivo because of the crowded conditions of the cell where weak contacts can more readily compete with specific binding interactions.

Novel crystallization conditions for tandem variant R67 DHFR yield a wild-type crystal structure

Trimethoprim is an antibiotic that targets bacterial dihydrofolate reductase (DHFR). A plasmid-encoded DHFR known as R67 DHFR provides resistance to trimethoprim in bacteria. To better understand the mechanism of this homotetrameric enzyme, a tandem dimer construct was created that linked two monomeric R67 DHFR subunits together and mutated the sequence of residues 66-69 of the first subunit from VQIY to INSF. Using a modified crystallization protocol for this enzyme that included in situ proteolysis using chymotrypsin, the tandem dimer was crystallized and the structure was solved at 1.4 Å resolution. Surprisingly, only wild-type protomers were incorporated into the crystal. Further experiments demonstrated that the variant protomer was selectively degraded by chymotrypsin, although no canonical chymotrypsin cleavage site had been introduced by these mutations.

Thermodynamics and solvent effects on substrate and cofactor binding in Escherichia coli chromosomal dihydrofolate reductase

Chromosomal dihydrofolate reductase from Escherichia coli catalyzes the reduction of dihydrofolate to tetrahydrofolate using NADPH as a cofactor. The thermodynamics of ligand binding were examined using an isothermal titration calorimetry approach. Using buffers with different heats of ionization, zero to a small, fractional proton release was observed for dihydrofolate binding, while a proton was released upon NADP(+) binding. The role of water in binding was additionally monitored using a number of different osmolytes. Binding of NADP(+) is accompanied by the net release of ∼5-24 water molecules, with a dependence on the identity of the osmolyte. In contrast, binding of dihydrofolate is weakened in the presence of osmolytes, consistent with "water uptake". Different effects are observed depending on the identity of the osmolyte. The net uptake of water upon dihydrofolate binding was previously observed in the nonhomologous R67-encoded dihydrofolate reductase (dfrB or type II enzyme) [Chopra, S., et al. (2008) J. Biol. Chem. 283, 4690-4698]. As R67 dihydrofolate reductase possesses a nonhomologous sequence and forms a tetrameric structure with a single active site pore, the observation of weaker DHF binding in the presence of osmolytes in both enzymes implicates cosolvent effects on free dihydrofolate. Consistent with this analysis, stopped flow experiments find betaine mostly affects DHF binding via changes in k(on), while betaine mostly affects NADPH binding via changes in k(off). Finally, nonadditive enthalpy terms when binary and ternary cofactor binding events are compared suggest the presence of long-lived conformational transitions that are not included in a simple thermodynamic cycle.

Radical redesign of a tandem array of four R67 dihydrofolate reductase genes yields a functional, folded protein possessing 45 substitutions

R67 dihydrofolate reductase (DHFR) is a plasmid-encoded, type II enzyme. Four monomers (78 amino acids long) assemble into a homotetramer possessing 222 symmetry. In previous studies, a tandem array of four R67 DHFR gene copies was fused in frame to generate a functional monomer named Quad1. This protein possessed the essential tertiary structure of the R67 "parent". To facilitate mutagenesis reactions, restriction enzyme sites were introduced in the tandem gene array. S59A and H362L mutations were also added to minimize possible folding topologies; this protein product, named Quad3, possesses 10 substitutions and is functional. Since R67 DHFR possesses a stable scaffold, a large jump in sequence space was performed by the further addition of 45 amino acid substitutions. The mutational design utilized alternate sequences from other type II DHFRs. In addition, most of the mutations were positioned on the surface of the protein as well as in the disordered N-terminal sequence, which serves as the linker between the fused domains. The resulting Quad4 protein is quite functional; however, it is less stable than Quad1, suffering a DeltaDeltaG loss of 5 kcal/mol at pH 5. One unexpected result was formation of Quad4 dimers and higher order oligomers at pH 8. R67 DHFR, and its derivative Quad proteins, possesses a robust scaffold, capable of withstanding introduction of >or=55 substitutions.

A bacterial antirepressor with SH3 domain topology mimics operator DNA in sequestering the repressor DNA recognition helix

Direct targeting of critical DNA-binding elements of a repressor by its cognate antirepressor is an effective means to sequester the repressor and remove a transcription initiation block. Structural descriptions for this, though often proposed for bacterial and phage repressor–antirepressor systems, are unavailable. Here, we describe the structural and functional basis of how the Myxococcus xanthus CarS antirepressor recognizes and neutralizes its cognate repressors to turn on a photo-inducible promoter. CarA and CarH repress the carB operon in the dark. CarS, produced in the light, physically interacts with the MerR-type winged-helix DNA-binding domain of these repressors leading to activation of carB. The NMR structure of CarS1, a functional CarS variant, reveals a five-stranded, antiparallel β-sheet fold resembling SH3 domains, protein–protein interaction modules prevalent in eukaryotes but rare in prokaryotes. NMR studies and analysis of site-directed mutants in vivo and in vitro unveil a solvent-exposed hydrophobic pocket lined by acidic residues in CarS, where the CarA DNA recognition helix docks with high affinity in an atypical ligand-recognition mode for SH3 domains. Our findings uncover an unprecedented use of the SH3 domain-like fold for protein–protein recognition whereby an antirepressor mimics operator DNA in sequestering the repressor DNA recognition helix to activate transcription.

Crystal structure of a type II dihydrofolate reductase catalytic ternary complex

Type II dihydrofolate reductase (DHFR) is a plasmid-encoded enzyme that confers resistance to bacterial DHFR-targeted antifolate drugs. It forms a symmetric homotetramer with a central pore which functions as the active site. Its unusual structure, which results in a promiscuous binding surface that accommodates either the dihydrofolate (DHF) substrate or the NADPH cofactor, has constituted a significant limitation to efforts to understand its substrate specificity and reaction mechanism. We describe here the first structure of a ternary R67 DHFR.DHF.NADP+ catalytic complex, resolved to 1.26 A. This structure provides the first clear picture of how this enzyme, which lacks the active site carboxyl residue that is ubiquitous in Type I DHFRs, is able to function. In the catalytic complex, the polar backbone atoms of two symmetry-related I68 residues provide recognition motifs that interact with the carboxamide on the nicotinamide ring, and the N3-O4 amide function on the pteridine ring. This set of interactions orients the aromatic rings of substrate and cofactor in a relative endo geometry in which the reactive centers are held in close proximity. Additionally, a central, hydrogen-bonded network consisting of two pairs of Y69-Q67-Q67'-Y69' residues provides an unusually tight interface, which appears to serve as a "molecular clamp" holding the substrates in place in an orientation conducive to hydride transfer. In addition to providing the first clear insight regarding how this extremely unusual enzyme is able to function, the structure of the ternary complex provides general insights into how a mutationally challenged enzyme, i.e., an enzyme whose evolution is restricted to four-residues-at-a-time active site mutations, overcomes this fundamental limitation.

A balancing act between net uptake of water during dihydrofolate binding and net release of water upon NADPH binding in R67 dihydrofolate reductase

R67 dihydrofolate reductase (DHFR) catalyzes the reduction of dihydrofolate (DHF) to tetrahydrofolate using NADPH as a cofactor. This enzyme is a homotetramer possessing 222 symmetry, and a single active site pore traverses the length of the protein. A promiscuous binding surface can accommodate either DHF or NADPH, thus two nonproductive complexes can form (2NADPH or 2DHF) as well as a productive complex (NADPH.DHF). The role of water in binding was monitored using a number of different osmolytes. From isothermal titration calorimetry (ITC) studies, binding of NADPH is accompanied by the net release of 38 water molecules. In contrast, from both steady state kinetics and ITC studies, binding of DHF is accompanied by the net uptake of water. Although different osmolytes have similar effects on NADPH binding, variable results are observed when DHF binding is probed. Sensitivity to water activity can also be probed by an in vivo selection using the antibacterial drug, trimethoprim, where the water content of the media is decreased by increasing concentrations of sorbitol. The ability of wild type and mutant clones of R67 DHFR to allow host Escherichia coli to grow in the presence of trimethoprim plus added sorbitol parallels the catalytic efficiency of the DHFR clones, indicating water content strongly correlates with the in vivo function of R67 DHFR.

Structure of the Q67H mutant of R67 dihydrofolate reductase-NADP+ complex reveals a novel cofactor binding mode

Plasmid-encoded bacterial R67 dihydrofolate reductase (DHFR) is a NADPH-dependent enzyme unrelated to chromosomal DHFR in amino acid sequence and structure. R67 DHFR is insensitive to the bacterial drug trimethoprim in contrast to chromosomal DHFR. The crystal structure of Q67H mutant of R67 DHFR bound to NADP(+) has been determined at 1.15 angstroms resolution. The cofactor assumes an extended conformation with the nicotinamide ring bound near the center of the active site pore, the ribose and pyrophosphate group (PP(i)) extending toward the outer pore. The ribonicotinamide exhibits anti conformation as in chromosomal DHFR complexes. The relative orientation between the PP(i) and the nicotinamide ribose differs from that observed in chromosomal DHFR-NADP(+) complexes. The coenzyme displays symmetrical binding mode with several water-mediated hydrogen bonds with the protein besides ionic, stacking, and van der Waals interactions. The structure provides a molecular basis for the observed stoichiometry and cooperativity in ligand binding. The ternary model based on the present structure and the previous R67 DHFR-folate complex provides insight into the catalytic mechanism and indicates that the relative orientation of the reactants in plasmid DHFR is different from that seen in chromosomal DHFRs.

How sensitive are nanosecond molecular dynamics simulations of proteins to changes in the force field?

The sensitivity of molecular dynamics simulations to variations in the force field has been examined in relation to a set of 36 structures corresponding to 31 proteins simulated by using different versions of the GROMOS force field. The three parameter sets used (43a1, 53a5, and 53a6) differ significantly in regard to the nonbonded parameters for polar functional groups and their ability to reproduce the correct solvation and partitioning behavior of small molecular analogues of the amino acid side chains. Despite the differences in the force field parameters no major differences could be detected in a wide range of structural properties such as the root-mean-square deviation from the experimental structure, radii of gyration, solvent accessible surface, secondary structure, or hydrogen bond propensities on a 5 to 10 ns time scale. The small differences that were observed correlated primarily with the presence of charged residues as opposed to residues that differed most between the parameter sets. The work highlights the variation that can be observed in nanosecond simulations of protein systems and implications of this for force field validation, as well as for the analysis of protein simulations in general.

Protein structure prediction by all-atom free-energy refinement

Background

The reliable prediction of protein tertiary structure from the amino acid sequence remains challenging even for small proteins. We have developed an all-atom free-energy protein forcefield (PFF01) that we could use to fold several small proteins from completely extended conformations. Because the computational cost of de-novo folding studies rises steeply with system size, this approach is unsuitable for structure prediction purposes. We therefore investigate here a low-cost free-energy relaxation protocol for protein structure prediction that combines heuristic methods for model generation with all-atom free-energy relaxation in PFF01.

Results

We use PFF01 to rank and cluster the conformations for 32 proteins generated by ROSETTA. For 22/10 high-quality/low quality decoy sets we select near-native conformations with an average C_α root mean square deviation of 3.03 Å/6.04 Å. The protocol incorporates an inherent reliability indicator that succeeds for 78% of the decoy sets. In over 90% of these cases near-native conformations are selected from the decoy set. This success rate is rationalized by the quality of the decoys and the selectivity of the PFF01 forcefield, which ranks near-native conformations an average 3.06 standard deviations below that of the relaxed decoys (Z-score).

Conclusion

All-atom free-energy relaxation with PFF01 emerges as a powerful low-cost approach toward generic de-novo protein structure prediction. The approach can be applied to large all-atom decoy sets of any origin and requires no preexisting structural information to identify the native conformation. The study provides evidence that a large class of proteins may be foldable by PFF01.

Integron-sequestered dihydrofolate reductase: a recently redeployed enzyme

The introduction and wide use of antibacterial drugs has resulted in the emergence of resistant organisms. DfrB dihydrofolate reductase (DHFR) is a bacterial enzyme that is uniquely associated with mobile gene cassettes within integrons, and confers resistance to the drug trimethoprim. This enzyme has intrigued microbiologists since it was discovered more than thirty years ago because of its simple structure, enzymatic inefficiency and its virtual insensitivity to trimethoprim. Here, for the first time, a comprehensive discussion of genetic, evolutionary, structural and functional studies of this enzyme is presented together. This information supports the ideas that DfrB DHFR is a poorly adapted catalyst and has recently been recruited to perform a novel enzymatic activity in response to selective pressure.

Multiple ligand-binding modes in bacterial R67 dihydrofolate reductase

R67 dihydrofolate reductase (DHFR), a bacterial plasmid-encoded enzyme associated with resistance to the drug trimethoprim, shows neither sequence nor structural homology with the chromosomal DHFR. It presents a highly symmetrical toroidal structure, where four identical monomers contribute to the unique central active-site pore. Two reactants (dihydrofolate, DHF), two cofactors (NADPH) or one of each (R67*DHF*NADPH) can be found simultaneously within the active site, the last one being the reactive ternary complex. As the positioning of the ligands has proven elusive to empirical determination, we addressed the problem from a theoretical perspective. Several potential structures of the ternary complex were generated using the docking programs AutoDock and FlexX. The variability among the final poses, many of which conformed to experimental data, prompted us to perform a comparative scoring analysis and molecular dynamics simulations to assess the stability of the complexes. Analysis of ligand-ligand and ligand-protein interactions along the 4 ns trajectories of eight different structures allowed us to identify important inter-ligand contacts and key protein residues. Our results, combined with published empirical data, clearly suggest that multipe binding modes of the ligands are possible within R67 DHFR. While the pterin ring of DHF and the nicotinamide ring of NADPH assume a stacked endo-conformation at the centre of the pore, probably assisted by V66, Q67 and I68, the tails of the molecules extend towards opposite ends of the cavity, adopting multiple configurations in a solvent rich-environment where hydrogen-bond interactions with K32 and Y69 may play important roles.

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.

Combinatorial exploration of the catalytic site of a drug-resistant dihydrofolate reductase: creating alternative functional configurations

We have applied a global approach to enzyme active site exploration, where multiple mutations were introduced combinatorially at the active site of Type II R67 dihydrofolate reductase (R67 DHFR), creating numerous new active site environments within a constant framework. By this approach, we combinatorially modified all 16 principal amino acids that constitute the active site of this enzyme. This approach is fundamentally different from active site point mutation in that the native active site context is no longer accounted for. Among the 1536 combinatorially mutated active site variants of R67 DHFR we created, we selected and kinetically characterized three variants with highly altered active site compositions. We determined that they are of high fitness, as defined by a complex function consisting jointly of catalytic activity and resistance to trimethoprim. The k(cat) and K(M) values were similar to those for the native enzyme. The favourable Delta(DeltaG) values obtained (ranging from -0.72 to -1.08 kcal/mol) suggest that, despite their complex mutational pattern, no fundamental change in the catalytic mechanism has occurred. We illustrate that combinatorial active site mutagenesis can allow for the creation of compensatory mutations that could not be predicted and thus provides a route for more extensive exploration of functional sequence space than is allowed by point mutation.

"Catch 222," the effects of symmetry on ligand binding and catalysis in R67 dihydrofolate reductase as determined by mutations at Tyr-69

R67 dihydrofolate reductase (R67 DHFR) catalyzes the transfer of a hydride ion from NADPH to dihydrofolate, generating tetrahydrofolate. The homotetrameric enzyme provides a unique environment for catalysis as both ligands bind within a single active site pore possessing 222 symmetry. Mutation of one active site residue results in concurrent mutation of three additional symmetry-related residues, causing large effects on binding of both ligands as well as catalysis. For example, mutation of symmetry-related tyrosine 69 residues to phenylalanine (Y69F), results in large increases in Km values for both ligands and a 2-fold rise in the kcat value for the reaction (Strader, M. B., Smiley, R. D., Stinnett, L. G., VerBerkmoes, N. C., and Howell, E. E. (2001) Biochemistry 40, 11344-11352). To understand the interactions between specific Tyr-69 residues and each ligand, asymmetric Y69F mutants were generated that contain one to four Y69F mutations. A general trend observed from isothermal titration calorimetry and steady-state kinetic studies of these asymmetric mutants is that increasing the number of Y69F mutations results in an increase in the Kd and Km values. In addition, a comparison of steady-state kinetic values suggests that two Tyr-69 residues in one half of the active site pore are necessary for NADPH to exhibit a wild-type Km value. A tyrosine 69 to leucine mutant was also generated to approach the type(s) of interaction(s) occurring between Tyr-69 residues and the ligands. These studies suggest that the hydroxyl group of Tyr-69 is important for interactions with NADPH, whereas both the hydroxyl group and hydrophobic ring atoms of the Tyr-69 residues are necessary for proper interactions with dihydrofolate.

Role of Lys-32 residues in R67 dihydrofolate reductase probed by asymmetric mutations

R67 dihydrofolate reductase (R67 DHFR) is a novel protein encoded by an R-plasmid that confers resistance to the antibiotic, trimethoprim. This homotetrameric enzyme possesses 222 symmetry, which imposes numerous constraints on the single active site pore, including a "one-site-fits-both" strategy for binding its ligands, dihydrofolate (DHF) and NADPH. Previous studies uncovered salt effects on binding and catalysis (Hicks, S. N., Smiley, R. D., Hamilton, J. B., and Howell, E. E. (2003) Biochemistry 42, 10569-10578), however the one or more residues that participate in ionic contacts with the negatively charged tail of DHF as well as the phosphate groups in NADPH were not identified. Several studies predict that Lys-32 residues were involved, however mutations at this residue destabilize the R67 DHFR homotetramer. To study the role of Lys-32 in binding and catalysis, asymmetric K32M mutations have been utilized. To create asymmetry, individual mutations were added to a tandem array of four in-frame gene copies. These studies show one K32M mutation is tolerated quite well, whereas addition of two mutations has variable effects. Two double mutants, K32M:1+2 and K32M: 1+4, which place the mutations on opposite sides of the pore, reduce kcat. However a third double mutant, K32M: 1+3, that places two mutations on the same half pore, enhances kcat 4- to 5-fold compared with the parent enzyme, albeit at the expense of weaker binding of ligands. Because the kcat/Km values for this double mutant series are similar, these mutations appear to have uncovered some degree of non-productive binding. This non-productive binding mode likely arises from formation of an ionic interaction that must be broken to allow access to the transition state. The K32M:1+3 mutant data suggest this interaction is an ionic interaction between Lys-32 and the charged tail of dihydrofolate. This unusual catalytic scenario arises from the 222 symmetry imposed on the single active site pore.

Effect of multiple symmetries on the association of R67 DHFR subunits bearing interfacial complementing mutations

It was shown previously that complementation could be a powerful mean to probe protein-protein interactions in the normally tetrameric R67 DHFR. Indeed, mixing complementing inactive dimeric mutants produced active heterotetramers. This approach turned a homo-oligomer into a hetero-oligomer and thus allowed the use of combinatorial assays, a subtle analysis of the association forces, and a precise determination of the equilibrium dissociation constants (K(D)) by titrimetry. However, for some of the complementing pairs, the experimental data implied multiple equilibria involving heterodimers, although no monomers could be detected. Thus, the reactions involved had to be identified to elaborate a suitable model to determine the K(D) of those pairs correctly. That model suggested that homodimers associated rapidly before the protomers could be redistributed in a multiple equilibrium system. Kinetic data confirmed that view. The association data at equilibrium were analyzed by multiple curve fitting with all plausible combinations of parameters. This gave a confidence interval for K(D) that is safer than the usual 67% or 90% confidence interval. Finally, the K(D) of one specific reaction, the dissociation of a heterotetramer with the relevant symmetry into two homodimers could be determined with the relevant model for each complementing pair, although multiple equilibria were present. These K(D) can thus be used as a set of references data to test and improve theoretical methods such as association free energy calculations.

Solution structure of folic acid. Molecular mechanics and NMR investigation

The structure of folic acid in solution was investigated by nuclear magnetic resonance (NMR) and theoretical calculations. Dynamical information and geometrical constraints were obtained by carbon-13 relaxation study, homo-nuclear NOESY spectra and hetero-nuclear 1H-13C NOE experiments. This set of experimental data was used for the molecular mechanics and molecular dynamic calculations. The accuracy of the final structure was established by the R(NMR) factor, which was calculated comparing the experimental NOESY cross-peaks intensities and the corresponding values simulated by using the complete relaxation matrix analysis (CORMA) approach.

Dramatic acceleration of protein folding by stabilization of a nonnative backbone conformation

Through a mutagenic investigation of Gly-48, a highly conserved position in the Src homology 3 domain, we have discovered a series of amino acid substitutions that are highly destabilizing, yet dramatically accelerate protein folding, some up to 10-fold compared with the wild-type rate. The unique folding properties of these mutants allowed for accurate measurement of their folding and unfolding rates in water with no denaturant by using an NMR spin relaxation dispersion technique. A strong correlation was found between beta-sheet propensity and the folding rates of the Gly-48 mutants, even though Gly-48 lies in an unusual non-beta-strand backbone conformation in the native state. This finding indicates that the accelerated folding rates of the Gly-48 mutants are the result of stabilization of a nonnative beta-strand conformation in the transition-state structure at this position, thus providing the first, to our knowledge, experimentally elucidated example of a mechanism by which folding can occur fastest through a nonnative conformation. We also demonstrate that residues that are most stabilizing in the transition-state structure are most destabilizing in the native state, and also cause the greatest reductions in in vitro functional activity. These data indicate that the unusual native conformation of the Gly-48 position is important for function, and that evolutionary selection for function can result in a domain that folds at a rate far below the maximum possible.

Refinement of homology-based protein structures by molecular dynamics simulation techniques

The use of classical molecular dynamics simulations, performed in explicit water, for the refinement of structural models of proteins generated ab initio or based on homology has been investigated. The study involved a test set of 15 proteins that were previously used by Baker and coworkers to assess the efficiency of the ROSETTA method for ab initio protein structure prediction. For each protein, four models generated using the ROSETTA procedure were simulated for periods of between 5 and 400 nsec in explicit solvent, under identical conditions. In addition, the experimentally determined structure and the experimentally derived structure in which the side chains of all residues had been deleted and then regenerated using the WHATIF program were simulated and used as controls. A significant improvement in the deviation of the model structures from the experimentally determined structures was observed in several cases. In addition, it was found that in certain cases in which the experimental structure deviated rapidly from the initial structure in the simulations, indicating internal strain, the structures were more stable after regenerating the side-chain positions. Overall, the results indicate that molecular dynamics simulations on a tens to hundreds of nanoseconds time scale are useful for the refinement of homology or ab initio models of small to medium-size proteins.

Functional evolution of two subtly different (similar) folds

BACKGROUND

The function of proteins is a direct consequence of their three-dimensional structure. The structural classification of proteins describes the ways of folding patterns all proteins could adopt. Although, the protein folds were described in many ways the functional properties of individual folds were not studied.

RESULTS

We have analyzed two beta-barrel folds generally adopted by small proteins to be looking similar but have different topology. On the basis of the topology they could be divided into two different folds named SH3-fold and OB-fold. There was no sequence homology between any of the proteins considered. The sequence diversity and loop variability was found to be important for various binding functions.

CONCLUSIONS

The function of Oligonucleotide/oligosaccharide-binding (OB) fold proteins was restricted to either DNA/RNA binding or sugar binding whereas the Src homology 3 (SH3) domain like proteins bind to a variety of ligands through loop modulations. A question was raised whether the evolution of these two folds was through DNA shuffling.

Relative stability of protein structures determined by X-ray crystallography or NMR spectroscopy: a molecular dynamics simulation study

The relative stability of protein structures determined by either X-ray crystallography or nuclear magnetic resonance (NMR) spectroscopy has been investigated by using molecular dynamics simulation techniques. Published structures of 34 proteins containing between 50 and 100 residues have been evaluated. The proteins selected represent a mixture of secondary structure types including all alpha, all beta, and alpha/beta. The proteins selected do not contain cysteine-cysteine bridges. In addition, any crystallographic waters, metal ions, cofactors, or bound ligands were removed before the systems were simulated. The stability of the structures was evaluated by simulating, under identical conditions, each of the proteins for at least 5 ns in explicit solvent. It is found that not only do NMR-derived structures have, on average, higher internal strain than structures determined by X-ray crystallography but that a significant proportion of the structures are unstable and rapidly diverge in simulations.

SycE allows secretion of YopE-DHFR hybrids by the Yersinia enterocolitica type III Ysc system

The Ysc type III secretion system allows Yersinia enterocolitica to translocate virulence proteins, called Yop effectors, into the cytosol of eukaryotic cells. Some of the Yop effectors possess an individual chaperone called a Syc protein. The first 15 amino acids of the YopE effector constitute a secretion signal that is sufficient to promote secretion of several reporter proteins. Residues 15-50 of YopE comprise the minimal binding domain for the SycE chaperone. In this study, we investigated the secretion by the Ysc system of several YopE-DHFR hybrid proteins with different folding properties, and evaluated the role of SycE, the cognate chaperone of YopE, in this context. We have analysed the secretion of hybrids containing 16 (YopE16), 52 (YopE52) and 80 (the complete region covered by the chaperone, YopE80) amino acids of YopE or full-length YopE (YopEFL) with wild-type DHFR and two mutants with altered folding properties. The hybrids containing DHFR delta77, the mutant whose folding properties are the most drastically affected, could be secreted in all the conditions tested, even in the absence of the chaperone SycE. In contrast, DHFRwt could only be secreted fused to the first 52 amino acids of YopE, and its secretion was strictly dependent on SycE. The hybrids YopE80-DHFRwt and YopEFL-DHFRwt were not secreted. YopEFL-DHFRwt completely jammed the channel in an SycE-dependent fashion. Our experiments indicate that, in order to be secreted, proteins must be unfolded or only partially folded, and that TSS chaperones could keep their substrates in a secretion-competent conformation, probably by preventing their folding. In addition, they show that the secretion apparatus can reject folded proteins if they are not deeply engaged into the injectisome.

Intact protein analysis for site-directed mutagenesis overexpression products: plasmid-encoded R67 dihydrofolate reductase

Mass spectrometry is currently the method of choice for the analysis of recombinant protein expression products. By combining proteolytic digestion with peptide mapping and tandem mass spectrometry techniques, verification of site-directed mutagenesis products can be obtained. The proteolytic digestion step converts a purified recombinant protein into a mixture that must be reseparated, thus greatly increasing the analysis time associated with the confirmation of site-directed mutagenesis products. Ion/ion reaction chemistry combined with quadrupole ion trap mass spectrometry provides a fast and efficient way to analyze intact proteins for the correct site-directed mutagenesis products, without heavy reliance on the proteolytic digestion step. Analysis of a series of protein variants (I68M, I68Q, Y69F, and Q67Y) from plasmid-encoded R67 dihydrofolate reductase using ion/ion reaction chemistry confirmed the presence of the correct site-directed mutagenesis products. For the I68M mutant, ion/ion separations detected the presence of extensive degradation from the N-terminal end of the protein. In the case of the Q67Y mutant, a mixture of Q67Y and Q67C species was detected by employing tandem mass spectrometry combined with ion/ion reactions. The ion/ion reaction technique was also performed on a partially purified lysate of the Q67Y/C mixture and successfully screened for the presence of both components in a complex mixture. The ion/ion reaction approach achieved the same results as the proteolytic-digestion-based methodology in a much shorter analysis time.

An Src homology 3-like domain is responsible for dimerization of the repressor protein KorB encoded by the promiscuous IncP plasmid RP4

KorB is a regulatory protein encoded by the conjugative plasmid RP4 and a member of the ParB family of bacterial partitioning proteins. The protein regulates the expression of plasmid genes whose products are involved in replication, transfer, and stable inheritance of RP4 by binding to palindromic 13-bp DNA sequences (5'-TTTAGC(G/C)GCTAAA-3') present 12 times in the 60-kb plasmid. Here we report the crystal structure of KorB-C, the C-terminal domain of KorB comprising residues 297-358. The structure of KorB-C was solved in two crystal forms. Quite unexpectedly, we find that KorB-C shows a fold closely resembling the Src homology 3 (SH3) domain, a fold well known from proteins involved in eukaryotic signal transduction. From the arrangement of molecules in the asymmetric unit, it is concluded that two molecules form a functionally relevant dimer. The detailed analysis of the dimer interface and a chemical cross-linking study suggest that the C-terminal domain is responsible for stabilizing the dimeric form of KorB in solution to facilitate binding to the palindromic operator sequence. The KorB-C crystal structure extends the range of protein-protein interactions known to be promoted by SH3 and SH3-like domains.

Vibrational structure of dihydrofolate bound to R67 dihydrofolate reductase

R67 is a Type II dihydrofolate reductase (DHFR) that catalyzes the reduction of dihydrofolate (DHF) to tetrahydrofolate by facilitating the addition of a proton to N5 of DHF and the transfer of a hydride ion from NADPH to C6. Because this enzyme is a plasmid-encoded DHFR from trimethoprim-resistant bacteria, extensive studies on R67 with various methods have been performed to elucidate its reaction mechanism. Here, Raman difference measurements, conducted on the ternary complex of R67.NADP(+).DHF believed to be an accurate mimic of the productive DHFR.NADPH.DHF complex, show that the pK(a) of N5 in the complex is less than 4. This is in clear contrast to the behavior observed in Escherichia coli DHFR, a substantially more efficient enzyme, where the pK(a) of bound DHF at N5 is increased to 6.5 compared with its solution value of 2.6. A comparison of the ternary complexes in R67 and E. coli DHFRs suggests that enzymic raising of the pK(a) at N5 can significantly increase the catalytic efficiency of the hydride transfer step. However, R67 shows that even without such a strategy an effective DHFR can still be designed.

One site fits both: a model for the ternary complex of folate + NADPH in R67 dihydrofolate reductase, a D2 symmetric enzyme

R67 dihydrofolate reductase (DHFR) is a novel enzyme that confers resistance to the antibiotic trimethoprim. The crystal structure of R67 DHFR displays a toroidal structure with a central active-site pore. This homotetrameric protein exhibits 222 symmetry, with only a few residues from each chain contributing to the active site, so related sites must be used to bind both substrate (dihydrofolate) and cofactor (NADPH) in the productive R67 DHFR.NADPH.dihydrofolate complex. Whereas the site of folate binding has been partially resolved crystallographically, an interesting question remains: how can the highly symmetrical active site also bind and orient NADPH for catalysis? To model this ternary complex, we employed DOCK and SLIDE, two methods for docking flexible ligands into proteins using quite different algorithms. The bound pteridine ring of folate (Fol I) from the crystal structure of R67 DHFR was used as the basis for docking the nicotinamide-ribose-Pi (NMN) moiety of NADPH. NMN was positioned by both DOCK and SLIDE on the opposite side of the pore from Fol I, where it interacts with Fol I at the pore's center. Numerous residues serve dual roles in binding. For example, Gln 67 from both the B and D subunits has several contacts with the pteridine ring, while the same residue from the A and C subunits has several contacts with the nicotinamide ring. The residues involved in dual roles are generally amphipathic, allowing them to make both hydrophobic and hydrophilic contacts with the ligands. The result is a 'hot spot' binding surface allowing the same residues to co-optimize the binding of two ligands, and orient them for catalysis.