Effects of a nonionic surfactant on the behavior of Bacillus amyloliquefaciens α-amylase in the hydrolysis of malto-oligosaccharide

The effects of a nonionic surfactant, polyoxyethylene mono-N-dodecyl-ether (Brij 35; number of units of ethylene oxide moieties, 23), on the behavior of Bacillus amyloliquefaciens α-amylase in the hydrolysis of maltoheptaose were studied at 25°C and pH 7.0. The addition of the surfactant to the enzyme-substrate system increased the amount of reducing sugars during the hydrolysis. A high-performance liquid chromatography technique was used for the determination of the hydrolysis products. During the hydrolysis, the concentration of maltose was increased by the addition of the surfactant, whereas the concentrations of other products were unchanged. Therefore, the increase in the amount of reducing sugars was thought to be due to the increase of the maltose fragment. Because the transglycosylation of hydrolyzed maltose takes place during the hydrolysis, the suppression of transglycosylation accelerates the production of total reducing sugars. The interaction between the maltose fragment and micelle surface might induce the release of this fragment from the enzyme active site and then suppress the transglycosylation. The bindings of the enzyme, the substrate, and/or the products to the micelle surface may affect the hydrolytic behavior including both the hydrolysis and transglycosylation catalytic conversions.     

Molecular Dynamics Simulations of Bromodomains Reveal Binding-Site Flexibility and Multiple Binding Modes of the Natural Ligand Acetyl-Lysine

Experimental protein structures provide spatial information at the atomic level. A further dimension, time, is supplemented by molecular dynamics. Since the pioneering work on the 58-residue inhibitor of bovine pancreatic trypsin in the group of Martin Karplus in the seventies, molecular dynamics simulations have shown that the intrinsic flexibility of proteins is essential for their function. Here, we review simulation studies of bromodomains. These protein modules are involved in the recognition of acetylated lysine side chains, a post-translational modification frequently observed in histone tails. The molecular dynamics simulations have unmasked: (i) the large plasticity of the loops lining the acetyl-lysine binding site (coupled to its self-occlusion), and (ii) multiple binding modes of acetyl-lysine. These simulation results suggest that recognition of histone tails by bromodomains is modulated by their intrinsic flexibility, and further corroborate the utility of molecular dynamics in understanding (macro)molecular recognition.     

Loop Protein Engineering for Improved Transglycosylation Activity of a β‐N‐Acetylhexosaminidase

Certain enzymes of the glycoside hydrolase family 20 (GH20) exert transglycosylation activity and catalyze the transfer of β‐N‐acetylglucosamine (GlcNAc) from a chitobiose donor to lactose to produce lacto‐N‐triose II (LNT2), a key human milk oligosaccharide backbone moiety. The present work is aimed at increasing the transglycosylation activity of two selected hexosaminidases, HEX1 and HEX2, to synthesize LNT2 from lactose and chitobiose. Peptide pattern recognition analysis was used to categorize all GH20 proteins in subgroups. On this basis, we identified a series of proteins related to HEX1 and HEX2. By sequence alignment, four additional loop sequences were identified that were not present in HEX1 and HEX2. Insertion of these loop sequences into the wild‐type sequences induced increased transglycosylation activity for three out of eight mutants. The best mutant, HEX1_GTEPG, had a transglycosylation yield of LNT2 on the donor that was nine times higher than that of the wild‐type enzyme. Homology modeling of the enzymes revealed that the loop insertion produced a more shielded substrate‐binding pocket. This shielding is suggested to explain the reduced hydrolytic activity, which in turn resulted in the increased transglycosylation activity of HEX1_GTEPG.     

Engineering a substrate-specific cold-adapted subtilisin

One region predicted to be highly flexible for a psychrophilic enzyme, TA39 subtilisin (S39), was transferred in silico to the mesophilic subtilisin, savinase (EC 3.4.21.62), from Bacillus lentus (clausii). The engineered hybrid and savinase were initially investigated by molecular dynamic simulations at 300 K to show binding region and global flexibility. The predicted S39 region consists of 12 residues, which due to homology between the subtilisins, results in a total change of eight residues. By site-directed modifications, the region was transferred to the binding region of savinase, thus a savinase-S39 hybrid, named H5, was constructed. The designed hybrid showed the same temperature optimum and pH profile as savinase, but H5 had higher specific activity on the synthetic substrate N-succinyl-L-Ala-L-Ala-L-Pro-L-Phe-p-nitroanilide (AAPF) at all temperatures measured and, at the same time, H5 showed a decrease in thermostability. The H5 hybrid showed broader substrate specificity, measured at room temperature, due to an increase in catalytic efficiency on AAPF, AAPA and FAAF compared with savinase (N-succinyl-XXXX-pNA; XXXX = AAPF, AAPA and FAAF). The H5 hybrid showed increased activity at low temperature, increased binding region and global flexibility, as investigated by molecular dynamic simulations, and global destabilization from differential scanning calorimetry measurements. These psychrophilic characteristics indicated an increase in binding site flexibility, probably due to the modifications P129S, S130G, P131E, and thus we show that it is possible to increase low temperature activity and global flexibility by engineered flexibility in the binding region.     

Molecular Mechanism of the Hydration of Candida antarctica Lipase B in the Gas Phase: Water Adsorption Isotherms and Molecular Dynamics Simulations

Hydration is a major determinant of activity and selectivity of enzymes in organic solvents or in gas phase. The molecular mechanism of the hydration of Candida antarctica lipase B (CALB) and its dependence on the thermodynamic activity of water (a_w) was studied by molecular dynamics simulations and compared to experimentally determined water sorption isotherms. Hydration occurred in two phases. At low water activity, single water molecules bound to specific water binding sites at the protein surface. As the water activity increased, water networks gradually developed. The number of protein‐bound water molecules increased linearly with a_w, until at a_w=0.5 a spanning water network was formed consisting of 311 water molecules, which covered the hydrophilic surface of CALB, with the exception of the hydrophobic substrate‐binding site. At higher water activity, the thickness of the hydration shell increased up to 10 Å close to a_w=1. Above a limit of 1600 protein‐bound water molecules the hydration shell becomes unstable and the formation of pure water droplets occurs in these oversaturated simulation conditions. While the structure and the overall flexibility of CALB was independent of the hydration state, the flexibility of individual loops was sensitive to hydration: some loops, such as those part of the substrate‐binding site, became more flexible, while other parts of the protein became more rigid upon hydration. However, the molecular mechanism of how flexibility is related to activity and selectivity is still elusive.     

Molecular Dynamics Simulations of Antibiotic Ceftaroline at the Allosteric Site of Penicillin-Binding Protein 2a (PBP2a)

Methicillin-resistant Staphylococcus aureus (MRSA) tolerates β-lactam antibiotics by carrying out cell wall synthesis with the transpeptidase Penicillin-binding protein 2a (PBP2a), which cannot be inhibited by β-lactams. It has been proposed that PBP2a's active site is protected by two loops to reduce the probability of it binding with β-lactams. Previous crystallographic studies suggested that this protected active site opens for reaction once a native substrate binds at an allosteric domain of PBP2a. This opening was proposed for the new β-lactam ceftaroline's mechanism in successfully treating MRSA infections, i. e. by it binding to the allosteric site, thereby opening the active site to inhibition. In this work, we investigate the binding of ceftaroline at this proposed allosteric site using molecular dynamics simulations. Unstable binding was observed using the major force fields CHARMM36 and Amber ff14SB, and free energy calculations were unable to confirm a strong allosteric effect. Our study suggests that the allosteric effect induced by ceftaroline is weak at best.     

Structure and Dynamics of Meprin β in Complex with a Hydroxamate-Based Inhibitor

The astacin protease Meprin β represents an emerging target for drug development due to its potential involvement in disorders such as acute and chronic kidney injury and fibrosis. Here, we elaborate on the structural basis of inhibition by a specific Meprin β inhibitor. Our analysis of the crystal structure suggests different binding modes of the inhibitor to the active site. This flexibility is caused, at least in part, by movement of the C-terminal region of the protease domain (CTD). The CTD movement narrows the active site cleft upon inhibitor binding. Compared with other astacin proteases, among these the highly homologous isoenzyme Meprin α, differences in the subsites account for the unique selectivity of the inhibitor. Although the inhibitor shows substantial flexibility in orientation within the active site, the structural data as well as binding analyses, including molecular dynamics simulations, support a contribution of electrostatic interactions, presumably by arginine residues, to binding and specificity. Collectively, the results presented here and previously support an induced fit and substantial movement of the CTD upon ligand binding and, possibly, during catalysis. To the best of our knowledge, we here present the first structure of a Meprin β holoenzyme containing a zinc ion and a specific inhibitor bound to the active site. The structural data will guide rational drug design and the discovery of highly potent Meprin inhibitors.     

氧气和一氧化碳在人血红蛋白迁移过程研究

为了揭示CO和O₂竞争性结合人血红蛋白血红素位点的机制及其与人血红蛋白结构转换之间的关系，本文采用全原子分子动力学模拟（MD）结合马尔科夫状态模型（MSMs）研究氧气（O₂）和一氧化碳（CO）分子从水溶液迁移进入人血红蛋白四聚体α链和β链的全过程。分子动力学模拟揭示了O₂和CO结合α链和β链的稳态结合位点和瞬态结合位点、迁移通道以及α链的结构变化。结果显示，分子模拟不仅仅能够再现全部实验中所观察到的离散Xe结合位点和分子扩散通道，而且揭示了实验中无法观测的瞬态结合位点和多重气体迁移途径。上述结果表明人血红蛋白因其结构柔性所形成的瞬态通道对于气体分子迁移过程的重要性。除此之外，利用MSM和过渡路径理论（TPT）构建了人血红蛋白α链结构变化与气体分子迁移之间的关系，阐释了血红蛋白中影响气体迁移的关键结构及其微观机制。     

In Silico Design of MDM2-Targeting Peptides from a Naturally Occurring Constrained Peptide

Naturally occurring constrained peptides are frequently used as scaffolds for bioactive peptide grating due to their high stability. Here, we used in silico methods to design several constrained peptides comprising a scorpion toxin scaffold, a MDM2 binding epitope, and a cluster of positively charged residues. The designed peptides displayed varied binding affinity to MDM2 despite differing by only one or two residues. One of the peptides, SC426, had nanomolar binding affinity (K_D=6.6±2.6 nm ) to MDM2, and exhibited stronger inhibitory activity on the proliferation of HCT116 cells (p53-wild type) and SW480 cells (p53-mutant) than that of nutlin-3a. Binding mode analysis of the designed peptide at MDM2 suggests that the conserved “FWL” epitope was buried in the hydrophobic binding pocket, and the residues located at the periphery of the binding site contributed to the high binding affinity of SC426. Overall, in silico design of miniproteins with therapeutic potential through epitope grafting to the naturally occurring constrained peptide is an effective strategy.     

Loop flexibility and solvent dynamics as determinants for the selective inhibition of cyclin-dependent kinase 4: comparative molecular dynamics simulation studies of CDK2 and CDK4

The design and discovery of selective cyclin-dependent kinase 4 (CDK4) inhibitors have been actively pursued in order to develop therapeutic cancer treatments. By means of a consecutive computational protocol involving homology modeling, docking experiments, and molecular dynamics simulations, we examine the characteristic structural and dynamic properties that distinguish CDK4 from CDK2 in its complexation with selective inhibitors. The results for all three CDK4-selective inhibitors under investigation show that the large-amplitude motion of a disordered loop of CDK4 is damped out in the presence of the inhibitors whereas their binding in the CDK2 active site has little effect on the loop flexibility. It is also found that the binding preference of CDK4- selective inhibitors for CDK4 over CDK2 stems from the reduced solvent accessibility in the active site of the former due to the formation of a stable hydrogen-bond triad by the Asp99, Arg101, and Thr102 side chains at the top of the active-site gorge. Besides the differences in loop flexibility and solvent accessibility, the dynamic stabilities of the hydrogen bonds between the inhibitors and the side chain of the lysine residue at the bottom of the active site also correlate well with the relative binding affinities of the inhibitors for the two CDKs. These results highlight the usefulness of this computational approach in evaluating the selectivity of a CDK inhibitor, and demonstrate the necessity of considering protein flexibility and solvent effects in designing new selective CDK4-selective inhibitors.     

A Novel Fluorescent Probe for Retaining Galactosyltransferases

Glycosyltransferases (GTs) are a large class of carbohydrate‐active enzymes that are involved, in both pro‐ and eukaryotic organisms, in numerous important biological processes, from cellular adhesion to carcinogenesis. GTs have enormous potential as molecular targets for chemical biology and drug discovery. For the full realisation of this potential, operationally simple and generally applicable GT bioassays, especially for inhibitor screening, are indispensable tools. In order to facilitate the development of GT high‐throughput screening assays for the identification of GT inhibitors, we have developed novel, fluorescent derivatives of UDP‐galactose (UDP‐Gal) that are recognised as donor analogues by several different retaining galactosyltransferases (GalTs). We demonstrate for one of these derivatives that fluorescence emission is quenched upon specific binding to individual GalTs, and that this effect can be used as the read‐out in ligand‐displacement experiments. The novel fluorophore acts as an excellent sensor for several different enzymes and is suitable for the development of a new type of GalT bioassay, whose modular nature and operational simplicity will significantly facilitate inhibitor screening. Importantly, the structural differences between the natural donor UDP‐Gal and the new fluorescent derivatives are minimal, and the general assay principle described herein may therefore also be applicable to other GalTs and/or proteins that use nucleotides or nucleotide conjugates as their cofactor.     

Protein Contacts and Ligand Binding in the Inward‐Facing Model of Human P‐Glycoprotein

The primary aim of this work was to analyze the contacts between residues in the nucleotide binding domains (NBDs) and at the interface between the transmembrane domains (TMDs) and the NBDs in the inward‐open homology model of human P‐glycoprotein (P‐gp). The analysis revealed communication nets through hydrogen bonding in the NBD and at the NBD–TMD interface of each half involving residues from the adenosine triphosphate (ATP) motifs and the coupling helices of the intracellular loops. Similar networks have been identified in P‐gp conformations generated by molecular dynamics simulation. Differences have been recorded in the networking between both halves of P‐gp. Many of the residue contacts have also been observed in the X‐ray crystal structures of other ATP binding cassette (ABC) transporters, which confirms their validity. Next, possible binding pockets involving residues of importance for the TMD–NBD communication were identified. By studying these pockets, binding sites were suggested for rhodamine 123 (R‐site) and prazosin (regulatory site) at the NBD–TMD interface that agreed with the experimental data on their location. Additionally, one more R‐site in the protein cavity was proposed, in accordance with the available biochemical data. Together with the previously suggested Hoechst 33342 site (H‐site), all sites were interpreted with respect to their effects on the protein ATPase activity, in correspondence with the experimental observations. Several residues involved in key contacts in the P‐gp NBDs were proposed for further targeted mutagenesis experiments.     

Negative Cooperativity in NAD(P)H Quinone Oxidoreductase 1 (NQO1)

NAD(P)H quinone oxidoreductase-1 (NQO1) is a homodimeric protein that acts as a detoxifying enzyme or as a chaperone protein. Dicourmarol interacts with NQO1 at the NAD(P)H binding site and can both inhibit enzyme activity and modulate the interaction of NQO1 with other proteins. We show that the binding of dicoumarol and related compounds to NQO1 generates negative cooperativity between the monomers. This does not occur in the presence of the reducing cofactor, NAD(P)H, alone. Alteration of Gly150 (but not Gly149 or Gly174) abolished the dicoumarol-induced negative cooperativity. Analysis of the dynamics of NQO1 with the Gaussian network model indicates a high degree of collective motion by monomers and domains within NQO1. Ligand binding is predicted to alter NQO1 dynamics both proximal to the ligand binding site and remotely, close to the second binding site. Thus, drug-induced modulation of protein motion might contribute to the biological effects of putative inhibitors of NQO1.     

Molecular Dynamics Simulation of the Crystallizable Fragment of IgG1—Insights for the Design of Fcabs

An interesting format in the development of therapeutic monoclonal antibodies uses the crystallizable fragment of IgG1 as starting scaffold. Engineering of its structural loops allows generation of an antigen binding site. However, this might impair the molecule’s conformational stability, which can be overcome by introducing stabilizing point mutations in the CH3 domains. These point mutations often affect the stability and unfolding behavior of both the CH2 and CH3 domains. In order to understand this cross-talk, molecular dynamics simulations of the domains of the Fc fragment of human IgG1 are reported. The structure of human IgG1-Fc obtained from X-ray crystallography is used as a starting point for simulations of the wild-type protein at two different pH values. The stabilizing effect of a single point mutation in the CH3 domain as well as the impact of the hinge region and the glycan tree structure connected to the CH2 domains is investigated. Regions of high local flexibility were identified as potential sites for engineering antigen binding sites. Obtained data are discussed with respect to the available X-ray structure of IgG1-Fc, directed evolution approaches that screen for stability and use of the scaffold IgG1-Fc in the design of antigen binding Fc proteins.     

Understanding the Key Factors that Control the Inhibition of Type II Dehydroquinase by (2R)‐2‐Benzyl‐3‐dehydroquinic Acids

The binding mode of several substrate analogues, (2R)‐2‐benzyl‐3‐dehydroquinic acids 4 , which are potent reversible competitive inhibitors of type II dehydroquinase (DHQ2), the third enzyme of the shikimic acid pathway, has been investigated by structural and computational studies. The crystal structures of Mycobacterium tuberculosis and Helicobacter pylori DHQ2 in complex with one of the most potent inhibitor, p‐methoxybenzyl derivative 4 a , have been solved at 2.40 Å and 2.75 Å, respectively. This has allowed the resolution of the M. tuberculosis DHQ2 loop containing residues 20–25 for the first time. These structures show the key interactions of the aromatic ring in the active site of both enzymes and additionally reveal an important change in the conformation and flexibility of the loop that closes over substrate binding. The loop conformation and the binding mode of compounds 4 b – d has been also studied by molecular dynamics simulations, which suggest that the benzyl group of inhibitors 4 prevent appropriate orientation of the catalytic tyrosine of the loop for proton abstraction and disrupts its basicity.     

Effect of Ions and Sequence Variants on the Antagonist Binding Properties of the Histamine H1 Receptor

The histamine H₁ receptor (H₁R) is a G protein-coupled receptor (GPCR) and represents a main target in the treatment of allergic reactions as well as inflammatory reactions and depressions. Although the overall effect of antagonists on H₁ function has been extensively investigated, rather little is known about the potential modulatory effect of ions or sequence variants on antagonist binding. We investigated the dynamics of a phosphate ion present in the crystal structure and of a sodium ion, for which we determined the position in the allosteric pocket by metadynamics simulations. Both types of ions exhibit significant dynamics within their binding site; however, some key contacts remain stable over the simulation time, which might be exploited to develop more potent drugs targeting these sites. The dynamics of the ions is almost unaffected by the presence or absence of doxepin, as also reflected in their small effect (less than 1 kcal·mol⁻¹) on doxepin binding affinity. We also examined the effect of four H₁R sequence variants observed in the human population on doxepin binding. These variants cause a reduction in doxepin affinity of up to 2.5 kcal·mol⁻¹, indicating that personalized medical treatments that take into account individual mutation patterns could increase precision in the dosage of GPCR-targeting drugs.     

Chemoenzymatic synthesis of branched oligo- and polysaccharides as potential substrates for starch active enzymes

Oligo- and polysaccharides embodying the alpha-maltotriosyl-6(II)-maltotetraosyl structure were readily synthesized by transglycosylation of maltosyl fluoride onto panose and pullulan catalysed by the bacterial transglycosylase cyclodextrin glycosyltransferase (CGTase). The two products obtained proved useful for increasing the knowledge of substrate binding and processing at the active site of barley limit dextrinase that is involved in the metabolism of amylopectin by acting upon its branch points.     

A Second,Druggable Binding Site in UDP‐Galactopyranose Mutase from Mycobacterium tuberculosis?

UDP‐galactopyranose mutase (UGM), a key enzyme in the biosynthesis of mycobacterial cell walls, is a potential target for the treatment of tuberculosis. In this work, we investigate binding models of a non‐substrate‐like inhibitor, MS‐208, with M. tuberculosis UGM. Initial saturation transfer difference (STD) NMR experiments indicated a lack of direct competition between MS‐208 and the enzyme substrate, and subsequent kinetic assays showed mixed inhibition. We thus hypothesized that MS‐208 binds at an allosteric binding site (A‐site) instead of the enzyme active site (S‐site). A candidate A‐site was identified in a subsequent computational study, and the overall hypothesis was supported by ensuing mutagenesis studies of the A‐site. Further molecular dynamics studies led us to propose that MS‐208 inhibition occurs by preventing complete closure of an active site mobile loop that is necessary for productive substrate binding. The results suggest the presence of an A‐site with potential druggability, opening up new opportunities for the development of novel drug candidates against tuberculosis.     

A mutant D-amino acid aminotransferase with broad substrate specificity: construction by replacement of the interdomain loop Pro119- Arg120-Pro121 by Gly-Gly-Gly

D-amino acid aminotransferase (EC 2.6.1.21) catalyzes the interconversion of various D-amino acids and 2-oxo acids. Each homodimer subunit consists of two domains, which are connected by a single loop, Asn118-Pro119-Arg120-Pro121. The loop has no direct contact with the active site region or the cofactor, pyridoxal 5'- phosphate. We attempted to increase the conformational flexibility of this loop through a triple glycine substitution. The resultant mutant P119G-R120G-P121G has features clearly different from the wild-type enzyme under overall as well as half-reaction conditions. The pre- steady-state kinetic analyses of half reactions showed that the mutant enzyme has kmax values higher than the wild-type enzyme towards most D- amino acids examined. A concomitant decrease in substrate affinity (1/Kd), particularly for acidic amino acids, was also observed. A putative binding site for the distal carboxyl group of acidic amino acids in the wild-type enzyme was incidentally displaced by the loop mutation, indicating a functional linkage between the interdomain loop and the active site region. This study has exemplified the usefulness of engineering relatively distant loops as a means to modify substrate specificity of an enzyme.      

Specificity of Donor Structures for endo‐β‐N‐Acetylglucosaminidase‐Catalyzed Transglycosylation Reactions

To demonstrate the structural specificity of the glycosyl donor for the transglycosylation reaction by using endo‐β‐N‐acetylglucosaminidase from Mucor hiemalis (endo‐M), a series of tetrasaccharide oxazoline derivatives was synthesized. These derivatives correspond to the core structure of an asparagine‐linked glycoprotein glycan with a β‐mannose unit of a non‐natural‐type monosaccharide, including β‐glucose, β‐galactose, and β‐talose in place of the β‐mannose moiety. The transglycosylation activity of wildtype (WT) endo‐M and two mutants, N175Q and N175A, was examined by using these tetrasaccharide donors with p‐nitrophenyl N‐acetylglucosaminide (GlcNAc‐pNp). The essential configuration of the hydroxy group for the transglycosylation reaction was determined. On the basis of these results, the transglycosylation reaction was investigated by using chemically modified donors, and transglycosylated products were successfully obtained.