Peptide nucleic acid Hoogsteen strand linker design for major groove recognition of DNA thymine bases.

Sequence-specific targeting of double-stranded DNA and non-coding RNA via triple-helix-forming peptide nucleic acids (PNAs) has attracted considerable attention in therapeutic, diagnostic and nanotechnological fields. An E-base (3-oxo-2,3-dihydropyridazine), attached to the polyamide backbone of a PNA Hoogsteen strand by a side-chain linker molecule, is typically used in the hydrogen bond recognition of the 4-oxo group of thymine and uracil nucleic acid bases in the major groove. We report on the application of quantum chemical computational methods, in conjunction with spatial constraints derived from the experimental structure of a homopyrimidine PNA·DNA-PNA hetero-triplex, to investigate the influence of linker flexibility on binding interactions of the E-base with thymine and uracil bases in geometry-optimised model systems. Hydrogen bond formation between the N2 E-base atom and target pyrimidine base 4-oxo groups in model systems containing a ß-alanine linker (J Am Chem Soc 119:11116, 1997) was found to incur significant internal strain energy and the potential disruption of intra-stand aromatic base stacking interactions in an oligomeric context. In geometry-optimised model systems containing a 3-trans olefin linker (Bioorg Med Chem Lett 14:1551, 2004) the E-base swung out away from the target pyrimidine bases into the solvent. These findings are in qualitative agreement with calorimetric measurements in hybridisation experiments at T-A and U-A inversion sites. In contrast, calculations on a novel 2-cis olefin linker design indicate that it could permit simultaneous E-base hydrogen bonding with the thymine 4-oxo group, circumvention and solvent screening of the thymine 5-methyl group, and maintenance of triplex intra-stand base stacking interactions.

Construction of a synthetic pathway for the production of 1,3-propanediol from glucose.

In this work, we describe the construction of a synthetic metabolic pathway enabling direct biosynthesis of 1,3-propanediol (PDO) from glucose via the Krebs cycle intermediate malate. This non-natural pathway extends a previously published synthetic pathway for the synthesis of (L)-2,4-dihydroxybutyrate (L-DHB) from malate by three additional reaction steps catalyzed respectively, by a DHB dehydrogenase, a 2-keto-4-hydroxybutyrate (OHB) dehydrogenase and a PDO oxidoreductase. Screening and structure-guided protein engineering provided a (L)-DHB dehydrogenase from the membrane-associated (L)-lactate dehydrogenase of E. coli and OHB decarboxylase variants derived from the branched-chain keto-acid decarboxylase encoded by kdcA from Lactococcus lactis or pyruvate decarboxylase from Zymomonas mobilis. The simultaneous overexpression of the genes encoding these enzymes together with the endogenous ydhD-encoded aldehyde reductase enabled PDO biosynthesis from (L)-DHB. While the simultaneous expression of the six enzymatic activities in a single engineered E. coli strain resulted in a low production of 0.1 mM PDO from 110 mM glucose, a 40-fold increased PDO titer was obtained by co-cultivation of an E. coli strain expressing the malate-DHB pathway with another strain harboring the DHB-to-PDO pathway.

Rational engineering of a malate dehydrogenase for microbial production of 2,4-dihydroxybutyric acid via homoserine pathway.

A synthetic pathway for the production of 2,4-dihydroxybutyric acid from homoserine (HMS), composed of two consecutive enzymatic reaction steps has been recently reported. An important step in this pathway consists in the reduction in 2-keto-4-hydroxybutyrate (OHB) into (l)-dihydroxybutyrate (DHB), by an enzyme with OHB reductase activity. In the present study, we used a rational approach to engineer an OHB reductase by using the cytosolic (l)-malate dehydrogenase from Escherichia coli (Ec-Mdh) as the template enzyme. Structural analysis of (l)-malate dehydrogenase and (l)-lactate dehydrogenase enzymes acting on sterically cognate substrates revealed key residues in the substrate and co-substrate-binding sites responsible for substrate discrimination. Accordingly, amino acid changes were introduced in a stepwise manner into these regions of the protein. This rational engineering led to the production of an Ec-Mdh-5E variant (I12V/R81A/M85E/G179D/D86S) with a turnover number (kcat) on OHB that was increased by more than 2000-fold (from 0.03 up to 65.0 s-1), which turned out to be 7-fold higher than that on its natural substrate oxaloacetate. Further kinetic analysis revealed the engineered enzyme to possess comparable catalytic efficiencies (kcat/Km) between natural and synthetic OHB substrates (84 and 31 s-1 mM-1, respectively). Shake-flask cultivation of a HMS-overproducing E. coli strain expressing this improved OHB reductase together with a transaminase encoded by aspC able to convert HMS to OHB resulted in 89% increased DHB production as compared with our previous report using a E. coli host strain expressing an OHB reductase derived from the lactate dehydrogenase A of Lactococcus lactis.

A generic HTS assay for kinase screening: Validation for the isolation of an engineered malate kinase.

An end-point ADP/NAD+ acid/alkali assay procedure, directly applicable to library screening of any type of ATP-utilising/ADP producing enzyme activity, was implemented. Typically, ADP production is coupled to NAD+ co-enzyme formation by the conventional addition of pyruvate kinase and lactate dehydrogenase. Transformation of enzymatically generated NAD+ into a photometrically active alkali derivative product is then achieved through the successive application of acidic/alkali treatment steps. The assay was successfully miniaturized to search for malate kinase activity in a structurally-guided library of LysC aspartate kinase variants comprising 6,700 clones. The screening procedure enabled the isolation of nine positive variants showing novel kinase activity on (L)-malate, the best mutant, LysC V115A:E119S:E434V exhibited strong substrate selectivity for (L)-malate compared to (L)-aspartate with a (kcat/Km)malate/(kcat/Km)aspartate ratio of 86. Double mutants V115A:E119S, V115A:E119C and E119S:E434V were constructed to further probe the origins of stabilising substrate binding energy gains for (L)-malate due to mutation. The introduction of less sterically hindering side-chains in engineered enzymes carrying E119S and V115A mutations increases the effective volume available for substrate binding in the catalytic pocket. Improved binding of the (L)-malate substrate may be assisted by less hindered movement of the Phe184 aromatic side-chain. Additional favourable long-range electostatic effects on binding arising from the E434V surface mutation are conditionally dependent upon the presence of the V115A mutation close to Phe184 in the active-site.

Construction of a synthetic metabolic pathway for biosynthesis of the non-natural methionine precursor 2,4-dihydroxybutyric acid.

2,4-Dihydroxybutyric acid (DHB) is a molecule with considerable potential as a versatile chemical synthon. Notably, it may serve as a precursor for chemical synthesis of the methionine analogue 2-hydroxy-4-(methylthio)butyrate, thus, targeting a considerable market in animal nutrition. However, no natural metabolic pathway exists for the biosynthesis of DHB. Here we have therefore conceived a three-step metabolic pathway for the synthesis of DHB starting from the natural metabolite malate. The pathway employs previously unreported malate kinase, malate semialdehyde dehydrogenase and malate semialdehyde reductase activities. The kinase and semialdehyde dehydrogenase activities were obtained by rational design based on structural and mechanistic knowledge of candidate enzymes acting on sterically cognate substrates. Malate semialdehyde reductase activity was identified from an initial screening of several natural enzymes, and was further improved by rational design. The pathway was expressed in a minimally engineered Escherichia coli strain and produces 1.8 g l-1 DHB with a molar yield of 0.15.

An Atomistic Statistically Effective Energy Function for Computational Protein Design.

Shortcomings in the definition of effective free-energy surfaces of proteins are recognized to be a major contributory factor responsible for the low success rates of existing automated methods for computational protein design (CPD). The formulation of an atomistic statistically effective energy function (SEEF) suitable for a wide range of CPD applications and its derivation from structural data extracted from protein domains and protein-ligand complexes are described here. The proposed energy function comprises nonlocal atom-based and local residue-based SEEFs, which are coupled using a novel atom connectivity number factor to scale short-range, pairwise, nonbonded atomic interaction energies and a surface-area-dependent cavity energy term. This energy function was used to derive additional SEEFs describing the unfolded-state ensemble of any given residue sequence based on computed average energies for partially or fully solvent-exposed fragments in regions of irregular structure in native proteins. Relative thermal stabilities of 97 T4 bacteriophage lysozyme mutants were predicted from calculated energy differences for folded and unfolded states with an average unsigned error (AUE) of 0.84 kcal mol(-1) when compared to experiment. To demonstrate the utility of the energy function for CPD, further validation was carried out in tests of its capacity to recover cognate protein sequences and to discriminate native and near-native protein folds, loop conformers, and small-molecule ligand binding poses from non-native benchmark decoys. Experimental ligand binding free energies for a diverse set of 80 protein complexes could be predicted with an AUE of 2.4 kcal mol(-1) using an additional energy term to account for the loss in ligand configurational entropy upon binding. The atomistic SEEF is expected to improve the accuracy of residue-based coarse-grained SEEFs currently used in CPD and to extend the range of applications of extant atom-based protein statistical potentials.

Tri-peptide reference structures for the calculation of relative solvent accessible surface area in protein amino acid residues.

Relative amino acid residue solvent accessibility values allow the quantitative comparison of atomic solvent-accessible surface areas in different residue types and physical environments in proteins and in protein structural alignments. Geometry-optimised tri-peptide structures in extended solvent-exposed reference conformations have been obtained for 43 amino acid residue types at a high level of quantum chemical theory. Significant increases in side-chain solvent accessibility, offset by reductions in main-chain atom solvent exposure, were observed for standard residue types in partially geometry-optimised structures when compared to non-minimised models built from identical sets of proper dihedral angles abstracted from the literature. Optimisation of proper dihedral angles led most notably to marked increases of up to 54% in proline main-chain atom solvent accessibility compared to literature values. Similar effects were observed for fully-optimised tri-peptides in implicit solvent. The relief of internal strain energy was associated with systematic variation in N, C(α) and C(ß) atom solvent accessibility across all standard residue types. The results underline the importance of optimisation of 'hard' degrees of freedom (bond lengths and valence bond angles) and improper dihedral angle values from force field or other context-independent reference values, and impact on the use of standardised fixed internal co-ordinate geometry in sampling approaches to the determination of absolute values of protein amino acid residue solvent accessibility. Quantum chemical methods provide a useful and accurate alternative to molecular mechanics methods to perform energy minimisation of peptides containing non-standard (chemically modified) amino acid residues frequently present in experimental protein structure data sets, for which force field parameters may not be available. Reference tri-peptide atomic co-ordinate sets including hydrogen atoms are made freely available.

Probing impact of active site residue mutations on stability and activity of Neisseria polysaccharea amylosucrase.

The amylosucrase from Neisseria polysaccharea is a transglucosidase from the GH13 family of glycoside-hydrolases that naturally catalyzes the synthesis of α-glucans from the widely available donor sucrose. Interestingly, natural molecular evolution has modeled a dense hydrogen bond network at subsite -1 responsible for the specific recognition of sucrose and conversely, it has loosened interactions at the subsite +1 creating a highly promiscuous subsite +1. The residues forming these subsites are considered to be likely involved in the activity as well as the overall stability of the enzyme. To assess their role, a structure-based approach was followed to reshape the subsite -1. A strategy based on stability change predictions, using the FoldX algorithm, was considered to identify the best candidates for site-directed mutagenesis and guide the construction of a small targeted library. A miniaturized purification protocol was developed and both mutant stability and substrate promiscuity were explored. A range of 8 °C between extreme melting temperature values was observed and some variants were able to synthesize series of oligosaccharides with distributions differing from that of the parental enzyme. The crucial role of subsite -1 was thus highlighted and the biocatalysts generated can now be considered as starting points for further engineering purposes.

Role of glycoside phosphorylases in mannose foraging by human gut bacteria.

To metabolize both dietary fiber constituent carbohydrates and host glycans lining the intestinal epithelium, gut bacteria produce a wide range of carbohydrate-active enzymes, of which glycoside hydrolases are the main components. In this study, we describe the ability of phosphorylases to participate in the breakdown of human N-glycans, from an analysis of the substrate specificity of UhgbMP, a mannoside phosphorylase of the GH130 protein family discovered by functional metagenomics. UhgbMP is found to phosphorolyze ß-D-Manp-1,4-ß-D-GlcpNAc-1,4-D-GlcpNAc and is also a highly efficient enzyme to catalyze the synthesis of this precious N-glycan core oligosaccharide by reverse phosphorolysis. Analysis of sequence conservation within family GH130, mapped on a three-dimensional model of UhgbMP and supported by site-directed mutagenesis results, revealed two GH130 subfamilies and allowed the identification of key residues responsible for catalysis and substrate specificity. The analysis of the genomic context of 65 known GH130 sequences belonging to human gut bacteria indicates that the enzymes of the GH130_1 subfamily would be involved in mannan catabolism, whereas the enzymes belonging to the GH130_2 subfamily would rather work in synergy with glycoside hydrolases of the GH92 and GH18 families in the breakdown of N-glycans. The use of GH130 inhibitors as therapeutic agents or functional foods could thus be considered as an innovative strategy to inhibit N-glycan degradation, with the ultimate goal of protecting, or restoring, the epithelial barrier.

Adaptive Smith-Waterman residue match seeding for protein structural alignment.

The POLYFIT rigid-body algorithm for automated global pairwise and multiple protein structural alignment is presented. Smith-Waterman local alignment is used to establish a set of seed equivalences that are extended using Needleman-Wunsch dynamic programming techniques. Structural and functional interaction constraints provided by evolution are encoded as one-dimensional residue physical environment strings for alignment of highly structurally overlapped protein pairs. Local structure alignment of more distantly related pairs is carried out using rigid-body conformational matching of 15-residue fragments, with allowance made for less stringent conformational matching of metal-ion and small molecule ligand-contact, disulphide bridge, and cis-peptide correspondences. Protein structural plasticity is accommodated through the stepped adjustment of a single empirical distance parameter value in the calculation of the Smith-Waterman dynamic programming matrix. Structural overlap is used both as a measure of similarity and to assess alignment quality. Pairwise alignment accuracy has been benchmarked against that of 10 widely used aligners on the Sippl and Wiederstein set of difficult pairwise structure alignment problems, and more extensively against that of Matt, SALIGN, and MUSTANG in pairwise and multiple structural alignments of protein domains with low shared sequence identity in the SCOP-ASTRAL 40% compendium. The results demonstrate the advantages of POLYFIT over other aligners in the efficient and robust identification of matching seed residue positions in distantly related protein targets and in the generation of longer structurally overlapped alignment lengths. Superposition-based application areas include comparative modeling and protein and ligand design. POLYFIT is available on the Web server at http://polyfit.insa-toulouse.fr.

Orientation preferences of backbone secondary amide functional groups in peptide nucleic acid complexes: quantum chemical calculations reveal an intrinsic preference of cationic D-amino acid-based chiral PNA analogues for the P-form.

Geometric descriptions of nonideal interresidue hydrogen bonding and backbone-base water bridging in the minor groove are established in terms of polyamide backbone carbonyl group orientation from analyses of residue junction conformers in experimentally determined peptide nucleic acid (PNA) complexes. Two types of interresidue hydrogen bonding are identified in PNA conformers in heteroduplexes with nucleic acids that adopt A-like basepair stacking. Quantum chemical calculations on the binding of a water molecule to an O2 base atom in glycine-based PNA thymine dimers indicate that junctions modeled with P-form backbone conformations are lower in energy than a dimer comprising the predominant conformation observed in A-like helices. It is further shown in model systems that PNA analogs based on D-lysine are better able to preorganize in a conformation exclusive to P-form helices than is glycine-based PNA. An intrinsic preference for this conformation is also exhibited by positively charged chiral PNA dimers carrying 3-amino-D-alanine or 4-aza-D-leucine residue units that provide for additional rigidity by side-chain hydrogen bonding to the backbone carbonyl oxygen. Structural modifications stabilizing P-form helices may obviate the need for large heterocycles to target DNA pyrimidine bases via PNA.DNA-PNA triplex formation. Quantum chemical modeling methods are used to propose candidate PNA Hoogsteen strand designs.