Effects of Distal Mutations on the Structure, Dynamics and Catalysis of Human Monoacylglycerol Lipase

Recent advances in structure-based enzyme engineering for functional reconstruction

Structural information can help engineer enzymes. Usually, specific amino acids in particular regions are targeted for functional reconstruction to enhance the catalytic performance, including activity, stereoselectivity, and thermostability. Appropriate selection of target sites is the key to structure-based design, which requires elucidation of the structure-function relationships. Here, we summarize the mutations of residues in different specific regions, including active center, access tunnels, and flexible loops, on fine-tuning the catalytic performance of enzymes, and discuss the effects of altering the local structural environment on the functions. In addition, we keep up with the recent progress of structure-based approaches for enzyme engineering, aiming to provide some guidance on how to take advantage of the structural information.

Structural dynamics and mechanistic action guided engineering of lipolytic enzymes

Lipases have been established as important biocatalysts in several industrial applications, owing to their diverse substrate specificity. The availability of data on three-dimensional crystal structures for various lipases offers an opportunity for modulating their structural and functional aspects to design and engineer better versions of lipases. With the aim of investigating the structural components governing the extremophilic behavior of lipases, structural analysis of microbial lipases was performed using advanced bioinformatics and molecular dynamics simulation approaches. In sequences and functionally distinct alkaliphilic and thermophilic lipases were investigated for their functional properties to understand the distinguishing features of their structures. The alkaliphilic lipase from Bacillus subtilis (LipA) showed conformational changes in the loop region Ala132-Met137, subsequently, the active site residue His156 shows two conformations, toward the active site nucleophilic residues Ser77 and away from the Ser77. Interestingly, the active site of LipA is more solvent-exposed and can be correlated with the adoption of an open conformation which might extend and expose the active site region to solvents during the catalysis process. Furthermore, the MD simulation of thermophilic lipase from marine Streptomyces (MAS1) revealed the role of N- and C-terminal regions with disulfide bridges and identified a metal ion binding site that facilitates the enzyme stability. The novel thermo-alkaliphilic lipase can be designed to integrate the stability features of MAS1 into the alkaliphilic LipA. These structural-level intrinsic characteristics can be used for lipase engineering to amend the lipase activity and stability as per the requirements of the industrial processes.

Non-active Site Residue in Loop L4 Alters Substrate Capture and Product Release in d-Arginine Dehydrogenase

Numerous studies demonstrate that enzymes undergo multiple conformational changes during catalysis. The malleability of enzymes forms the basis for allosteric regulation: residues located far from the active site can exert long-range dynamical effects on the active site residues to modulate catalysis. The structure of Pseudomonas aeruginosa d-arginine dehydrogenase (PaDADH) shows four loops (L1, L2, L3, and L4) that span the substrate and the FAD-binding domains. Loop L4 comprises residues 329-336, spanning over the flavin cofactor. The I335 residue on loop L4 is ∼10 Å away from the active site and ∼3.8 Å from N(1)-C(2)═O atoms of the flavin. In this study, we used molecular dynamics and biochemical techniques to investigate the effect of the mutation of I335 to histidine on the catalytic function of PaDADH. Molecular dynamics showed that the conformational dynamics of PaDADH are shifted to a more closed conformation in the I335H variant. In agreement with an enzyme that samples more in a closed conformation, the kinetic data of the I335H variant showed a 40-fold decrease in the rate constant of substrate association (k₁), a 340-fold reduction in the rate constant of substrate dissociation from the enzyme-substrate complex (k₂), and a 24-fold decrease in the rate constant of product release (k₅), compared to that of the wild-type. Surprisingly, the kinetic data are consistent with the mutation having a negligible effect on the reactivity of the flavin. Altogether, the data indicate that the residue at position 335 has a long-range dynamical effect on the catalytic function in PaDADH.

Structural features, temperature adaptation and industrial applications of microbial lipases from psychrophilic, mesophilic and thermophilic origins

Microbial lipases are very prominent biocatalysts because of their ability to catalyze a wide variety of reactions in aqueous and non-aqueous media. Here microbial lipases from different origins (psychrophiles, mesophiles, and thermophiles) have been reviewed. This review emphasizes an update of structural diversity in temperature adaptation and industrial applications, of psychrophilic, mesophilic, and thermophilic lipases. The microbial origins of lipases are logically dynamic, proficient, and also have an extensive range of industrial uses with the manufacturing of altered molecules. It is therefore of interest to understand the molecular mechanisms of adaptation to temperature in occurring lipases. However, lipases from extremophiles (psychrophiles, and thermophiles) are widely used to design biotransformation reactions with higher yields, fewer byproducts, or useful side products and have been predicted to catalyze those reactions also, which otherwise are not possible with the mesophilic lipases. Lipases as a multipurpose biological catalyst have given a favorable vision in meeting the needs of several industries such as biodiesel, foods, and drinks, leather, textile, detergents, pharmaceuticals, and medicals.

Leveraging intrinsic flexibility to engineer enhanced enzyme catalytic activity

Dynamic motions of enzymes occurring on a broad range of timescales play a pivotal role in all steps of the reaction pathway, including substrate binding, catalysis, and product release. However, it is unknown whether structural information related to conformational flexibility can be exploited for the directed evolution of enzymes with higher catalytic activity. Here, we show that mutagenesis of residues exclusively located at flexible regions distal to the active site of Homo sapiens kynureninase (HsKYNase) resulted in the isolation of a variant (BF-HsKYNase) in which the rate of the chemical step toward kynurenine was increased by 45-fold. Mechanistic pre–steady-state kinetic analysis of the wild type and the evolved enzyme shed light on the underlying effects of distal mutations (>10 Å from the active site) on the rate-limiting step of the catalytic cycle. Hydrogen-deuterium exchange coupled to mass spectrometry and molecular dynamics simulations revealed that the amino acid substitutions in BF-HsKYNase allosterically affect the flexibility of the pyridoxal-5′-phosphate (PLP) binding pocket, thereby impacting the rate of chemistry, presumably by altering the conformational ensemble and sampling states more favorable to the catalyzed reaction.

Characterizing and explaining the impact of disease-associated mutations in proteins without known structures or structural homologs

Mutations in human proteins lead to diseases. The structure of these proteins can help understand the mechanism of such diseases and develop therapeutics against them. With improved deep learning techniques, such as RoseTTAFold and AlphaFold, we can predict the structure of proteins even in the absence of structural homologs. We modeled and extracted the domains from 553 disease-associated human proteins without known protein structures or close homologs in the Protein Databank. We noticed that the model quality was higher and the Root mean square deviation (RMSD) lower between AlphaFold and RoseTTAFold models for domains that could be assigned to CATH families as compared to those which could only be assigned to Pfam families of unknown structure or could not be assigned to either. We predicted ligand-binding sites, protein–protein interfaces and conserved residues in these predicted structures. We then explored whether the disease-associated missense mutations were in the proximity of these predicted functional sites, whether they destabilized the protein structure based on ddG calculations or whether they were predicted to be pathogenic. We could explain 80% of these disease-associated mutations based on proximity to functional sites, structural destabilization or pathogenicity. When compared to polymorphisms, a larger percentage of disease-associated missense mutations were buried, closer to predicted functional sites, predicted as destabilizing and pathogenic. Usage of models from the two state-of-the-art techniques provide better confidence in our predictions, and we explain 93 additional mutations based on RoseTTAFold models which could not be explained based solely on AlphaFold models.

Mechanistic Modeling of Monoglyceride Lipase Covalent Modification Elucidates the Role of Leaving Group Expulsion and Discriminates Inhibitors with High and Low Potency

Inhibition of monoglyceride lipase (MGL), also known as monoacylglycerol lipase (MAGL), has emerged as a promising approach for treating neurological diseases. To gain useful insights in the design of agents with balanced potency and reactivity, we investigated the mechanism of MGL carbamoylation by the reference triazole urea SAR629 (IC₅₀ = 0.2 nM) and two recently described inhibitors featuring a pyrazole (IC₅₀ = 1800 nM) or a 4-cyanopyrazole (IC₅₀ = 8 nM) leaving group (LG), using a hybrid quantum mechanics/molecular mechanics (QM/MM) approach. Opposite to what was found for substrate 2-arachidonoyl-sn-glycerol (2-AG), covalent modification of MGL by azole ureas is controlled by LG expulsion. Simulations indicated that changes in the electronic structure of the LG greatly affect reaction energetics with triazole and 4-cyanopyrazole inhibitors following a more accessible carbamoylation path compared to the unsubstituted pyrazole derivative. The computational protocol provided reaction barriers able to discriminate between MGL inhibitors with different potencies. These results highlight how QM/MM simulations can contribute to elucidating structure–activity relationships and provide insights for the design of covalent inhibitors.

Inhibition of triple negative breast cancer-associated inflammation, tumor growth and brain colonization by targeting monoacylglycerol lipase

While the prevalence of breast cancer metastasis in the brain is significantly higher in triple negative breast cancers (TNBCs), there is a lack of novel and/or improved therapies for these patients. Monoacylglycerol lipase (MAGL) is a hydrolase involved in lipid metabolism that catalyzes the degradation of 2-arachidonoylglycerol (2-AG) linked to generation of pro- and anti-inflammatory molecules. Here, we targeted MAGL in TNBCs, using a potent carbamate-based inhibitor AM9928 (hMAGL IC₅₀ = 9 nM) with prolonged pharmacodynamic effects (46 h of target residence time). AM9928 blocked TNBC cell adhesion and transmigration across human brain microvascular endothelial cells (HBMECs) in 3D co-cultures. In addition, AM9928 inhibited the secretion of IL-6, IL-8, and VEGF-A from TNBC cells. TNBC-derived exosomes activated HBMECs resulting in secretion of elevated levels of IL-8 and VEGF, which were inhibited by AM9928. Using in vivo studies of syngeneic GFP-4T1-BrM5 mammary tumor cells, AM9928 inhibited tumor growth in the mammary fat pads and attenuated blood brain barrier (BBB) permeability changes, resulting in reduced TNBC colonization in brain. Together, these results support the potential clinical application of MAGL inhibitors as novel treatments for TNBC.

A review on structurally diversified synthesized molecules as monoacylglycerol lipase inhibitors and their therapeutic uses

Monoacylglycerol is a metabolic key serine hydrolase, engaged in the regulation of signalling network system of endocannabinoids, which is associated with various physiological processes like pain, inflammation, feeding cognition and neurodegenerative diseases like Alzheimer, Parkinson's disease. The monoacylglycerol also found to act as a regulator and the free fatty acid provider in the proliferation of cancer cells, numerous aggressive tumours such as colorectal cancer, neuroblastoma and nasopharyngeal carcinoma. It also played an important role in increasing the concentration of specific lipids derived from free fatty acids like phosphatidic acid, lysophosphatidic acid, sphingosine-1-phosphate and prostaglandin E2. These signalling lipids are associated with cell proliferation, survival, tumour cell migration, contributing to tumour development, maturation and metastases. In the present study here, we are presenting a review on structurally diverse MAGL inhibitors, their development and their evaluation for different pharmacological activities.

Importance of protein intrinsic conformational dynamics and transient nature of non-covalent interactions in ligand binding affinity

We have recently identified BEN1 as a protein interactor of seryl-tRNA synthetase (SerRS) from model plant Arabidopsis thaliana. BEN1 contains an NADP⁺ binding domain and possesses acidic N-terminal extension essential for interaction with A. thaliana SerRS. This extension, specific for BEN1 homologues from Brassicaceae family, is solvent-exposed and distant to the nucleotide-binding site. We prepared a truncated BEN1 variant ΔN17BEN1 lacking the first 17 amino acid of this N-terminal extension as well as full-length BEN1 to investigate how the truncation affects the binding affinity towards coenzyme NADP⁺. By performing microscale thermophoresis (MST) experiments we have shown that both BEN1 variants bind the NADP⁺ cofactor, however, truncated BEN1 showed 34-fold higher affinity towards NADP⁺ indicating that its core protein structure is not just preserved but it binds NADP⁺ even stronger. To further corroborate the obtained results, we opted for a computational approach based on classical molecular dynamics simulations of both complexes. Our results have shown that both truncated and intact BEN1 variants form the same number of interactions with the NADP⁺ cofactor; however, it was the interaction occupancy that was affected. Namely, three independent MD simulations showed that the ΔN17BEN1 variant in complex with NADP⁺ has significantly higher interaction occupancy thus binds NADP⁺ with more than one order of magnitude higher affinity. Contrary to our expectations, the truncation of this distant region that does not communicate with the nucleotide-binding site didn't result in the gain of interaction but affected the intrinsic conformational dynamics which in turn fine-tuned the binding affinity by increasing the interaction occupancy and strength of the key conserved cation-π interaction between Arg69 and adenine of NADP⁺ and hydrogen bond between Ser244 and phosphate of NADP⁺.

Tyrp1 Mutant Variants Associated with OCA3: Computational Characterization of Protein Stability and Ligand Binding

Oculocutaneous albinism type 3 (OCA3) is an autosomal recessive disorder caused by mutations in the TYRP1 gene. Tyrosinase-related protein 1 (Tyrp1) is involved in eumelanin synthesis, catalyzing the oxidation of 5,6-dihydroxyindole-2-carboxylic acid oxidase (DHICA) to 5,6-indolequinone-2-carboxylic acid (IQCA). Here, for the first time, four OCA3-causing mutations of Tyrp1, C30R, H215Y, D308N, and R326H, were investigated computationally to understand Tyrp1 protein stability and catalytic activity. Using the Tyrp1 crystal structure (PDB:5M8L), global mutagenesis was conducted to evaluate mutant protein stability. Consistent with the foldability parameter, C30R and H215Y should exhibit greater instability, and two other mutants, D308N and R326H, are expected to keep a native conformation. SDS-PAGE and Western blot analysis of the purified recombinant proteins confirmed that the foldability parameter correctly predicted the effect of mutations critical for protein stability. Further, the mutant variant structures were built and simulated for 100 ns to generate free energy landscapes and perform docking experiments. Free energy landscapes formed by Y362, N378, and T391 indicate that the binding clefts of C30R and H215Y mutants are larger than the wild-type Tyrp1. In docking simulations, the hydrogen bond and salt bridge interactions that stabilize DHICA in the active site remain similar among Tyrp1, D308N, and R326H. However, the strengths of these interactions and stability of the docked ligand may decrease proportionally to mutation severity due to the larger and less well-defined natures of the binding clefts in mutants. Mutational perturbations in mutants that are not unfolded may result in allosteric alterations to the active site, reducing the stability of protein-ligand interactions.

Characterization of Temperature-Dependent Kinetics of Oculocutaneous Albinism-Causing Mutants of Tyrosinase

Human tyrosinase (Tyr) is a glycoenzyme that catalyzes the first and rate-limiting step in melanin production, and its gene (TYR) is mutated in many cases of oculocutaneous albinism type 1 (OCA1). The mechanisms by which individual mutations contribute to the diverse pigmentation phenotype in patients with OCA1 have only began to be examined and remain to be delineated. Here, we analyze the temperature-dependent kinetics of wild-type Tyr (WT) and two OCA1B mutant variants (R422Q and P406L) using Michaelis–Menten and Van’t Hoff analyses. Recombinant truncated human Tyr proteins (residues 19–469) were produced in the whole insect Trichoplusia Ni larvae. Proteins were purified by a combination of affinity and size-exclusion chromatography. The temperature dependence of diphenol oxidase protein activities and kinetic parameters were measured by dopachrome absorption. Using the same experimental conditions, computational simulations were performed to assess the temperature-dependent association of L-DOPA and Tyr. Our results revealed, for the first time, that the association of L-DOPA with R422Q and P406L followed by dopachrome formation is a complex reaction supported by enthalpy and entropy forces. We show that the WT has a higher turnover number as compared with both R422Q and P406L. Elucidating the kinetics and thermodynamics of mutant variants of Tyr in OCA1B helps to understand the mechanisms by which they lower Tyr catalytic activity and to discover novel therapies for patients.

Methemoglobin formation in mutant hemoglobin α chains: electron transfer parameters and rates

Hemoglobin-mediated transport of dioxygen (O₂) critically depends on the stability of the reduced (Fe²⁺) form of the heme cofactors. Some protein mutations stabilize the oxidized (Fe³⁺) state (methemoglobin, Hb M), causing methemoglobinemia, and can be lethal above 30%. The majority of the analyses of factors influencing Hb oxidation are retrospective and give insights only for inner-sphere mutations of heme (His58, His87). Herein, we report the first all-atom molecular dynamics simulations on both redox states and calculations of the Marcus electron transfer (ET) parameters for the α chain Hb oxidation and reduction rates for Hb M. The Hb wild-type (WT) and most of the studied α chain variants maintain globin structure except the Hb M Iwate (H87Y). The mutants forming Hb M tend to have lower redox potentials and thus stabilize the oxidized (Fe³⁺) state (in particular, the Hb Miyagi variant with K61E mutation). Solvent reorganization (λ_solv 73-96%) makes major contributions to reorganization free energy, whereas protein reorganization (λ_prot) accounts for 27-30% except for the Miyagi and J-Buda variants (λ_prot ∼4%). Analysis of heme-solvent H-bonding interactions among variants provide insights into the role of Lys61 residue in stabilizing the Fe²⁺ state. Semiclassical Marcus ET theory-based calculations predict experimental k_ET for the Cyt b5-Hb complex and provide insights into relative reduction rates for Hb M in Hb variants. Thus, our methodology provides a rationale for the effect of mutations on the structure, stability, and Hb oxidation reduction rates and has potential for identification of mutations that result in methemoglobinemia.

D3DistalMutation: a Database to Explore the Effect of Distal Mutations on Enzyme Activity

Enzyme activity is affected by amino acid mutations, particularly mutations near the active site. Increasing evidence has shown that distal mutations more than 10 Å away from the active site may significantly affect enzyme activity. However, it is difficult to study the enzyme regulation mechanism of distal mutations due to the lack of a systematic collection of three-dimensional (3D) structures, highlighting distal mutation site and the corresponding enzyme activity change. Therefore, we constructed a distal mutation database, namely, D3DistalMutation, which relates the distal mutation to enzyme activity. As a result, we observed that approximately 80% of distal mutations could affect enzyme activity and 72.7% of distal mutations would decrease or abolish enzyme activity in D3DistalMutation. Only 6.6% of distal mutations in D3DistalMutation could increase enzyme activity, which have great potential to the industrial field. Among these mutations, the Y to F, S to D, and T to D mutations are most likely to increase enzyme activity, which sheds some light on industrial catalysis. Distal mutations decreasing enzyme activity in the allosteric pocket play an indispensable role in allosteric drug design. In addition, the pockets in the enzyme structures are provided to explore the enzyme regulation mechanism of distal mutations. D3DistalMutation is accessible free of charge at https://www.d3pharma.com/D3DistalMutation/index.php.

Conformational gating, dynamics and allostery in human monoacylglycerol lipase

Inhibition of human Monoacylglycerol Lipase (hMGL) offers a novel approach for treating neurological diseases. The design of inhibitors, targeting active-inactive conformational transitions of the enzyme, can be aided by understanding the interplay between structure and dynamics. Here, we report the effects of mutations within the catalytic triad on structure, conformational gating and dynamics of hMGL by combining kinetics, NMR, and HDX-MS data with metadynamics simulations. We found that point mutations alter delicate conformational equilibria between active and inactive states. HDX-MS reveals regions of the hMGL that become substantially more dynamic upon substitution of catalytic acid Asp-239 by alanine. These regions, located far from the catalytic triad, include not only loops but also rigid α-helixes and β-strands, suggesting their involvement in allosteric regulation as channels for long-range signal transmission. The results identify the existence of a preorganized global communication network comprising of tertiary (residue-residue contacts) and quaternary (rigid-body contacts) networks that mediate robust, rapid intraprotein signal transmission. Catalytic Asp-239 controls hMGL allosteric communications and may be considered as an essential residue for the integration and transmission of information to enzymes' remote regions, in addition to its well-known role to facilitate Ser-122 activation. Our findings may assist in the identification of new druggable sites in hMGL.

Modeling the Role of a Flexible Loop and Active Site Side Chains in Hydride Transfer Catalyzed by Glycerol-3-phosphate Dehydrogenase

Glycerol-3-phosphate dehydrogenase is a biomedically important enzyme that plays a crucial role in lipid biosynthesis. It is activated by a ligand-gated conformational change that is necessary for the enzyme to reach a catalytically competent conformation capable of efficient transition-state stabilization. While the human form (hlGPDH) has been the subject of extensive structural and biochemical studies, corresponding computational studies to support and extend experimental observations have been lacking. We perform here detailed empirical valence bond and Hamiltonian replica exchange molecular dynamics simulations of wild-type hlGPDH and its variants, as well as providing a crystal structure of the binary hlGPDH·NAD R269A variant where the enzyme is present in the open conformation. We estimated the activation free energies for the hydride transfer reaction in wild-type and substituted hlGPDH and investigated the effect of mutations on catalysis from a detailed structural study. In particular, the K120A and R269A variants increase both the volume and solvent exposure of the active site, with concomitant loss of catalytic activity. In addition, the R269 side chain interacts with both the Q295 side chain on the catalytic loop, and the substrate phosphodianion. Our structural data and simulations illustrate the critical role of this side chain in facilitating the closure of hlGPDH into a catalytically competent conformation, through modulating the flexibility of a key catalytic loop (292-LNGQKL-297). This, in turn, rationalizes a tremendous 41,000 fold decrease experimentally in the turnover number, k_cat, upon truncating this residue, as loop closure is essential for both correct positioning of key catalytic residues in the active site, as well as sequestering the active site from the solvent. Taken together, our data highlight the importance of this ligand-gated conformational change in catalysis, a feature that can be exploited both for protein engineering and for the design of allosteric inhibitors targeting this biomedically important enzyme.

Effects of Distal Mutations on Prolyl-Adenylate Formation of Escherichia coli Prolyl-tRNA Synthetase

Enzymes play important roles in many biological processes. Amino acid residues in the active site pocket of an enzyme, which are in direct contact with the substrate(s), are generally believed to be critical for substrate recognition and catalysis. Identifying and understanding how these "catalytic" residues help enzymes achieve enormous rate enhancement has been the focus of many structural and biochemical studies over the past several decades. Recent studies have shown that enzymes are intrinsically dynamic and dynamic coupling between distant structural elements is essential for effective catalysis in modular enzymes. Therefore, distal residues are expected to have impact on enzyme function. However, few studies have investigated the role of distal residues on enzymatic catalysis. In the present study, the effects of distal residue mutations on the catalytic function of an aminoacyl-tRNA synthetase, namely, prolyl-tRNA synthase, were investigated. The present study demonstrates that distal residues significantly contribute to catalysis of the modular Escherichia coli prolyl-tRNA synthetase by maintaining intrinsic protein flexibility.

A beginner's guide to molecular dynamics simulations and the identification of cross-correlation networks for enzyme engineering

The functional properties of proteins are decided not only by their relatively rigid overall structures, but even more importantly, by their dynamic properties. In a protein, some regions of structure exhibit highly correlated or anti-correlated motions with others, some are highly dynamic but uncorrelated, while other regions are relatively static. The residues with correlated or anti-correlated motions can form a so-called dynamic cross-correlation network, through which information can be transmitted. Such networks have been shown to be critical to allosteric transitions, and ligand binding, and have also been shown to be able to mediate epistatic interactions between mutations. As a result, they are likely to play a significant role in the development of new enzyme engineering strategies. In this chapter, protocols are provided for the assessment of dynamic cross-correlation networks, and for their application in protein engineering. Transketolase from E. coli is used as a model and the software GROMACS is applied for carrying out MD simulations to generate trajectories containing structural ensembles. The trajectory is then used for a dynamic cross correlation analysis using the R package, Bio3D. A matrix of all atom-wise cross-correlation coefficients is finally obtained, which can be displayed in a graphical representation termed a dynamical cross-correlation matrix.

A Thermolabile Phospholipase B from Talaromyces marneffei GD-0079: Biochemical Characterization and Structure Dynamics Study

Phospholipase B (EC 3.1.1.5) are a distinctive group of enzymes that catalyzes the hydrolysis of fatty acids esterified at the sn-1 and sn-2 positions forming free fatty acids and lysophospholipids. The structural information and catalytic mechanism of phospholipase B are still not clear. Herein, we reported a putative phospholipase B (TmPLB1) from Talaromyces marneffei GD-0079 synthesized by genome mining library. The gene (TmPlb1) was expressed and the TmPLB1 was purified using E. coli shuffle T7 expression system. The putative TmPLB1 was purified by affinity chromatography with a yield of 13.5%. The TmPLB1 showed optimum activity at 35 °C and pH 7.0. The TmPLB1 showed enzymatic activity using Lecithin (soybean > 98% pure), and the hydrolysis of TmPLB1 by ³¹P NMR showed phosphatidylcholine (PC) as a major phospholipid along with lyso-phospholipids (1-LPC and 2-LPC) and some minor phospholipids. The molecular modeling studies indicate that its active site pocket contains Ser125, Asp183 and His215 as the catalytic triad. The structure dynamics and simulations results explained the conformational changes associated with different environmental conditions. This is the first report on biochemical characterization and structure dynamics of TmPLB1 enzyme. The present study could be helpful to utilize TmPLB1 in food industry for the determination of food components containing phosphorus. Additionally, such enzyme could also be useful in Industry for the modifications of phospholipids.

Structural and free energy landscape of novel mutations in ribosomal protein S1 (rpsA) associated with pyrazinamide resistance

Resistance to key first-line drugs is a major hurdle to achieve the global end tuberculosis (TB) targets. A prodrug, pyrazinamide (PZA) is the only drug, effective in latent TB, recommended in drug resistance and susceptible Mycobacterium tuberculosis (MTB) isolates. The prodrug conversion into active form, pyrazinoic acid (POA), required the activity of pncA gene encoded pyrazinamidase (PZase). Although pncA mutations have been commonly associated with PZA resistance but a small number of resistance cases have been associated with mutationss in RpsA protein. Here in this study a total of 69 PZA resistance isolates have been sequenced for pncA mutations. However, samples that were found PZA resistant but pncA wild type (pncAWT), have been sequenced for rpsA and panD genes mutation. We repeated a drug susceptibility testing according to the WHO guidelines on 18 pncAWT MTB isolates. The rpsA and panD genes were sequenced. Out of total 69 PZA resistant isolates, 51 harbored 36 mutations in pncA gene (GeneBank Accession No. MH46111) while, fifteen different mutations including seven novel, were detected in the fourth S1 domain of RpsA known as C-terminal (MtRpsACTD) end. We did not detect any mutations in panD gene. Among the rpsA mutations, we investigated the molecular mechanism of resistance behind mutations, D342N, D343N, A344P, and I351F, present in the MtRpsACTD through molecular dynamic simulations (MD). WT showed a good drug binding affinity as compared to mutants (MTs), D342N, D343N, A344P, and I351F. Binding pocket volume, stability, and fluctuations have been altered whereas the total energy, protein folding, and geometric shape analysis further explored a significant variation between WT and MTs. In conclusion, mutations in MtRpsACTD might be involved to alter the RpsA activity, resulting in drug resistance. Such molecular mechanism behind resistance may provide a better insight into the resistance mechanism to achieve the global TB control targets.

Protein engineering: the potential of remote mutations

Engineered proteins, especially enzymes, are now commonly used in many industries owing to their catalytic power, specific binding of ligands, and properties as materials and food additives. As the number of potential uses for engineered proteins has increased, the interest in engineering or designing proteins to have greater stability, activity and specificity has increased in turn. With any rational engineering or design pursuit, the success of these endeavours relies on our fundamental understanding of the systems themselves; in the case of proteins, their structure-dynamics-function relationships. Proteins are most commonly rationally engineered by targeting the residues that we understand to be functionally important, such as enzyme active sites or ligand-binding sites. This means that the majority of the protein, i.e. regions remote from the active- or ligand-binding site, is often ignored. However, there is a growing body of literature that reports on, and rationalises, the successful engineering of proteins at remote sites. This minireview will discuss the current state of the art in protein engineering, with a particular focus on engineering regions that are remote from active- or ligand-binding sites. As the use of protein technologies expands, exploiting the potential improvements made possible through modifying remote regions will become vital if we are to realise the full potential of protein engineering and design.

Structural Perspective on Revealing and Altering Molecular Functions of Genetic Variants Linked with Diseases

Structural information of biological macromolecules is crucial and necessary to deliver predictions about the effects of mutations-whether polymorphic or deleterious (i.e., disease causing), wherein, thermodynamic parameters, namely, folding and binding free energies potentially serve as effective biomarkers. It may be emphasized that the effect of a mutation depends on various factors, including the type of protein (globular, membrane or intrinsically disordered protein) and the structural context in which it occurs. Such information may positively aid drug-design. Furthermore, due to the intrinsic plasticity of proteins, even mutations involving radical change of the structural and physico⁻chemical properties of the amino acids (native vs. mutant) can still have minimal effects on protein thermodynamics. However, if a mutation causes significant perturbation by either folding or binding free energies, it is quite likely to be deleterious. Mitigating such effects is a promising alternative to the traditional approaches of designing inhibitors. This can be done by structure-based in silico screening of small molecules for which binding to the dysfunctional protein restores its wild type thermodynamics. In this review we emphasize the effects of mutations on two important biophysical properties, stability and binding affinity, and how structures can be used for structure-based drug design to mitigate the effects of disease-causing variants on the above biophysical properties.