Dynamics of proteins in different solvent systems: analysis of essential motion in lipases

Solvent-dependent activity of Candida antarctica lipase B and its correlation with a regioselective mono aza-Michael addition - experimental and molecular dynamics simulation studies

With the aim of gaining understanding of the molecular basis of commercially available Candida antarctica lipase B (CALB) immobilized on polyacrylic resin catalyzed regioselective mono aza-Michael addition of Benzhydrazide to Diethyl maleate we decided to carry out molecular dynamics (MD) simulation studies in parallel with our experimental study. We found a correlation between the activity of CALB and the choice of solvent. Our study showed that solvent affects the performance of the enzyme due to the binding of solvent molecules to the enzyme active site region, and the solvation energy of substrates in the different solvents. We also found that CALB is only active in nonpolar solvent (i.e. Hexane), and therefore we investigated the influence of Hexane on the catalytic activity of CALB for the reaction. The results of this study and related experimental validation from our studies have been discussed here.

The role of plant cell wall encapsulation and porosity in regulating lipolysis during the digestion of almond seeds

Intact cell walls of almond prevent lipase penetration thus hindering lipid digestion.

Needle-Free Dermal Delivery of a Diphtheria Toxin CRM197 Mutant on Potassium-Doped Hydroxyapatite Microparticles

Injections with a hypodermic needle and syringe (HNS) are the current standard of care globally, but the use of needles is not without limitation. While a plethora of needle-free injection devices exist, vaccine reformulation is costly and presents a barrier to their widespread clinical application. To provide a simple, needle-free, and broad-spectrum protein antigen delivery platform, we developed novel potassium-doped hydroxyapatite (K-Hap) microparticles with improved protein loading capabilities that can provide sustained local antigen presentation and release. K-Hap showed increased protein adsorption over regular hydroxyapatite (P < 0.001), good structural retention of the model antigen (CRM197) with 1% decrease in α-helix content and no change in β-sheet content upon adsorption, and sustained release in vitro. Needle-free intradermal powder inoculation with K-Hap-CRM197 induced significantly higher IgG1 geometric mean titers (GMTs) than IgG2a GMTs in a BALB/c mouse model (P < 0.001) and induced IgG titer levels that were not different from the current clinical standard (P > 0.05), namely, alum-adsorbed CRM197 by intramuscular (i.m.) delivery. The presented results suggest that K-Hap microparticles may be used as a novel needle-free delivery vehicle for some protein antigens.

Sodium hexadecyl sulfate as an interfacial substance adjusting the adsorption of a protein on carbon nanotubes

Carbon nanotubes (CNTs) were functionalized with sodium hexadecyl sulfate (SHS). The lysozyme adsorbed on the SHS-CNTs exhibited a higher activity than that immobilized on the nonfunctionalized CNTs. To explain the experimental results and explore the mechanism of lysozyme adsorption, large-scale molecular dynamics simulations have been performed for a four-component system, including lysozyme, SHS, CNTs in explicit water. It has been found that the assembled SHS molecules form a soft layer on the surface of CNTs. The interactions between lysozyme and SHS induce the rearrangement of SHS molecules, forming a saddle-like structure on the CNT surface. The saddle-like structure fits the shape of the lysozyme, and the active-site cleft of the lysozyme is exposed to the water phase. Whereas, for the lysozyme adsorbed on the nonfunctionalized CNT, due to the hydrophobic interactions, the active-site cleft of the enzyme tends to face the wall of the CNT. The results of this work demonstrate that the SHS molecules as the interfacial substance have a function of adjusting the lysozyme with an appropriate orientation, which is favorable for the lysozyme having a higher activity.

Catalytic behavior of lipase immobilized onto congo red and PEG-decorated particles

Poly(ethylene glycol) (PEG)-decorated polystyrene (PS) nanoparticles with mean hydrodynamic diameter (D) and zeta–potential (ζ) of (286 ± 15) nm and (−50 ± 5) mV, respectively, were modified by the adsorption of Congo red (CR). The PS/PEG/CR particles presented D and ζ values of (290 ± 19) nm and (−36 ± 5) mV, respectively. The adsorption of lipase onto PS/PEG or PS/PEG/CR particles at (24 ± 1) °C and pH 7 changed the mean D value to (380 ± 20) and (405 ± 11) nm, respectively, and ζ value to (−32 ± 4) mV and (−25 ± 2) mV, respectively. The kinetic parameters of the hydrolysis of p-nitrophenyl butyrate were determined for free lipase, lipase immobilized onto PS/PEG and PS/PEG/CR particles. Lipase on PS/PEG/CR presented the largest Michaelis-Menten constant (K_M), but also the highest V_max and k_cat values. Moreover, it could be recycled seven times, losing a maximum 10% or 30% of the original enzymatic activity at 40 °C or 25 °C, respectively. Although lipases immobilized onto PS/PEG particles presented the smallest K_M values, the reactions were comparatively the slowest and recycling was not possible. Hydrolysis reactions performed in the temperature range of 25 °C to 60 °C with free lipases and lipases immobilized onto PS/PEG/CR particles presented an optimal temperature at 40 °C. At 60 °C free lipases and lipases immobilized onto PS/PEG/CR presented ~80% and ~50% of the activity measured at 40 °C, indicating good thermal stability. Bioconjugation effects between CR and lipase were evidenced by circular dichroism spectroscopy and spectrophotometry. CR molecules mediate the open state conformation of the lipase lid and favor the substrate approaching.

New insights in the activation of human cholesterol esterase to design potent anti-cholesterol drugs

Primary hypercholesterolemia is the root cause for major health issues like coronary heart disease and atherosclerosis. Regulating plasma cholesterol level, which is the product of biosynthesis as well as dietary intake, has become one of the major therapeutic strategies to effectively control these diseases. Human cholesterol esterase (hCEase) is an interesting target involved in the regulation of plasma cholesterol level and thus inhibition of this enzyme is highly effective in the treatment of hypercholesterolemia. This study was designed to understand the activation mechanism that enables the enzyme to accommodate long chain fatty acids and to identify the structural elements for the successful catalysis. Primarily the activation efficiencies of three different bile salts were studied and compared using molecular dynamics simulations. Based on the conformations of major surface loops, hydrogen bond interactions, and distance analyses, taurocholate was concluded as the preferred activator of the enzyme. Furthermore, the importance of two bile salt binding sites (proximal and remote) and the crucial role of 7α-OH group of the bile salts in the activation of hCEase was examined and evidenced. The results of our study explain the structural insights of the activation mechanism and show the key features of the bile salts responsible for the enzyme activation which are very useful in hypolipidemic drug designing strategies.

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.

Study of Thermomyces lanuginosa lipase in the presence of tributyrylglycerol and water

The Thermomyces lanuginosa lipase has been extensively studied in industrial and biotechnological research because of its potential for triacylglycerol transformation. This protein is known to catalyze both hydrolysis at high water contents and transesterification in quasi-anhydrous conditions. Here, we investigated the Thermomyces lanuginosa lipase structure in solution in the presence of a tributyrin aggregate using 30 ns molecular-dynamics simulations. The water content of the active-site groove was modified between the runs to focus on the protein-water molecule interactions and their implications for protein structure and protein-lipid interactions. The simulations confirmed the high plasticity of the lid fragment and showed that lipid molecules also bind to a secondary pocket beside the lid. Together, these results strongly suggest that the lid plays a role in the anchoring of the protein to the aggregate. The simulations also revealed the existence of a polar channel that connects the active-site groove to the outside solvent. At the inner extremity of this channel, a tyrosine makes hydrogen bonds with residues interacting with the catalytic triad. This system could function as a pipe (polar channel) controlled by a valve (the tyrosine) that could regulate the water content of the active site.

Modeling of solvent-dependent conformational transitions in Burkholderia cepacia lipase

BACKGROUND

The characteristic of most lipases is the interfacial activation at a lipid interface or in non-polar solvents. Interfacial activation is linked to a large conformational change of a lid, from a closed to an open conformation which makes the active site accessible for substrates. While for many lipases crystal structures of the closed and open conformation have been determined, the pathway of the conformational transition and possible bottlenecks are unknown. Therefore, molecular dynamics simulations of a closed homology model and an open crystal structure of Burkholderia cepacia lipase in water and toluene were performed to investigate the influence of solvents on structure, dynamics, and the conformational transition of the lid.

RESULTS

The conformational transition of B. cepacia lipase was dependent on the solvent. In simulations of closed B. cepacia lipase in water no conformational transition was observed, while in three independent simulations of the closed lipase in toluene the lid gradually opened during the first 10-15 ns. The pathway of conformational transition was accessible and a barrier was identified, where a helix prevented the lid from opening to the completely open conformation. The open structure in toluene was stabilized by the formation of hydrogen bonds.In simulations of open lipase in water, the lid closed slowly during 30 ns nearly reaching its position in the closed crystal structure, while a further lid opening compared to the crystal structure was observed in toluene. While the helical structure of the lid was intact during opening in toluene, it partially unfolded upon closing in water. The closing of the lid in water was also observed, when with eight intermediate structures between the closed and the open conformation as derived from the simulations in toluene were taken as starting structures. A hydrophobic beta-hairpin was moving away from the lid in all simulations in water, which was not observed in simulations in toluene. The conformational transition of the lid was not correlated to the motions of the beta-hairpin structure.

CONCLUSION

Conformational transitions between the experimentally observed closed and open conformation of the lid were observed by multiple molecular dynamics simulations of B. cepacia lipase. Transitions in both directions occurred without applying restraints or external forces. The opening and closing were driven by the solvent and independent of a bound substrate molecule.

Modeling hydration mechanisms of enzymes in nonpolar and polar organic solvents

A comprehensive study of the hydration mechanism of an enzyme in nonaqueous media was done using molecular dynamics simulations in five organic solvents with different polarities, namely, hexane, 3-pentanone, diisopropyl ether, ethanol, and acetonitrile. In these solvents, the serine protease cutinase from Fusarium solani pisi was increasingly hydrated with 12 different hydration levels ranging from 5% to 100% (w/w) (weight of water/weight of protein). The ability of organic solvents to 'strip off' water from the enzyme surface was clearly dependent on the nature of the organic solvent. The rmsd of the enzyme from the crystal structure was shown to be lower at specific hydration levels, depending on the organic solvent used. It was also shown that organic solvents determine the structure and dynamics of water at the enzyme surface. Nonpolar solvents enhance the formation of large clusters of water that are tightly bound to the enzyme, whereas water in polar organic solvents is fragmented in small clusters loosely bound to the enzyme surface. Ions seem to play an important role in the stabilization of exposed charged residues, mainly at low hydration levels. A common feature is found for the preferential localization of water molecules at particular regions of the enzyme surface in all organic solvents: water seems to be localized at equivalent regions of the enzyme surface independently of the organic solvent employed.

Structure and dynamics of Candida rugosa lipase: the role of organic solvent

The effect of organic solvent on the structure and dynamics of proteins was investigated by multiple molecular dynamics simulations (1 ns each) of Candida rugosa lipase in water and in carbon tetrachloride. The choice of solvent had only a minor structural effect. For both solvents the open and the closed conformation of the lipase were near to their experimental X-ray structures (C(alpha) rms deviation 1-1.3 A). However, the solvents had a highly specific effect on the flexibility of solvent-exposed side chains: polar side chains were more flexible in water, but less flexible in organic solvent. In contrast, hydrophobic residues were more flexible in organic solvent, but less flexible in water. As a major effect solvent changed the dynamics of the lid, a mobile element involved in activation of the lipase, which fluctuated as a rigid body about its average position. While in water the deviations were about 1.6 A, organic solvent reduced flexibility to 0.9 A. This increase rigidity was caused by two salt bridges (Lys85-Asp284, Lys75-Asp79) and a stable hydrogen bond (Lys75-Asn 292) in organic solvent. Thus, organic solvents stabilize the lid but render the side chains in the hydrophobic substrate-binding site more mobile. [figure: see text]. Superimposition of open (black, PDB entry 1CRL) and closed (gray, PDB entry 1TRH) conformers of C. rugosa lipase. The mobile lid is indicated.

Protein structure and dynamics in nonaqueous solvents: insights from molecular dynamics simulation studies

Protein structure and dynamics in nonaqueous solvents are here investigated using molecular dynamics simulation studies, by considering two model proteins (ubiquitin and cutinase) in hexane, under varying hydration conditions. Ionization of the protein groups is treated assuming "pH memory," i.e., using the ionization states characteristic of aqueous solution. Neutralization of charged groups by counterions is done by considering a counterion for each charged group that cannot be made neutral by establishing a salt bridge with another charged group; this treatment is more physically reasonable for the nonaqueous situation, contrasting with the usual procedures. Our studies show that hydration has a profound effect on protein stability and flexibility in nonaqueous solvents. The structure becomes more nativelike with increasing values of hydration, up to a certain point, when further increases render it unstable and unfolding starts to occur. There is an optimal amount of water, approximately 10% (w/w), where the protein structure and flexibility are closer to the ones found in aqueous solution. This behavior can explain the experimentally known bell-shaped dependence of enzyme catalysis on hydration, and the molecular reasons for it are examined here. Water and counterions play a fundamental and dynamic role on protein stabilization, but they also seem to be important for protein unfolding at high percentages of bound water.

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

As an aid to understanding the influence of dynamic fluctuations during esterolytic catalysis, we follow protein flexibility at three different steps along the catalytic pathway from substrate binding to product clearance via a covalently attached inhibitor, which represents a transition-state mimic. We have applied a classical approach, using molecular dynamics simulations to monitor protein dynamics in the nanosecond regime. We filter out small amplitude fluctuations and focus on the anharmonic contributions to the overall dynamics. This 'essential dynamics' analysis reveals different modes of response along the pathway suggesting that binding, catalysis and product clearance occur along different energy surfaces. Motions in the enzyme with a covalently attached ligand are more complex and occur along several eigenvectors. The magnitudes of the fluctuations in these individual subspaces are significantly smaller than those observed for the substrate and product molecules, indicating that the energy surface is shallow and that a relatively large number of conformational substates are accessible. On the other hand, substrate binding and product release occur at distinct modes of the protein flexibility suggesting that these processes occur along rough energy surfaces with only a few minima. Detailed energetic analyses along the trajectories indicated that in all cases binding is dominated by van der Waals interactions. The carboxylate form of the product is stabilized by a tight hydrogen bond network involving in particular Ser82, which may be a potential cause of product inhibition. Considerations such as these should aid the understanding of mechanisms of substrate, inhibitor or product recognition and could become of importance in the design of new substrates or inhibitors for enzymes.

Orientation and conformation of a lipase at an interface studied by molecular dynamics simulations

Electron density profiles calculated from molecular dynamics trajectories are used to deduce the orientation and conformation of Thermomyces lanuginosa lipase and a mutant adsorbed at an air-water interface. It is demonstrated that the profiles display distinct fine structures, which uniquely characterize enzyme orientation and conformation. The density profiles are, on the nanosecond timescale, determined by the average enzyme conformation. We outline a computational scheme that from a single molecular dynamics trajectory allows for extraction of electron density profiles referring to different orientations of the lipase relative to an implicit interface. Profiles calculated for the inactive and active conformations of the lipase are compared with experimental electron density profiles measured by x-ray reflectivity for the lipase adsorbed at an air-water interface. The experimental profiles contain less fine structural information than the calculated profiles because the resolution of the experiment is limited by the intrinsic surface roughness of water. Least squares fits of the calculated profiles to the experimental profiles provide areas per adsorbed enzyme and suggest that Thermomyces lanuginosa lipase adsorbs to the air-water interface in a semiopen conformation with the lid oriented away from the interface.

Revisiting the structural flexibility of the complex p21(ras)-GTP: the catalytic conformation of the molecular switch II

The hydrolysis of GTP in p21(ras) triggers conformational changes that regulate the ras/ERK signaling pathway. An important active site residue is Gln61, which has been found to be mutated in 30% of human tumors. The dynamics of the active site conformation is studied by using molecular dynamics simulation of two independent structures of the GTP-bound uncomplexed enzyme. Two distinct conformations of the enzyme are observed, in which the side-chain residue Gln61 is in different orientations. Essential dynamics analysis is used to describe the essential motions in the transition between the two conformations. Results are compared with earlier simulations of p21(ras) and its complex with GTPase activating protein p21-GAP.

Influence of a lipid interface on protein dynamics in a fungal lipase

Lipases catalyze lipolytic reactions and for optimal activity they require a lipid interface. To study the effect of a lipid aggregate on the behavior of the enzyme at the interfacial plane and how the aggregate influences an attached substrate or product molecule in time and space, we have performed molecular dynamics simulations. The simulations were performed over 1 to 2 ns using explicit SPC water. The interaction energies between protein and lipid are mainly due to van der Waals contributions reflecting the hydrophobic nature of the lipid molecules. Estimations of the protonation state of titratable residues indicated that the negative charge on the fatty acid is stabilized by interactions with the titratable residues Tyr-28, His-143, and His-257. In the presence of a lipid patch, the active site lid opens wider than observed in the corresponding simulations in an aqueous environment. In that lid conformation, the hydrophobic residues Ile-85, Ile-89, and Leu-92 are embedded in the lipid patch. The behavior of the substrate or product molecule is sensitive to the environment. Entering and leaving of substrate molecules could be observed in presence of the lipid patch, whereas the product forms strong hydrogen bonds with Ser-82, Ser-144, and Trp-88, suggesting that the formation of hydrogen bonds may be an important contribution to the mechanism by which product inhibition might take place.

A model of the pressure dependence of the enantioselectivity of Candida rugosalipase towards (+/-)-menthol

Transesterification of (+/-)-menthol using propionic acid anhydride and Candida rugosa lipase was performed in chloroform and water at different pressures (1, 10, 50, and 100 bar) to study the pressure dependence of enantioselectivity E. As a result, E significantly decreased with increasing pressure from E = 55 (1 bar) to E = 47 (10 bar), E = 37 (50 bar), and E = 9 (100 bar). To rationalize the experimental findings, molecular dynamics simulations of Candida rugosa lipase were carried out. Analyzing the lipase geometry at 1, 10, 50, and 100 bar revealed a cavity in the Candida rugosa lipase. The cavity leads from a position on the surface distinct from the substrate binding site to the core towards the active site, and is limited by F415 and the catalytic H449. In the crystal structure of the Candida rugosa lipase, this cavity is filled with six water molecules. The number of water molecules in this cavity gradually increased with increasing pressure: six molecules in the simulation at 1 bar, 10 molecules at 10 bar, 12 molecules at 50 bar, and 13 molecules at 100 bar. Likewise, the volume of the cavity progressively increased from about 1864 A(3) in the simulation at 1 bar to 2529 A(3) at 10 bar, 2526 A(3) at 50 bar, and 2617 A(3) at 100 bar. At 100 bar, one water molecule slipped between F415 and H449, displacing the catalytic histidine side chain and thus opening the cavity to form a continuous water channel. The rotation of the side chain leads to a decreased distance between the H449-N epsilon and the (+)-menthyl-oxygen (nonpreferred enantiomer) in the acyl enzyme intermediate, a factor determining the enantioselectivity of the lipase. Although the geometry of the preferred enantiomer is similar in all simulations, the geometry of the nonpreferred enantiomer gets gradually more reactive. This observation correlates with the gradually decreasing enantioselectivity E.

Dynamics of the substrate binding pocket in the presence of an inhibitor covalently attached to a fungal lipase

To gain insight into the mobility of the occupied ligand-binding pocket of the Rhizomucor miehei lipase we have conducted a rigorous molecular dynamics analysis. The covalently attached inhibitor, ethylhexylphosphonate, was employed as a mimic of the putative tetrahedral intermediate in the esterolytic reaction. Our results show that in this lipase, ligand recognition is influenced by the flexibility of the binding pocket, a feature that is common to many other enzymes. Several regions around the active site were found to move significantly to adapt to the inhibitor. These motions are correlated to the flexibility of the inhibitor. In particular, the hexyl chain of the inhibitor shows considerable mobility, and adjacent residues in the binding cleft accommodate to this flexibility. Pronounced fluctuations in the binding pocket induced by the flexibility of the inhibitor are observed in the hinge region F79-S82, the active site loop region W88-V95 and the protein regions P209-F215/H257-Y260. The flexibility in the regions F79-S82 and H257-Y260, where the shorter ethyl chain is located, indicates that additional space in this binding cleft region is available for accommodating a larger moiety. Fluctuations in the region W88-V95 and P209-F215 are due to the relatively short flexible hexyl carbon chain. This part of the binding pocket could be stiffened by the presence of a longer carbon chain. Though the inhibitor is covalently attached through the phosphonate moiety, interaction of the remainder of the molecule and the enzyme are determined by hydrophobic interactions, where the Van der Waals energies are approximately 25% lower than the electrostatic contributions.

Concerted motions in copper plastocyanin and azurin: an essential dynamics study

Essential dynamics analysis of molecular dynamics simulation trajectories (1.1 ns) of two copper containing electron transfer proteins, plastocyanin and azurin, has been performed. The protein essential modes have been analysed in order to identify large concerted motions which could be relevant for the electron transfer function exerted by these proteins. The analysis, conducted for temporal windows of different lengths along the protein trajectories, shows a rapid convergence and indicates that for both the proteins the predominant internal motions occur in a subspace of only a few degrees of freedom. Moreover, it is found that for both the proteins the likely binding sites (i.e. the hydrophobic and negative patches) with the reaction partners move in a concerted fashion with a few structural regions far from the active site. Such results are discussed in connection with the possible involvement of large concerted motions in the recognition and binding interaction with physiological electron transfer partners.

Computational studies of essential dynamics of Pseudomonas cepacia lipase

In order to investigate the interfacial activation of a lipase from Pseudomonas cepacia (PcL), molecular dynamics (MD) simulations and essential dynamics (ED) analysis were performed in different solvent environments: vacuum and explicit water solvents. Starting from the active (open) structure of PcL, the essential dynamics analysis of the simulations revealed large correlated motions, which may be responsible for the activation of the enzyme. Fluctuations in the U1 (active-site lid) and U2 domains are found to be important in the activation of PcL. In contrast, the catalytic triad exhibits very little displacement. These results are consistent with the previous X-ray structural determination. A combined analysis of the trajectories showed some differences for the simulations in different solvent environments. It was found that the region around the helix alpha5 showed larger displacements in the water simulations. It can be concluded that the open structure of PcL becomes unstable in water solvents, leading to the closing of the so-called 'lid' region. The simulations and ED analysis on the closed structure of PgL provided additional information concerning the structural changes involved in the activation of the lipases. It was found that structural changes for PcL and PgL, which are responsible for the essential motions of the protein, showed contrasting behavior in the different solvent environments.

Directed evolution of an enantioselective lipase

BACKGROUND

The biocatalytic production of enantiopure compounds is of steadily increasing importance to the chemical and biotechnological industry. In most cases, however, it is impossible to identify an enzyme that possesses the desired enantioselectivity. Therefore, there is a strong need to create by molecular biological methods novel enzymes which display high enantioselectivity.

RESULTS

A bacterial lipase from Pseudomonas aeruginosa (PAL) was evolved to catalyze with high enantioselectivity the hydrolysis of the chiral model substrate 2-methyldecanoic acid p-nitrophenyl ester. Successive rounds of random mutagenesis by ep-PCR and saturation mutagenesis resulted in an increase in enantioselectivity from E=1.1 for the wild-type enzyme to E=25.8 for the best variant which carried five amino acid substitutions. The recently solved three-dimensional structure of PAL allowed us to analyze the structural consequences of these substitutions.

CONCLUSIONS

A highly enantioselective lipase was created by increasing the flexibility of distinct loops of the enzyme. Our results demonstrate that enantioselective enzymes can be created by directed evolution, thereby opening up a large area of novel applications in biotechnology.

Theoretical studies of the response of a protein structure to cavity-creating mutations

We have investigated the response of a protein structure to cavity-creating mutations by molecular dynamics (MD) simulations for the wild-type and the five mutants of phage T4 lysozyme. Essential dynamics (ED) analysis and the methods for calculating different components of local interaction energies are used to examine the structural and energetic characteristics associated with the mutations. In agreement with the x-ray results, it is found that the structural changes due to the replacements of a bulky side chain such as Leu or Phe with Ala within the hydrophobic core can be characterized as slight adjustments rather than substantial reorganization of the protein. The relative stability of different mutant structures can be related with the extent of structural readjustments in response to the mutation. The destabilization of the mutant Leu-->Ala proteins relative to the wild-type is closely related with the loss of van der Waals contacts due to the cavity-creating mutations.

Conformational substates in different crystal forms of the photoactive yellow protein--correlation with theoretical and experimental flexibility

The conformational changes during the photocycle of the photoactive yellow protein have been the subject of many recent studies. Spectroscopic measurements have shown that the photocycle also occurs in a crystalline environment, and this has been the basis for subsequent Laue diffraction and cryocrystallographic studies. These studies have shown that conformational changes during the photocycle are limited to the chromophore and its immediate environment. However, spectroscopic studies suggest the presence of large conformational changes in the protein. Here, we address this apparent discrepancy in two ways. First, we obtain a description of large concerted motions in the ground state of the yellow protein from NMR data and theoretical calculations. Second, we describe the high-resolution structure of the yellow protein crystallized in a different space group. The structure of the yellow protein differs significantly between the two crystal forms. We show that these differences can be used to obtain a description of the flexibility of the protein that is consistent with the motions observed in solution.

Computational analysis of chain flexibility and fluctuations in Rhizomucor miehei lipase

We have performed molecular dynamics simulation of Rhizomucor miehei lipase (Rml) with explicit water molecules present. The simulation was carried out in periodic boundary conditions and conducted for 1. 2 ns in order to determine the concerted protein dynamics and to examine how well the essential motions are preserved along the trajectory. Protein motions are extracted by means of the essential dynamics analysis method for different lengths of the trajectory. Motions described by eigenvector 1 converge after approximately 200 ps and only small changes are observed with increasing simulation time. Protein dynamics along eigenvectors with larger indices, however, change with simulation time and generally, with increasing eigenvector index, longer simulation times are required for observing similar protein motions (along a particular eigenvector). Several regions in the protein show relatively large fluctuations and in particular motions in the active site lid and the segments Thr57-Asn63 and the active site hinge region Pro101-Gly104 are seen along several eigenvectors. These motions are generally associated with glycine residues, while no direct correlations are observed between these fluctuations and the positioning of prolines in the protein structure. The partial opening/closing of the lid is an example of induced fit mechanisms seen in other enzymes and could be a general mechanism for the activation of Rml.

Molecular dynamics simulations of protein-tyrosine phosphatase 1B. I. ligand-induced changes in the protein motions

Activity of enzymes, such as protein tyrosine phosphatases (PTPs), is often associated with structural changes in the enzyme, resulting in selective and stereospecific reactions with the substrate. To investigate the effect of a substrate on the motions occurring in PTPs, we have performed molecular dynamics simulations of PTP1B and PTP1B complexed with a high-affinity peptide DADEpYL, where pY stands for phosphorylated tyrosine. The peptide sequence is derived from the epidermal growth factor receptor (EGFR988-993). Simulations were performed in water for 1 ns, and the concerted motions in the protein were analyzed using the essential dynamics technique. Our results indicate that the predominately internal motions in PTP1B occur in a subspace of only a few degrees of freedom. Upon substrate binding, the flexibility of the protein is reduced by approximately 10%. The largest effect is found in the protein region, where the N-terminal of the substrate is located, and in the loop region Val198-Gly209. Displacements in the latter loop are associated with the motions in the WPD loop, which contains a catalytically important aspartic acid. Estimation of the pKa of the active-site cysteine along the trajectory indicates that structural inhomogeneity causes the pKa to vary by approximately +/-1 pKa unit. In agreement with experimental observations, the active-site cysteine is negatively charged at physiological pH.

A solvent model for simulations of peptides in bilayers. I. Membrane-promoting alpha-helix formation

We describe an efficient solvation model for proteins. In this model atomic solvation parameters imitating the hydrocarbon core of a membrane, water, and weak polar solvent (octanol) were developed. An optimal number of solvation parameters was chosen based on analysis of atomic hydrophobicities and fitting experimental free energies of gas-cyclohexane, gas-water, and octanol-water transfer for amino acids. The solvation energy term incorporated into the ECEPP/2 potential energy function was tested in Monte Carlo simulations of a number of small peptides with known energies of bilayer-water and octanol-water transfer. The calculated properties were shown to agree reasonably well with the experimental data. Furthermore, the solvation model was used to assess membrane-promoting alpha-helix formation. To accomplish this, all-atom models of 20-residue homopolypeptides-poly-Leu, poly-Val, poly-Ile, and poly-Gly in initial random coil conformation-were subjected to nonrestrained Monte Carlo conformational search in vacuo and with the solvation terms mimicking the water and hydrophobic parts of the bilayer. All the peptides demonstrated their largest helix-forming tendencies in a nonpolar environment, where the lowest-energy conformers of poly-Leu, Val, Ile revealed 100, 95, and 80% of alpha-helical content, respectively. Energetic and conformational properties of Gly in all environments were shown to be different from those observed for residues with hydrophobic side chains. Applications of the solvation model to simulations of peptides and proteins in the presence of membrane, along with limitations of the approach, are discussed.

Analysis of the dynamics of rhizomucor miehei lipase at different temperatures

The dynamics of Rhizomucor miehei lipase has been studied by molecular dynamics simulations at temperatures ranging from 200-500K. Simulations carried out in periodic boundary conditions and using explicit water molecules were performed for 400 ps at each temperature. Our results indicate that conformational changes and internal motions in the protein are significantly influenced by the temperature increase. With increasing temperature, the number of internal hydrogen bonds decreases, while surface accessibility, radius of gyration and the number of residues in random coil conformation increase. In the temperature range studied, the motions can be described in a low dimensional subspace, whose dimensionality decreases with increasing temperature. Approximately 80% of the total motion is described by the first (i) 80 eigenvectors at T=200K, (ii) 30 eigenvectors at T=300K and (iii) 10 eigenvectors at T=400K. At high temperature, the alpha-helix covering the active site in the native Rhizomucor miehei lipase, the helix at which end the active site is located, and in particular, the loop (Gly35-Lys50) show extensive flexibility.

Assignment of side-chain conformation using adiabatic energy mapping, free energy perturbation, and molecular dynamic simulations

NMR spectroscopic analysis of the C-terminal Kunitz domain fragment (alpha3(VI)) from the human alpha3-chain of type VI collagen has revealed that the side chain of Trp21 exists in two unequally populated conformations. The major conformation (M) is identical to the conformation observed in the X-ray crystallographic structure, while the minor conformation (m) cannot structurally be resolved in detail by NMR due to insufficient NOE data. In the present study, we have applied: (1) rigid and adiabatic mapping, (2) free energy simulations, and (3) molecular dynamic simulations to elucidate the structure of the m conformer and to provide a possible pathway of the Trp21 side chain between the two conformers. Adiabatic energy mapping of conformations of the Trp21 side chain obtained by energy minimization identified two energy minima: One corresponding to the conformation of Trp21 observed in the X-ray crystallographic structure and solution structure of alpha3(VI) (the M conformation) and the second corresponding to the m conformation predicted by NMR spectroscopy. A transition pathway between the M and m conformation is suggested. The free-energy difference between the two conformers obtained by the thermodynamic integration method is calculated to 1.77+/-0.7 kcal/mol in favor of the M form, which is in good agreement with NMR results. Structural and dynamic properties of the major and minor conformers of the alpha3(VI) molecule were investigated by molecular dynamic. Essential dynamics analysis of the two resulting 800 ps trajectories reveals that when going from the M to the m conformation only small, localized changes in the protein structure are induced. However, notable differences are observed in the mobility of the binding loop (residues Thr13-Ile18), which is more flexible in the m conformation than in the M conformation. This suggests that the reorientation of Trp2 might influence the inhibitory activity against trypsin, despite the relative large distance between the binding loop and Trp21.

Conformational change in the activation of lipase: an analysis in terms of low-frequency normal modes

The interfacial activation of Rhizomucor miehei lipase (RmL) involves the motion of an alpha-helical region (residues 82-96) which acts as a "lid" over the active site of the enzyme, undergoing a displacement from a "closed" to an "open" conformation upon binding of substrate. Normal mode analyses performed in both low and high dielectric media reveal that low-frequency vibrational modes contribute significantly to the conformational transition between the closed and open conformations. In these modes, the lid displacement is coupled to local motions of active site loops as well as global breathing motions. Atomic fluctuations of the first hinge of the lid (residues 83-84) are substantially larger in the low dielectric medium than in the high dielectric medium. Our results also suggest that electrostatic interactions of Arg86 play an important role in terms of both the intrinsic stability of the lid and its displacement, through enhancement of hinge mobility in a high dielectric medium. Additional calculations demonstrate that the observed patterns of atomic fluctuations are an intrinsic feature of the protein structure and not dependent on the nature of specific energy minima.

Enzymes in low water systems

Water is fundamental for enzyme action and for formation of the three-dimensional structure of proteins. Hence, it may be assumed that studies on the interplay between water and enzymes can yield insight into enzyme function and formation. This has proven correct, because the numerous studies that have been made on the behavior of water-soluble and membrane enzymes in systems with a low water content (reverse micelles or enzymes suspended in nonpolar organic solvents) have revealed properties of enzymes that are not easily appreciated in aqueous solutions. In the low water systems, it has been possible to probe the relation between solvent and enzyme kinetics, as well as some of the factors that affect enzyme thermostability and catalysis. Furthermore, the studies show that low water environments can be used to stabilize conformers that exhibit unsuspected catalytic properties, as well as intermediates of enzyme function and formation that in aqueous media have relatively short life-times. The structure of enzymes in these unnatural conditions is actively being explored.