Molecular Motions and Free-Energy Landscape of Serine Proteinase K in Relation to Its Cold-Adaptation: A Comparative Molecular Dynamics Simulation Study and the Underlying Mechanisms

Peng Sang,Xing Du,Li-Quan Yang,Zhao-Hui Meng,Shu-Qun Liu

DOI: https://doi.org/10.1039/c6ra23230b

IF: 4.036

2017-01-01

RSC Advances

Abstract:The physicochemical bases for enzyme cold-adaptation remain elusive.

What problem does this paper attempt to address?

ns and free-energy landscape of serine proteinase K in relation to its cold-adaptation : a comparative molecular dynamics simulation study and the underlying mechanisms †

Peng Sang,Xing Du,Li-Quan Yang,Zhao-Hui Meng,Shu-Qun Liu

2017-01-01

Abstract:Laboratory of Molecular Cardiology, Depa Hospital of Kunming Medical University, K aliyun.com; shuqunliu@ynu.edu.cn College of Agriculture and Biological Scienc Laboratory for Conservation and Utilizati Kunming, P. R. China Department of Biochemistry and Molecul Ningxia, P. R. China † Electronic supplementary informa 10.1039/c6ra23230b ‡ These authors contributed equally to th Cite this: RSC Adv., 2017, 7, 28580
Insights Derived from Molecular Dynamics Simulation into the Molecular Motions of Serine Protease Proteinase K

Shu-Qun Liu,Zhao-Hui Meng,Yun-Xin Fu,Ke-Qin Zhang

DOI: https://doi.org/10.1007/s00894-009-0518-x

IF: 2.172

2009-01-01

Journal of Molecular Modeling

Abstract:Serine protease proteinase K, a member of the subtilisin family of enzymes, is of significant industrial, agricultural and biotechnological importance. Despite the wealth of structural information about proteinase K provided by static X-ray structures, a full understanding of the enzymatic mechanism requires further insight into the dynamic properties of this enzyme. Molecular dynamics simulations and essential dynamics (ED) analysis were performed to investigate the molecular motions in proteinase K. The results indicate that the internal core of proteinase K is relatively rigid, whereas the surface-exposed loops, most notably the substrate-binding regions, exhibit considerable conformational fluctuations. Further ED analysis reveals that the large concerted motions in the substrate-binding regions cause opening/closing of the substrate-binding pockets, thus supporting the proposed induced-fit mechanism of substrate binding. The distinct electrostatic/hydrogen-bonding interactions between Asp39 and His69 and between His69 and Ser224 within the catalytic triad lead to different thermal motions and orientations of these three catalytic residues, which can be related to their different functional roles in the catalytic process. Statistical analyses of the geometrical/functional properties as well as evolutionary conservation of the glycines in proteinase K-like proteins reveal that glycines may play an important role in determining the folding architecture and structural flexibility of this class of enzymes. Our simulation study complements the biochemical and structural studies and provides new insights into the dynamic structural basis of the functional properties of this class of enzymes.
Effect of the Solvent Temperatures on Dynamics of Serine Protease Proteinase K

Peng Sang,Qiong Yang,Xing Du,Nan Yang,Li-Quan Yang,Xing-Lai Ji,Yun-Xin Fu,Zhao-Hui Meng,Shu-Qun Liu

DOI: https://doi.org/10.3390/ijms17020254

IF: 5.6

2016-02-19

International Journal of Molecular Sciences

Abstract:To obtain detailed information about the effect of the solvent temperatures on protein dynamics, multiple long molecular dynamics (MD) simulations of serine protease proteinase K with the solute and solvent coupled to different temperatures (either 300 or 180 K) have been performed. Comparative analyses demonstrate that the internal flexibility and mobility of proteinase K are strongly dependent on the solvent temperatures but weakly on the protein temperatures. The constructed free energy landscapes (FELs) at the high solvent temperatures exhibit a more rugged surface, broader spanning range, and higher minimum free energy level than do those at the low solvent temperatures. Comparison between the dynamic hydrogen bond (HB) numbers reveals that the high solvent temperatures intensify the competitive HB interactions between water molecules and protein surface atoms, and this in turn exacerbates the competitive HB interactions between protein internal atoms, thus enhancing the conformational flexibility and facilitating the collective motions of the protein. A refined FEL model was proposed to explain the role of the solvent mobility in facilitating the cascade amplification of microscopic motions of atoms and atomic groups into the global collective motions of the protein.

biochemistry & molecular biology,chemistry, multidisciplinary
Exploring the molecular mechanism of cold-adaption of an alkaline protease mutant by molecular dynamics simulations and residue interaction network

Ailan Huang,Fuping Lu,Fufeng Liu

DOI: https://doi.org/10.1002/pro.4837

Abstract:Psychrophilic proteases have attracted enormous attention in past decades, due to their high catalytic activity at low temperatures in a wide range of industrial processes, especially in the detergent and leather industries. Among them, H5 is an alkaline protease mutant, which featuring psychrophilic-like behavior, but the reasons that H5 with higher activity at low temperatures are still poorly understood. Herein, the molecular dynamics (MD) simulations combined with residue interaction network (RIN) were utilized to investigate the mechanisms of the cold-adaption of mutant H5. The results demonstrated that two loops involved in the substrate binding G100-S104 and S125-S129 in H5 had higher mobility, and the distance enlargement between the two loops modulated the substrate's accessibility compared with wild type counterpart. Besides, H5 enhanced conformational flexibility by weakening salt bridges and increasing interaction with the solvent. In particular, the absence of Lys251-Asp197-Arg247 salt bridge network may contribute to the structural mobility. Based on the free energy landscape and molecular mechanics Poisson-Boltzmann surface area of the wild type and H5, it was elucidated that H5 possessed a large population of interconvertible conformations, resulting in the weaker substrate binding free energy. The calculated RIN topology parameters such as the average degree, average cluster coefficient, and average path length further verified that the mutant H5 attenuated residue-to-residue interactions. The investigation of the mechanisms by which how the residue mutation affects the stability and activity of enzymes provides a theoretical basis for the development of cold-adapted protease.
Insight Derived from Molecular Dynamics Simulation into Substrate-Induced Changes in Protein Motions of Proteinase K.

Yan Tao,Zi-He Rao,Shu-Qun Liu

DOI: https://doi.org/10.1080/073911010010524953

2010-01-01

Journal of Biomolecular Structure and Dynamics

Abstract:Because of the significant industrial, agricultural and biotechnological importance of serine protease proteinase K, it has been extensively investigated using experimental approaches such as X-ray crystallography, site-directed mutagenesis and kinetic measurement. However, detailed aspects of enzymatic mechanism such as substrate binding, release and relevant regulation remain unstudied. Molecular dynamics (MD) simulations of the proteinase K alone and in complex with the peptide substrate AAPA were performed to investigate the effect of substrate binding on the dynamics/molecular motions of proteinase K. The results indicate that during simulations the substrate-complexed proteinase K adopt a more compact and stable conformation than the substrate-free form. Further essential dynamics (ED) analysis reveals that the major internal motions are confined within a subspace of very small dimension. Upon substrate binding, the overall flexibility of the protease is reduced; and the noticeable displacements are observed not only in substrate-binding regions but also in regions opposite the substrate-binding groove/pockets. The dynamic pockets caused by the large concerted motions are proposed to be linked to the substrate recognition, binding, orientation and product release; and the significant displacements in regions opposite the binding groove/pockets are considered to play a role in modulating the dynamics of enzyme-substrate interaction. Our simulation results complement the biochemical and structural studies, highlighting the dynamic mechanism of the functional properties of proteinase K.
Computer simulations explain the anomalous temperature optimum in a cold-adapted enzyme

Jaka Sočan,Miha Purg,Johan Åqvist

DOI: https://doi.org/10.1038/s41467-020-16341-2

IF: 16.6

2020-05-26

Nature Communications

Abstract:Abstract Cold-adapted enzymes from psychrophilic species show the general characteristics of being more heat labile, and having a different balance between enthalpic and entropic contributions to free energy barrier of the catalyzed reaction compared to mesophilic orthologs. Among cold-adapted enzymes, there are also examples that show an enigmatic inactivation at higher temperatures before unfolding of the protein occurs. Here, we analyze these phenomena by extensive computer simulations of the catalytic reactions of psychrophilic and mesophilic α-amylases. The calculations yield temperature dependent reaction rates in good agreement with experiment, and also elicit the anomalous rate optimum for the cold-adapted enzyme, which occurs about 15 °C below the melting point. This result allows us to examine the structural basis of thermal inactivation, which turns out to be caused by breaking of a specific enzyme-substrate interaction. This type of behaviour is also likely to be relevant for other enzymes displaying such anomalous temperature optima.

multidisciplinary sciences
Temperature effects on structure and dynamics of the psychrophilic protease subtilisin S41 and its thermostable mutants in solution

Ronny Martinez,Ulrich Schwaneberg,Danilo Roccatano

DOI: https://doi.org/10.1093/protein/gzr014

Abstract:The psychrophilic protease subtilisin S41 from the Antarctic bacillus TA41, and two variants with two and seven amino acid substitutions were studied using molecular dynamics simulation at 283 and 363 K. The analysis of protein dynamics revealed that the average global flexibility of both variants was slightly higher than wild type at both 283 and 363 K. Essential dynamics analysis evidenced that the most relevant collective motions, especially at 363 K, differ in distribution and intensity for each protein variant. At high temperature and for the thermo labile wild type, an amplification of a subset of the low-temperature largest collective motions was observed. On the other hand, the two thermostable variants showed a rather different pattern of essential motions at 363 K from those at 283 K. These results support the hypothesis that the introduced amino acid substitutions, rather than improving the global stability of the variants by increasing its rigidity, lead to a change on the principal fluxional modes allowing the protein to explore a different subset of conformations. A better understanding of this process can open alternative strategies to increase the enzyme stability in addition to increasing the rigidity of the protein scaffold.
Structure and in silico simulations of a cold-active esterase reveals its prime cold-adaptation mechanism

Nehad Noby,Husam Sabah Auhim,Samuel Winter,Harley L Worthy,Amira M Embaby,Hesham Saeed,Ahmed Hussein,Christopher R Pudney,Pierre J Rizkallah,Stephen A Wells,D Dafydd Jones

DOI: https://doi.org/10.1098/rsob.210182

Abstract:Here we determined the structure of a cold active family IV esterase (EstN7) cloned from Bacillus cohnii strain N1. EstN7 is a dimer with a classical α/β hydrolase fold. It has an acidic surface that is thought to play a role in cold-adaption by retaining solvation under changed water solvent entropy at lower temperatures. The conformation of the functionally important cap region is significantly different to EstN7's closest relatives, forming a bridge-like structure with reduced helical content providing greater access to the active site through more than one substrate access tunnel. However, dynamics do not appear to play a major role in cold adaption. Molecular dynamics at different temperatures, rigidity analysis, normal mode analysis and geometric simulations of motion confirm the flexibility of the cap region but suggest that the rest of the protein is largely rigid. Rigidity analysis indicates the distribution of hydrophobic tethers is appropriate to colder conditions, where the hydrophobic effect is weaker than in mesophilic conditions due to reduced water entropy. Thus, it is likely that increased substrate accessibility and tolerance to changes in water entropy are important for of EstN7's cold adaptation rather than changes in dynamics.
Mapping transiently formed and sparsely populated conformations on a complex energy landscape

Yong Wang,Elena Papaleo,Kresten Lindorff-Larsen

DOI: https://doi.org/10.7554/elife.17505

IF: 7.7

2016-01-01

eLife

Abstract:Determining the structures, kinetics, thermodynamics and mechanisms that underlie conformational exchange processes in proteins remains extremely difficult. Only in favourable cases is it possible to provide atomic-level descriptions of sparsely populated and transiently formed alternative conformations. Here we benchmark the ability of enhanced-sampling molecular dynamics simulations to determine the free energy landscape of the L99A cavity mutant of T4 lysozyme. We find that the simulations capture key properties previously measured by NMR relaxation dispersion methods including the structure of a minor conformation, the kinetics and thermodynamics of conformational exchange, and the effect of mutations. We discover a new tunnel that involves the transient exposure towards the solvent of an internal cavity, and show it to be relevant for ligand escape. Together, our results provide a comprehensive view of the structural landscape of a protein, and point forward to studies of conformational exchange in systems that are less characterized experimentally.
Parallel molecular mechanisms for enzyme temperature adaptation

Margaux M Pinney,Daniel A Mokhtari,Eyal Akiva,Filip Yabukarski,David M Sanchez,Ruibin Liang,Tzanko Doukov,Todd J Martinez,Patricia C Babbitt,Daniel Herschlag,Margaux M. Pinney,Daniel A. Mokhtari,David M. Sanchez,Todd J. Martinez,Patricia C. Babbitt

DOI: https://doi.org/10.1126/science.aay2784

IF: 56.9

2021-03-05

Science

Abstract:Some like it hot, others not Enzymes strike a delicate balance between features that enhance chemical reactivity and those that contribute to stable structure. Both features are important and can be unrelated or antagonistic. Pinney et al. combined rich experimental work on thermophilic and mesophilic variants of the enzyme ketosteroid isomerase (KSI) with bioinformatic data from a diverse set of bacterial enzymes to reveal the molecular determinants of thermal adaptation in enzymes. For KSI, they observed a trade-off between activity and thermal stability that comes down to a single active-site residue. With their larger dataset, they identified patterns of individual amino acid substitutions that are favored at higher temperatures, and also consider how networks of stabilizing interactions develop. Science , this issue p. eaay2784

multidisciplinary sciences
Comparative analysis of thermal adaptations of extremophilic prolyl oligopeptidases

Elizabeth M Diessner,Gemma R Takahashi,Carter T Butts,Rachel W Martin

DOI: https://doi.org/10.1016/j.bpj.2024.07.013

2024-09-17

Abstract:Prolyl oligopeptidases from psychrophilic, mesophilic, and thermophilic organisms found in a range of natural environments are studied using a combination of protein structure prediction, atomistic molecular dynamics, and trajectory analysis to determine how the S9 protease family adapts to extreme thermal conditions. We compare our results with hypotheses from the literature regarding structural adaptations that allow proteins to maintain structure and function at extreme temperatures, and we find that, in the case of prolyl oligopeptidases, only a subset of proposed adaptations are employed for maintaining stability. The catalytic and propeller domains are highly structured, limiting the range of mutations that can be made to enhance hydrophobicity or form disulfide bonds without disrupting the formation of necessary secondary structure. Rather, we observe a pattern in which overall prevalence of bound interactions (salt bridges and hydrogen bonds) is conserved by using increasing numbers of increasingly short-lived interactions as temperature increases. This suggests a role for an entropic rather than energetic strategy for thermal adaptation in this protein family.
Adaptations of protein structure and function to temperature: there is more than one way to ‘skin a cat’

Peter A. Fields,Yunwei Dong,Xianliang Meng,George N. Somero,Jason E. Podrabsky,Jonathon H. Stillman,Lars Tomanek

DOI: https://doi.org/10.1242/jeb.114298

IF: 3.308

2015-06-01

Journal of Experimental Biology

Abstract:ABSTRACT Sensitivity to temperature helps determine the success of organisms in all habitats, and is caused by the susceptibility of biochemical processes, including enzyme function, to temperature change. A series of studies using two structurally and catalytically related enzymes, A4-lactate dehydrogenase (A4-LDH) and cytosolic malate dehydrogenase (cMDH) have been especially valuable in determining the functional attributes of enzymes most sensitive to temperature, and identifying amino acid substitutions that lead to changes in those attributes. The results of these efforts indicate that ligand binding affinity and catalytic rate are key targets during temperature adaptation: ligand affinity decreases during cold adaptation to allow more rapid catalysis. Structural changes causing these functional shifts often comprise only a single amino acid substitution in an enzyme subunit containing approximately 330 residues; they occur on the surface of the protein in or near regions of the enzyme that move during catalysis, but not in the active site; and they decrease stability in cold-adapted orthologs by altering intra-molecular hydrogen bonding patterns or interactions with the solvent. Despite these structure–function insights, we currently are unable to predict a priori how a particular substitution alters enzyme function in relation to temperature. A predictive ability of this nature might allow a proteome-wide survey of adaptation to temperature and reveal what fraction of the proteome may need to adapt to temperature changes of the order predicted by global warming models. Approaches employing algorithms that calculate changes in protein stability in response to a mutation have the potential to help predict temperature adaptation in enzymes; however, using examples of temperature-adaptive mutations in A4-LDH and cMDH, we find that the algorithms we tested currently lack the sensitivity to detect the small changes in flexibility that are central to enzyme adaptation to temperature.

biology
Biomolecular conformational changes and ligand binding: from kinetics to thermodynamics.

Yong Wang,João Miguel Martins,Kresten Lindorff-Larsen

DOI: https://doi.org/10.1039/c7sc01627a

IF: 8.4

2017-01-01

Chemical Science

Abstract:The behaviour of biomolecular systems is governed by their thermodynamic and kinetic properties. It is thus important to be able to calculate, for example, both the affinity and rate of binding and dissociation of a protein-ligand complex, or the populations and exchange rates between distinct conformational states. Because these are typically rare events, calculating these properties from long molecular dynamics simulations remains extremely difficult. Instead, one often adopts a divide-and-conquer strategy in which equilibrium free-energy differences and the fastest state-to-state transition (e.g. ligand association or minor-to-major state conversion) are combined to estimate the slow rate (e.g. ligand dissociation) using a two-state assumption. Here we instead address these problems by using a previously developed method to calculate both the forward and backward rates directly from simulations. We then estimate the thermodynamics from the rates, and validate these values by independent means. We applied the approach to three systems of increasing complexity, including the association and dissociation of benzene to a fully buried cavity inside the L99A mutant variant of T4 lysozyme. In particular, we were able to determine both millisecond association and dissociation rates, and the affinity, of the protein-ligand system by directly observing dozens of rare events in atomic detail. Our approach both sheds light on the precision of methods for calculating kinetics and further provides a generally useful test for the internal consistency of kinetics and thermodynamics. We also expect our route to be useful for obtaining both the kinetics and thermodynamics at the same time in more challenging cases.
Molecular determinants for cold adaptation in an Antarctic Na+/K+-ATPase

Gaddiel Galarza-Muñoz,Sonia I Soto-Morales,Song Jiao,Miguel Holmgren,Joshua J C Rosenthal

DOI: https://doi.org/10.1073/pnas.2301207120

2023-10-10

Abstract:Enzymes from ectotherms living in chronically cold environments have evolved structural innovations to overcome the effects of temperature on catalysis. Cold adaptation of soluble enzymes is driven by changes within their primary structure or the aqueous milieu. For membrane-embedded enzymes, like the Na+/K+-ATPase, the situation is different because changes to the lipid bilayer in which they operate may also be relevant. Although much attention has been focused on thermal adaptation within lipid bilayers, relatively little is known about the contribution of structural changes within membrane-bound enzymes themselves. The identification of specific mutations that confer temperature compensation is complicated by the presence of neutral mutations, which can be more numerous. In the present study, we identified specific amino acids in a Na+/K+-ATPase from an Antarctic octopus that underlie cold resistance. Our approach was to generate chimeras between an Antarctic clone and a temperate ortholog and then study their temperature sensitivities in Xenopus oocytes using an electrophysiological approach. We identified 12 positions in the Antarctic Na+/K+-ATPase that, when transferred to the temperate ortholog, were sufficient to confer cold tolerance. Furthermore, although all 12 Antarctic mutations were required for the full phenotype, a single leucine in the third transmembrane segment (M3) imparted most of it. Mutations that confer cold resistance are mostly in transmembrane segments, at positions that face the lipid bilayer. We propose that the interface between a transmembrane enzyme and the lipid bilayer is a critical determinant of temperature sensitivity and, accordingly, has been a prime evolutionary target for thermal adaptation.
Exploration of Multi-State Conformational Dynamics and Underlying Global Functional Landscape of Maltose Binding Protein

Yong Wang,Chun Tang,Erkang Wang,Jin Wang

DOI: https://doi.org/10.1371/journal.pcbi.1002471

2012-01-01

PLoS Computational Biology

Abstract:An increasing number of biological machines have been revealed to have more than two macroscopic states. Quantifying the underlying multiple-basin functional landscape is essential for understanding their functions. However, the present models seem to be insufficient to describe such multiple-state systems. To meet this challenge, we have developed a coarse grained triple-basin structure-based model with implicit ligand. Based on our model, the constructed functional landscape is sufficiently sampled by the brute-force molecular dynamics simulation. We explored maltose-binding protein (MBP) which undergoes large-scale domain motion between open, apo-closed (partially closed) and holo-closed (fully closed) states responding to ligand binding. We revealed an underlying mechanism whereby major induced fit and minor population shift pathways co-exist by quantitative flux analysis. We found that the hinge regions play an important role in the functional dynamics as well as that increases in its flexibility promote population shifts. This finding provides a theoretical explanation of the mechanistic discrepancies in PBP protein family. We also found a functional “backtracking” behavior that favors conformational change. We further explored the underlying folding landscape in response to ligand binding. Consistent with earlier experimental findings, the presence of ligand increases the cooperativity and stability of MBP. This work provides the first study to explore the folding dynamics and functional dynamics under the same theoretical framework using our triple-basin functional model.
Computational design of the temperature optimum of an enzyme reaction

Florian van der Ent,Susann Skagseth,Bjarte A. Lund,Jaka Sočan,Julia J. Griese,Bjørn O. Brandsdal,Johan Åqvist,Jaka Sǒan

DOI: https://doi.org/10.1126/sciadv.adi0963

IF: 13.6

2023-06-29

Science Advances

Abstract:Cold-adapted enzymes are characterized both by a higher catalytic activity at low temperatures and by having their temperature optimum down-shifted, compared to mesophilic orthologs. In several cases, the optimum does not coincide with the onset of protein melting but reflects some other type of inactivation. In the psychrophilic α-amylase from an Antarctic bacterium, the inactivation is thought to originate from a specific enzyme-substrate interaction that breaks around room temperature. Here, we report a computational redesign of this enzyme aimed at shifting its temperature optimum upward. A set of mutations designed to stabilize the enzyme-substrate interaction were predicted by computer simulations of the catalytic reaction at different temperatures. The predictions were verified by kinetic experiments and crystal structures of the redesigned α-amylase, showing that the temperature optimum is indeed markedly shifted upward and that the critical surface loop controlling the temperature dependence approaches the target conformation observed in a mesophilic ortholog.

multidisciplinary sciences
Hidden Conformational States and Strange Temperature Optima in Enzyme Catalysis

Johan Åqvist,Jaka Sočan,Miha Purg,Jaka Sočan

DOI: https://doi.org/10.1021/acs.biochem.0c00705

IF: 3.321

2020-09-25

Biochemistry

Abstract:The existence of temperature optima in enzyme catalysis that occur before protein melting sets in can be described by different types of kinetic models. Such optima cause distinctly curved Arrhenius plots and have, for example, been observed in several cold-adapted enzymes from psychrophilic species. The two main explanations proposed for this behavior either invoke conformational equilibria with inactive substrate-bound states or postulate differences in heat capacity between the reactant and transition states. Herein, we analyze the implications of the different types of kinetic models in terms of apparent activation enthalpies, entropies, and heat capacities, using the catalytic reaction of a cold-adapted α-amylase as a prototypic example. We show that the behavior of these thermodynamic activation parameters is fundamentally different between equilibrium and heat capacity models, and in the α-amylase case, computer simulations have shown the former model to be correct. A few other enzyme-catalyzed reactions are also discussed in this context.This article has not yet been cited by other publications.

biochemistry & molecular biology
Molecular dynamics method on thermostability of cold active Lipase 5

LV Jian-Ping,WEI Dong-Qing,WANG Yong-Hua,Xu Qin

DOI: https://doi.org/10.3969/j.issn.1000-0364.2016.02.022

2016-01-01

Abstract:The thermostability of naturally cold -adapted lipases is often relatively poor , limiting their long cata-lytic effect and long-term preservation.In the current study , one cold-adapted lipase Lipase 5 from Candida albicans was investigated by means of molecular dynamics method in order to improve its thermal stability.First the three-dimensional model of Lipase 5 was constructed by homology modeling.Then 18 ns molecular dynam-ics simulations were utilized to find the unstable region in Lipase 5 out of the protein active center with high con-formational fluctuations , and one possible mutant site 279-Gly was detected.It is finally proved that mutation of flexible 279-Gly to rigid Proline endows Lipase 5 enhances thermostability.Our studies may provide theoretical direction and support to design experiments of reforming thermal stability of cold active lipases.
Enhance Cold Adaptation of Bidomain Amylases via High-throughput Computational Engineering

Ning Ding,Yaoyukun Jiang,Robbie Ge,Xinchun Ran,Zhongyue J. Yang

DOI: https://doi.org/10.26434/chemrxiv-2024-rstbz

2024-08-20

Abstract:Cold-adapted bidomain enzymes are vital for transforming modern industries by decreasing energy consumption, delivering economic benefits, and fostering sustainability through reduced greenhouse gas emissions. Yet, the design strategies guiding their acquisition of cold adaptation remain unknown. Here, we developed an integrated computational-experimental strategy to engineer bidomain enzymes for enhanced cold-adaptation. Using five model amylase variants exhibiting different degrees of cold adaptation, we identified a descriptor from molecular dynamics simulations, namely domain separation index (DSI), which positively correlates with bidomain amylases’ relative activity at 0°C. The bidomain amylase variants with a longer distance between its catalytic domain and carbohydrate-binding module (i.e., a high DSI) were observed to demonstrate cold adaptation. Guided by DSI, we developed a high-throughput molecular modeling protocol to convert the thermophilic Pseudomonas saccharophila amylase (psA) into a cold-adapted bidomain enzyme, virtually screening 120 psA variants with different linkers. Two psA variants with a greater DSI value were selected and experimentally confirmed to be cold-adapted, with the psA121 variant achieving a 12-fold increase in relative activity at 0°C from 2.4% (specific activity: 14 U/mg) to 30.5% (specific activity: 219 U/mg). Conformational analyses reveal that compared to non-cold-adapted counterparts, cold-adapted variants leverage its linker to induce domain separation and enhance flexibility of active-site and binding loop via dynamic allostery, thereby promoting substrate recruiting, binding, and catalysis at lower temperatures. Statistical analyses of 120 variants demonstrate that helical motifs within linkers drive interdomain separation. Overall, our study offers strategies for engineering cold-adapted bidomain enzymes and suggests the molecular basis of cold adaptation in bidomain amylases.

Chemistry
Comparative molecular dynamics simulations provided insights into the mechanisms of cold-adaption of alginate lyases from the PL7 family

Fei Bian,Xiao-Yue Liang,Meng Wang,Zhong-Zhi Sun,Bin-Bin Xie

DOI: https://doi.org/10.1007/s00792-024-01340-8

2024-04-11

Extremophiles

Abstract:Alginate is an important polysaccharide that is abundant in the marine environments, including the Polar Regions, and bacterial alginate lyases play key roles in its degradation. Many reported alginate lyases show characteristics of cold-adapted enzymes, including relatively low temperature optimum of activities (T opt ) and low thermal stabilities. However, the cold-adaption mechanisms of alginate lyases remain unclear. Here, we studied the cold-adaptation mechanisms of alginate lyases by comparing four members of the PL7 family from different environments: AlyC3 from the Arctic ocean ( Psychromonas sp. C-3), AlyA1 from the temperate ocean ( Zobellia galactanivorans ), PA1167 from the human pathogen ( Pseudomonas aeruginosa PAO1), and AlyQ from the tropic ocean ( Persicobacter sp. CCB-QB2). Sequence comparison and comparative molecular dynamics (MD) simulations revealed two main strategies of cold adaptation. First, the Arctic AlyC3 and temperate AlyA1 increased the flexibility of the loops close to the catalytic center by introducing insertions at these loops. Second, the Arctic AlyC3 increased the electrostatic attractions with the negatively charged substrate by introducing a high portion of positively charged lysine at three of the insertions mentioned above. Furthermore, our study also revealed that the root mean square fluctuation (RMSF) increased greatly when the temperature was increased to T opt or higher, suggesting the RMSF increase temperature as a potential indicator of the cold adaptation level of the PL7 family. This study provided new insights into the cold-adaptation mechanisms of bacterial alginate lyases and the marine carbon cycling at low temperatures.

microbiology,biochemistry & molecular biology

Molecular Motions and Free-Energy Landscape of Serine Proteinase K in Relation to Its Cold-Adaptation: A Comparative Molecular Dynamics Simulation Study and the Underlying Mechanisms

ns and free-energy landscape of serine proteinase K in relation to its cold-adaptation : a comparative molecular dynamics simulation study and the underlying mechanisms †

Insights Derived from Molecular Dynamics Simulation into the Molecular Motions of Serine Protease Proteinase K

Effect of the Solvent Temperatures on Dynamics of Serine Protease Proteinase K

Exploring the molecular mechanism of cold-adaption of an alkaline protease mutant by molecular dynamics simulations and residue interaction network

Insight Derived from Molecular Dynamics Simulation into Substrate-Induced Changes in Protein Motions of Proteinase K.

Computer simulations explain the anomalous temperature optimum in a cold-adapted enzyme

Temperature effects on structure and dynamics of the psychrophilic protease subtilisin S41 and its thermostable mutants in solution

Structure and in silico simulations of a cold-active esterase reveals its prime cold-adaptation mechanism

Mapping transiently formed and sparsely populated conformations on a complex energy landscape

Parallel molecular mechanisms for enzyme temperature adaptation

Comparative analysis of thermal adaptations of extremophilic prolyl oligopeptidases

Adaptations of protein structure and function to temperature: there is more than one way to ‘skin a cat’

Biomolecular conformational changes and ligand binding: from kinetics to thermodynamics.

Molecular determinants for cold adaptation in an Antarctic Na+/K+-ATPase

Exploration of Multi-State Conformational Dynamics and Underlying Global Functional Landscape of Maltose Binding Protein

Computational design of the temperature optimum of an enzyme reaction

Hidden Conformational States and Strange Temperature Optima in Enzyme Catalysis

Molecular dynamics method on thermostability of cold active Lipase 5

Enhance Cold Adaptation of Bidomain Amylases via High-throughput Computational Engineering

Comparative molecular dynamics simulations provided insights into the mechanisms of cold-adaption of alginate lyases from the PL7 family