Computational studies on the catalytic potential of the double active site for enzyme engineering.

Proteins possessing double active sites have the potential to revolutionise enzyme design strategies. This study extensively explored an enzyme that contains both a natural active site (NAS) and an engineered active site (EAS), focusing on understanding its structural and functional properties. Metadynamics simulations were employed to investigate how substrates interacted with their respective active sites. The results revealed that both the NAS and EAS exhibited similar minimum energy states, indicating comparable binding affinities. However, it became apparent that the EAS had a weaker binding site for the substrate due to its smaller pocket and constrained conformation. Interestingly, the EAS also displayed dynamic behaviour, with the substrate observed to move outside the pocket, suggesting the possibility of substrate translocation. To gain further insights, steered molecular dynamics (SMD) simulations were conducted to study the conformational changes of the substrate and its interactions with catalytic residues. Notably, the substrate adopted distinct conformations, including near-attack conformations, in both the EAS and NAS. Nevertheless, the NAS demonstrated superior binding minima for the substrate compared to the EAS, reinforcing the observation that the engineered active site was less favourable for substrate binding due to its limitations. The QM/MM (Quantum mechanics and molecular mechanics) analyses highlight the energy disparity between NAS and EAS. Specifically, EAS exhibited elevated energy levels due to its engineered active site being located on the surface. This positioning exposes the substrate to solvents and water molecules, adding to the energy challenge. Consequently, the engineered enzyme did not provide a significant advantage in substrate binding over the single active site protein. Further, the investigation of internal channels and tunnels within the protein shed light on the pathways facilitating transport between the two active sites. By unravelling the complex dynamics and functional characteristics of this double-active site protein, this study offers valuable insights into novel strategies of enzyme engineering. These findings establish a solid foundation for future research endeavours aimed at harnessing the potential of double-active site proteins in diverse biotechnological applications.

Towards universal synthetic heterotrophy using a metabolic coordinator.

Engineering the utilization of non-native substrates, or synthetic heterotrophy, in proven industrial microbes such as Saccharomyces cerevisiae represents an opportunity to valorize plentiful and renewable sources of carbon and energy as inputs to bioprocesses. We previously demonstrated that activation of the galactose (GAL) regulon, a regulatory structure used by this yeast to coordinate substrate utilization with biomass formation during growth on galactose, during growth on the non-native substrate xylose results in a vastly altered gene expression profile and faster growth compared with constitutive overexpression of the same heterologous catabolic pathway. However, this effort involved the creation of a xylose-inducible variant of Gal3p (Gal3pSyn4.1), the sensor protein of the GAL regulon, preventing this semi-synthetic regulon approach from being easily adapted to additional non-native substrates. Here, we report the construction of a variant Gal3pMC (metabolic coordinator) that exhibits robust GAL regulon activation in the presence of structurally diverse substrates and recapitulates the dynamics of the native system. Multiple molecular modeling studies suggest that Gal3pMC occupies conformational states corresponding to galactose-bound Gal3p in an inducer-independent manner. Using Gal3pMC to test a regulon approach to the assimilation of the non-native lignocellulosic sugars xylose, arabinose, and cellobiose yields higher growth rates and final cell densities when compared with a constitutive overexpression of the same set of catabolic genes. The subsequent demonstration of rapid and complete co-utilization of all three non-native substrates suggests that Gal3pMC-mediated dynamic global gene expression changes by GAL regulon activation may be universally beneficial for engineering synthetic heterotrophy.

In-depth Sequence-Function Characterization Reveals Multiple Pathways to Enhance Enzymatic Activity.

Deep mutational scanning (DMS) has recently emerged as a powerful method to study protein sequence-function relationships but it has not been well-explored as a guide to enzyme engineering and identifying pathways by which their catalytic cycle may be improved. We report such a demonstration in this work using a Phenylalanine ammonia-lyase (PAL), which deaminates L-phenylalanine to trans-cinnamic acid and has widespread application in chemo-enzymatic synthesis, agriculture, and medicine. In particular, the PAL from Anabaena variabilis (AvPAL*) has garnered significant attention as the active ingredient in Pegvaliase®, the only FDA-approved drug treating classical Phenylketonuria (PKU). Although an extensive body of literature exists on the structure, substrate-specificity, and catalytic cycle, protein-wide sequence determinants of function remain unknown, as do intermediate reaction steps that limit turnover frequency, all of which has hindered rational engineering of these enzymes. Here, we created a detailed sequence-function landscape of AvPAL* by performing DMS and revealed 112 mutations at 79 functionally relevant sites that affect a positive change in enzyme fitness. Using fitness values and structure-function analysis, we picked a subset of positions for comprehensive single- and multi-site saturation mutagenesis and identified combinations of mutations that led to improved reaction kinetics in cell-free and cellular contexts. We then performed QM/MM and MD to understand the mechanistic role of the most beneficial mutations and observed that different mutants confer improvements via different mechanisms, including stabilizing transition and intermediate states, improving substrate diffusion into the active site, and decreasing product inhibition. This work demonstrates how DMS can be combined with computational analysis to effectively identify significant mutations that enhance enzyme activity along with the underlying mechanisms by which these mutations confer their benefit.

Connecting the dots in the mechanism of action of Cucurbitacin E (CurE) - path analysis and steered molecular dynamics reveal the precise site of entry and the passage of CurE in filamentous actin.

Cucurbitacin E (CurE) modulates actin cytoskeleton by forming an irreversible covalent bond with Cys257 of actin. The reported binding conformation of CurE is deeply buried in the subdomain 4 of actin and is closely situated to the ATP-binding site. The entry and the path taken by CurE to reach this buried site remain a mystery. In this study, steered molecular dynamics (SMD) simulations were conducted to delineate the diffusion of CurE to its binding site. SMD simulations reveal that the distinctive entry site of CurE found in subdomain 4 is by itself a closed and compact region of two loops lying beside each other like a closed door and CurE induces it to open. From this point, CurE moves toward its binding site through a path facilitated by Thr188, Leu261, Ile267, Ile309, Tyr306. This study is also an insight into how CurE distinctly differentiates its tunnel to the binding site from the ATP-binding site. The conformational changes of CurE along the path to the binding site are surprisingly very minimal and closely resemble the attack conformation at the end of simulation. The study reveals that the little energy spent by the molecule is compensated by the enthalpic contribution to binding-free energy barrier making it undoubtedly the most preferred path of CurE. This study is the first of its kind in which the SMD was used to derive the complete and continuous translocation of one of the most potent phytochemicals, CurE through the binding site gorge of actin. [Formula: see text] Communicated by Ramaswamy H. Sarma.

Azadirachtin(A) distinctively modulates subdomain 2 of actin - novel mechanism to induce depolymerization revealed by molecular dynamics study.

Azadirachtin(A) (AZA), a potential insecticide from neem, binds to actin and induces depolymerization in Drosophila. AZA binds to the pocket same as that of Latrunculin A (LAT), but LAT inhibits actin polymerization by stiffening the actin structure and affects the ADP-ATP exchange. The mechanism by which AZA induces actin depolymerization is not clearly understood. Therefore, different computational experiments were conducted to delineate the precise mechanism of AZA-induced actin depolymerization. Molecular dynamics studies showed that AZA strongly interacted with subdomain 2 and destabilized the interactions between subdomain 2 of one actin and subdomains 1 and 4 of the adjacent actin, causing the separation of actin subunits. The separation was observed between subdomain 3 of subunit n and subdomain 4 of subunit n + 2. However, the specific triggering point for the separation of the subunits was the destabilization of direct interactions between subdomain 2 of subunit n (Arg39, Val45, Gly46 and Arg62) and subdomain 4 of subunit n + 2 (Asp286, Ile287, Asp288, Ile289, Asp244 and Lys291). These results reveal a unique mechanism of an actin filament modulator that induces depolymerization. This mechanism of AZA can be used to design similar molecules against mammalian actins for cancer therapy.

In silico approach of azadirachtin binding with actins.

In silico docking analysis reported here suggests that insect cellular cytoskeletal beta-actin could be the target of Azadirachtin A (Aza-the principle bioactive compound of neem seeds). The best docking energy of -40.09 kcal/mol at 8.73 A RMSD and predicted hydrogen bond between Arg210 and carboxymethyl group of Aza accompanied with seven hydrophobic interactions in the proposed binding site strongly support this hypothesis. This is of specific interest due to the non-affinity of Aza to mammalian beta-actins under the same set of conditions, although beta-actins across the species are highly conserved. Our results show that few scattered amino acid changes have caused significant steric hindrance in the binding pocket for mammalian beta-actin, causing a reverse orientation of Aza. These results suggest a model to support the recently observed biological effects caused by Aza in Drosophila cytoskeletal elements and explain why Aza is highly specific to insects than mammals.