Design

To use BL as a substrate in an efficient way, we designed a metabolic network involving organelle engineering. To achieve this, we divided the project into three main parts:

• Selection of synthetic metabolic pathway. This pathway converts vanillic acid, a LDIC, into muconic acid, a nylon precursor.
• Organelle engineering. Protein engineering to encapsulate enzymes and enhance permeability of a bacterial microcompartment.
• Genetic circuit design. Would provide slightly regulation and optimization variables, as BMC assembly, growth and dynamic metabolic control.

Selection of the metabolic pathway.

To select the pathway, we established three objectives:

• The final product may be higher value added product that BL burning to produce energy.
• It must be valuable in the market
• The final product can also be extracted from petroleum derivatives.

BL composition analysis

The LDICs are very diverse because their structure depends on the composition of lignin from de organic-vegetal source. In this case, we analyzed the source reported from kraft pulp, generated by the paper industry. The components that are present in higher quantities are guaiacol, catechol, vanillin, vanillic acid and isoeugenol. (1,2) Aromatic phenolic compounds are hydrophobically predominant. Most of them can cross lipid membranes easily, which causes their specific toxicity. Their IC50 varies from 0.2g/L until 2 g/L for vanillin and vanillic acid. These factors were important when designing the system. However, the selection of the first substrate would impact the overall system.

Synthetic metabolic pathways for revaluation purposes

We investigated several ways to build or implement a synthetic route to use LDICs as a primary substrate, with the main goal of using them as raw material to synthesize other high-value products. As the LDICs are hard to work with; since there are not a lot of metabolic systems to implement new synthetic pathways. We first used software like RetroPath2.0 (3), a program that is useful to discover and predict new metabolic pathways; however, time and resources were limited. To understand and model metabolic systems, we eventually decided to select a constructed pathway reported before. We chose the synthetic pathway, which has a set of known enzymes and previously expressed in E. coli(4). The synthetic pathway uses as primary metabolite vanillic acid, as they oxidize to cis-cis muconic acid. Figure 2 (on the left side) shows information about the synthetic pathway. Cis-cis-muconic acid may be transformed in adipic acid, and caprolactam, which are used in nylon production. We design metabolic analysis experiments to detect metabolites and their intermediates in culture media. Using TLC as easy control and further quantitative data with HPLC-UV/Vis. (4),(5)

Protein engineering design of synthetic organelles

Spatial organization systems can develop a lot of advantages, one of them is compartmentalize chemical reactions, sequestering toxic intermediate and enhancing pathways flux (6) (Figure 3). As LDICs and the pathway selected has toxic intermediates, we implement the idea of encapsulation and isolation of metabolic systems using bacterial microcompartments (BMCs) (7–9). We explain first the basic functioning of BMCs, then which parts we design to use as a nanoreactor and organelle.

BMCs organelles, a useful and straightforward bioengineering chassis

BMCs are protein-based organelles. They are composed of multiple protein subunits, which self-assemble into an icosahedron complex between 100-200 nm of diameter. This structure can be compared to some viral capsids, and they can direct proteins inside through an encapsulation-tag peptide. When the proteins encapsulate enzymes and pathways, they are called metabolosomes, and a pore permits the metabolite diffusion through the different subunits. Their genomic organization is an operon of several genes and can be used to make easy organelle engineering, compared to lipid-based ones (10,11).

Function in nature

On the right side, the representation of Pdu system from Salmonella enterica. The shell is conformed by Pdu units which most are self-assembly protein hexamers (represented as PduA), and they have a pore with selectivity for substrates. Wild type Pdu has selectivity to propanediol, and inside the shell, it encapsulates enzymes and cofactors necessary to produce 1-propanol, an energetic pathway. Protein encapsulation is possible by an encapsulation peptide (EP) that redirects protein inside(Left side). We use these principles to create a new metabolosome.

Pathways and different types of BMCs

Different compartmentalization systems have been described in the past. We could have used Propanediol Utilization System (Pdu), ethanolamine utilization (Eut) and carboxysome (Ccm). We chose Pdu form Salmonella enterica, which is the best described in the literature.

Subunits

Pdu has several different units with the same structural and predominant domain called Pfam0936 (9) (Figure 7) to form empty shell complex, characteristic of BMCs systems. Pfam0936 permits interaction with other subunits in a zipper-like interaction and its highly stable. All units have convex and concave faces, corresponding with outer and inner shell faces respectively. This particularity has an impact on pore selectivity and shell assembly stability. In our system, there are only seven necessary units: PduA, PduB, PduU, PduJ, PduK, and PduN…pduABJKNU (14).

Shell proteins can be classified into different types as the roles in BMCs. BMC-H, BMC-P depending on their structure(11). BMC-H: PduA and PduJ, a hexameric protein which can form pores, PduU is a permuted pore with structural co-localization function in the cytoskeleton. BMC-T: PduB, a trimeric protein with pseudohexameric symmetry, which can form pores, and BMC-P: PduN, a pentameric protein with particular structural function.

Structural analysis for viable pore and protein engineering

BMC permeability engineering is one of the most challenge task actually in synthetic BMCs as the pore selectivity has been determined by two important factors, size and electrostatic surface interaction with metabolites(10,15). To determine the permeability of our desired metabolites we need to analyze different pores and engineering them.

PduA, PduB and PduJ (PduABJ) are the most abundant shell protein, with similar function and structure. They had a central pore with similar width that permit the pass of metabolites through them, which naturally is propanediol, however these are the biggest pores in BMCs systems reported. Indeed, this characteristic may not be deterministic way to know if vanillate, catechol, protocatechuate and muconic acid, and other cofactors may be permeable enough.We use CaverSuite(16–18) and Molecular Dynamics Simulations in GROMACS software to determine these parameters important in BMC engineering.

In a structural analysis, it is reported that for all shell monomer-monomer and subunit-subunit interactions, is essential the Pfam0936 domain (represented in orange, figure 8). An α Helix (denoted in white in concave face, figure 8) are motifs with interaction of encapsulation peptide. In PduABJ the subunits are almost full composed by this domain. In the structural view of other subunits, we found particularities in PduU subunit. Has an additional domain (not officially reported as a pfam) and distinct folding which makes possibly enhance permeability. PduU has a β barrel (denoted in white in convex face in PduU, figure 8) which has been computational determined that interact with pduV, at the same time this complex associate with cytoskeleton to attach in specific place in the cytoplasm. However, the shell can be formed without PduV, which essentially the β barrel an active function reported in shell. (19)

Analysis of cavity, electrostatic surfaces and tunnel generation with possible trajectories were made with CaverAnalyst software, reporting the real pore size in PduA are approximately 2.8 Å, almost the size of their selective metabolite propanediol. Residues T38,G39, S40, G41 that conform a little loop in the center pore, has almost all their selective characteristics. Polar groups and positively charged pore have an affinity to negatively-charged molecules, like propanediol.

Fig 10.Structure analysis of A. PduA and B. PduU. Showed sequence, secondary structures, and pfam domains (the only one is pfam 0936 denoted in the light blue bar). Red rectangles marked the pore residues in similar ubication of the domain, however there are structural differences in residues corresponding to the pore. PduU shows a larger sequence and the yellow rectangle is denoted the β barrel, in N terminal, of 15 residues long.

Sequence showed that almost all the protein is formed by their main domain, pfam0936, which is difficult to make a rational protein engineering to enhance the pore permeability, because mutations in secondary structures may make an unstable unit and unfolding of the protein. Free energy analysis of protein allows only loops to mutate them without unstabilize the protein, but it would-be high-risk mutations. However, PduU subunit its a candidate to engineering permeability, because has the biggest pore, and it is possible the engineering in β barrel that permute the pore.

Modification of PduU subunit to engineer pore permeability of Pdu

We made a deletion of 15 residues in the N-terminal sequence of PduU, creating PduUmd variant. β barrel was deleted, and the pore with bigger size and nonpolar partially charged pore may augment BMC permeability.

Fig. 12 Structural differences in PduA/PduUmd (Left) and PduU/Pdumd. Light blue PduUmd and violet PduA. Modelling demosntrate wide pores in the loops of center pore, however, their structures do not align exactly. (Right) Alignment between PduU (red) and PduUmd (blue), the pore is permuted by the β barrel which do not interacts with pfam0936 domain.

We use phyre2 and I-TASSER software to predict tertiary structures, from the deleted sequence. However, has the biggest deletion do not clarifies the stability of the protein. We proceed to make stability analysis by computational methods described in simulations-modeling page, as a part of design process. In resume, we use moleculr dynamics and docking softare, the analysis confirms that the complex can be formed and being functional.

It is showed 3 proteins, after their solvation (except for PduU) in GROMACS, which represent a realistic view of the pore. On the left of every set of images are results of CaverAnalyst, which are tunnel surfaces. Every possible tunnel is showed in different colors and get joined in the pore. On the right side, the surface model of every protein. A. PduA B. PduU C PduUmd.

As its showed, PduUmd gets a bigger pore, compared to PduA and it is not permuted like PduU wild type unit.

We can use experimental wet lab work to demonstrate the different units assembled. There are BMC extraction protocols with easy steps, further we can use SDS-PAGE and Atomic Force Microscopy to analyze the complex formed, encapsulated enzymes and subunits present. Permeability differences can be measured by titers in culture media of muconic acid and vanillic acid, also proven toxicity assays. We expect higher titers than reported, by the flux sequestration and toxic diminishing, however, energy demand of expression systems may slow bacterial metabolism. (8),(20)

Encapsulation peptides

Proteins can be encapsulated by signal peptide, named as Encapsulation Peptide (EP).(14,21) EPs may vary between wild type proteins as PduD, PduP and other enzimes, which is 18 to 20 residues in alpha helix conformation. It is highly conserved and may interact with several units in the concave face of Pdu shell units, directing a protein aggregate inside the BMC. (22)Also, these EP can interact with other EPs, augmenting the encapsulation rate (8,20). The interactions with N terminal of PduP and C-terminal of PduA was determined, indeed, we extract EP from PduP with a 19 long residues. The linker sequence may be required and important in the attachment of new proteins. then a short linker G4S was added.

Enzymes with EP

We add the selected EP (See add sequence structures) to the N- teminal and C- terminal sequences, depending their structural conformation using Phyre2 to modelling and predict tridimensional structures. CatA and AroY seems that N-terminal added EP (name) do not interact with catalityc sites, as LigM the signal peptide interact between monomers and get “buried” between them, blocking catalytic sites. EP was then added in C-terminal in a way to prevent this conformation.

Assembly

BMCs have the property of self-assembling. Analysis of the residues in the interaction zone between units has a “zipper” mechanism involving polar interactions(12). This residues are highly conserved, which makes possible hybrid BMCs. However, there are several variables involving the right assembling of BMCs, some are the protein-encapsulated concentration and BMCs subunits ratio. As the assembly depends in equilibrium between EP and subunits interactions, the size and stability depend in the expression rate of proteins. This is a problem afforded in the construction of genetic circuit in which is the regulation of expression can be slightly controlled. However, attempts reported in encapsulated protein vary between 1-5% of total protein expressed.(23) Our circuit constructed may be a great way to optimize the encapsulation and BMC stability, preventing formation of other types of structures, like nanotubes or inclusion bodies.(24)

Expession and regulation systems: Genetic circuits

Importance of expression regulation in our system:

• Assembly and encapsulation control(25)
• Dynamic control of metabolic systems(26)
• Strategies for metabolic optimization(27)

References

Links

Esta es una liga, al picarle te manda a la página que va dentro del código