Luis Fernando Garcia-Alles

Deep learning has revolutionized computational protein design, enabling the generation of sequences that fold onto target backbones with unprecedented accuracy. However, state-of-the-art inverse folding tools largely rely on auto-regressive sampling. While powerful, this paradigm is increasingly recognized for its inability to “think ahead”, a crucial capacity to reliably create the complex, long-range inter-residue dependencies essential for most biological functions. To overcome these fundamental limitations, we introduce EffieDes, a generative neuro-symbolic AI framework that synergizes the predictive capabilities of deep learning with the logical precision of automated reasoning. EffieDes leverages deep learning to encode the target backbone’s fitness landscape into Effie— a fully decomposable probabilistic graphical model (Potts model). This landscape is then rigorously explored by an automated reasoning prover to identify sequences that simultaneously satisfy complex design constraints and optimize backbone fitness. We validated this neuro-symbolic approach through the design of orthogonal sequence pairs that adopt identical folds but exhibit selective self-assembly, as well as the design of a de novo selective nanobody with nanomolar affinity for an immune-evasive SARS-CoV-2 variant. EffieDes provides a robust architecture for precisely dissecting learned fitness landscapes, offering a new path toward proteins with highly optimized performances and sophisticated functional objectives.

Bacterial microcompartments (BMC) are submicrometric reactors that encapsulate dedicated metabolic activities. BMC-H hexamers, the most abundant components of BMC shells, play major roles for shell plasticity and permeability. In part, chemical exchanges between the BMC lumen and the cellular cytosol will be defined by the disposition of amino acids lining the central BMC-H pores. Current models attribute to BMC-H a homo-oligomeric nature. The hexagonal symmetry of corresponding pores, however, would break down if hetero-hexamers formed, a possibility suggested by the frequent presence of multiple paralogs within BMC operons. Here, we gauged the degree of structural promiscuity between the 11 BMC-H paralogs from Klebsiella pneumoniae, a potential human pathogen endowed with the capacity to express three different BMC types. Concomitant activation of transcription of several BMC operons was first shown to be possible. By leveraging an adapted tripartite GFP technology, all possible BMC-H pair combinations were screened in E. coli. Multiple structural cross-compatibilities were pinpointed between homologs arising not only from the same BMC operon, but also from different BMC types, results supported by Alphafold and ESMFold predictions. The structural stability and assembly propensity of selected hetero-associations was established by biochemical means. In light of these results, we reinterpreted published lysine cross-linking mass spectrometry data to demonstrate that one of these hetero-hexamers, involving PduA and PduJ, was already detected to form in the shell of a recombinantly-expressed 1,2-propanediol utilization compartment from Salmonella enterica. Altogether, this study points to the need to embrace an augmented structural complexity in BMC shells.

Cytochrome P450 1A2, as many isoenzymes, can generate multiple metabolites from a single substrate. A loose coupling between substrate binding and oxygen activation makes possible substrate reorientations at the active site prior to catalysis. In the present work, caffeine oxidation to alternative bioactive compounds was used to decipher this pluripotency. A model involving two interacting subsites capable of sequentially accommodating one or two caffeine molecules was considered. Molecular dynamics was used to characterize subsite interactions and feed a dedicated geometric encoding of trajectories that was coupled to dimensional reductions and differential machine learning. The two subsites differentially control caffeine orientations and can exchange substrate through a phenylalanine gated mechanism. This exchange can be locked by the presence of a second bound molecule. Complementary roles of subsites in progressively determining the caffeine orientation during its approach to active oxygen were examined. Interestingly, substrate face flipping becomes impaired upon entry into the rather flat active site. This makes the mechanisms that define the orientation of caffeine relative to active oxygen dependent on the substrate face oriented toward heme. Globally, this evidenced that P450 1A2 regioselectivity results from local determinants combined with subsite interactions and caffeine face preselection at a longer distance.

Bacterial microcompartments (BMC) are submicrometric reactors that encapsulate dedicated metabolic activities. BMC-H hexamers, the most abundant components of BMC shells, play major roles for shell plasticity and permeability. In part, chemical exchanges between the BMC lumen and the cellular cytosol will be defined by the disposition of amino acids lining the central BMC-H pores. Current models attribute to BMC-H a homo-oligomeric nature. The hexagonal symmetry of corresponding pores, however, would break down if hetero-hexamers formed, a possibility suggested by the frequent presence of multiple paralogs within BMC operons. Here, we gauged the degree of structural promiscuity between the 11 BMC-H paralogs from Klebsiella pneumoniae , a potential human pathogen endowed with the capacity to express three different BMC types. Concomitant activation of transcription of several BMC operons was first shown to be possible. By leveraging an adapted tripartite GFP technology, all possible BMC-H pair combinations were screened in E. coli . Multiple structural cross-compatibilities were pinpointed between homologs arising not only from the same BMC operon, but also from different BMC types, results supported by Alphafold and ESMFold predictions. The structural stability and assembly propensity of selected hetero-associations was established by biochemical means. In light of these results, we reinterpreted published lysine cross-linking mass spectrometry data to demonstrate that one of these hetero-hexamers, involving PduA and PduJ, was already detected to form in the shell of a recombinantly-expressed 1,2-propanediol utilization compartment from Salmonella enterica . Altogether, this study points to the need to embrace an augmented structural complexity in BMC shells.

Proteins, the fundamental building blocks of biological function, orchestrate complex cellular processes by assembling into intricate structures through meticulous interactions. The design of specific protein-protein interfaces to create customized protein assemblies holds immense potential for various biotechnological applications. To address the current limitations in designing heteromeric interactions for multi-component assemblies, we developed a hybrid generative AI design approach. This method combines deep learning and automated reasoning, explicitly considering both positive and negative interaction states to favor heteromeric desired over undesired interactions. The approach leverages Effie, a deep-learned pairwise decomposable scoring function, and an advanced reasoning tool extended for multicriteria optimization of this function. Here, we tested the ability of this hybrid AI method to redesign homomeric interfaces of bacterial microcompartment components (BMC-H) into heteromeric assemblies. We benchmarked its performance against ProteinMPNN, a sequence design autoregressive model. Our in silico assessment, complemented by experimental validation, highlights the outperformance of the hybrid AI generative design approach, and its potential to unlock the engineering of complex multi-component self-assembling protein entities.

Terpenes are high-value chemicals which can be produced by engineered cyanobacteria from sustainable resources, solar energy, water and CO₂. We previously reported that the euryhaline unicellular cyanobacteria Synechocystis sp. PCC 6803 (S.6803) and Synechococcus sp. PCC 7002 (S.7002) produce farnesene and limonene, respectively, more efficiently than other terpenes. In the present study, we attempted to enhance farnesene production in S.6803 and limonene production in S.7002. Practically, we tested the influence of key cyanobacterial enzymes acting in carbon fixation (RubisCO, PRK, CcmK3 and CcmK4), utilization (CrtE, CrtR and CruF) and storage (PhaA and PhaB) on terpene production in S.6803, and we compared some of the findings with the data obtained in S.7002. We report that the overproduction of RubisCO from S.7002 and PRK from Cyanothece sp. PCC 7425 increased farnesene production in S.6803, but not limonene production in S.7002. The overexpression of the crtE genes (synthesis of terpene precursors) from S.6803 or S.7002 did not increase farnesene production in S.6803. In contrast, the overexpression of the crtE gene from S.6803, but not S.7002, increased farnesene production in S.7002, emphasizing the physiological difference between these two model cyanobacteria. Furthermore, the deletion of the crtR and cruF genes (carotenoid synthesis) and phaAB genes (carbon storage) did not increase the production of farnesene in S.6803. Finally, as a containment strategy of genetically modified strains of S.6803, we report that the deletion of the ccmK3K4 genes (carboxysome for CO₂ fixation) did not affect the production of limonene, but decreased the production of farnesene in S.6803.

Enzymatic parameters are classically determined in vitro , under conditions that are far from those encountered in cells, casting doubt on their physiological relevance. We developed a generic approach combining tools from synthetic and systems biology to measure enzymatic parameters in vivo . In the context of a synthetic carotenoid pathway in Saccharomyces cerevisiae , we focused on a phytoene synthase and three phytoene desaturases, which are difficult to study in vitro . We designed, built, and analyzed a collection of yeast strains mimicking substantial variations in substrate concentration by strategically manipulating the expression of geranyl-geranyl pyrophosphate (GGPP) synthase. We successfully determined in vivo Michaelis-Menten parameters ( K _M , V _max and k _cat ) for GGPP-converting phytoene synthase from absolute metabolomics, fluxomics and proteomics data, highlighting differences between in vivo and in vitro parameters. Leveraging the versatility of the same set of strains, we then extracted enzymatic parameters for two of the three phytoene desaturases. Our approach demonstrates the feasibility of assessing enzymatic parameters directly in vivo , providing a novel perspective on the kinetic characteristics of enzymes in real cellular conditions.

Bacterial micro-compartments (BMC) are complex macromolecular assemblies that participate in varied metabolic processes in about 20% of bacterial species. Most of these organisms carry BMC genetic information organized in operons that often include several paralog genes coding for components of the compartment shell. BMC shell constituents can be classified depending on their oligomerization state as hexamers (BMC-H), pentamers (BMC-P) or trimers (BMC-T). Formation of hetero-oligomers combining different protein homologs is theoretically feasible, something that could ultimately modify BMC shell rigidity or permeability, for instance. Despite that, it remains largely unknown whether hetero-oligomerization is a widespread phenomenon. Here, we demonstrated that the tripartite GFP (tGFP) reporter technology is an appropriate tool that might be exploited for such purposes. Thus, after optimizing parameters such as the size of linkers connecting investigated proteins to GFP10 or GFP11 peptides, the type and strength of promoters, or the impact of placing coding cassettes in the same or different plasmids, homo-oligomerization processes could be successfully monitored for any of the three BMC shell classes. Moreover, the screen perfectly reproduced published data on hetero-association between couples of CcmK homologues from Syn. sp. PCC6803, which were obtained following a different approach. This study paves the way for mid/high throughput screens to characterize the extent of hetero-oligomerization occurrence in BMC-possessing bacteria, and most especially in organisms endowed with several BMC types and carrying numerous shell paralogs. On the other hand, our study also unveiled technology limitations deriving from the low solubility of one of the components of this modified split-GFP approach, the GFP1-9.

Terpenes are high-value chemicals which can be produced by engineered cyanobacteria from sustainable resources, solar energy, water and CO2. We previously reported that the euryhaline unicellular cyanobacteria Synechocystis PCC 6803 (S.6803) and Synechococcus PCC 7002 (S.7002) produce farnesene and limonene, respectively, more efficiently than other terpenes. In the present study, we attempted at enhancing farnesene production in S.6803 and limonene production in S.7002. Practically, we tested the influence of key cyanobacterial enzymes acting in carbon fixation (RubisCO, PRK, CcmK3 and CcmK4), utilization (CrtE, CrtR and CruF) and storage (PhaA and PhaB) on terpene production in S.6803, and we compared some of the findings with the data obtained in S.7002. We report that the overproduction of RubisCO from S.7002 and PRK from Cyanothece PCC 7425 increased farnesene production in S.6803, but not limonene production in S.7002. The overexpression of the crtE genes (synthesis or terpene precursors) from S.6803 or S.7002 did not increased farnesene production in S.6803. In contrast, the overexpression of the crtE gene from S.6803, but not S.7002, increased farnesene production in S.7002, emphasizing the physiological difference between these two model cyanobacteria. Furthermore, the deletion of the crtR and cruF genes (carotenoid synthesis) and phaAB genes (carbon storage) did not increase the production of farnesene in S.6803. Finally, as a containment strategy of genetically modified (GMO) strains of S.6803 we report that the deletion of the ccmK3K4 genes (carboxysome for CO2 fixation) did not affect the production of limonene but decreased the production of farnesene in S.6803.

Bacterial microcompartments (BMC) are complex macromolecular assemblies that participate in varied chemical processes in about one fourth of bacterial species. BMC-encapsulated enzymatic activities are segregated from other cell contents by means of semipermeable shells, justifying why BMC are viewed as prototype nano-reactors for biotechnological applications. Herein, we undertook a comparative study of bending propensities of BMC hexamers (BMC-H), the most abundant shell constituents. Published data show that some BMC-H, like β-carboxysomal CcmK, tend to assemble flat whereas other BMC-H often build curved objects. Inspection of available crystal structures presenting BMC-H in tiled arrangements permitted us to identify two major assembly modes with a striking connection with experimental trends. All-atom molecular dynamics (MD) supported that BMC-H bending is triggered robustly only from the arrangement adopted in crystals by BMC-H that experimentally form curved objects, leading to very similar arrangements to those found in structures of recomposed BMC shells. Simulations on triplets of planar-behaving hexamers, which were previously reconfigured to comply with such organization, confirmed that bending propensity is mostly defined by the precise lateral positioning of hexamers, rather than by BMC-H identity. Finally, an interfacial lysine was pinpointed as the most decisive residue in controlling PduA spontaneous curvature. Globally, results presented herein should contribute to improve our understanding of the variable mechanisms of biogenesis characterized for BMC, and of possible strategies to regulate BMC size and shape.