Increasing diterpene yield with a modular metabolic engineering system in E. coli: comparison of MEV and MEP isoprenoid precursor pathway engineering

Engineering biosynthetic pathways in heterologous microbial host organisms offers an elegant approach to pathway elucidation via the incorporation of putative biosynthetic enzymes and characterization of resulting novel metabolites. Our previous work in Escherichia coli demonstrated the feasibility of a facile modular approach to engineering the production of labdane-related diterpene (20 carbon) natural products. However, yield was limited (<0.1 mg/L), presumably due to reliance on endogenous production of the isoprenoid precursors dimethylallyl diphosphate and isopentenyl diphosphate. Here, we report incorporation of either a heterologous mevalonate pathway (MEV) or enhancement of the endogenous methyl erythritol phosphate pathway (MEP) with our modular metabolic engineering system. With MEP pathway enhancement, it was found that pyruvate supplementation of rich media and simultaneous overexpression of three genes (idi, dxs, and dxr) resulted in the greatest increase in diterpene yield, indicating distributed metabolic control within this pathway. Incorporation of a heterologous MEV pathway in bioreactor grown cultures resulted in significantly higher yields than MEP pathway enhancement. We have established suitable growth conditions for diterpene production levels ranging from 10 to >100 mg/L of E. coli culture. These amounts are sufficient for nuclear magnetic resonance analyses, enabling characterization of enzymatic products and hence, pathway elucidation. Furthermore, these results represent an up to >1,000-fold improvement in diterpene production from our facile, modular platform, with MEP pathway enhancement offering a cost effective alternative with reasonable yield. Finally, we reiterate here that this modular approach is expandable and should be easily adaptable to the production of any terpenoid natural product

Four small puzzles that Rosetta doesn't solve

A complete macromolecule modeling package must be able to solve the simplest structure prediction problems. Despite recent successes in high resolution structure modeling and design, the Rosetta software suite fares poorly on deceptively small protein and RNA puzzles, some as small as four residues. To illustrate these problems, this manuscript presents extensive Rosetta results for four well-defined test cases: the 20-residue mini-protein Trp cage, an even smaller disulfide-stabilized conotoxin, the reactive loop of a serine protease inhibitor, and a UUCG RNA tetraloop. In contrast to previous Rosetta studies, several lines of evidence indicate that conformational sampling is not the major bottleneck in modeling these small systems. Instead, approximations and omissions in the Rosetta all-atom energy function currently preclude discriminating experimentally observed conformations from de novo models at atomic resolution. These molecular "puzzles" should serve as useful model systems for developers wishing to make foundational improvements to this powerful modeling suite.Comment: Published in PLoS One as a manuscript for the RosettaCon 2010 Special Collectio

Assessment of the Performance Potential of Advanced Subsonic Transport Concepts for NASA's Environmentally Responsible Aviation Project

NASA's Environmentally Responsible Aviation (ERA) project has matured technologies to enable simultaneous reductions in fuel burn, noise, and nitrogen oxide (NOx) emissions for future subsonic commercial transport aircraft. The fuel burn reduction target was a 50% reduction in block fuel burn (relative to a 2005 best-in-class baseline aircraft), utilizing technologies with an estimated Technology Readiness Level (TRL) of 4-6 by 2020. Progress towards this fuel burn reduction target was measured through the conceptual design and analysis of advanced subsonic commercial transport concepts spanning vehicle size classes from regional jet (98 passengers) to very large twin aisle size (400 passengers). Both conventional tube-and-wing (T+W) concepts and unconventional (over-wing-nacelle (OWN), hybrid wing body (HWB), mid-fuselage nacelle (MFN)) concepts were developed. A set of propulsion and airframe technologies were defined and integrated onto these advanced concepts which were then sized to meet the baseline mission requirements. Block fuel burn performance was then estimated, resulting in reductions relative to the 2005 best-in-class baseline performance ranging from 39% to 49%. The advanced single-aisle and large twin aisle T+W concepts had reductions of 43% and 41%, respectively, relative to the 737-800 and 777-200LR aircraft. The single-aisle OWN concept and the large twin aisle class HWB concept had reductions of 45% and 47%, respectively. In addition to their estimated fuel burn reduction performance, these unconventional concepts have the potential to provide significant noise reductions due, in part, to engine shielding provided by the airframe. Finally, all of the advanced concepts also have the potential for significant NOx emissions reductions due to the use of advanced combustor technology. Noise and NOx emissions reduction estimates were also generated for these concepts as part of the ERA project

Conceptual Design and Structural Optimization of NASA Environmentally Responsible Aviation (ERA) Hybrid Wing Body Aircraft

Simultaneously achieving the fuel consumption and noise reduction goals set forth by NASA's Environmentally Responsible Aviation (ERA) project requires innovative and unconventional aircraft concepts. In response, advanced hybrid wing body (HWB) aircraft concepts have been proposed and analyzed as a means of meeting these objectives. For the current study, several HWB concepts were analyzed using the Hybrid wing body Conceptual Design and structural optimization (HCDstruct) analysis code. HCDstruct is a medium-fidelity finite element based conceptual design and structural optimization tool developed to fill the critical analysis gap existing between lower order structural sizing approaches and detailed, often finite element based sizing methods for HWB aircraft concepts. Whereas prior versions of the tool used a half-model approach in building the representative finite element model, a full wing-tip-to-wing-tip modeling capability was recently added to HCDstruct, which alleviated the symmetry constraints at the model centerline in place of a free-flying model and allowed for more realistic center body, aft body, and wing loading and trim response. The latest version of HCDstruct was applied to two ERA reference cases, including the Boeing Open Rotor Engine Integration On an HWB (OREIO) concept and the Boeing ERA-0009H1 concept, and results agreed favorably with detailed Boeing design data and related Flight Optimization System (FLOPS) analyses. Following these benchmark cases, HCDstruct was used to size NASA's ERA HWB concepts and to perform a related scaling study

HOXB5 Cooperates with NKX2-1 in the Transcription of Human RET

The enteric nervous system (ENS) regulates peristaltic movement of the gut, and abnormal ENS causes Hirschsprung's disease (HSCR) in newborns. HSCR is a congenital complex genetic disorder characterised by a lack of enteric ganglia along a variable length of the intestine. The receptor tyrosine kinase gene (RET) is the major HSCR gene and its expression is crucial for ENS development. We have previously reported that (i) HOXB5 transcription factor mediates RET expression, and (ii) mouse with defective HOXB5 activity develop HSCR phenotype. In this study, we (i) elucidate the underlying mechanisms that HOXB5 mediate RET expression, and (ii) examine the interactions between HOXB5 and other transcription factors implicated in RET expression. We show that human HOXB5 binds to the promoter region 5′ upstream of the binding site of NKX2-1 and regulates RET expression. HOXB5 and NKX2-1 form a protein complex and mediate RET expression in a synergistic manner. HSCR associated SNPs at the NKX2-1 binding site (-5G>A rs10900296; -1A>C rs10900297), which reduce NKX2-1 binding, abolish the synergistic trans-activation of RET by HOXB5 and NKX2-1. In contrast to the synergistic activation of RET with NKX2-1, HOXB5 cooperates in an additive manner with SOX10, PAX3 and PHOX2B in trans-activation of RET promoter. Taken together, our data suggests that HOXB5 in coordination with other transcription factors mediates RET expression. Therefore, defects in cis- or trans-regulation of RET by HOXB5 could lead to reduction of RET expression and contribute to the manifestation of the HSCR phenotype

Surface free energy of polyurethane coatings with improved hydrophobicity

The polarity of polyurethane coats was studied on the basis of the goniometric method for determination of wetting angle values, on the basis of calculated surface free energy (SFE) values by the van Oss–Good and Owens–Wendt methods, and on the basis of polarity measurements with the use of the 1H NMR spectra. Test polyurethanes were synthesised in the reaction of methylene diphenyl 4,4′-diisocyanate (MDI) or 3-izocyanatomethyl –3,5,5- trimethylcyclohexyl isocyanate (IPDI) and polyoxyethylene glycols or polyesters poly(ε-caprolactone) diols and poly(ethyleneadipate) diol with different molecular weights, and some diols as chain extenders, in dioxane. The type of raw material was found to significantly affect the phase structure of the obtained polyurethane elastomers and to control physical interactions within those structures, thus influencing the SFE values. Fundamental reduction in the SFE value of a coating below 28 mJ/m2 was achieved by the use of 2,2,3,3-tetrafluoro-1,4-butanediol as the urethane prepolymer chain extender

Effect of the integration method on the accuracy and computational efficiency of free energy calculations using thermodynamic integration

Although calculations of free energy using molecular dynamics simulations have gained significant importance in the chemical and biochemical fields, they still remain quite computationally intensive. Furthermore, when using thermodynamic integration, numerical evaluation of the integral of the Hamiltonian with respect to the coupling parameter may introduce unwanted errors in the free energy. In this paper, we compare the performance of two numerical integration techniques-the trapezoidal and Simpson's rules and propose a new method, based on the analytic integration of physically based fitting functions that are able to accurately describe the behavior of the data. We develop and test our methodology by performing detailed studies on two prototype systems, hydrated methane and hydrated methanol, and treat Lennard-Jones and electrostatic contributions separately. We conclude that the widely used trapezoidal rule may introduce systematic errors in the calculation, but these errors are reduced if Simpson's rule is employed, at least for the electrostatic component. Furthermore, by fitting thermodynamic integration data, we are able to obtain precise free energy estimates using significantly fewer data points (5 intermediate states for the electrostatic component and 11 for the Lennard-Jones term), thus significantly decreasing the associated computational cost. Our method and improved protocol were successfully validated by computing the free energy of more complex systems hydration of 2-methylbutanol and of 4-nitrophenol-thus paving the way for widespread use in solvation free energy calculations of drug molecules

Molecular dynamics simulation of biomolecular systems

The group for computer-aided chemistry at the ETH Zurich focuses its research on the development of methodology to simulate the behavior of biomolecular systems and the use of simulation techniques to analyze and understand biomolecular processes at the atomic level. Here, the current research directions are briefly reviewed and illustrated with a few examples

Expression of Fraser syndrome genes in normal and polycystic murine kidneys

BACKGROUND: Fraser syndrome (FS) features renal agenesis and cystic kidneys. Mutations of FRAS1 (Fraser syndrome 1)and FREM2 (FRAS1-related extracellular matrix protein 2)cause FS. They code for basement membrane proteins expressed in metanephric epithelia where they mediate epithelial/mesenchymal signalling. Little is known about whether and where these molecules are expressed in more mature kidneys. METHODS: In healthy and congenital polycystic kidney (cpk)mouse kidneys we sought Frem2 expression using a LacZ reporter gene and quantified Fras family transcripts. Fras1 immunohistochemistry was undertaken in cystic kidneys from cpk mice and PCK (Pkhd1 mutant) rats (models of autosomal recessive polycystic kidney disease) and in wildtype metanephroi rendered cystic by dexamethasone. RESULTS: Nascent nephrons transiently expressed Frem2 in both tubule and podocyte epithelia. Maturing and adult collecting ducts also expressed Frem2. Frem2 was expressed in cpk cystic epithelia although Frem2 haploinsufficiency did not significantly modify cystogenesis in vivo. Fras1 transcripts were significantly upregulated, and Frem3 downregulated, in polycystic kidneys versus the non-cystic kidneys of littermates. Fras1 was immunodetected in cpk, PCK and dexamethasone-induced cystepithelia. CONCLUSIONS: These descriptive results are consistent with the hypothesis that Fras family molecules play diverse roles in kidney epithelia. In future, this should be tested by conditional deletion of FS genes in nephron segments and collecting ducts

Effects of Restrained Sampling Space and Nonplanar Amino Groups on Free-Energy Predictions for RNA with Imino and Sheared Tandem GA Base Pairs Flanked by GC, CG, iGiC or iCiG Base Pairs

Guanine-adenine (GA) base pairs play important roles in determining the structure, dynamics, and stability of RNA. In RNA internal loops, GA base pairs often occur in tandem arrangements and their structure is context and sequence dependent. Calculations reported here test the thermodynamic integration (TI) approach with the amber99 force field by comparing computational predictions of free energy differences with the free energy differences expected on the basis of NMR determined structures of the RNA motifs (5′-GCGGACGC-3′)2, (5′-GCiGGAiCGC-3′)2, (5′-GGCGAGCC-3′)2, and (5′-GGiCGAiGCC-3′)2. Here, iG and iC denote isoguanosine and isocytidine, which have amino and carbonyl groups transposed relative to guanosine and cytidine. The NMR structures show that the GA base pairs adopt either imino (cis Watson−Crick/Watson−Crick A-G) or sheared (trans Hoogsteen/Sugar edge A-G) conformations depending on the identity and orientation of the adjacent base pair. A new mixing function for the TI method is developed that allows alchemical transitions in which atoms can disappear in both the initial and final states. Unrestrained calculations gave ΔG° values 2−4 kcal/mol different from expectations based on NMR data. Restraining the structures with hydrogen bond restraints did not improve the predictions. Agreement with NMR data was improved by 0.7 to 1.5 kcal/mol, however, when structures were restrained with weak positional restraints to sample around the experimentally determined NMR structures. The amber99 force field was modified to partially include pyramidalization effects of the unpaired amino group of guanosine in imino GA base pairs. This provided little or no improvement in comparisons with experiment. The marginal improvement is observed when the structure has potential cross-strand out-of-plane hydrogen bonding with the G amino group. The calculations using positional restraints and a nonplanar amino group reproduce the signs of ΔG° from the experimental results and are, thus, capable of providing useful qualitative insights complementing the NMR experiments. Decomposition of the terms in the calculations reveals that the dominant terms are from electrostatic and interstrand interactions other than hydrogen bonds in the base pairs. The results suggest that a better description of the backbone is key to reproducing the experimental free energy results with computational free energy predictions