Pressure-induced conformational switch of an interfacial protein

The importance of helical structures to the overall activity and structural stability of a lipase from Pseudomonas aeruginosa PAO1 in n-hexane

Bacterial lipases are versatile extracellular enzymes with a catalytic triad at the active site and a flexible 'lid' that modulates catalytic accessibility. We combined computational modeling with preliminary in vitro testing to assess the structural stability and activity of the Pseudomonas aeruginosa PAO1 lipase (PAL). We evaluated several systems consisting of the native and mutant forms of the lipase in n-hexane using molecular dynamics simulations. Structural stability was assessed by calculating the B-factor for each system. We measured the gorge radius of the catalytic channel and the RMSD of the catalytic triad to approximate enzymatic activity. Based on the correlation of these metrics, mutant forms were selected for their potential activity and stability. Selected mutant forms were expressed in E. coli BL21, mass-produced, and validated through a lipase-catalyzed esterification assay in n-hexane. Several helices outside the 'lid' region were found to influence lid conformational switching. Moreover, our preliminary experimental results show promise in validating our in silico predictions. Our integrated in silico and in vitro pipeline offers a promising approach for designing and producing industrially relevant lipases.

High hydrostatic pressure reduces inflammation induced by litchi thaumatin-like protein via altering active domain

Thaumatin-like proteins (TLP), existing in various fruits, have allergenic and pro-inflammatory activities. The current research attempts to reduce the pro-inflammatory activity of litchi TLP (LcTLP) through high hydrostatic pressure (HHP). This study demonstrated that HHP (250-500 MPa, 5-10 min) was a potential technique to reduce the pro-inflammatory activity of LcTLP, which was attributed to the irreversible destruction of the active domain, ie., V-cleft. SDS-PAGE showed no change in the protein profile. Continuous HHP treatment promoted LcTLP unfolding and then reassembling (400 MPa was the transition pressure), and the content of β-sheets decreased from 35.4% to 31.1%. HHP treatment could mitigate inflammatory responses of LcTLP, as confirmed by ELISA and western blot. Molecular dynamics simulations showed significant changes in the residue network under HHP, thereby affecting the V-cleft. These findings provide a theoretical explanation and structural insights into the HHP-induced reduction of pro-inflammatory activity of LcTLP.

Isothermal Compressibility Perturbation as a Protein Design Principle for T1 Lipase Stability-Activity Trade-Off Counteracting

Given the widely existing stability-activity trade-off in enzyme evolution, it is still a goal to obtain enzymes embracing both high activity and stability. Herein, we employed an isothermal compressibility (β_T) perturbation engineering (ICPE) strategy to comprehensively understand the stability-activity seesaw-like mechanism. The stability and activity of mutants derived from ICPE uncovered a high Pearson correlation (r = 0.93) in a prototypical enzyme T1 lipase. The best variant A186L/L188M/A190Y exhibited a high T_m value up to 78.70 °C, catalytic activity of 474.04 U/mg, and a 73.33% increase in dimethyl sulfoxide resistance compared to the wild type, one of the highest comprehensive performances reported to date. The elastic activation mechanism mediated by conformational change with a Δβ_T range of -6.81 × 10^-6 to -1.90 × 10^-6 bar^-1 may account for the balancing of stability and activity to achieve better performing enzymes. The ICPE strategy deepens our understanding of stability-activity trade-off and boosts its applications in enzyme engineering.

Rebuilding the lid region from conformational and dynamic features to engineering applications of lipase in foods: Current status and future prospects

The applications of lipases in esterification, amidation, and transesterification have broadened their potential in the production of fine compounds with high cumulative values. Mostly, the catalytic triad of lipases is covered by either one or two mobile peptides called the "lid" that control the substrate channel to the catalytic center. The lid holds unique conformational allostery via interfacial activation to regulate the dynamics and catalytic functions of lipases, thereby highlighting its importance in redesigning these enzymes for industrial applications. The structural characteristic of lipase, the dynamics of lids, and the roles of lid in lipase catalysis were summarized, providing opportunities for rebuilding lid region by biotechniques (e.g., metagenomic technology and protein engineering) and enzyme immobilization. The review focused on the advantages and disadvantages of strategies rebuilding the lid region. The main shortcomings of biotechnologies on lid rebuilding were discussed such as negative effects on lipase (e.g., a decrease of activity). Additionally, the main shortcomings (e.g., enzyme desorption at high temperatre) in immobilization on hydrophobic supports via interfacial action were presented. Solutions to the mentioned problems were proposed by combinations of computational design with biotechnologies, and improvements of lipase immobilization (e.g., immobilization protocols and support design). Finally, the review provides future perspectives about designing hyperfunctional lipases as biocatalysts in the food industry based on lid conformation and dynamics.

Increased activity of alcohol oxidase at high hydrostatic pressure

Alcohol oxidase (AOx) from P. pastoris has potential applications in the production of carbonyl compounds and for the detection and quantification of alcohols. However, AOx's poor stability and low activity have hindered its practical application. There are two fractions of AOx in P. pastoris with different thermal stability. High hydrostatic pressure (HHP) increased the activity of the labile (L) + resistant (R) combined fractions but not of the R fraction alone. The activity of the L + R fractions increased 2.4-fold at 160 MPa and 30 °C compared to the activity at 0.1 MPa. At higher temperatures, the increase in activity with pressure was greater due to the combined stabilization and activation effects. The reaction rate of the R fraction at 50 °C was 17.9 ± 3.6 or 17.7 ± 0.8 μM min^-1 at 80 or 160 MPa, respectively, and was not significantly different from the activity of the L + R fractions under the same conditions (18.4 ± 2.7 μM min^-1). The activation energy of the R fraction was not significantly different between 80 MPa (41.5 ± 10.5 kJ mol^-1) and 160 MPa (43.8 ± 7.8 kJ mol^-1). The combined increase in the stability of the R fraction at HHP enables the use of the enzyme at 50 °C with little loss of activity and an increased catalytic rate.

The effect of ligand affinity to the contact dynamics of the ligand binding domain of thyroid hormone receptor - retinoid X receptor

Ligand-based allostery has been gaining attention for its importance in protein regulation and implication in drug design. One of the interesting cases of protein allostery is the thyroid hormone receptor - retinoid x receptor (TR:RXR), which regulates the gene expression of important physiological processes, such as development and metabolism. It is regulated by the TR native ligand triiodothyronine (T3), which displays anticooperative behavior to the RXR ligand 9-cis retinoic acid (9C). In contrast to this anticooperative behavior, 9C has been shown to increase the activity of TR:RXR. Here we probed the influence of the affinity and the interactions of the TR ligand to the allostery of the TR:RXR through contact dynamics and residue networks. The TR ligand analogs were designed to have higher (G2) and lower (N1) binding energies than T3 when docked to the TR:RXR(9C) complex. The aqueous TR(N1/T3/G2):RXR(9C) complexes were subjected to 30 ns all-atom simulations using theNAMD. The program CAMERRA was used to capture the subtle perturbations of TR:RXR by mapping the residue contact dynamics. Various parts of the TR ligands; including the hydrophilic head, the iodine substituents, and the ligand tail; have been probed for their significance in ligand affinity. The results on the T3 and G2 complexes suggest that ligand affinity can be utilized as a predictor for anticooperative systems on which ligand is more likely to dissociate or remain bound. All 3 complexes also display distinct contact networks for cross-dimer signalling and ligand communication. Understanding ligand-based allostery could potentially unveil secrets of ligand-regulated protein dynamics, a foundation for the design of better and more efficient allosteric drugs.

In situ and real-time insight into Rhizopus chinensis lipase under high pressure and temperature: Conformational traits and biobehavioural analysis

An in situ and real-time investigation was performed using an optical cell system and in-silico analysis to reveal the impacts of pressure and temperature on the conformational state and behaviours of Rhizopus chinensis lipase (RCL). The fluorescence intensity (FI) of RCL increased remarkably under high pressure, and part of this increase was recovered after depressurization. This result displayed the partially reversible conformational change of RCL, which may be associated with the local change of Trp224 near the catalytic centre. High temperature caused a significant loss of secondary structure, whereas the α-helical segments including the lid were preserved by high pressure even at temperatures over 60 °C. The parameters of enzymatic reaction monitored by UV showed that the hydrolysis rate was remarkably enhanced by the pressure of 200 MPa. In the pressure range of 0.1-200 MPa, the active volume measured by the in situ system decreased from -2.85 to -6.73 mL/mol with the temperature increasing from 20 °C to 40 °C. The high catalytic capacity of the lipase under high pressure and high temperature was primarily attributed to pressure protection on RCL.

Dynamics, a Powerful Component of Current and Future in Silico Approaches for Protein Design and Engineering

Computational prediction has become an indispensable aid in the processes of engineering and designing proteins for various biotechnological applications. With the tremendous progress in more powerful computer hardware and more efficient algorithms, some of in silico tools and methods have started to apply the more realistic description of proteins as their conformational ensembles, making protein dynamics an integral part of their prediction workflows. To help protein engineers to harness benefits of considering dynamics in their designs, we surveyed new tools developed for analyses of conformational ensembles in order to select engineering hotspots and design mutations. Next, we discussed the collective evolution towards more flexible protein design methods, including ensemble-based approaches, knowledge-assisted methods, and provable algorithms. Finally, we highlighted apparent challenges that current approaches are facing and provided our perspectives on their further development.

The Membrane-Integrated Steric Chaperone Lif Facilitates Active Site Opening of Pseudomonas aeruginosa Lipase A

Lipases are essential and widely used biocatalysts. Hence, the production of lipases requires a detailed understanding of the molecular mechanism of its folding and secretion. Lipase A from Pseudomonas aeruginosa, PaLipA, constitutes a prominent example that has additional relevance because of its role as a virulence factor in many diseases. PaLipA requires the assistance of a membrane-integrated steric chaperone, the lipase-specific foldase Lif, to achieve its enzymatically active state. However, the molecular mechanism of how Lif activates its cognate lipase has remained elusive. Here, we show by molecular dynamics simulations at the atomistic level and potential of mean force computations that Lif catalyzes the activation process of PaLipA by structurally stabilizing an intermediate PaLipA conformation, particularly a β-sheet in the region of residues 17-30, such that the opening of PaLipA's lid domain is facilitated. This opening allows substrate access to PaLipA's catalytic site. A surprising and so far not fully understood aspect of our study is that the open state of PaLipA is unstable compared to the closed one according to our computational and in vitro biochemical results. We thus speculate that further interactions of PaLipA with the Xcp secretion machinery and/or components of the extracellular matrix contribute to the remaining activity of secreted PaLipA. © 2019 Wiley Periodicals, Inc.

Characterizing the 3D structure and dynamics of chromosomes and proteins in a common contact matrix framework

Conformational ensembles of biopolymers, whether proteins or chromosomes, can be described using contact matrices. Principal component analysis (PCA) on the contact data has been used to interrogate both protein and chromosome structures and/or dynamics. However, as these fields have developed separately, variants of PCA have emerged. Previously, a variant we hereby term Implicit-PCA (I-PCA) has been applied to chromosome contact matrices and revealed the spatial segregation of active and inactive chromatin. Separately, Explicit-PCA (E-PCA) has previously been applied to proteins and characterized their correlated structure fluctuations. Here, we swapped analysis methods (I-PCA and E-PCA), applying each to a different biopolymer type (chromosome or protein) than the one for which they were initially developed. We find that applying E-PCA to chromosome distance matrices derived from microscopy data can reveal the dominant motion (concerted fluctuation) of these chromosomes. Further, by applying E-PCA to Hi-C data across the human blood cell lineage, we isolated the aspects of chromosome structure that most strongly differentiate cell types. Conversely, when we applied I-PCA to simulation snapshots of proteins, the major component reported the consensus features of the structure, making this a promising approach for future analysis of semi-structured proteins.

Molecular dynamics simulation of proteins under high pressure: Structure, function and thermodynamics

BACKGROUND

Molecular dynamics (MD) simulation is well-recognized as a powerful tool to investigate protein structure, function, and thermodynamics. MD simulation is also used to investigate high pressure effects on proteins. For conducting better MD simulation under high pressure, the main issues to be addressed are: (i) protein force fields and water models were originally developed to reproduce experimental properties obtained at ambient pressure; and (ii) the timescale to observe the pressure effect is often much longer than that of conventional MD simulations.

SCOPE OF REVIEW

First, we describe recent developments in MD simulation methodologies for studying the high-pressure structure and dynamics of protein molecules. These developments include force fields for proteins and water molecules, and enhanced simulation techniques. Then, we summarize recent studies of MD simulations of proteins in water under high pressure.

MAJOR CONCLUSIONS

Recent MD simulations of proteins in solution under pressure have reproduced various phenomena identified by experiments using high pressure, such as hydration, water penetration, conformational change, helix stabilization, and molecular stiffening.

GENERAL SIGNIFICANCE

MD simulations demonstrate differences in the properties of proteins and water molecules between ambient and high-pressure conditions. Comparing the results obtained by MD calculations with those obtained experimentally could reveal the mechanism by which biological molecular machines work well in collaboration with water molecules.

Effector-Binding-Directed Dimerization and Dynamic Communication between Allosteric Sites of Ribonucleotide Reductase

Proteins forming dimers or larger complexes can be strongly influenced by their effector-binding status. We investigated how the effector-binding event is coupled with interface formation via computer simulations, and we quantified the correlation of two types of contact interactions: between the effector and its binding pocket and between protein monomers. This was achieved by connecting the protein dynamics at the monomeric level with the oligomer interface information. We applied this method to ribonucleotide reductase (RNR), an essential enzyme for de novo DNA synthesis. RNR contains two important allosteric sites, the s-site (specificity site) and the a-site (activity site), which bind different effectors. We studied these different binding states with atomistic simulation and used their coarse-grained contact information to analyze the protein dynamics. The results reveal that the effector-protein dynamics at the s-site and dimer interface formation are positively coupled. We further quantify the resonance level between these two events, which can be applied to other similar systems. At the a-site, different effector-binding states (ATP vs dATP) drastically alter the protein dynamics and affect the activity of the enzyme. On the basis of these results, we propose a new mechanism of how the a-site regulates enzyme activation.

Molecular Dynamics Simulation for Mechanism Elucidation of Food Processing and Safety: State of the Art

Molecular dynamics (MD) simulation is a useful technique to study the interaction between molecules and how they are affected by various processes and processing conditions. This review summarizes the application of MD simulations in food processing and safety, with an emphasis on the effects that emerging nonthermal technologies (for example, high hydrostatic pressure, pulsed electric field) have on the molecular and structural characteristics of foods and biomaterials. The advances and potential projection of MD simulations in the science and engineering aspects of food materials are discussed and focused on research work conducted to study the effects of emerging technologies on food components. It is expected by showing key case studies that it will stir novel developments as a valuable tool to study the effects of emerging food technologies on biomaterials. This review is useful to food researchers and the food industry, as well as researchers and practitioners working on flavor and nutraceutical encapsulations, dietary carbohydrate product developments, modified starches, protein engineering, and other novel food applications.

Characterizing protein conformations by correlation analysis of coarse-grained contact matrices

We have developed a method to capture the essential conformational dynamics of folded biopolymers using statistical analysis of coarse-grained segment-segment contacts. Previously, the residue-residue contact analysis of simulation trajectories was successfully applied to the detection of conformational switching motions in biomolecular complexes. However, the application to large protein systems (larger than 1000 amino acid residues) is challenging using the description of residue contacts. Also, the residue-based method cannot be used to compare proteins with different sequences. To expand the scope of the method, we have tested several coarse-graining schemes that group a collection of consecutive residues into a segment. The definition of these segments may be derived from structural and sequence information, while the interaction strength of the coarse-grained segment-segment contacts is a function of the residue-residue contacts. We then perform covariance calculations on these coarse-grained contact matrices. We monitored how well the principal components of the contact matrices is preserved using various rendering functions. The new method was demonstrated to assist the reduction of the degrees of freedom for describing the conformation space, and it potentially allows for the analysis of a system that is approximately tenfold larger compared with the corresponding residue contact-based method. This method can also render a family of similar proteins into the same conformational space, and thus can be used to compare the structures of proteins with different sequences.

CAMERRA: An analysis tool for the computation of conformational dynamics by evaluating residue-residue associations

A computational method which extracts the dominant motions from an ensemble of biomolecular conformations via a correlation analysis of residue-residue contacts is presented. The algorithm first renders the structural information into contact matrices, then constructs the collective modes based on the correlated dynamics of a selected set of dynamic contacts. Associated programs can bridge the results for further visualization using graphics software. The aim of this method is to provide an analysis of conformations of biopolymers from the contact viewpoint. It may assist a systematical uncovering of conformational switching mechanisms existing in proteins and biopolymer systems in general by statistical analysis of simulation snapshots. In contrast to conventional correlation analyses of Cartesian coordinates (such as distance covariance analysis and Cartesian principal component analysis), this program also provides an alternative way to locate essential collective motions in general. Herein, we detail the algorithm in a stepwise manner and comment on the importance of the method as applied to decoding allosteric mechanisms. © 2018 Wiley Periodicals, Inc.

Coupled effects of salt and pressure on catalytic ability of Rhizopus chinensis lipase

BACKGROUND

Both high pressure and environmental factors could influence the catalytic abilities of enzymes. This work investigated coupled effects of pressure and salts on Rhizopus chinensis lipase (RCL) to provide significant information for its further applications.

RESULTS

The maximum activity of RCL was observed under 200 MPa at 40 °C. The highest activity was achieved at concentrations of 0.06-0.1 mol L^-1 for tested salts. The effect of monovalent cations on RCL activity followed the Hofmeister series (K⁺ > Na⁺ > Li⁺ ) at 0.1 MPa but the order of Na⁺ and K⁺ was changed under 200 MPa. Meanwhile, the effects of anions did not follow the Hofmeister series. KCl slightly improved the thermostability of RCL at moderate concentration. At 60 °C, LiCl only stabilised RCL at 0.1 mol L^-1 . The pre-transition unfolding point was shifted from 4.5 to 3.5 mol L^-1 with pressure increasing from 0.1 to 600 MPa. In addition, KCl could not change the lipase's extrinsic fluorescence evolution versus pressure.

CONCLUSION

Pressure and salts could improve catalytic ability and stability of RCL under appropriate conditions. The effect of high pressure on RCL was influenced by salts. Meanwhile salts cannot prevent high pressure-induced damage to RCL. © 2017 Society of Chemical Industry.

DMSO enhanced conformational switch of an interfacial enzyme

Interfacial proteins function in unique heterogeneous solvent environments, such as water-oil interfaces. One important example is microbial lipase, which is activated in an oil-water emulsion phase and has many important enzymatic functions. A unique aprotic dipolar organic solvent, dimethyl sulfoxide (DMSO), has been shown to increase the activity of lipases, but the mechanism behind this enhancement is still unknown. Here, all-atom molecular dynamics simulations of lipase in a binary solution were performed to examine the effects of DMSO on the dynamics of the gating mechanism. The amphiphilic α5 region of the lipase was a focal point for the analysis, since the structural ordering of α5 has been shown to be important for gating under other perturbations. Compared to the closed-gorge ensemble in an aqueous environment, the conformational ensemble shifts towards open-gorge structures in the presence of DMSO solvents. Increased width of the access channel is particularly prevalent in 45% and 60% DMSO concentrations (w/w). As the amount of DMSO increases, the α5 region of the lipase becomes more α-helical, as we previously observed in studies that address water-oil interfacial and high pressure activation. We believe that the structural ordering of α5 plays an essential role on gating and lipase activity.