Macromolecular modeling and design in Rosetta: recent methods and frameworks.

The Rosetta software for macromolecular modeling, docking and design is extensively used in laboratories worldwide. During two decades of development by a community of laboratories at more than 60 institutions, Rosetta has been continuously refactored and extended. Its advantages are its performance and interoperability between broad modeling capabilities. Here we review tools developed in the last 5 years, including over 80 methods. We discuss improvements to the score function, user interfaces and usability. Rosetta is available at http://www.rosettacommons.org.

Metal Dependence of the Xylose Isomerase from Piromyces sp. E2 Explored by Activity Profiling and Protein Crystallography.

Xylose isomerase from Piromyces sp. E2 (PirXI) can be used to equip Saccharomyces cerevisiae with the capacity to ferment xylose to ethanol. The biochemical properties and structure of the enzyme have not been described even though its metal content, catalytic parameters, and expression level are critical for rapid xylose utilization. We have isolated the enzyme after high-level expression in Escherichia coli, analyzed the metal dependence of its catalytic properties, and determined 12 crystal structures in the presence of different metals, substrates, and substrate analogues. The activity assays revealed that various bivalent metals can activate PirXI for xylose isomerization. Among these metals, Mn2+ is the most favorable for catalytic activity. Furthermore, the enzyme shows the highest affinity for Mn2+, which was established by measuring the activation constants (Kact) for different metals. Metal analysis of the purified enzyme showed that in vivo the enzyme binds a mixture of metals that is determined by metal availability as well as affinity, indicating that the native metal composition can influence activity. The crystal structures show the presence of an active site similar to that of other xylose isomerases, with a d-xylose binding site containing two tryptophans and a catalytic histidine, as well as two metal binding sites that are formed by carboxylate groups of conserved aspartates and glutamates. The binding positions and conformations of the metal-coordinating residues varied slightly for different metals, which is hypothesized to contribute to the observed metal dependence of the isomerase activity.

Coupling Protein Side-Chain and Backbone Flexibility Improves the Re-design of Protein-Ligand Specificity.

Interactions between small molecules and proteins play critical roles in regulating and facilitating diverse biological functions, yet our ability to accurately re-engineer the specificity of these interactions using computational approaches has been limited. One main difficulty, in addition to inaccuracies in energy functions, is the exquisite sensitivity of protein-ligand interactions to subtle conformational changes, coupled with the computational problem of sampling the large conformational search space of degrees of freedom of ligands, amino acid side chains, and the protein backbone. Here, we describe two benchmarks for evaluating the accuracy of computational approaches for re-engineering protein-ligand interactions: (i) prediction of enzyme specificity altering mutations and (ii) prediction of sequence tolerance in ligand binding sites. After finding that current state-of-the-art "fixed backbone" design methods perform poorly on these tests, we develop a new "coupled moves" design method in the program Rosetta that couples changes to protein sequence with alterations in both protein side-chain and protein backbone conformations, and allows for changes in ligand rigid-body and torsion degrees of freedom. We show significantly increased accuracy in both predicting ligand specificity altering mutations and binding site sequences. These methodological improvements should be useful for many applications of protein-ligand design. The approach also provides insights into the role of subtle conformational adjustments that enable functional changes not only in engineering applications but also in natural protein evolution.

Semi-rational engineering of cytochrome CYP153A from Marinobacter aquaeolei for improved ω-hydroxylation activity towards oleic acid.

ω-Hydroxy oleic acid is an important intermediate for the synthesis of certain polyesters and polyamides. In this study, a functional CYP153A/putidaredoxin (Pdx)/putidaredoxin reductase (Pdr) hybrid system was engineered for improved ω-hydroxylation activity towards oleic acid. By the combination of site-directed saturation mutagenesis (SDSM) and iterative saturation mutagenesis (ISM), a best mutant (Variant II) was obtained with mutations at two sites (S120 and P165) at the Pdx interaction interface with CYP153A, and one site (S453) in the substrate binding pocket. The in vitro-reconstituted activity of Variant II with purified Pdx and Pdr was 2.7-fold that of the template, while the whole cell transformation activity was 2.0-fold that of the template. A 96-well format-based screening scheme for CYP153A was also developed, which should be useful for engineering of other P450s with low activity. Kinetic analyses indicated that the activity improvement for CYP153A variants largely resulted from enhanced electron transfer. This further demonstrates the importance of the electron transfer between P450s and the non-native redox partners for the overall performance of hybrid P450 systems.

Valencene synthase from the heartwood of Nootka cypress (Callitropsis nootkatensis) for biotechnological production of valencene.

Nootkatone is one of the major terpenes in the heartwood of the Nootka cypress Callitropsis nootkatensis. It is an oxidized sesquiterpene, which has been postulated to be derived from valencene. Both valencene and nootkatone are used for flavouring citrus beverages and are considered among the most valuable terpenes used at commercial scale. Functional evaluation of putative terpene synthase genes sourced by large-scale EST sequencing from Nootka cypress wood revealed a valencene synthase gene (CnVS). CnVS expression in different tissues from the tree correlates well with nootkatone content, suggesting that CnVS represents the first dedicated gene in the nootkatone biosynthetic pathway in C. nootkatensis The gene belongs to the gymnosperm-specific TPS-d subfamily of terpenes synthases and its protein sequence has low similarity to known citrus valencene synthases. In vitro, CnVS displays high robustness under different pH and temperature regimes, potentially beneficial properties for application in different host and physiological conditions. Biotechnological production of sesquiterpenes has been shown to be feasible, but productivity of microbial strains expressing valencene synthase from Citrus is low, indicating that optimization of valencene synthase activity is needed. Indeed, expression of CnVS in Saccharomyces cerevisiae indicated potential for higher yields. In an optimized Rhodobacter sphaeroides strain, expression of CnVS increased valencene yields 14-fold to 352 mg/L, bringing production to levels with industrial potential.

Structural and time-resolved mechanistic investigations of protein hydrolysis by the acidic proline-specific endoprotease from Aspergillus niger.

Proline-specific endoproteases have been successfully used in, for example, the in-situ degradation of gluten, the hydrolysis of bitter peptides, the reduction of haze during beer production, and the generation of peptides for mass spectroscopy and proteomics applications. Here we present the crystal structure of the extracellular proline-specific endoprotease from Aspergillus niger (AnPEP), a member of the S28 peptidase family with rarely observed true proline-specific endoprotease activity. Family S28 proteases have a conventional Ser-Asp-His catalytic triad, but their oxyanion-stabilizing hole shows a glutamic acid, an amino acid not previously observed in this role. Since these enzymes have an acidic pH optimum, the presence of a glutamic acid in the oxyanion hole may confine their activity to an acidic pH. Yet, considering the presence of the conventional catalytic triad, it is remarkable that the A. niger enzyme remains active down to pH 1.5. The determination of the primary cleavage site of cytochrome c along with molecular dynamics-assisted docking studies indicate that the active site pocket of AnPEP can accommodate a reverse turn of approximately 12 amino acids with proline at the S1 specificity pocket. Comparison with the structures of two S28-proline-specific exopeptidases reveals not only a more spacious active site cavity but also the absence of any putative binding sites for amino- and carboxyl-terminal residues as observed in the exopeptidases, explaining AnPEP's observed endoprotease activity.

Introducing an Artificial Deazaflavin Cofactor in Escherichia coli and Saccharomyces cerevisiae.

Deazaflavin-dependent whole-cell conversions in well-studied and industrially relevant microorganisms such as Escherichia coli and Saccharomyces cerevisiae have high potential for the biocatalytic production of valuable compounds. The artificial deazaflavin FOP (FO-5'-phosphate) can functionally substitute the natural deazaflavin F420 and can be synthesized in fewer steps, offering a solution to the limited availability of the latter due to its complex (bio)synthesis. Herein we set out to produce FOP in vivo as a scalable FOP production method and as a means for FOP-mediated whole-cell conversions. Heterologous expression of the riboflavin kinase from Schizosaccharomyces pombe enabled in vivo phosphorylation of FO, which was supplied by either organic synthesis ex vivo, or by a coexpressed FO synthase in vivo, producing FOP in E. coli as well as in S. cerevisiae. Through combined approaches of enzyme engineering as well as optimization of expression systems and growth media, we further improved the in vivo FOP production in both organisms. The improved FOP production yield in E. coli is comparable to the F420 yield of native F420-producing organisms such as Mycobacterium smegmatis, but the former can be achieved in a significantly shorter time frame. Our E. coli expression system has an estimated production rate of 0.078 µmol L-1 h-1 and results in an intracellular FOP concentration of about 40 µM, which is high enough to support catalysis. In fact, we demonstrate the successful FOP-mediated whole-cell conversion of ketoisophorone using E. coli cells. In S. cerevisiae, in vivo FOP production by SpRFK using supplied FO was improved through media optimization and enzyme engineering. Through structure-guided enzyme engineering, a SpRFK variant with 7-fold increased catalytic efficiency compared to the wild type was discovered. By using this variant in optimized media conditions, FOP production yield in S. cerevisiae was 20-fold increased compared to the very low initial yield of 0.24 ± 0.04 nmol per g dry biomass. The results show that bacterial and eukaryotic hosts can be engineered to produce the functional deazaflavin cofactor mimic FOP.

Computation-Aided Engineering of Cytochrome P450 for the Production of Pravastatin.

CYP105AS1 is a cytochrome P450 from Amycolatopsis orientalis that catalyzes monooxygenation of compactin to 6-epi-pravastatin. For fermentative production of the cholesterol-lowering drug pravastatin, the stereoselectivity of the enzyme needs to be inverted, which has been partially achieved by error-prone PCR mutagenesis and screening. In the current study, we report further optimization of the stereoselectivity by a computationally aided approach. Using the CoupledMoves protocol of Rosetta, a virtual library of mutants was designed to bind compactin in a pro-pravastatin orientation. By examining the frequency of occurrence of beneficial substitutions and rational inspection of their interactions, a small set of eight mutants was predicted to show the desired selectivity and these variants were tested experimentally. The best CYP105AS1 variant gave >99% stereoselective hydroxylation of compactin to pravastatin, with complete elimination of the unwanted 6-epi-pravastatin diastereomer. The enzyme-substrate complexes were also examined by ultrashort molecular dynamics simulations of 50 × 100 ps and 5 × 22 ns, which revealed that the frequency of occurrence of near-attack conformations agreed with the experimentally observed stereoselectivity. These results show that a combination of computational methods and rational inspection could improve CYP105AS1 stereoselectivity beyond what was obtained by directed evolution. Moreover, the work lays out a general in silico framework for specificity engineering of enzymes of known structure.

Rhodoxanthin synthase from honeysuckle; a membrane diiron enzyme catalyzes the multistep conversation of ß-carotene to rhodoxanthin.

Rhodoxanthin is a vibrant red carotenoid found across the plant kingdom and in certain birds and fish. It is a member of the atypical retro class of carotenoids, which contain an additional double bond and a concerted shift of the conjugated double bonds relative to the more widely occurring carotenoid pigments, and whose biosynthetic origins have long remained elusive. Here, we identify LHRS (Lonicera hydroxylase rhodoxanthin synthase), a variant ß-carotene hydroxylase (BCH)-type integral membrane diiron enzyme that mediates the conversion of ß-carotene into rhodoxanthin. We identify residues that are critical to rhodoxanthin formation by LHRS. Substitution of only three residues converts a typical BCH into a multifunctional enzyme that mediates a multistep pathway from ß-carotene to rhodoxanthin via a series of distinct oxidation steps in which the product of each step becomes the substrate for the next catalytic cycle. We propose a biosynthetic pathway from ß-carotene to rhodoxanthin.

Saturation mutagenesis of Asn152 reveals a substrate selectivity switch in P99 cephalosporinase.

In class C beta-lactamases, the strictly conserved Asn152 forms part of an extended active-site hydrogen-bonding network. To probe its role in catalysis, all 19 mutants of Enterobacter cloacae P99 cephalosporinase Asn152 were simultaneously constructed and screened in Escherichia coli for their in vivo activity. The screen identified the previously uncharacterized mutants Asn152Ser, Asn152Thr, and Asn152Gly, which possess significant activity and altered substrate selectivity. In vitro measurement of Michaelis-Menten kinetic constants revealed that the Asn152Ser mutation causes a selectivity switch for penicillin G versus cefoxitin. Asn152Thr showed a 63-fold increase in k (cat) for oxacillin, a slow substrate for wild-type cephalosporinase. The results contribute to a growing body of data showing that mutation of highly conserved residues in the active site can result in substrate selectivity changes. The library screening method presented here would be applicable to substrate selectivity determination in other readily screenable enzymes.

Structure and mechanism of bacterial dehalogenases: different ways to cleave a carbon-halogen bond.

The dehalogenases make use of fundamentally different strategies to cleave carbon-halogen bonds. The structurally characterized haloalkane dehalogenases, haloacid dehalogenases and 4-chlorobenzoate-coenzyme A dehalogenases use substitution mechanisms that proceed via a covalent aspartyl intermediate. Recent X-ray crystallographic analysis of a haloalcohol dehalogenase and a trans-3-chloroacrylic acid dehalogenase has provided detailed insight into a different intramolecular substitution mechanism and a hydratase-like mechanism, respectively. The available information on the various dehalogenases supports different views on the possible evolutionary origins of their activities.

CorNet: Assigning function to networks of co-evolving residues by automated literature mining.

CorNet is a web-based tool for the analysis of co-evolving residue positions in protein super-family sequence alignments. CorNet projects external information such as mutation data extracted from literature on interactively displayed groups of co-evolving residue positions to shed light on the functions associated with these groups and the residues in them. We used CorNet to analyse six enzyme super-families and found that groups of strongly co-evolving residues tend to consist of residues involved in a same function such as activity, specificity, co-factor binding, or enantioselectivity. This finding allows to assign a function to residues for which no data is available yet in the literature. A mutant library was designed to mutate residues observed in a group of co-evolving residues predicted to be involved in enantioselectivity, but for which no literature data is available yet. The resulting set of mutations indeed showed many instances of increased enantioselectivity.

An engineered cryptic Hxt11 sugar transporter facilitates glucose-xylose co-consumption in Saccharomyces cerevisiae.

BACKGROUND: The yeast Saccharomyces cerevisiae is unable to ferment pentose sugars like d-xylose. Through the introduction of the respective metabolic pathway, S. cerevisiae is able to ferment xylose but first utilizes d-glucose before the d-xylose can be transported and metabolized. Low affinity d-xylose uptake occurs through the endogenous hexose (Hxt) transporters. For a more robust sugar fermentation, co-consumption of d-glucose and d-xylose is desired as d-xylose fermentation is in particular prone to inhibition by compounds present in pretreated lignocellulosic feedstocks. RESULTS: Evolutionary engineering of a d-xylose-fermenting S. cerevisiae strain lacking the major transporter HXT1-7 and GAL2 genes yielded a derivative that shows improved growth on xylose because of the expression of a normally cryptic HXT11 gene. Hxt11 also supported improved growth on d-xylose by the wild-type strain. Further selection for glucose-insensitive growth on d-xylose employing a quadruple hexokinase deletion yielded mutations at N366 of Hxt11 that reversed the transporter specificity for d-glucose into d-xylose while maintaining high d-xylose transport rates. The Hxt11 mutant enabled the efficient co-fermentation of xylose and glucose at industrially relevant sugar concentrations when expressed in a strain lacking the HXT1-7 and GAL2 genes. CONCLUSIONS: Hxt11 is a cryptic sugar transporter of S. cerevisiae that previously has not been associated with effective d-xylose transport. Mutagenesis of Hxt11 yielded transporters that show a better affinity for d-xylose as compared to d-glucose while maintaining high transport rates. d-glucose and d-xylose co-consumption is due to a redistribution of the sugar transport flux while maintaining the total sugar conversion rate into ethanol. This method provides a single transporter solution for effective fermentation on lignocellulosic feedstocks.

Engineering of an endogenous hexose transporter into a specific D-xylose transporter facilitates glucose-xylose co-consumption in Saccharomyces cerevisiae.

BACKGROUND: Engineering of Saccharomyces cerevisiae for the simultaneous utilization of hexose and pentose sugars is vital for cost-efficient cellulosic bioethanol production. This yeast lacks specific pentose transporters and depends on endogenous hexose transporters for low affinity pentose uptake. Consequently, engineered xylose-fermenting yeast strains first utilize D-glucose before D-xylose can be transported and metabolized. RESULTS: We have used an evolutionary engineering approach that depends on a quadruple hexokinase deletion xylose-fermenting S. cerevisiae strain to select for growth on D-xylose in the presence of high D-glucose concentrations. This resulted in D-glucose-tolerant growth of the yeast of D-xylose. This could be attributed to mutations at N367 in the endogenous chimeric Hxt36 transporter, causing a defect in D-glucose transport while still allowing specific uptake of D-xylose. The Hxt36-N367A variant transports D-xylose with a high rate and improved affinity, enabling the efficient co-consumption of D-glucose and D-xylose. CONCLUSIONS: Engineering of yeast endogenous hexose transporters provides an effective strategy to construct glucose-insensitive xylose transporters that are well integrated in the carbon metabolism regulatory network, and that can be used for efficient lignocellulosic bioethanol production.

Valencene oxidase CYP706M1 from Alaska cedar (Callitropsis nootkatensis).

(+)-Nootkatone is a natural sesquiterpene ketone used in grapefruit and citrus flavour compositions. It occurs in small amounts in grapefruit and is a major component of Alaska cedar (Callitropsis nootkatensis) heartwood essential oil. Upon co-expression of candidate cytochrome P450 enzymes from Alaska cedar in yeast with a valencene synthase, a C. nootkatensis valencene oxidase (CnVO) was identified to produce trans-nootkatol and (+)-nootkatone. Formation of (+)-nootkatone was detected at 144±10µg/L yeast culture. CnVO belongs to a new subfamily of the CYP706 family of cytochrome P450 oxidases.

Crystal structures of native and inactivated cis-3-chloroacrylic acid dehalogenase. Structural basis for substrate specificity and inactivation by (R)-oxirane-2-carboxylate.

The bacterial degradation pathways for the nematocide 1,3-dichloropropene rely on hydrolytic dehalogenation reactions catalyzed by cis- and trans-3-chloroacrylic acid dehalogenases (cis-CaaD and CaaD, respectively). X-ray crystal structures of native cis-CaaD and cis-CaaD inactivated by (R)-oxirane-2-carboxylate were elucidated. They locate four known catalytic residues (Pro-1, Arg-70, Arg-73, and Glu-114) and two previously unknown, potential catalytic residues (His-28 and Tyr-103'). The Y103F and H28A mutants of these latter two residues displayed reductions in cis-CaaD activity confirming their importance in catalysis. The structure of the inactivated enzyme shows covalent modification of the Pro-1 nitrogen atom by (R)-2-hydroxypropanoate at the C3 position. The interactions in the complex implicate Arg-70 or a water molecule bound to Arg-70 as the proton donor for the epoxide ring-opening reaction and Arg-73 and His-28 as primary binding contacts for the carboxylate group. This proposed binding mode places the (R)-enantiomer, but not the (S)-enantiomer, in position to covalently modify Pro-1. The absence of His-28 (or an equivalent) in CaaD could account for the fact that CaaD is not inactivated by either enantiomer. The cis-CaaD structures support a mechanism in which Glu-114 and Tyr-103' activate a water molecule for addition to C3 of the substrate and His-28, Arg-70, and Arg-73 interact with the C1 carboxylate group to assist in substrate binding and polarization. Pro-1 provides a proton at C2. The involvement of His-28 and Tyr-103' distinguishes the cis-CaaD mechanism from the otherwise parallel CaaD mechanism. The two mechanisms probably evolved independently as the result of an early gene duplication of a common ancestor.

The X-ray structure of the haloalcohol dehalogenase HheA from Arthrobacter sp. strain AD2: insight into enantioselectivity and halide binding in the haloalcohol dehalogenase family.

Haloalcohol dehalogenases are bacterial enzymes that cleave the carbon-halogen bond in short aliphatic vicinal haloalcohols, like 1-chloro-2,3-propanediol, some of which are recalcitrant environmental pollutants. They use a conserved Ser-Tyr-Arg catalytic triad to deprotonate the haloalcohol oxygen, which attacks the halogen-bearing carbon atom, producing an epoxide and a halide ion. Here, we present the X-ray structure of the haloalcohol dehalogenase HheA(AD2) from Arthrobacter sp. strain AD2 at 2.0-A resolution. Comparison with the previously reported structure of the 34% identical enantioselective haloalcohol dehalogenase HheC from Agrobacterium radiobacter AD1 shows that HheA(AD2) has a similar quaternary and tertiary structure but a much more open substrate-binding pocket. Docking experiments reveal that HheA(AD2) can bind both enantiomers of the haloalcohol substrate 1-p-nitrophenyl-2-chloroethanol in a productive way, which explains the low enantiopreference of HheA(AD2). Other differences are found in the halide-binding site, where the side chain amino group of Asn182 is in a position to stabilize the halogen atom or halide ion in HheA(AD2), in contrast to HheC, where a water molecule has taken over this role. These results broaden the insight into the structural determinants that govern reactivity and selectivity in the haloalcohol dehalogenase family.

Structural basis for the enantioselectivity of an epoxide ring opening reaction catalyzed by halo alcohol dehalogenase HheC.

Halo alcohol dehalogenase HheC catalyzes the highly enantioselective dehalogenation of vicinal halo alcohols to epoxides, as well as the reverse reaction, the enantioselective and beta-regioselective nucleophilic ring opening of epoxides by pseudo-halides such as azide and cyanide. To investigate this latter reaction, we determined X-ray structures of complexes of HheC with the favored and unfavored enantiomers of para-nitrostyrene oxide. The aromatic parts of the two enantiomers bind in a very similar way, but the epoxide ring of the unfavored (S)-enantiomer binds in a nonproductive inverted manner, with the epoxide oxygen and Cbeta atom positions interchanged with respect to those of the favored (R)-enantiomer. The calculated difference in relative Gibbs binding energy is in agreement with the observed loss of a single hydrogen bond in the S bound state with respect to the R bound state. Our results indicate that it is the nonproductive binding of the unfavored (S)-enantiomer, rather than the difference in affinity for the two enantiomers, that allows HheC to catalyze the azide-mediated ring opening of para-nitrostyrene oxide with high enantioselectivity. This work represents a rare opportunity to explain the enantioselectivity of an enzymatic reaction by comparison of crystallographic data on the binding of both the favored and unfavored enantiomers.

Improved catalytic properties of halohydrin dehalogenase by modification of the halide-binding site.

Halohydrin dehalogenase (HheC) from Agrobacterium radiobacter AD1 catalyzes the dehalogenation of vicinal haloalcohols by an intramolecular substitution reaction, resulting in the formation of the corresponding epoxide, a halide ion, and a proton. Halide release is rate-limiting during the catalytic cycle of the conversion of (R)-p-nitro-2-bromo-1-phenylethanol by the enzyme. The recent elucidation of the X-ray structure of HheC showed that hydrogen bonds between the OH group of Tyr187 and between the Odelta1 atom of Asn176 and Nepsilon1 atom of Trp249 could play a role in stabilizing the conformation of the halide-binding site. The possibility that these hydrogen bonds are important for halide binding and release was studied using site-directed mutagenesis. Steady-state kinetic studies revealed that mutant Y187F, which has lost both hydrogen bonds, has a higher catalytic activity (k(cat)) with two of the three tested substrates compared to the wild-type enzyme. Mutant W249F also shows an enhanced k(cat) value with these two substrates, as well as a remarkable increase in enantiopreference for (R)-p-nitro-2-bromo-1-phenylethanol. In case of a mutation at position 176 (N176A and N176D), a 1000-fold lower catalytic efficiency (k(cat)/K(m)) was obtained, which is mainly due to an increase of the K(m) value of the enzyme. Pre-steady-state kinetic studies showed that a burst of product formation precedes the steady state, indicating that halide release is still rate-limiting for mutants Y187F and W249F. Stopped-flow fluorescence experiments revealed that the rate of halide release is 5.6-fold higher for the Y187F mutant than for the wild-type enzyme and even higher for the W249F enzyme. Taken together, these results show that the disruption of two hydrogen bonds around the halide-binding site increases the rate of halide release and can enhance the overall catalytic activity of HheC.

Escherichia coli MltA: MAD phasing and refinement of a tetartohedrally twinned protein crystal structure.

Crystals were grown of a mutant form of the bacterial cell-wall maintenance protein MltA that diffracted to 2.15 A resolution. When phasing with molecular replacement using the native structure failed, selenium MAD was used to obtain initial phases. However, after MAD phasing the crystals were found to be tetartohedrally twinned, hampering correct space-group determination and refinement. A refinement protocol was designed to take tetartohedral twinning into account and was successfully applied to refine the structure. The refinement protocol is described and the reasons for the failure of molecular replacement and the success of MAD are discussed in terms of the effects of the tetartohedral twinning.