Coenzyme models: 46. “Remote control” of flavin reactivities by an intramolecular crown ring serving as a metal binding site: Relationship between spectral properties and dissociation of the 8-sulfonamide group

Crown ether flavin mimics (CrSO₂NHFl and NCrSO₂NHFl) which have a flavin moiety serving as a catalytic site and a crown ether moiety serving as a metal binding site at the two sides of a sulfonamide group were synthesized. We have found that the absorption spectra of these flavins are very sensitive to solvent effects; that is, they are yellow to orange in nonpolar solvents like “regular” flavins but imparted a red color to polar solvents characteristic of the intramolecular charge transfer like roseoflavin. This is due to the dissociation of the 8-sulfonamide group in polar solvents. The fluorescence spectra were also sensitive to solvent effects: the quantum yields of neutral NCrSO₂NHFl increased with decreasing solvent polarity. In acetonitrile, Ca²⁺ ion bound to the crown ether cavity in NCrSO₂NHFl facilitated deprotonation of the sulfonamide group to give a new absorption maximum at 452 nm. Correspondingly, the quantum yields for photooxidation of benzyl alcohol by NCrSO₂NHFl increased with increasing Ca²⁺ concentration. These findings indicate that Ca²⁺ ion can control the catalytic activity of NCrSO₂NHFl through the interaction with the sulfonamide group serving as a cap for Ca²⁺ bound to the crown cavity. The changes in the absorption spectra and the quantum yields were not observed for CrSO₂NHFl and a reference flavin, PhSO₂NHFl. Therefore, NCrSO₂NHFl acts as a new model system relevant to allosteric enzymes in which binding of an effector to a remote, allosteric site induces activity changes in the active sites.     

In silico approaches towards understanding CALB using molecular dynamics simulation and docking

Candida antarctica lipase B (CALB), a serine protease, is involved in the hydrolysis of substrates at the aqueous lipid interface. There is a significant role played by the helices in serine proteases including acting as a flap covering the active site region. The α5 and α10 helices in the path to the active site of CALB appear to play an important role in the region. This study investigates these helices by mutational studies, docking and molecular dynamics simulations. The mutations were selected based on their proximity to the active site and their presence at the α10-helix in the path of the active site. Molecular dynamics studies reveal the flexibility, stability and hydrogen bonding ability of the α5 helix. The radius of gyration (R _g) clearly showed the compactness of the structure. Docking studies show the changes occurring at the protein's binding site before and after 15 ns of simulation. Results from the study demonstrate the importance of the two helices α5 and α10 in the stability of CALB.     

Esterification degree of fructose laurate exerted by Candida antarctica lipase B in organic solvents

Sugar esters of fatty acids have many applications as biocompatible and biodegradable emulsifiers, which are determined by their degrees of esterification (DE). Direct esterification of fructose with lauric acid in organic media used commercial immobilized Candida antarctica lipase B (CALB) was investigated for DE. Significant difference of DE was observed between 2-methyl-2-butanol (2M2B) and methyl ethyl ketone (MEK), as di-ester/mono-ester molar ratio of 1.05:1 in 2M2B and 2.79:1 in MEK. Fourier transform infrared (FTIR) spectra showed that the secondary structure of the enzyme binding mono-ester presented distinct difference in 2M2B and MEK. Contents of β-turn and antiparallel β-sheet of CALB in 2M2B were 26.9% and 16.2%, respectively, but 19.1% and 13.2% in MEK. To understand the relationship between the conformational changes and differences of DE, mono-ester and fatty acid were directly employed for synthesis of di-ester. The maximum initial velocity of di-ester synthesis in MEK was 0.59 mmol g (enzyme)⁻¹ h⁻¹, which was 2.19-fold as greater as that in 2M2B, indicating that CALB conformation in MEK was preferred for the synthesis of di-ester. These results demonstrated that the conformation of CALB binding mono-ester affected by organic solvents essentially determined DE.     

Functional motions of Candida antarctica lipase B: a survey through open-close conformations

Candida antarctica lipase B (CALB) belongs to psychrophilic lipases which hydrolyze carboxyl ester bonds at low temperatures. There have been some features reported about cold-activity of the enzyme through experimental methods, whereas there is no detailed information on its mechanism of action at molecular level. Herein, a comparative molecular dynamics simulation and essential dynamics analysis have been carried out at three temperatures (5, 35 and 50 °C) to trace the dominant factors in the psychrophilic properties of CALB under cold condition. The results clearly describe the effect of temperature on CALB with meaningful differences in the flexibility of the lid region (α5 helix), covering residues 141-147. Open- closed conformations have been obtained from different sets of long-term simulations (60 ns) at 5 °C gave two reproducible distinct forms of CALB. The starting open conformation became closed immediately at 35 and 50 °C during 60 ns of simulation, while a sequential open-closed form was observed at 5 °C. These structural alterations were resulted from α5 helical movements, where the closed conformation of active site cleft was formed by displacement of both helix and its side chains. Analysis of normal mode showed concerted motions that are involved in the movement of both α5 and α10 helices. It is suggested that the functional motions needed for lypolytic activity of CALB is constructed from short-range movement of α5, accompanied by long-range movement of the domains connected to the lid region.     

Hydration of enzyme in nonaqueous media is consistent with solvent dependence of its activity

Water plays an important role in enzyme structure and function in aqueous media. That role becomes even more important when one focuses on enzymes in low water media. Here we present results from molecular dynamics simulations of surfactant-solubilized subtilisin BPN' in three organic solvents (octane, tetrahydrofuran, and acetonitrile) and in pure water. Trajectories from simulations are analyzed with a focus on enzyme structure, flexibility, and the details of enzyme hydration. The overall enzyme and backbone structures, as well as individual residue flexibility, do not show significant differences between water and the three organic solvents over a timescale of several nanoseconds currently accessible to large-scale molecular dynamics simulations. The key factor that distinguishes molecular-level details in different media is the partitioning of hydration water between the enzyme and the bulk solvent. The enzyme surface and the active site region are well hydrated in aqueous medium, whereas with increasing polarity of the organic solvent (octane --> tetrahydrofuran --> acetonitrile) the hydration water is stripped from the enzyme surface. Water stripping is accompanied by the penetration of tetrahydrofuran and acetonitrile molecules into crevices on the enzyme surface and especially into the active site. More polar organic solvents (tetrahydrofuran and acetonitrile) replace mobile and weakly bound water molecules in the active site and leave primarily the tightly bound water in that region. In contrast, the lack of water stripping in octane allows efficient hydration of the active site uniformly by mobile and weakly bound water and some structural water similar to that in aqueous solution. These differences in the active site hydration are consistent with the inverse dependence of enzymatic activity on organic solvent polarity and indicate that the behavior of hydration water on the enzyme surface and in the active site is an important determinant of biological function especially in low water media.     

Insights into diterpene cyclization from structure of bifunctional abietadiene synthase from Abies grandis

Abietadiene synthase from Abies grandis (AgAS) is a model system for diterpene synthase activity, catalyzing class I (ionization-initiated) and class II (protonation-initiated) cyclization reactions. Reported here is the crystal structure of AgAS at 2.3 Å resolution and molecular dynamics simulations of that structure with and without active site ligands. AgAS has three domains (α, β, and γ). The class I active site is within the C-terminal α domain, and the class II active site is between the N-terminal γ and β domains. The domain organization resembles that of monofunctional diterpene synthases and is consistent with proposed evolutionary origins of terpene synthases. Molecular dynamics simulations were carried out to determine the effect of substrate binding on enzymatic structure. Although such studies of the class I active site do lead to an enclosed substrate-Mg²⁺ complex similar to that observed in crystal structures of related plant enzymes, it does not enforce a single substrate conformation consistent with the known product stereochemistry. Simulations of the class II active site were more informative, with observation of a well ordered external loop migration. This “loop-in” conformation not only limits solvent access but also greatly increases the number of conformational states accessible to the substrate while destabilizing the nonproductive substrate conformation present in the “loop-out” conformation. Moreover, these conformational changes at the class II active site drive the substrate toward the proposed transition state. Docked substrate complexes were further assessed with regard to the effects of site-directed mutations on class I and II activities.     

Effects of organic solvents and substrate binding on trypsin in acetonitrile and hexane media

In this work, we used molecular dynamic (MD) simulation to study trypsin with and without a six-amino-acid peptide bound in three different solvents (water, acetonitrile and hexane) in order to provide molecular information for well understanding the structure and function of enzymes in non-aqueous media. The results show that the enzyme is more compact and less native-like in hexane than in the other two polar solvents. The substrate could stabilize the native protein structure in the two polar media, but not in the non-polar hexane. There are no significant differences in the conformation of the S1 pocket upon the substrate binding in water and acetonitrile media while a reverse behavior is observed in hexane media, implying a possible induced fit binding mechanism in the non-polar media. The substrate binding enhances the stability of catalytic H-bond network since it could expel the solvent molecules from the active site. The enzyme and the substrate appear to be more appropriate to the reactive conformation in the organic solvents compared with aqueous solution. There is much greater substrate binding strength in hexane media than the water and acetonitrile ones since the polar solvent significantly weakens electrostatic interactions, which are observed to be the main driving force to the binding. In addition, some residues of the S1 pocket could remain favorable contribution to the binding despite the solvent change, but with differences in the contribution extent, the number and the type of residues between the three media.

Synthesis and conformational analysis of cyclo-TRI[l-valyl-d-hexahydromandelyl]

The synthesis of the cyclo-hexadepsipeptide [l-valyl-d-hexahydromandelyl]₃ is described. Examination of this macrocyclic compound by 220-MHz nuclear magnetic resonance spectroscopy shows that symmetrical conformations are stabilized in strongly polar solvents (trifluoroacetic acid, acetonitrile), whereas asymmetric conformations are preferred in nonpolar or slightly polar media such as carbon tetrachloride, chloroform, cyclohexane, and benzene.From analysis of the temperature dependence of the chemical shift and of the coupling constants, together with conformational energy calculations, a model is proposed for the preferred conformation of this molecule in nonpolar solvents.     

CLEAs of Candida antarctica lipase B (CALB) with a bovine serum albumin (BSA) cofeeder core: Study of their catalytic activity

Highly active CALB cross-linked enzyme aggregates (CLEAs) were synthesized using a layered methodology based on the synthesis of a cross-linked protein cofeeder core over which an external layer of lipase was later cross-linked. The layered CALB CLEAs were characterized in terms of their catalytic activity in three different test reactions: esterification of oleic acid and ethanol in absence of solvents, esterification of oleic acid and heptanol in organic medium, and hydrolysis of triolein in emulsioned medium. The impact of the cross-linker/protein mass ratio on CLEAs activity, and its evolution with storage time were evaluated in the solventless synthesis of ethyloleate. The amount of cross-linker used showed to be a key parameter for the evolution of the catalytic activity of CLEAs during storage. Under the best conditions found, hyperactivated CALB CLEAs with up to 188% of recovered activity in ethyl oleate synthesis were obtained. In terms of hydrolytic activity mature layered CALB CLEAs showed a retained activity of 68%. The assay of dried mature layered CALB CLEAs in heptyl oleate synthesis showed catalytic activities much higher than the one exhibited by free CALB, reaching 1 h-fatty acid conversions of 14% and 2%, respectively. The high catalytic activity shown by layered CALB CLEAs, suggests that they are an interesting alternative specially for the catalysis of fatty acid esterifications in both organic and solventless medium.     

The Molecular Determinants of the Increased Reduction Potential of the Rubredoxin Domain of Rubrerythrin Relative to Rubredoxin

Based on the crystal structures, three possible sequence determinants have been suggested as the cause of a 285 mV increase in reduction potential of the rubredoxin domain of rubrerythrin over rubredoxin by modulating the polar environment around the redox site. Here, electrostatic calculations of crystal structures of rubredoxin and rubrerythrin and molecular dynamics simulations of rubredoxin wild-type and mutants are used to elucidate the contributions to the increased reduction potential. Asn¹⁶⁰ and His¹⁷⁹ in rubrerythrin versus valines in rubredoxins are predicted to be the major contributors, as the polar side chains contribute significantly to the electrostatic potential in the redox site region. The mutant simulations show both side chains rotating on a nanosecond timescale between two conformations with different electrostatic contributions. Reduction also causes a change in the reduction energy that is consistent with a linear response due to the interesting mechanism of shifting the relative populations of the two conformations. In addition to this, a simulation of a triple mutant indicates the side-chain rotations are approximately anticorrelated so whereas one is in the high potential conformation, the other is in the low potential conformation. However, Ala¹⁷⁶ in rubrerythrin versus a leucine in rubredoxin is not predicted to be a large contributor, because the solvent accessibility increases only slightly in mutant simulations and because it is buried in the interface of the rubrerythrin homodimer.     

A Substrate-induced Biotin Binding Pocket in the Carboxyltransferase Domain of Pyruvate Carboxylase

Biotin-dependent enzymes catalyze carboxyl transfer reactions by efficiently coordinating multiple reactions between spatially distinct active sites. Pyruvate carboxylase (PC), a multifunctional biotin-dependent enzyme, catalyzes the bicarbonate- and MgATP-dependent carboxylation of pyruvate to oxaloacetate, an important anaplerotic reaction in mammalian tissues. To complete the overall reaction, the tethered biotin prosthetic group must first gain access to the biotin carboxylase domain and become carboxylated and then translocate to the carboxyltransferase domain, where the carboxyl group is transferred from biotin to pyruvate. Here, we report structural and kinetic evidence for the formation of a substrate-induced biotin binding pocket in the carboxyltransferase domain of PC from Rhizobium etli. Structures of the carboxyltransferase domain reveal that R. etli PC occupies a symmetrical conformation in the absence of the biotin carboxylase domain and that the carboxyltransferase domain active site is conformationally rearranged upon pyruvate binding. This conformational change is stabilized by the interaction of the conserved residues Asp⁵⁹⁰ and Tyr⁶²⁸ and results in the formation of the biotin binding pocket. Site-directed mutations at these residues reduce the rate of biotin-dependent reactions but have no effect on the rate of biotin-independent oxaloacetate decarboxylation. Given the conservation with carboxyltransferase domains in oxaloacetate decarboxylase and transcarboxylase, the structure-based mechanism described for PC may be applicable to the larger family of biotin-dependent enzymes.     

Minimal Zn Binding Site of Amyloid-β

Zinc-induced aggregation of amyloid-β peptide (Aβ) is a hallmark molecular feature of Alzheimer's disease. Here we provide direct thermodynamic evidence that elucidates the role of the Aβ region 6-14 as the minimal Zn²⁺ binding site wherein the ion is coordinated by His⁶, Glu¹¹, His¹³, and His¹⁴. With the help of isothermal titration calorimetry and quantum mechanics/molecular mechanics simulations, the region 11-14 was determined as the primary zinc recognition site and considered an important drug-target candidate to prevent Zn²⁺-induced aggregation of Aβ.     

1H nuclear magnetic resonance studies of N-acetyl-L-proline N-methylamide. Molecular conformations,hydrogen bondings,and thermodynamic quantities in various solvents

Concentration and temperature dependences of the ¹H nmr spectra of N-acetyl-L -proline N-methylamide were observed in various solvents [CCl₄, CDCl₃, (CD₃)₂CO, (CD₃)₂SO, H₂O, and D₂O]. The fraction of the cis isomer (with respect to the bond between the acetyl carbonyl carbon and prolyl nitrogen atoms) depends greatly on the solvent used; the fraction of the cis isomer is higher in polar solvents than in nonpolar solvents. It depends also on concentration and temperature in nonpolar solvents but not in polar solvents. In nonpolar solvents the trans isomer mostly exists in the γ-turn structure with an intramolecular hydrogen bond and the cis isomer tends to form molecular aggregates by intermolecular hydrogen bonds. In polar solvents both the cis and trans isomers exist in monomeric forms which interact with solvent molecules. The pH dependences of the N-methyl proton resonances indicate that the γ-turn structure of the trans isomer is present also in aqueous solution, though its population is difficult to determine. Apparent enthalpy and entropy changes for the conversion of the trans isomer to cis isomer are evaluated for various solvents. The results are discussed in terms of the intra- and intermolecular hydrogen bondings.     

Lipase in aqueous‐polar organic solvents: Activity,structure, and stability

Studying alterations in biophysical and biochemical behavior of enzymes in the presence of organic solvents and the underlying cause(s) has important implications in biotechnology. We investigated the effects of aqueous solutions of polar organic solvents on ester hydrolytic activity, structure and stability of a lipase. Relative activity of the lipase monotonically decreased with increasing concentration of acetone, acetonitrile, and DMF but increased at lower concentrations (upto ～20% v/v) of dimethylsulfoxide, isopropanol, and methanol. None of the organic solvents caused any appreciable structural change as evident from circular dichorism and NMR studies, thus do not support any significant role of enzyme denaturation in activity change. Change in 2D [15N, 1H]‐HSQC chemical shifts suggested that all the organic solvents preferentially localize to a hydrophobic patch in the active‐site vicinity and no chemical shift perturbation was observed for residues present in protein's core. This suggests that activity alteration might be directly linked to change in active site environment only. All organic solvents decreased the apparent binding of substrate to the enzyme (increased K_m); however significantly enhanced the k_cat. Melting temperature (T_m) of lipase, measured by circular dichroism and differential scanning calorimetry, altered in all solvents, albeit to a variable extent. Interestingly, although the effect of all organic solvents on various properties on lipase is qualitatively similar, our study suggest that magnitudes of effects do not appear to follow bulk solvent properties like polarity and the solvent effects are apparently dictated by specific and local interactions of solvent molecule(s) with the protein.     

Highly polar environments catalyze the unfolding of PrPC helix 1

The first α-helix (H1) likely plays an important role in the conversion of the cellular prion protein (PrP^C) into its pathogenic isoform (PrP^Sc). In this conversion, H1 may either have to unfold or may represent a site of intermolecular contact. A recent molecular dynamics simulation suggested that H1 can unfold if it is detached from the protein core (Hirschberger et al. in Biophys J 90:3908, 2006). It has been hypothesized that the high dielectric constant ε _S of the bulk water environment facilitates the unfolding of H1. To check this hypothesis, we performed a number of replica exchange molecular dynamics simulations of an H1 peptide in solvents of different ε _S . We found that the equilibrium helix fraction in water is less than 40%, in agreement with previous experimental findings, and that the helix unfolds much faster in water than in less polar solvents. The kinetically stabilizing effect of the organic solvents is largely unspecific and correlates well with their dielectric constant ε _S .     

Crystal structure of Escherichia coli diaminopropionate ammonia-lyase reveals mechanism of enzyme activation and catalysis

Pyridoxal 5′-phosphate (PLP)-dependent enzymes utilize the unique chemistry of a pyridine ring to carry out diverse reactions involving amino acids. Diaminopropionate (DAP) ammonia-lyase (DAPAL) is a prokaryotic PLP-dependent enzyme that catalyzes the degradation of d- and l-forms of DAP to pyruvate and ammonia. Here, we report the first crystal structure of DAPAL from Escherichia coli (EcDAPAL) in tetragonal and monoclinic forms at 2.0 and 2.2 Å resolutions, respectively. Structures of EcDAPAL soaked with substrates were also determined. EcDAPAL has a typical fold type II PLP-dependent enzyme topology consisting of a large and a small domain with the active site at the interface of the two domains. The enzyme is a homodimer with a unique biological interface not observed earlier. Structure of the enzyme in the tetragonal form had PLP bound at the active site, whereas the monoclinic structure was in the apo-form. Analysis of the apo and holo structures revealed that the region around the active site undergoes transition from a disordered to ordered state and assumes a conformation suitable for catalysis only upon PLP binding. A novel disulfide was found to occur near a channel that is likely to regulate entry of ligands to the active site. EcDAPAL soaked with dl-DAP revealed density at the active site appropriate for the reaction intermediate aminoacrylate, which is consistent with the observation that EcDAPAL has low activity under crystallization conditions. Based on the analysis of the structure and results of site-directed mutagenesis, a two-base mechanism of catalysis involving Asp¹²⁰ and Lys⁷⁷ is suggested.     

Functional interactions in bacteriorhodopsin: a theoretical analysis of retinal hydrogen bonding with water.

The light-driven proton pump, bacteriorhodopsin (bR) contains a retinal molecule with a Schiff base moiety that can participate in hydrogen-bonding interactions in an internal, water-containing channel. Here we combine quantum chemistry and molecular mechanics techniques to determine the geometries and energetics of retinal Schiff base-water interactions. Ab initio molecular orbital calculations are used to determine potential surfaces for water-Schiff base hydrogen-bonding and to characterize the energetics of rotation of the C-C single bond distal and adjacent to the Schiff base NH group. The ab initio results are combined with semiempirical quantum chemistry calculations to produce a data set used for the parameterization of a molecular mechanics energy function for retinal. Using the molecular mechanics force field the hydrated retinal and associated bR protein environment are energy-minimized and the resulting geometries examined. Two distinct sites are found in which water molecules can have hydrogen-bonding interactions with the Schiff base: one near the NH group of the Schiff base in a polar region directed towards the extracellular side, and the other near a retinal CH group in a relatively nonpolar region, directed towards the cytoplasmic side.     

Stability and structure of Penicillium chrysogenum lipase in the presence of organic solvents

Abstract

The present work describes the enzymatic properties of Penicillium chrysogenum lipase and its behavior in the presence of organic solvents. The temperature and pH optima of the purified lipase was found to be 55?°C and pH 8.0 respectively. The lipase displayed remarkable stability in both polar and non-polar solvents upto 50% (v/v) concentrations for 72?h. A structural perspective of the purified lipase in different organic solvents was gained by using circular dichroism and intrinsic fluorescence spectroscopy. The native lipase consisted of a predominant α-helix structure which was maintained in both polar and non-polar solvents with the exception of ethyl butyrate where the activity was decreased and the structure was disrupted. The quenching of fluorescence intensity in the presence of organic solvents indicated the transformation of the lipase microenviroment P. chrysogenum lipase offers an interesting system for understanding the solvent stability mechanisms which could be used for rationale designing of engineered lipase biocatalysts for application in organic synthesis in non-aqueous media.     

Molecular dynamics simulations of an enzyme surrounded by vacuum, water, or a hydrophobic solvent.

We report on molecular dynamics simulations of a medium-sized protein, a lipase from Rhizomucor miehei, in vacuum, in water, and in a nonpolar solvent, methyl hexanoate. Depending on force field and solvent, the molecular dynamics structures obtained as averages over 150 ps had root-mean-square deviations in the range of 1.9 to 3.6 A from the crystal structure. The largest differences between the structures were in hydrogen bonding and exposed surface areas of the protein. The surface area increased in both solvents and became smaller in vacuum. The change of surface exposure varied greatly between different residues and occurred in accordance with the hydrophobicity of the residue and the nature of the solvent. The fluctuations of the atoms were largest in the external loops and agreed well with crystallographic temperature factors. Root-mean-square fluctuations were significantly smaller in the nonpolar solvents than they were in water, which is in accordance with the notion that proteins become more rigid in nonpolar solvents. In methyl hexanoate a partial opening of the lid covering the active site occurred, letting a methyl hexanoate molecule approach the active site.     

X-ray Crystallographic and MD Simulation Studies on the Mechanism of Interfacial Activation of a Family I.3 Lipase with Two Lids

The interfacial activation mechanism of family I.3 lipase from Pseudomonas sp. MIS38 (PML), which has two α-helical lids (lid1 and lid2), was investigated using a combination of X-ray crystallography and molecular dynamics (MD) simulation. The crystal structure of PML in an open conformation was determined at 2.1 Å resolution in the presence of Ca²⁺ and Triton X-100. Comparison of this structure with that in the closed conformation indicates that both lids greatly change their positions and lid1 is anchored by the calcium ion (Ca1) in the open conformation. This structure was not seriously changed even when the protein was dialyzed extensively against the Ca²⁺-free buffer containing Triton X-100 before crystallization, indicating that the open conformation is fairly stable unless a micellar substance is removed. The crystal structure of the PML derivative, in which the active site serine residue (Ser207) is diethylphosphorylated by soaking the crystal of PML in the open conformation in a solution containing diethyl p-nitrophenyl phosphate, was also determined. This structure greatly resembles that in the open conformation, indicating that PML structure in the open conformation represents that in the active form. MD simulation of PML in the open conformation in the absence of micelles showed that lid2 closes first, while lid1 maintains its open conformation. Likewise, MD simulation of PML in the closed conformation in the absence of Ca²⁺ and in the presence of octane or trilaurin micelles showed that lid1 opens, while lid2 remains closed. These results suggest that Ca1 functions as a hook for stabilization of a fully opened conformation of lid1 and for initiation of subsequent opening of lid2.