Structure of a ternary complex of proteinase K,mercury, and a substrate-analogue hexa-peptide at 2.2 Å resolution

The crystal structure of a ternary complex of proteinase K, Hg(II) and a hexapeptide N-Ac-Pro-Ala-Pro-Phe-Pro-Ala-NH₂ has been determined at 2.2 Å resolution and refined to an R factor of 0.172 for 12,910 reflections. The mercury atom occupies two alternate sites, each of which was assigned an occupancy of 0.45. These two sites are bridged by Cys-73 Sγ which forms covalent bonds to both. Both mercury sites form regular polyhedrons involving atoms from residues Asp-39, His-69, Cys-73, His-72, Met-225, and Wat-324. The complex formation with mercury seems to disturb the stereochemistry of the residues of the catalytic triad Asp-39, His-69, and Ser-224 appreciably, thus reducing the enzymatic activity of proteinase K to 15%. The electron density in the difference Fourier map shows that the hexapeptide occupies the S1 subsite predominantly and the standard recognition site constituted by Ser-132 to Gly-136 and Gly-100 to Tyr-104 segments is virtually empty. The hexapeptide is held firmly through a series of hydrogen bonds involving protein atoms and water molecules. As a result of complex formation, Asp-39, His-69, Met-225, Ile-220, Ser-219, Thr-223, and Ser-224 residues move appreciably to accommodate the mercury atoms and the hexapeptide. The largest movement is observed for Met-225 which is involved in multiple interactions with both mercury and the hexapeptide. The activity results indicate an inhibition rate of 95%, as a result of the combined effect of mercury and hexapeptide. © 1996 Wiley-Liss, Inc.     

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.     

Intramembrane proton binding site linked to activation of bacterial pentameric ion channel

Prokaryotic orthologs of eukaryotic Cys-loop receptor channels recently emerged as structural and mechanistic surrogates to investigate this superfamily of intercellular signaling proteins. Here, we examine proton activation of the prokaryotic ortholog GLIC using patch clamp electrophysiology, mutagenesis, and molecular dynamics (MD) simulations. Whole-cell current recordings from human embryonic kidney (HEK) 293 cells expressing GLIC show half-maximal activation at pH 6, close to the pK(a) of histidine, implicating the three native His residues in proton sensing linked to activation. The mutation H235F abolishes proton activation, H277Y is without effect, and all nine mutations of His-127 prevent expression on the cell surface. In the GLIC crystal structure, His-235 on transmembrane (TM) α-helix 2, hydrogen bonds to the main chain carbonyl oxygen of Ile-259 on TM α-helix 3. MD simulations show that when His-235 is protonated, the hydrogen bond persists, and the channel remains in the open conformation, whereas when His-235 is deprotonated, the hydrogen bond dissociates, and the channel closes. Mutations of the proximal Tyr-263, which also links TM α-helices 2 and 3 via a hydrogen bond, alter proton sensitivity over a 1.5 pH unit range. MD simulations show that mutations of Tyr-263 alter the hydrogen bonding capacity of His-235. The overall findings show that His-235 in the TM region of GLIC is a novel proton binding site linked to channel activation.     

Active site dynamics of acyl-chymotrypsin

The motions of water molecules, the acyl moiety, the catalytic triad, and the oxyanion binding site of acyl-chymotrypsin were studied by means of a stochastic boundary molecular dynamics simulation. A water molecule that could provide the nucleophilic OH^? for the deacylation stage of the catalysis was found to be trapped between the imidazole ring of His-57 and the carbonyl carbon of the acyl group. It makes a hydrogen bond with the N^ε2 of His-57 and is heldin place through a network of hydrogen-bonded water molecules in theactive site. The water molecule was found as close as 2.8 Å to the carbonyl carbon. This appears to be due to the constraints imposed by nonbonded interaction in the active site. Configurations were found in which one hydrogen of the trapped water shared a bifurcated hydrogen bond with His-57-N^ε2 and Ser-195-0^γ with the water oxygen very close to the carbonyl carbon. The existence of such a water molecule suggests that large movement of the His-57 imidazole ring between positions suitable for providing general-base catalyzed assistance and for providing general-acid catalyzed assistance may notbe required during the reaction. The simulation indicates that the side chains of residues involved in catalysis (i.e., His-57, Ser-195, and Asp-102) are significantly less flexible than other side chains in the protein. The 40% reduction in rms fluctuations is consistent with a comparable reduction calculated from the temperature factors obtained in the X-ray crystal-lographic data of γ-chymotrypsin. The greater rigidity of active site residues seems to result from interconnected hydrogen bonding networks among the residues and between the residues and the solvent water in the active site. © Wiley-Liss, Inc.     

Crystal structure of human gastric lipase and model of lysosomal acid lipase, two lipolytic enzymes of medical interest.

Fat digestion in humans requires not only the classical pancreatic lipase but also gastric lipase, which is stable and active despite the highly acidic stomach environment. We report here the structure of recombinant human gastric lipase at 3.0-A resolution, the first structure to be described within the mammalian acid lipase family. This globular enzyme (379 residues) consists of a core domain belonging to the alpha/beta hydrolase-fold family and a "cap" domain, which is analogous to that present in serine carboxypeptidases. It possesses a classical catalytic triad (Ser-153, His-353, Asp-324) and an oxyanion hole (NH groups of Gln-154 and Leu-67). Four N-glycosylation sites were identified on the electron density maps. The catalytic serine is deeply buried under a segment consisting of 30 residues, which can be defined as a lid and belonging to the cap domain. The displacement of the lid is necessary for the substrates to have access to Ser-153. A phosphonate inhibitor was positioned in the active site that clearly suggests the location of the hydrophobic substrate binding site. The lysosomal acid lipase was modeled by homology, and possible explanations for some previously reported mutations leading to the cholesterol ester storage disease are given based on the present model.     

Monte Carlo simulations of the peptide recognition at the consensus binding site of the constant fragment of human immunoglobulin G: the energy landscape analysis of a hot spot at the intermolecular interface

Monte Carlo simulations of molecular recognition at the consensus binding site of the constant fragment (Fc) of human immunoglobulin G (Ig) protein have been performed to analyze structural and thermodynamic aspects of binding for the 13-residue cyclic peptide DCAWHLGELVWCT. The energy landscape analysis of a hot spot at the intermolecular interface using alanine scanning and equilibrium-simulated tempering dynamics with the simplified, knowledge-based energy function has enabled the role of the protein hot spot residues in providing the thermodynamic stability of the native structure to be determined. We have found that hydrophobic interactions between the peptide and the Met-252, Ile-253, His-433, and His-435 protein residues are critical to guarantee the thermodynamic stability of the crystallographic binding mode of the complex. Binding free energy calculations, using a molecular mechanics force field and a solvation energy model, combined with alanine scanning have been conducted to determine the energetic contribution of the protein hot spot residues in binding affinity. The conserved Asn-434, Ser-254, and Tyr-436 protein residues contribute significantly to the binding affinity of the peptide-protein complex, serving as an energetic hot spot at the intermolecular interface. The results suggest that evolutionary conserved hot spot protein residues at the intermolecular interface may be partitioned in fulfilling thermodynamic stability of the native binding mode and contributing to the binding affinity of the complex.     

Theoretical investigation of the dynamics of the active site lid in Rhizomucor miehei lipase.

Interfacial activation of Rhizomucor miehei lipase is accompanied by a hinge-type motion of a single helix (residues 83-94) that acts as a lid over the active site. Activation of the enzyme involves the displacement of the lid to expose the active site, suggesting that the dynamics of the lid could be of mechanistic and kinetic importance. To investigate possible activation pathways and to elucidate the effect of a hydrophobic environment (as would be provided by a lipid membrane) on the lid opening, we have applied molecular dynamics and Brownian dynamics techniques. Our results indicate that the lipase activation is enhanced in a hydrophobic environment. In nonpolar low-dielectric surroundings, the lid opens in approximately 100 ns in the BD simulations. In polar high-dielectric (aqueous) surroundings, the lid does not always open up in simulations of up to 900 ns duration, but it does exhibit some gating motion, suggesting that the enzyme molecule may exist in a partially active form before the catalytic reaction. The activation is controlled by the charged residues ARG86 and ASP91. In the inactive conformation, ASP91 experiences repulsive forces and pushes the lid toward the open conformation. Upon activation ARG86 approaches ASP61, and in the active conformation, these residues form a salt bridge that stabilizes the open conformation.     

Hydrogen-bonding propensities of sphingomyelin in solution and in a bilayer assembly: a molecular dynamics study

Sphingomyelin is enriched within lipid microdomains of the cell membrane termed lipid rafts. These microdomains play a part in regulating a variety of cellular events. Computer simulations of the hydrogen-bonding properties of sphingolipids, believed to be central to the organization of these domains, can delineate the possible molecular interactions that underlie this lipid structure. We have therefore used molecular dynamics simulations to unravel the hydrogen-bonding behavior of palmitoylsphingomyelin (PSM). A series of eight simulations of 3 ns each of a single PSM molecule in water showed that the sphingosine OH and NH groups can form hydrogen bonds with the phosphate oxygens of their own polar head, in agreement with NMR data. Simulations of PSM in a bilayer assembly were carried out for 8 ns with three different force field parameterizations. The major physico-chemical parameters of the simulated bilayer agree with those established experimentally. The sphingosine OH group was mainly involved in intramolecular hydrogen bonds, in contrast to the almost exclusive intermolecular hydrogen bonds formed by the amide NH moiety. During the bilayer simulations the intermolecular hydrogen bonds among lipids formed a dynamic network characterized by the presence of hydrogen-bonded lipid clusters of up to nine PSM molecules.     

Detailed mechanism for AmtB conducting NH4+/NH3: molecular dynamics simulations

The mechanism by which the ammonium transporter, AmtB, conducts NH4+/NH3 into the cytoplasm was investigated using conventional molecular dynamics (MD) simulations. These simulations revealed that the neutral molecule, NH3, passes automatically through the channel upon its arrival at the Am2 site and that the function of the channel as a one-way valve for passage of NH3 is not determined by the cytoplasmic exit gate but, rather, by the periplasmic entrance gate of the channel. The NH3, produced by deprotonation of NH4+ at the entrance gate, is spontaneously conveyed to the central region of the channel via a hydrogen-bond network comprising His-168, His-318, Tyr-32, and the NH3 molecule. Finally, the NH3 molecule exits the channel through the exit gate formed by Phe-31, Ile-266, Val-314, and His-318. In addition, Ser-263 is shown to play a critical role in the final stages, acting as a pivoting arm to shunt the NH3 molecule from the cytoplasmic exit gate of the channel out into the cytoplasm. This finding is further supported by another simulation which shows that NH3 fails to be translocated through the channel formed by the Ser-263-Ala mutation. Thus, this study casts new insights on the mechanism of AmtB-mediated passage of NH3 across cellular membranes.     

Refined crystal structure of cytoplasmic malate dehydrogenase at 2.5-A resolution

The molecular structure of cytoplasmic malate dehydrogenase from pig heart has been refined by alternating rounds of restrained least-squares methods and model readjustment on an interactive graphics system. The resulting structure contains 333 amino acids in each of the two subunits, 2 NAD molecules, 471 solvent molecules, and 2 large noncovalently bound molecules that are assumed to be sulfate ions. The crystallographic study was done on one entire dimer without symmetry restraints. Analysis of the relative position of the two subunits shows that the dimer does not obey exact 2-fold rotational symmetry; instead, the subunits are related by a 173 degrees rotation. The structure results in a R factor of 16.7% for diffraction data between 6.0 and 2.5 A, and the rms deviations from ideal bond lengths and angles are 0.017 A and 2.57 degrees, respectively. The bound coenzyme in addition to hydrophobic interactions makes numerous hydrogen bonds that either are directly between NAD and the enzyme or are with solvent molecules, some of which in turn are hydrogen bonded to the enzyme. The carboxamide group of NAD is hydrogen bonded to the side chain of Asn-130 and via a water molecule to the backbone nitrogens of Leu-157 and Asp-158 and to the carbonyl oxygen of Leu-154. Asn-130 is one of the corner residues in a beta-turn that contains the lone cis peptide bond in cytoplasmic malate dehydrogenase, situated between Asn-130 and Pro-131. The active site histidine, His-186, is hydrogen bonded from nitrogen ND1 to the carboxylate of Asp-158 and from its nitrogen NE2 to the sulfate ion bound in the putative substrate binding site. In addition to interacting with the active site histidine, this sulfate ion is also hydrogen bonded to the guanidinium group of Arg-161, to the carboxamide group of Asn-140, and to the hydroxyl group of Ser-241. It is speculated that the substrate, malate or oxaloacetate, is bound in the sulfate binding site with the substrate 1-carboxyl hydrogen bonded to the guanidinium group of Arg-161.     

Crystal structure of the binary complex of ribulose-1,5-bisphosphate carboxylase and its product, 3-phospho-D-glycerate

The crystal structure of the binary complex of non-activated ribulose-1,5-bisphosphate carboxylase/oxygenase from Rhodospirillum rubrum and its product 3-phospho-D-glycerate has been determined to 2.9-A resolution. This structure determination confirms the proposed location of the active site (Schneider, G., Lindqvist, Y., Br?ndén, C.-I., and Lorimer, G. (1986) EMBO J. 5, 3409-3415) at the carboxyl end of the beta-strands of the alpha/beta-barrel in the carboxyl-terminal domain. One molecule of 3-phosphoglycerate is bound per active site. All oxygen atoms of 3-phosphoglycerate form hydrogen bonds to groups of the enzyme. The phosphate group interacts with the sidechains of residues Arg-288, His-321, and Ser-368, which are conserved between enzymes from different species as well as with the main chain nitrogens from residues Thr-322 and Gly-323. These amino acid residues constitute one of the two phosphate binding sites of the active site. The carboxyl group interacts with the side chains of His-287, Lys-191, and Asn-111. Implications of the activation process for the binding of 3-phosphoglycerate are discussed.     

The interaction of phospholipase A2 with a phospholipid bilayer: coarse-grained molecular dynamics simulations

A number of membrane-active enzymes act in a complex environment formed by the interface between a lipid bilayer and bulk water. Although x-ray diffraction studies yield structures of isolated enzyme molecules, a detailed characterization of their interactions with the interface requires a measure of how deeply such a membrane-associated protein penetrates into a lipid bilayer. Here, we apply coarse-grained (CG) molecular dynamics (MD) simulations to probe the interaction of porcine pancreatic phospholipase A2 (PLA2) with a lipid bilayer containing palmitoyl-oleoyl-phosphatidyl choline and palmitoyl-oleoyl-phosphatidyl glycerol molecules. We also used a configuration from a CG-MD trajectory to initiate two atomistic (AT) MD simulations. The results of the CG and AT simulations are evaluated by comparison with available experimental data. The membrane-binding surface of PLA2 consists of a patch of hydrophobic residues surrounded by polar and basic residues. We show this proposed footprint interacts preferentially with the anionic headgroups of the palmitoyl-oleoyl-phosphatidyl glycerol molecules. Thus, both electrostatic and hydrophobic interactions determine the location of PLA2 relative to the bilayer. From a general perspective, this study demonstrates that CG-MD simulations may be used to reveal the orientation and location of a membrane-surface-bound protein relative to a lipid bilayer, which may subsequently be refined by AT-MD simulations to probe more detailed interactions.     

Crystal structure of the Y52F/Y73F double mutant of phospholipase A2: increased hydrophobic interactions of the phenyl groups compensate for the disrupted hydrogen bonds of the tyrosines.

The enzyme phospholipase A2 (PLA2) catalyzes the hydrolysis of the sn-2 ester bond of membrane phospholipids. The highly conserved Tyr residues 52 and 73 in the enzyme form hydrogen bonds to the carboxylate group of the catalytic Asp-99. These hydrogen bonds were initially regarded as essential for the interfacial recognition and the stability of the overall catalytic network. The elimination of the hydrogen bonds involving the phenolic hydroxyl groups of the Tyr-52 and -73 by changing them to Phe lowered the stability but did not significantly affect the catalytic activity of the enzyme. The X-ray crystal structure of the double mutant Y52F/Y73F has been determined at 1.93 A resolution to study the effect of the mutation on the structure. The crystals are trigonal, space group P3(1)21, with cell parameters a = b = 46.3 A and c = 102.95 A. Intensity data were collected on a Siemens area detector, 8,024 reflections were unique with an R(sym) of 4.5% out of a total of 27,203. The structure was refined using all the unique reflections by XPLOR to a final R-factor of 18.6% for 955 protein atoms, 91 water molecules, and 1 calcium ion. The root mean square deviation for the alpha-carbon atoms between the double mutant and wild type was 0.56 A. The crystal structure revealed that four hydrogen bonds were lost in the catalytic network; three involving the tyrosines and one involving Pro-68. However, the hydrogen bonds of the catalytic triad, His-48, Asp-99, and the catalytic water, are retained. There is no additional solvent molecule at the active site to replace the missing hydroxyl groups; instead, the replacement of the phenolic OH groups by H atoms draws the Phe residues closer to the neighboring residues compared to wild type; Phe-52 moves toward His-48 and Asp-99 of the catalytic diad, and Phe-73 moves toward Met-8, both by about 0.5 A. The closing of the voids left by the OH groups increases the hydrophobic interactions compensating for the lost hydrogen bonds. The conservation of the triad hydrogen bonds and the stabilization of the active site by the increased hydrophobic interactions could explain why the double mutant has activity similar to wild type. The results indicate that the aspartyl carboxylate group of the catalytic triad can function alone without additional support from the hydrogen bonds of the two Tyr residues.     

Human lipoprotein lipase. Analysis of the catalytic triad by site-directed mutagenesis of Ser-132, Asp-156, and His-241.

Lipoprotein lipase (LPL) plays a central role in normal lipid metabolism as the key enzyme involved in the hydrolysis of triglycerides present in chylomicrons and very low density lipoproteins. LPL is a member of a family of hydrolytic enzymes that include hepatic lipase and pancreatic lipase. Based on primary sequence homology of LPL to pancreatic lipase, Ser-132, Asp-156, and His-241 have been proposed to be part of a domain required for normal enzymic activity. We have analyzed the role of these potential catalytic residues by site-directed mutagenesis and expression of the mutant LPL in human embryonic kidney-293 cells. Substitution of Ser-132, Asp-156, and His-241 by several different residues resulted in the expression of an enzyme that lacked both triolein and tributyrin esterase activities. Mutation of other conserved residues, including Ser-97, Ser-307, Asp-78, Asp-371, Asp-440, His-93, and His-439 resulted in the expression of active enzymes. Despite their effect on LPL activity, substitutions of Ser-132, Asp-156, and His-241 did not change either the heparin affinity or lipid binding properties of the mutant LPL. In summary, mutation of Ser-132, Asp-156, and His-241 specifically abolishes total hydrolytic activity without disrupting other important functional domains of LPL. These combined results strongly support the conclusion that Ser-132, Asp-156, and His-241 form the catalytic triad of LPL and are essential for LPL hydrolytic activity.     

Structural basis of phospholipase A2 inhibition for the synthesis of prostaglandins by the plant alkaloid aristolochic acid from a 1.7 A crystal structure

This is the first structural observation of a plant product showing high affinity for phospholipase A(2) and regulating the synthesis of arachidonic acid, an intermediate in the production of prostaglandins. The crystal structure of a complex formed between Vipera russelli phospholipase A(2) and a plant alkaloid aristolochic acid has been determined and refined to 1.7 A resolution. The structure contains two crystallographically independent molecules of phospholipase A(2) in the form of an asymmetric dimer with one molecule of aristolochic acid bound to one of them specifically. The most significant differences introduced by asymmetric molecular association in the structures of two molecules pertain to the conformations of their calcium binding loops, beta-wings, and the C-terminal regions. These differences are associated with a unique conformational behavior of Trp(31). Trp(31) is located at the entrance of the characteristic hydrophobic channel which works as a passage to the active site residues in the enzyme. In the case of molecule A, Trp(31) is found at the interface of two molecules and it forms a number of hydrophobic interactions with the residues of molecule B. Consequently, it is pulled outwardly, leaving the mouth of the hydrophobic channel wide open. On the other hand, Trp(31) in molecule B is exposed to the surface and moves inwardly due to the polar environment on the molecular surface, thus narrowing the opening of the hydrophobic channel. As a result, the aristolochic acid is bound to molecule A only while the binding site of molecule B is empty. It is noteworthy that the most critical interactions in the binding of aristolochic acid are provided by its OH group which forms two hydrogen bonds, one each with His(48) and Asp(49).     

Specificity and catalysis of uracil DNA glycosylase. A molecular dynamics study of reactant and product complexes with DNA.

The structure of uracil DNA glycosylase (UDG) in complex with a nonamer duplex DNA containing a uracil has been determined only in the product state. The reactant state was constructed by reattaching uracil to the deoxyribose, and both complexes were studied by molecular dynamics simulations. Significant changes in the positions of secondary structural elements in the enzyme are induced by the hydrolysis of the glycosidic bond. The simulations show that the specificity of the uracil pocket in the enzyme is largely retained in both complexes with the exception of Asn-204, which has been identified as a residue that contributes to discrimination between uracil and cytosine. The hydrogen bond between the amide group of Asn-204 and O(4) of uracil is disrupted by fluctuations of the side chain in the reactant state and is replaced by a hydrogen bond to water molecules trapped in the interior of the protein behind the uracil binding pocket. The role of two residues implicated by mutation experiments to be important in catalysis, His-268 and Asp-145, is clarified by the simulations. In the reactant state, His-268 is found 3.45 +/- 0.34 A from the uracil, allowing a water molecule to form a bridge to O(2). The environment in the enzyme raises the pK(a) value of His-268 to 7.1, establishing a protonated residue for assisting in the hydrolysis of the glycosidic bond. In agreement with the crystallographic structure, the DNA backbone retracts after the hydrolysis to allow His-268 to approach the O(2) of uracil with a concomitant release of the bridging water molecule and a reduction in the pK(a) to 5.5, which releases the proton to the product. The side chain of Asp-145 is fully solvated in the reactant state and H-bonded through a water molecule to the 3'-phosphate of uridine. Both the proximity of Asp-145 to the negatively charged phosphate and its pK(a) of 4.4 indicate that it cannot act as a general base catalyst. We propose a mechanism in which the bridging water between Asp-145 and the 3'-phosphate accepts a proton from another water to stabilize the bridge through a hydronium ion as well as to produce the hydroxide anion required for the hydrolytic step. The mechanism is consistent with known experimental data.     

Alkyl isocyanates as active site-specific reagents for serine proteases. Location of alkyl binding site in chymotrypsin by X-ray diffraction.

The structure of octylcarbamoyl-alpha-chymotrypsin to a resolution of 3.0 A is described. The n-octyl side chain of the active site directed irreversible inactivator octyl isocyanate is bound exclusively in the hydrophobic substrate binding pocket. The n-octyl isocyanate forms a planar urethane bond with the Ser-195 Ogamma and extends approximately 1 A deeper into the hydrophobic pocket than the indolyl group of indoleacryloyl-alpha-chymotrypsin (Henderson, R. (1970), J. Mol. Biol. 54, 341). All the structural changes are essentially identical with those observed in indoleacryloyl-alpha-chymotrypsin including the observation of a hydrogen bonded water molecule between the carbonyl oxygen of the octylcarbamoyl group and the imidazole group of His-57. The observed mode of n-octyl alkyl binding to chymotrypsin is consistent with the hypothesis proposed earlier (Brown, W. E. and Wold, F. (1973), Biochemistry 12, 828).     

Enhancement of catalytic efficiency of enzymes through exposure to anhydrous organic solvent at 70 degrees C. Three-dimensional structure of a treated serine proteinase at 2.2 A resolution

The enzyme behavior in anhydrous media has important applications in biotechnology. So far chemical modifications and protein engineering have been used to alter the catalytic power of the enzymes. For the first time, it is demonstrated that an exposure of enzyme to anhydrous organic solvents at optimized high temperature enhances its catalytic power through local changes at the binding region. Six enzymes: proteinase K, wheat germ acid phosphatase, alpha-amylase, beta-glucosidase, chymotrypsin and trypsin have been exposed to acetonitrile at 70 degrees C for three hours. The activities of these enzymes were found to be considerably enhanced. In order to understand the basis of this change in the activity of these enzymes, the structure of one of these treated enzymes, proteinase K has been analyzed in detail using X-ray diffraction method. The overall structure of the enzyme is similar to the native structure in aqueous environment. The hydrogen bonding system of the catalytic triad is intact after the treatment. However, the water structure in the substrate binding site undergoes some rearrangement as some of the water molecules are either displaced or completely absent. The most striking observation concerning the water structure pertains to the complete deletion of the water molecule which occupied the position at the so-called oxyanion hole in the active site of the native enzyme. Three acetonitrile molecules were found in the present structure. All the acetonitrile molecules are located in the recognition site. The sites occupied by acetonitrile molecules are independent of water molecules. The acetonitrile molecules are involved in extensive interactions with the protein atoms. All of them are interlinked through water molecules. The methyl group of one of the acetonitrile molecules (CCN1) interacts simultaneously with the hydrophobic side chains of Leu-96, Ile-107, and Leu-133. The development of such a hydrophobic environment at the recognition site introduces a striking conformation change in Ile-107 by rotating its side chain about C(alpha)--C(beta) bond by 180 degrees to bring about the delta-methyl group within the range of attractive van der Waals interactions with the methyl group of CCN1. A similar change has earlier been observed in proteinase K when it is complexed to a substrate analog lactoferrin fragment.     

Application of 3-Dimensional Homology Modeling of Cytochrome P450 2B1 for Interpretation of Site-Directed Mutagenesis Results

Abstract Three-dimensional structures of cytochrome P450 2B1 were modeled based on the crystallographic structure of P450(cam). The effect of the alignment, loop choice, and minimization with or without water was assessed. Although final models were similar in overall structure, the identity of active site residues depended upon the alignment. An example is Phe-206, which may or may not form part of the active site. The choice of the loop conformation had a lesser effect, while including water in the final minimization step was essential for preserving the shape and size of the active site. The best model (model 2) was in good agreement with the data from site-directed mutagenesis studies, and correctly predicted the effect of substitutions at 9 out of 10 amino acid positions. Thus, residues important for P450 2B1 activity, such as Ile- 114, Phe-206, Ile-290, Thr-302, Val-363, and Gly-478, constitute part of the active site and are able to interact with the substrate androstenedione through hydrophobic interactions. On the other hand, Ser-303, Ser-360 and Lys-473 are far from the active site and/or cannot interact with the substrate, in agreement with experimental data. The model indicates other residues likely to be important for enzyme function, such as Tyr- 111, Leu-209, Ile-477, and Ile- 480, which can be tested experimentally. The substrate may assume numerous binding orientations consistent with observed patterns of hydroxylation at C(5) and C(6). The replacement in the model of certain amino acid residues to mimic residue substitutions from site-directed mutagenesis studies and docking of the substrate into the modified active site allowed a plausible explanation for alterations in regio- and stereospecificities of some mutants of P450 2B1, such as Gly-478 → Ala or Val-363 Ala.     

Nuclear Overhauser enhancement spectroscopy cross-relaxation rates and ethanol distribution across membranes

Measurement of nuclear Overhauser enhancement spectroscopy cross-relaxation rates between ethanol and palmitoyloleoylphosphatidylcholine bilayers was combined with atomic-level molecular dynamics simulations. The molecular dynamics trajectories yielded autocorrelation functions of proton dipole-dipole interactions, and, consequently, relaxation times and cross-relaxation rates. These analyses allow the measured cross-relaxation rates to be interpreted in terms of relative interaction strengths with the various segments of the lipid molecule. We determined that cross-relaxation between ethanol and specific lipid resonances is primarily determined by the sites of interaction with some modulation due to lipid disorder and to local differences in intramolecular lipid dynamics. The rates scale linearly with the lifetime of temporary ethanol-lipid associations. Ethanol interacts with palmitoyloleoylphosphatidylcholine bilayers primarily via hydrophilic interactions, in particular the formation of hydrogen bonds to the lipid phosphate group. There is a weak contribution to binding from hydrophobic interaction with lipid chain segments near the glycerol. However, the strength of hydrophobic interactions is insufficient to compensate for the energetic loss of locating ethanol in an exclusively hydrophobic environment, resulting in a probability of locating ethanol in the bilayer center that is three orders of magnitude lower than locating ethanol at the lipid/water interface. The low cross-relaxation rates between terminal methyl protons of hydrocarbon chains and ethanol are as much the result of infrequent chain upturns as of brief excursions of ethanol into the region of lipid hydrocarbon chains near the glycerol. The combination of nuclear magnetic resonance measurements and molecular dynamics simulations offers a general pathway to study the interaction of small molecules with the lipid matrix at atomic resolution.