Molecular dynamics simulations of the unfolding of the starch binding domain from Aspergillus niger glucoamylase

The 600 ps molecular dynamics simulations to investigate the unfolding of the starch binding domain from Aspergillus niger glucoamylase were conducted in vacuum as well as in an external field with the dielectric constant of 80 with temperature jump technique. Electrostatic interactions play an important role in determining the stability of the beta-strands in this domain. The starch binding site 1 is less stable than site 2 since it is more exposed to the surface. The disulfide bond between C509 and C604 is unstable since these two residues are located near the flexible linker domain and in the mobile loop region between beta-strands 6 and 7, respectively. The melting temperature, at which the total residual beta-strand content is 50% that of the solution structure, is about 544K for the simulations with dielectric constant of 80, leading to the estimated unfolding timescale of 0.48 ms in vitro. In addition, the unfolding of the starch binding domain is proposed to initiate from the interior region by the lost of the integrity of the secondary structure.     

Molecular Dynamics Simulations of the Unfolding Mechanism of the Catalytic Domain from Aspergillus awamori var. X100 Glucoamylase

Abstract

In this study, 200 ps molecular dynamics simulations were conducted to investigate the unfolding mechanism of the catalytic domain of glucoamylase from Aspergillus awamori var. X100. The unfolding of this domain was suggested to follow a putative hierarchical manner, in which the heavily O-glycosylated belt region from residues T440 to A471 acted as the initiation site, followed by the a-helix secondary structure destruction, and then the collapse of the catalytic center pocket. The O-glycosylated belt region surrounded the surface of the catalytic domain in its native state at low temperature, whereas it was extended and is more suitable to be classified as part of the subsequent linker domain at high temperatures due to its high flexibility. The inner set helices of the (α/α)₆-barrel seemed to exhibit higher helical content than the outer set ones at all temperatures examined. The distances between the C_α of the three Cys residue pairs fluctuated rapidly at higher temperatures, indicating that these disulfide bonds have little effect on the structural stabilization. The melting temperature, at which the residual total helicity of the catalytic domain is 50%, is much lower than the critical temperature, at which the catalytic center pocket has lost its structural integrity.     

Biophysical insights into the binding characteristics of bovine serum albumin with dipyridamole and the influence of molecular interaction with β cyclodextrin

Abstract

The binding characteristic of anti-platelet drug dipyridamole has been investigated with a transport protein, serum albumin. A multi-spectroscopic approach has been employed, and the results were well supported by in silico molecular docking and simulation studies. The fluorescence quenching of serum albumin at three different temperatures revealed that the mechanism involved is static and the binding constant of the interaction was found to be of the order of 10⁴ M^?1. The reaction was found to be spontaneous and involved hydrophobic interactions. Synchronous, 3D fluorescence and CD spectroscopy indicated a change in conformation of bovine serum albumin (BSA) on interaction with DP. Using site-selective markers, the binding site of DP was found to be in subdomain IB. Molecular docking studies further corroborated these results. Molecular dynamic (MD) simulations showed lower RMSD values on interaction, suggesting the existence of a stable complex between DP and BSA. Furthermore, since β-Cyclodextrin (βCD) is used to improve the solubility of DP in ophthalmic solutions, therefore, the effect of (βCD) on the interaction of BSA and DP was also studied, and it was found that in the presence of βCD, the binding constant for BSA-DP interaction decreased. The present study is an attempt to characterize the transport of DP and to improve its bioavailability, consequently helping in dosage design to achieve optimum therapeutic levels.

Communicated by Ramaswamy H. Sarma     

Dissecting intrinsic and ligand-induced structural communication in the β3 headpiece of integrins

BackgroundGraph theory is widely used to dissect structural communication in biomolecular systems. Here, graph theory-based approaches were applied to the headpiece of integrins, adhesion cell-surface receptors that transmit signals across the plasma membranes.MethodsProtein Structure Network (PSN) analysis incorporating dynamic information either from molecular dynamics simulations or from Elastic Network Models was applied to the β3 domains from integrins α_Vβ₃ and α_IIbβ₃ in their apo and ligand-bound states.ResultsClosed and open states of the β headpiece are characterized by distinct allosteric communication pathways involving highly conserved amino acids at the two different α/β interfaces in the βI domain, the closed state being prompted to the closed-to-open transition. In the closed state, pure antagonism is associated with the establishment of communication pathways that start from the ligand, pass through the β1/α3,α4 interface, and end up in the hybrid domain by involving the Y110-Q82 link, which is weakened in the agonist-bound states.ConclusionsAllosteric communication in integrins relies on highly conserved and functionally relevant amino acid residues. The αβα-sandwich architecture of integrin βI domain dictates the structural communication between ligand binding site and hybrid domain. Differently from agonists, pure antagonists are directly involved in allosteric communication pathways and exert long-distance strengthening of the βI/hybrid interface. Release of the structure network in the ligand binding site is associated with the close-to-open transition accompanying the activation process.General significanceThe study strengthens the power of graph-based analyses to decipher allosteric communication intrinsic to protein folds and modified by functionally different ligands.     

Thermal unfolding molecular dynamics simulation of Escherichia coli dihydrofolate reductase: thermal stability of protein domains and unfolding pathway

Temperature induced unfolding of Escherichia coli dihydrofolate reductase was carried out by using molecular dynamic simulations. The simulations show that the unfolding generally involves an initial end-to-end collapse of the adenine binding domain into partially extended loops, followed by a gradual breakdown of the remaining beta sheet core structure. The core, which consists of beta strands 5-7, was observed to be the most resistant to thermal unfolding. This region, which is made up of part of the N terminus domain and part of the large domain of the E. coli dihydrofolate reductase, may constitute the nucleation site for protein folding and may be important for the eventual formation of both domains. The unfolding of different domains at different stages of the unfolding process suggests that protein domains vary in stability and that the rate at which they unfold can affect the overall outcome of the unfolding pathway. This observation is compared with the recently proposed hierarchical folding model. Finally, the results of the simulation were found to be consistent with a previous experimental study (Frieden, Proc Natl Acad Sci USA 1990;87:4413-4416) which showed that the folding process of E. coli dihydrofolate reductase involves sequential formation of the substrate binding sites.     

Computational prediction and experimental validation of the activator function of C2-β-D-glucopyranosyl-1,3,6,7-tetrahydroxyxanthone on pancreatic and hepatic hexokinase

Abstract

This study identifies and validates hexokinase type 4 (HK4), an isozyme of hexokinase in the liver and pancreas, as an important target of C2-β-D-glucopyranosyl-1,3,6,7-tetrahydroxyxanthone (βdGT), a xanthone glucoside suggested to have antidiabetic property. In the study, we applied the computational pipeline of molecular docking followed by the molecular dynamics simulations to shortlist potential βdGT protein targets. The analysis of protein dynamics and the binding free energy (ΔG) led us to the identification of HK4 as a key βdGT target, whereby the binding mode and domain dynamics suggested the activator function of βdGT. βdGT bound to the allosteric site of the isozyme ～13?Å away from the substrate (glucose)-binding site. The binding free energy of the ligand-protein complex was energetically feasible (ΔG, –41.61?kcal/mol) and the cleft angle deviation between the two (small and large) domains of HK4 revealed differential HK4 dynamics in response to βdGT binding. 3D structure analysis of the isozyme-ligand complex highlighted the role of Arg63, Glu67 and Lys458 in ligand stabilization and hydrophobic interactions mediated by Tyr214 and Met235. Experimental validation of the results of computational analysis confirmed the activator function of βdGT on HK4. The study has implication in diabetes as βdGT may be used to lower the blood glucose level by activating hepatic and pancreatic hexokinase without the risk of hypoglycemia.

Communicated by Ramaswamy H. Sarma     

Structural investigations into the binding mode of novel neolignans Cmp10 and Cmp19 microtubule stabilizers by in silico molecular docking,molecular dynamics,and binding free energy calculations

Microtubule stabilizers provide an important mode of treatment via mitotic cell arrest of cancer cells. Recently, we reported two novel neolignans derivatives Cmp10 and Cmp19 showing anticancer activity and working as microtubule stabilizers at micromolar concentrations. In this study, we have explored the binding site, mode of binding, and stabilization by two novel microtubule stabilizers Cmp10 and Cmp19 using in silico molecular docking, molecular dynamics (MD) simulation, and binding free energy calculations. Molecular docking studies were performed to explore the β-tubulin binding site of Cmp10 and Cmp19. Further, MD simulations were used to probe the β-tubulin stabilization mechanism by Cmp10 and Cmp19. Binding affinity was also compared for Cmp10 and Cmp19 using binding free energy calculations. Our docking results revealed that both the compounds bind at Ptxl binding site in β-tubulin. MD simulation studies showed that Cmp10 and Cmp19 binding stabilizes M-loop (Phe272-Val288) residues of β-tubulin and prevent its dynamics, leading to a better packing between α and β subunits from adjacent tubulin dimers. In addition, His229, Ser280 and Gln281, and Arg278, Thr276, and Ser232 were found to be the key amino acid residues forming H-bonds with Cmp10 and Cmp19, respectively. Consequently, binding free energy calculations indicated that Cmp10 (?113.655 kJ/mol) had better binding compared to Cmp19 (?95.216 kJ/mol). This study provides useful insight for better understanding of the binding mechanism of Cmp10 and Cmp19 and will be helpful in designing novel microtubule stabilizers.     

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 .     

Principal component analysis of chemical shift perturbation data of a multiple‐ligand‐binding system for elucidation of respective binding mechanism

Chemical shift perturbations (CSPs) in NMR spectra provide useful information about the interaction of a protein with its ligands. However, in a multiple‐ligand‐binding system, determining quantitative parameters such as a dissociation constant (K_d) is difficult. Here, we used a method we named CS‐PCA, a principal component analysis (PCA) of chemical shift (CS) data, to analyze the interaction between bovine β‐lactoglobulin (βLG) and 1‐anilinonaphthalene‐8‐sulfonate (ANS), which is a multiple‐ligand‐binding system. The CSP on the binding of ANS involved contributions from two distinct binding sites. PCA of the titration data successfully separated the CSP pattern into contributions from each site. Docking simulations based on the separated CSP patterns provided the structures of βLG–ANS complexes for each binding site. In addition, we determined the K_d values as 3.42 × 10⁻⁴M² and 2.51 × 10⁻³M for Sites 1 and 2, respectively. In contrast, it was difficult to obtain reliable K_d values for respective sites from the isothermal titration calorimetry experiments. Two ANS molecules were found to bind at Site 1 simultaneously, suggesting that the binding occurs cooperatively with a partial unfolding of the βLG structure. On the other hand, the binding of ANS to Site 2 was a simple attachment without a significant conformational change. From the present results, CS‐PCA was confirmed to provide not only the positions and the K_d values of binding sites but also information about the binding mechanism. Thus, it is anticipated to be a general method to investigate protein–ligand interactions. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Microscopic mechanisms that govern the titration response and pKa values of buried residues in staphylococcal nuclease mutants

To probe the microscopic mechanisms that govern the titration behavior of buried ionizable groups, microsecond explicit solvent molecular dynamics simulations are carried out for several mutants of Staphylococcal nuclease using both fixed charge and polarizable force fields. While the ionization of Asp 66, Glu 66, and Lys 125 lead to enhanced structural fluctuations and partial unfolding of adjacent α‐helical regions, the ionization of Lys 25 causes local unfolding of adjacent β sheets. Using the sampled conformational ensembles, good agreement with experimental pK_a values is obtained with Poisson–Boltzmann calculations using a protein dielectric constant of 2–4 for V66D/E; slightly larger dielectric constants are needed for Lys mutants especially L25K, suggesting that structural responses beyond microseconds are involved in ionization of Lys 25. Overall, the set of unbiased simulations provides insights into the spatial and temporal scales of protein and solvent motions that dictate the diverse titration behaviors of buried protein residues. Proteins 2017; 85:268–281. © 2016 Wiley Periodicals, Inc.     

Crystal structures of starch binding domain from Rhizopus oryzae glucoamylase in complex with isomaltooligosaccharide: Insights into polysaccharide binding mechanism of CBM21 family

Glucoamylases are responsible for hydrolysis of starch and polysaccharides to yield β‐d ‐glucose. Rhizopus oryzae glucoamylase (RoGA) is composed of an N‐terminal starch binding domain (SBD) and a C‐terminal catalytic domain connected by an O‐glycosylated linker. Two carbohydrate binding sites in RoSBD have been identified, site I is created by three highly conserved aromatic residues, Trp47, Tyr83, and Tyr94, and site II is built up by Tyr32 and Phe58. Here, the two crystal structures of RoSBD in complex with only α‐(1,6)‐linked isomaltotriose (RoSBD‐isoG3) and isomaltotetraose (RoSBD‐isoG4) have been determined at 1.2 and 1.3 Å, respectively. Interestingly, site II binding is observed in both complexes, while site I binding is only found in the RoSBD‐isoG4 complex. Hence, site II acts as the recognition binding site for carbohydrate and site I accommodates site II to bind isoG4. Site I participates in sugar binding only when the number of glucosyl units of oligosaccharides is more than three. Taken together, two carbohydrate binding sites in RoSBD cooperate to reinforce binding mode of glucoamylase with polysaccharides as well as the starch. Proteins 2014; 82:1079–1085. © 2013 Wiley Periodicals, Inc.     

Dynamics of starch granule biogenesis – the role of redox-regulated enzymes and low-affinity carbohydrate-binding modules

Abstract

The deposition and degradation of starch in plants is subject to extensive post-translational regulation. To permit degradation of B-type crystallites present in tuberous and leaf starch these starch types are phosphorylated by glucan, water dikinase (GWD). At the level of post-translational redox regulation, ADPglucose pyrophosphorylase, β-amylase (BAM1), limit dextrinase (LD), the starch phosphorylator GWD and the glucan phosphatase dual-specificity phosphatase 4 (DSP4), also named starch excess 4 (SEX4), are reductively activated in vitro. Redox screens now suggest the presence of a substantially more extensive and coordinated redox regulation involving a larger number of enzymes. Noticeably several of these enzymes contain a new type of low-affinity carbohydrate-binding module that we term a low-affinity starch-binding domain or LA-SBD. These are present in the CBM20, CBM45 and CBM53 families and can enable diurnal dynamics of starch–enzyme recognition. Such diurnal changes in starch binding have been indicated for the redox-regulated GWD and SEX4.     

The use of the AMBER force field in conformational analysis of carbohydrate molecules: Determination of the solution conformation of methyl α-lactoside by NMR spectroscopy,assisted by molecular mechanics and dynamics calculations

The solution conformation of methyl α-lactoside has been studied through nmr spectroscopy and molecular mechanics calculations using the assisted model building with energy refinement (AMBER) force field. The nmr data have included nuclear Overhauser effect (NOE) measurements hot It in the laboratory and rotating frames, longitudinal relaxation times, and homonuclear and heteronuclear coupling constants. The steady-state and transient NOEs have been interpreted in terms of an ensemble average distribution of conformers, making use of the complete relaxation matrix approach. The molecular mechanics calculations have been performed at two dielectric constants [ε = 1 * r and 80 Debyes (D)] in an exhaustive way, and have been complemented with specific calculations at intermediate ε values. Relaxed energy maps and adiabatic surfaces have been generated for the different dielectric constants. The probability distribution of conformers has been estimated from these steric energy maps. Molecular dynamics simulations in vacuo have also been performed. The experimental results indicate that the β(1 → 4)-glycosidic linkage shows some fluctuations among three low energy regions, although spends ca. 85% of its lime in the region close to the global minimum. It is shown that the over estimation of the electrostatic contributions in AMBER is responsible for the failure of this force field to explain the experimental results when Used at low dielectric constant (ε < 20 D). The matching between the expected and observed facts increases for ε > 40 D. Different conditions have been tested to perform temperature constant molecular dynamics simulations in vacuo, which have indicated that, when used without explicit solvent, this force field should only be employed in a qualitatively way when analyzing dynamical properties of oligosaccharides. © 1995 John Wiley & Sons, Inc.     

Novel stabilization pattern of thermolysin due to the binding substrate induced by electrostatic change in enzyme active site caused by the temperature elevation

During the synthesis of the dipeptide, N-benzyloxycarbonyl-l-phenylalanyl-l-phenylalanine methyl ester, from N-benzyloxycarbonyl-l-phenylalanine and l-phenylalanine methyl ester by thermolysin, the enzyme was stabilized by 20°C up to 110°C. The stabilization was caused by the interaction of the enzyme with Phe-OMe, a counterpart of the substrate, which was bound at the enzyme active site due to the drop in pH and dielectric constant following the temperature elevation of the medium. The binding of the enzyme to Phe-OMe suggested the induction of the transition state formation at around 80°C based on the UV spectra, resulting in the increase in the stability in the higher temperature region. The fluorescence second-order derivative spectra suggested that the binding Phe-OMe interacted with Trp 115 at the active site of the enzyme. The phenomenon was considered to be a novel stabilization pattern of the enzyme resulting from the conduction due to the chemical modification by the binding substrate.     

Interaction Between L-Aspartic Acid and L-Asparaginase from Escherichia Coli: Binding and Inhibition Studies

Abstract

Experiments using equilibrium dialysis and fluorescence quenching provided direct evidence that approximately four moles of L-aspartic acid were bound per mole of tetrameric L-asparaginase from Escherichia coli, with a dissociation constant on the order of 60-160 μM. In addition, a set of weaker binding sites with a dissociation constant in the millimolar range were detected. Kinetic studies also revealed that L-aspartic acid inhibited L-asparaginase competitively, with an inhibition constant of 80 μM at micromolar concentrations of L-asparagine; at millimolar concentrations of the amide, an increase in maximal velocity but a decrease in affinity for L-asparagine were observed. L-Aspartic acid at millimolar levels again displayed competitive inhibition. These and other observations suggest that L-aspartic acid binds not only to the active site but also a second site with lower intrinsic affinity for it. The observed “substrate activation” is most likely attributable to the binding of a second molecule of L-asparagine rather than negative cooperativity among the tight sites of the subunits of this tetrameric enzyme. Further support for L-aspartic acid binding to the active site comes from experiments in which the enzyme, when exposed to various group-specific reagents suffered parallel loss of catalytic activity and in its ability to bind L-aspartic acid. Different commercial preparations of Escherichia coli L-asparaginase were found to contain ～ 2-4 moles of L-aspartic acid; these were incompletely removed by dialysis, but could be removed by transamination or decarboxylation. Efficiency of dialysis increased with increasing pH. Taken together, this set of results is consistent with the existence of a covalent β-aspartyl enzyme intermediate.     

Force-Sensitive Autoinhibition of the von Willebrand Factor Is Mediated by Interdomain Interactions

Von Willebrand factor (VWF) plays a central role in hemostasis. Triggered by shear-stress, it adheres to platelets at sites of vascular injury. Inactivation of VWF has been associated to the shielding of its adhesion sites and proteolytic cleavage. However, the molecular nature of this shielding and its coupling to cleavage under shear-forces in flowing blood remain unknown. In this study, we describe, to our knowledge, a new force-sensory mechanism for VWF-platelet binding, which addresses these questions, based on a combination of molecular dynamics (MD) simulations, atomic force microscopy (AFM), and microfluidic experiments. Our MD simulations demonstrate that the VWF A2 domain targets a specific region at the VWF A1 domain, corresponding to the binding site of the platelet glycoprotein Ibα (GPIbα) receptor, thereby causing its blockage. This implies autoinhibition of the VWF for the binding of platelets mediated by the A1-A2 protein-protein interaction. During force-probe MD simulations, a stretching force dissociated the A1A2 complex, thereby unblocking the GPIbα binding site. Dissociation was found to be coupled to the unfolding of the A2 domain, with dissociation predominantly occurring before exposure of the cleavage site in A2, an observation that is supported by our AFM experiments. This suggests that the A2 domain prevents platelet binding in a force-dependent manner, ensuring that VWF initiates hemostasis before inactivation by proteolytic cleavage. Microfluidic experiments with an A2-deletion VWF mutant resulted in increased platelet binding, corroborating the key autoinhibitory role of the A2 domain within VWF multimers. Overall, autoinhibition of VWF mediated by force-dependent interdomain interactions offers the molecular basis for the shear-sensitive growth of VWF-platelet aggregates, and might be similarly involved in shear-induced VWF self-aggregation and other force-sensing functions in hemostasis.     

Unveiling the Unfolding Pathway of F5F8D Disorder-Associated D81H/V100D Mutant of MCFD2 via Multiple Molecular Dynamics Simulations

Abstract

Combined factor deficiency (F5F8D) is a rare autosomal recessive disorder caused by mutations in the LMAN1 or MCFD2 genes. It has been proposed that this pathogenic process occurs via a multi-step pathway involving metal loss, EF-hand-Ca²⁺ dissociation and assembly of misfolded MCFD2-LMAN1 complex. Here, we have investigated the solution conformations of the MCFD2^(D81H,V100D) protein mutant through extensive molecular dynamics (MD) simulations. The V100D, one of the many MCFD2 mutations known to be associated to F5F8D, is difficult to be reconciled with the pathway model because it is located far from the metal sites and the MCFD2/LMAN1 interface. Consequently, an inspection of all the steps involved in D81H/V100D MCFD2 misfolding is expected to provide hints in the understanding of the molecular basis of the disease. A comparison with parallel studies carried out for the Wild-Type (WT) MCFD2 pointed out that the mutation decreases the affinity of the protein for the Ca²⁺ ion. Multiple explicit solvents MD simulations (50 ns) performed on the two proteins revealed that in the WT protein, stable H-bond network and compact hydrophobic core region are created thus confirming a pivotal role of this region in driving the biophysical properties of the entire protein. In fact it is shown that the V100D mutation, although located far away the EF-hand domain, may induce subtle modification in the structural core of MCFD2 leading to the loosening of metal binding and to the formation of metastable intermediate states along the unfolding pathway. The native-like hydrophobic cluster formed near the V100 residue in the wild-type protein is disrupted by the negatively charged Asparagine residue. Furthermore, the presence of the D81H mutation in the EF-1 hand domain may also increase the protein unfolding rate and consequently prevent the formation of the MCFD2-LMAN1 complex. The detailed structural insights obtained from our large-scale simulations complement the clinical features and offer useful insights into the mechanism behind MCFD2 protein misfolding.     

Distinct Characteristics of Single Starch-Binding Domain SBD1 Derived from Tandem Domains SBD1-SBD2 of Halophilic <Emphasis Type="Italic">Kocuria varians</Emphasis> Alpha-Amylase

Kocuria varians alpha-amylase contains tandem starch-binding domains SBD1-SBD2 (SBD12) that possess typical halophilic characteristics. Recombinant tandem domains SBD12 and single domain SBD1, both with amino-terminal hexa-His tag, were expressed in and purified to homogeneity from Escherichia coli. The circular dichroism (CD) spectrum of His-SBD12 was characterized by a positive peak at 233 nm ascribed to the aromatic stacking. Although the signal occurred in the far UV region, it is an indication of tertiary structure folding. CD spectrum of single domain His-SBD1 exhibited the same peak position, signal intensity and spectral shape as those of His-SBD12, suggesting that the aromatic stacking must occur within the domain, and that two SBD domains in SBD12 and SBD1 has a similar folded structure. This structural observation was consistent with the biological activity that His-SBD1 showed binding activity against raw starch granules and amylose resin with 70–80% efficiency compared with binding of equimolar His-SBD12. Although the thermal unfolding rate of SBD12 and SBD1 were similar, the refolding rates of SBD12 and SBD1 from thermal melting were greatly different: His-SBD12 refolded slowly (T_1/2 = ~84 min), while refolding of single domain His-SBD1 was found to be 20-fold faster (T_1/2 = 4.2 min). The possible mechanism of this large difference in refolding rate was discussed. Maltose at 20 mM showed 5–6 °C increase in thermal melting of both His-SBD12 and His-SBD1, while its effects on the time course of unfolding and refolding were insignificant.     

Development of new α-amylases for raw starch hydrolysis

This paper describes the discovery of a new 4 domain α-amylase from Anoxybacillus contaminans which very efficiently hydrolyses raw starch granules. Compared to traditional starch liquefying α-amylases, this new 4 domain α-amylase contains a starch binding domain. The presence of this starch-binding domain enables the enzyme to efficiently hydrolyse starch at a temperature below the gelatinisation temperature. At a reaction temperature of 60°C and in combination with a glucoamylase from Aspergillusniger it was possible to liquefy 99% of the starch obtaining a DX value of 95%.

Furthermore, we describe how the current HFCS process can be turned into a low temperature simultaneous liquefaction and saccharification process by using this new 4 domain α-amylase in combination with a glucoamylase.