Essential motions in a fungal lipase with bound substrate, covalently attached inhibitor and product

Study on Lipase-Catalyzed Hydrolysis of Olive Oil at Oil-Water Interface

Olive oil was selected as the oil substrate and hydrolyzed by Candida sp. 99–125 lipase. The hydrolysis rate of olive oil was used as an indicator. Based on the single factor experiment, the effects of dosage of Candida sp. 99–125 lipase, reacting temperature, pH value and water-oil ratio were investigated. Box-Behnken center combination and response surface methodology were utilized to optimize the hydrolysis rate. The results showed that the significant differences of each single factor on lipase hydrolysis of olive oil on the oil-water interface were different. pH value is the first significance factor, and the significance of water oil ratio on lipase hydrolysis of olive oil is second only to pH value. Finally, the mechanism of Candida sp. 99–125 lipase hydrolyzing olive oil at the oil-water interface was discussed.

A tensegrity model for hydrogen bond networks in proteins

Hydrogen-bonding networks in proteins considered as structural tensile elements are in balance separately from any other stabilising interactions that may be in operation. The hydrogen bond arrangement in the network is reminiscent of tensegrity structures in architecture and sculpture. Tensegrity has been discussed before in cells and tissues and in proteins. In contrast to previous work only hydrogen bonds are studied here. The other interactions within proteins are either much stronger − covalent bonds connecting the atoms in the molecular skeleton or weaker forces like the so-called hydrophobic interactions. It has been demonstrated that the latter operate independently from hydrogen bonds. Each category of interaction must, if the protein is to have a stable structure, balance out. The hypothesis here is that the entire hydrogen bond network is in balance without any compensating contributions from other types of interaction. For sidechain-sidechain, sidechain-backbone and backbone-backbone hydrogen bonds in proteins, tensegrity balance (“closure”) is required over the entire length of the polypeptide chain that defines individually folding units in globular proteins (“domains”) as well as within the repeating elements in fibrous proteins that consist of extended chain structures. There is no closure to be found in extended structures that do not have repeating elements. This suggests an explanation as to why globular domains, as well as the repeat units in fibrous proteins, have to have a defined number of residues. Apart from networks of sidechain-sidechain hydrogen bonds there are certain key points at which this closure is achieved in the sidechain-backbone hydrogen bonds and these are associated with demarcation points at the start or end of stretches of secondary structure. Together, these three categories of hydrogen bond achieve the closure that is necessary for the stability of globular protein domains as well as repeating elements in fibrous proteins.

Effect of ligand binding on the dynamics of trypsin. Comparison of different approaches

The intramolecular signal transduction induced by the binding of ligands to trypsin was investigated by molecular dynamics simulations. Ligand binding changes the residue-residue interaction energies and suppresses the mobility of loops that are in direct contact with the ligand. The reduced mobility of these loops results in the altered flexibility of the nearby loops and thereby transmits the information from ligand binding site to the remote sites. The analysis of the flexibility of all residues confirmed the coupling between loops L1 (185-188) and L2 (221-224) and the residues in the active center. The significance of S1 pocket residues for the signal transduction from the active center to the substrate-binding site was confirmed by the dynamical network and covariance matrix analyses. Gaussian network model and principal component analysis demonstrated that the active center residues had zero amplitude in the slowest fluctuations acting as hinges or anchors. Overall, our results provide a new insight into protein-ligand interactions and show how the allosteric signaling may occur.

Accelerated simulation of unfolding and refolding of a large single chain globular protein

We have developed novel strategies for contracting simulation times in protein dynamics that enable us to study a complex protein with molecular weight in excess of 34 kDa. Starting from a crystal structure, we produce unfolded and then refolded states for the protein. We then compare these quantitatively using both established and new metrics for protein structure and quality checking. These include use of the programs Concoord and Darvols. Simulation of protein-folded structure well beyond the molten globule state and then recovery back to the folded state is itself new, and our results throw new light on the protein-folding process. We accomplish this using a novel cooling protocol developed for this work.

Inactivation and reactivation of ribonuclease A studied by computer simulation

The year 2011 marked the half-centenary of the publication of what came to be known as the Anfinsen postulate, that the tertiary structure of a folded protein is prescribed fully by the sequence of its constituent amino acid residues. This postulate has become established as a credo, and, indeed, no contradictions seem to have been found to date. However, the experiments that led to this postulate were conducted on only a single protein, bovine ribonuclease A (RNAse). We conduct molecular dynamics (MD) simulations on this protein with the aim of mimicking this experiment as well as making the methodology available for use with basically any protein. There have been many attempts to model denaturation and refolding processes of globular proteins in silico using MD, but only a few examples where disulphide-bond containing proteins were studied. We took the view that if the reductive deactivation and oxidative reactivation processes of RNAse could be modelled in silico, this would provide valuable insights into the workings of the classical Anfinsen experiment.

Combination of site-directed mutagenesis and yeast surface display enhances Rhizomucor miehei lipase esterification activity in organic solvent

To increase the activity of Rhizomucor miehei lipase (RML) in organic solvent, multiple sequence alignments and rational site-directed mutagenesis were used to create RML variants. The obtained proteins were surface-displayed on Pichia pastoris by fusion to Flo1p as an anchor protein. The synthetic activity of four variants showed from 1.1- to 5-fold the activity of native lipase in an esterification reaction in heptane with alcohol and caproic acid as substrates. The increase in esterification activity may be attributed to the four mutations changing the flexibility of RML or facilitating the reaction. In conclusion, this method demonstrated that multiple sequence alignments and rational site-directed mutagenesis combined with yeast display technology is a faster and more effective means of obtaining high-efficiency esterification lipase variants compared with previous similar methods.

The effect of Asp54 phosphorylation on the energetics and dynamics in the response regulator protein Spo0F studied by molecular dynamics

The response regulator protein Spo0F acts as an intermediate phospho-messenger in the signal transduction pathway that controls initiation of the differentiation process of sporculation in the bacterium Bacillus subtilis. The regulatory domain of Spo0F contains a triad of three conserved aspartate residues, whereof one aspartate (Asp54) is phosphorylated. Using molecular dynamics simulations, we have studied the changes in flexibility induced by phosphorylation and estimated the free energy cost of introducing a phosphate group at this position using alchemical free energy calculations. The deduced conformational flexibility compares well with experimental NMR results. We find that the apo-conformation of the protein explores a rough energy landscape resulting in a broad population of conformational substates. Phosphorylation of Spo0F reduces protein flexibility, and in particular, the so-called anchor and recognition regions exhibit lower mobility relative to the apo-conformation. Phosphorylation of Asp54 (P-Asp54), in which the apo-structure coordinates to the magnesium ion, results in extension of the sidechain, and depending on which carboxylate oxygen is phosphorylated, distinct interactions between P-Asp54 and magnesium ion as well as residues in its proximity are established. However, phosphorylation does not affect the coordination number of the magnesium ion yielding, within the statistical uncertainties, the same free energy change.

Molecular basis of phospholipase A2 activity toward phospholipids with sn-1 substitutions

We studied secretory phospholipase A(2) type IIA (sPLA(2)) activity toward phospholipids that are derivatized in the sn-1 position of the glycerol backbone. We explored what type of side group (small versus bulky groups, hydrophobic versus polar groups) can be introduced at the sn-1 position of the glycerol backbone of glycerophospholipids and at the same time be hydrolyzed by sPLA(2). The biophysical characterization revealed that the modified phospholipids can form multilamellar vesicles, and several of the synthesized sn-1 functionalized phospholipids were hydrolyzed by sPLA(2). Molecular dynamics simulations provided detailed insight on an atomic level that can explain the observed sPLA(2) activity toward the different phospholipid analogs. The simulations revealed that, depending on the nature of the side chain located at the sn-1 position, the group may interfere with an incoming water molecule that acts as the nucleophile in the enzymatic reaction. The simulation results are in agreement with the experimentally observed sPLA(2) activity toward the different phospholipid analogs.

Secretory phospholipase A2 hydrolysis of phospholipid analogues is dependent on water accessibility to the active site

A new and unnatural type of phospholipids with the head group attached to the 2-position of the glycerol backbone has been synthesized and shown to be a good substrate for secretory phospholipase A2 (sPLA2). To investigate the unexpected sPLA2 activity, we have compared three different phospholipids by using fluorescence techniques and HPLC, namely: (R)-1,2-dipalmitoyl-glycero-3-phosphocholine (hereafter referred to as 1R), (R)-1-O-hexadecyl-2-palmitoyl-glycero-3-phoshocholine (2R), and (S)-1-O-hexadecyl-3-palmitoyl-glycero-2-phosphocholine (3S). Furthermore, to understand the underlying mechanisms for the observed differences, we have performed molecular dynamics simulations to clarify on a structural level the substrate specificity of sPLA2 toward phospholipid analogues with their head groups in the 2-position of the glycerol backbone. We have studied the lipids above 1R, 2R, and 3S as well as their enantiomers 1S, 2S, and 3R. In the simulations of sPLA2-1S and sPLA2-3R, structural distortion in the binding cleft induced by the phospholipids showed that these are not substrates for sPLA2. In the case of the phospholipids 1R, 2R, and 3S, our simulations revealed that the difference observed experimentally in sPLA2 activity might be caused by reduced access of water molecules to the active site. We have monitored the number of water molecules that enter the active site region for the different sPLA2-phospholipid complexes and found that the probability of a water molecule reaching the correct position such that hydrolysis can occur is reduced for the unnatural lipids. The relative water count follows 1R > 2R > 3S. This is in good agreement with experimental data that indicate the same trend for sPLA2 activity: 1R > 2R > 3S.

Specificity of trypsin and chymotrypsin: loop-motion-controlled dynamic correlation as a determinant

Trypsin and chymotrypsin are both serine proteases with high sequence and structural similarities, but with different substrate specificity. Previous experiments have demonstrated the critical role of the two loops outside the binding pocket in controlling the specificity of the two enzymes. To understand the mechanism of such a control of specificity by distant loops, we have used the Gaussian network model to study the dynamic properties of trypsin and chymotrypsin and the roles played by the two loops. A clustering method was introduced to analyze the correlated motions of residues. We have found that trypsin and chymotrypsin have distinct dynamic signatures in the two loop regions, which are in turn highly correlated with motions of certain residues in the binding pockets. Interestingly, replacing the two loops of trypsin with those of chymotrypsin changes the motion style of trypsin to chymotrypsin-like, whereas the same experimental replacement was shown necessary to make trypsin have chymotrypsin's enzyme specificity and activity. These results suggest that the cooperative motions of the two loops and the substrate-binding sites contribute to the activity and substrate specificity of trypsin and chymotrypsin.

Specificity in lipases: a computational study of transesterification of sucrose

Computational conformational searches of putative transition states of the reaction of sucrose with vinyl laurate catalyzed by lipases from Candida antarctica B and Thermomyces lanuginosus have been carried out. The dielectric of the media have been varied to understand the role of protein plasticity in modulating the observed regioselective transesterification. The binding pocket of lipase from Candida adapts to the conformational variability of the various substates of the substrates by small, local adjustments within the binding pocket. In contrast, the more constrained pocket of the lipase from Thermomyces adapts by adjusting through concerted global motions between subdomains. This leads to the identification of one large pocket in Candida that accommodates both the sucrose and the lauroyl moieties of the transition state, whereas in Thermomyces the binding pocket is smaller, leading to the localization of the two moieties in two distinct pockets; this partly rationalizes the broader specificity of the former relative to the latter. Mutations have been suggested to exploit the differences towards changing the observed selectivities.