Surasak Chunsrivirot

Free energy of binding of a small molecule to an amorphous polymer in a solvent.

Crystallization is a commonly used purification process in industrial practice. It usually begins with heterogeneous nucleation on a foreign surface. The complicated mechanism of heterogeneous nucleation is not well understood, but we hypothesize that a possible correlation between binding affinity to a surface and nucleation enhancement might exist. Amorphous polymers have been used in controlling crystallization. However, to our knowledge, no attempt has been made to calculate the free energy of binding of a small molecule to an amorphous polymer in a solvent, and to characterize the binding sites/conformations of this system at a molecular level. We developed a two-step approach, first using Adsorption Locator to identify probable binding sites and molecular dynamics to screen for the best binding sites and then using the Blue-Moon Ensemble method to compute the free energy of binding. A system of ethylene glycol, polyvinyl alcohol (PVA), and heavy water (D(2)O) was used for validation, since experimental data exists on a related system. Looking at four independently constructed surfaces, we found that ethylene glycol binds to an indentation on the surface or in a hole beneath the surface. We focused on the indentation binding sites because they are easily accessible and do not have large free energy barriers. The closest system for which experimental data on binding energetics exists is ethylene glycol on PVA in aqueous solutions/gels, and the magnitudes of the free energy of binding to the three best indentation binding sites are close to the experimental value, 0.4-3.7 kcal/mol higher. Our approach offers a way to compute the free energy of binding and characterize the binding sites/conformations, and is general enough to apply to other small molecule/amorphous polymer/solvent systems.

Binding affinity of a small molecule to an amorphous polymer in a solvent. Part 1: free energy of binding to a binding site.

Crystallization is commonly used in a separation and purification process in the production of a wide range of materials in various industries. In industry, crystallization usually starts with heterogeneous nucleation on a foreign surface. The complicated mechanism of heterogeneous nucleation is not well understood; however, we hypothesize that there might be a possible correlation between binding affinity to a surface and enhancement of nucleation. Recent studies show that amorphous polymers can be used to control crystallization, selectively produce pharmaceutical polymorphs, and discover novel pharmaceutical polymorphs. To investigate the possible correlation between the binding affinity of one molecule to key binding sites (local binding) and heterogeneous nucleation activity as well as the possibility of using this binding affinity to help guide the selection of polymers that promote heterogeneous nucleation, we computed the free energy of binding of aspirin to four nonporous cross-linked polymers in an ethanol-water 38 v% mixture. These cross-linked polymers are poly(4-acryloylmorpholine) (PAM), poly(2-carboxyethyl acrylate) (PCEA), poly(4-hydroxylbutyl acrylate) (PHBA), and polystyrene (PS); all of them were cross-linked with divinylbenzene (DVB). These systems were used because their heterogeneous nucleation activities are available in literature, and the ranking is PAM > PCEA > PHBA ≈ PS. We generated three independent surfaces for each polymer and computed the free energy of binding of aspirin to the best binding site that we found on each surface. The average free energies of binding to the best sites of PAM, PCEA, PHBA, and PS are -20.4 ± 1.0, -16.7 ± 1.0, -14.4 ± 1.1, and -13.6 ± 1.1 kcal/mol, respectively. We found that the trend of the magnitudes of the average free energies of binding to the best sites is PAM > PCEA > PHBA ≈ PS. This trend is very similar to that of heterogeneous nucleation activity. Our results suggest the importance of the free energy of binding to key sites (local binding) and the possibility of using this quantity to help guide the selection of polymers that promote heterogeneous nucleation.

Design of peptides as inhibitors of human papillomavirus 16 transcriptional regulator E1‐E2

Human papillomavirus 16 (HPV 16) is a DNA virus that is capable of infecting humans and causing cervical cancer. HPV16 E2 plays an important role in viral gene regulation. This work aims to predict the binding conformations and interactions between the dodecapeptides and HPV16 E2 as well as to design novel peptide inhibitors that are capable of binding to HPV16 E2 and disrupt the transcriptional regulator E1–E2 complex formation, using computational protein design techniques. Based on previously reported peptide4 (TWFWPYPYPHLP), novel peptide inhibitors were designed and five peptides that showed lower binding energy to HPV16 E2 than that of peptide4, were selected for in vitro experiments. Enzyme‐linked immunosorbent (ELISA) assay showed that Y6R, W4H_Y6R, and W4H peptides bound to HPV16 E2 with higher affinity than peptide4 did. Moreover, Y6R, W4H_Y6R, and W4H peptides more effectively inhibited E1–E2 complex formation than peptide4. This work revealed important interactions between the peptides and E1–E2 complex, suggesting a strategy for development of more potent peptide inhibitors.

Replica exchange molecular dynamics simulations reveal the structural and molecular properties of levan-type fructo-oligosaccharides of various chain lengths

BackgroundLevan and levan-type fructo-oligosaccharides (LFOs) have various potential applications in pharmaceutical and food industries due to their beneficial properties such as their low intrinsic viscosity and high water solubility. Previous studies showed that they exhibited prebiotic effects, anti-inflammatory and anti-tumor activities against Sarcoma-180 tumor cells of human. Despite their various potential applications, the structural and molecular properties of LFOs of various chain lengths are not well understood.ResultsWe employed the replica-exchange molecular dynamics simulations method (REMD) in AMBER14 to elucidate structural and molecular properties of LFOs with chain lengths of 5 (LFO5), 10 (LFO10) and 15 (LFO15) residues in two models of generalized Born implicit solvent (GBHCT and GBOBC1). For LFO10 and LFO15, four distinct conformations (helix-like, partial helix, zig-zag and random structures) were characterized by their upper-middle and lower-middle torsions. For LFO5, two distinct conformations (partial helix and random structures) were characterized by their middle torsion and molecular angle of residues 1, 3 and 5. To determine hydrogen bonds important for the formation of helix-like structures of LFO10 and LFO15, occurrence frequencies of hydrogen bonds were analyzed, and the O6(i)--H3O(i+1) hydrogen bond was found with the highest frequency, suggesting its importance in helix formation. Among three dihedral angles between two fructosyl units [ϕ (O5’-C2’-O6-C6), ψ (C2’-O6-C6-C5) and ω (O6-C6-C5-C4)], dihedral angle distributions showed that ω was the most flexible dihedral angle and probably responsible for conformational differences of LFOs.ConclusionsOur study provides important insights into the structural and molecular properties of LFOs, which tend to form helical structures as the chain length increases from 5 to 15 residues. This information could be beneficial for the selection of LFOs with appropriate lengths and properties for pharmaceutical and biological applications.

Molecular dynamics provides insight into how N251A and N251Y mutations in the active site of Bacillus licheniformis RN-01 levansucrase disrupt production of long-chain levan

Produced by levansucrase, levan and levan oligosaccharides (GFn) have potential applications in food and pharmaceutical industries such as prebiotics, anti-tumor and anti-inflammatory agents. Previous study reported that Bacillus licheniformis RN-01 levansucrase could produce levan oligosaccharides and long-chain levan. However, its N251A and N251Y mutants could effectively produce short-chain oligosaccharides upto GF3, but they could not produce long-chain levan. We hypothesized that these mutations probably reduced GF3 binding affinity in levansucrase active site that contains fructosyl-Asp93 intermediate and caused GF3 to be in an unfavorable orientation for transfructosylation; therefore, levansucrase could not effectively extend GF3 by one fructosyl residue to produce GF4 and subsequently long-chain levan. However, these mutations probably did not significantly reduce binding affinity or drastically change orientation of GF2; therefore, levansucrase could still extend GF2 to produce GF3. Using this hypothesis, we employed molecular dynamics to investigate effects of these mutations on GF2/GF3 binding in levansucrase active site. Our results reasonably support this hypothesis as N251A and N251Y mutations did not significantly reduce GF2 binding affinity, as calculated by MM-GBSA technique and hydrogen bond occupations, or drastically change orientation of GF2 in levansucrase active site, as measured by distance between atoms necessary for transfructosylation. However, these mutations drastically decreased GF3 binding affinity and caused GF3 to be in an unfavorable orientation for transfructosylation. Furthermore, the free energy decomposition and hydrogen bond occupation results suggest the importance of Arg255 in GF2/GF3 binding in levansucrase active site. This study provides important and novel insight into the effects of N251A and N251Y mutations on GF2/GF3 binding in levansucrase active site and how they may disrupt production of long-chain levan. This knowledge could be beneficial in designing levansucrase to efficiently produce levan oligosaccharides with desired length.

Unraveling the structural and molecular properties of 34-residue levans with various branching degrees by replica exchange molecular dynamics simulations

Levan has various potential applications in the pharmaceutical and food industries, such as cholesterol-lowering agents and prebiotics, due to its beneficial properties, which depend on its length and branching degree. A previous study also found that the branching degree of levan affected anti-tumor activities against SNU-1 and HepG2 tumor cell lines. Despite its promising potential, the properties of levans with different branching degrees are not well understood at the molecular level. In two models of the generalized Born implicit solvent (GBHCT and GBOBC1), we employed replica-exchange molecular dynamics simulations to explore conformational spaces of 34-residue levans (L34) with branching degrees of zero (LFO34B0), one (LFO34B1), three (LFO34B3) and five (LFO34B5), as well as to elucidate their structural and molecular properties. To ensure a fair comparison of the effects of branching degree on these properties, we focused on analyzing the properties of the central 21-residue of the main chains of all systems. Our results show that all major representative conformations tend to form helix-like structures with kinks, where two-kink helix-like structures have the highest population. As branching degree increases, the population of helix-like structures with zero or one kink tends to increase slightly. As the number of kinks in the structures with the same branching degree increases, the average values of the lengths and angles among centers of masses of three consecutive turns of residue i, i+3, and i+6 tended to decrease. Due to its highest occurring frequencies, the O6 (i)—H3O (i+1) hydrogen bond could be important for helix-like structure formation. Moreover, hydrogen bonds forming among the branching residue (br), branching position (bp) and other residues of L34B1, L34B3 and L34B5 were identified. The O1(bp)—H3O(br), O1(br)—H3O(br) and O5(br)—H1O(br) hydrogen bonds were found in the first-, second- and third-highest occurrence frequencies, respectively. Our study provides novel and important insights into conformational spaces and the structural and molecular properties of 34-residue levans with various branching degrees, which tend to form helix-like structures with kinks.

Insight into the substrate specificity change caused by the Y227H mutation of α-glucosidase III from the European honeybee (Apis mellifera) through molecular dynamics simulations

Honey from the European honeybee, Apis mellifera, is produced by α-glucosidases (HBGases) and is widely used in food, pharmaceutical, and cosmetic industries. Categorized by their substrate specificities, HBGases have three isoforms: HBGase I, II and III. Previous experimental investigations showed that wild-type HBGase III from Apis mellifera (WT) preferred sucrose to maltose as a substrate, while the Y227H mutant (MT) preferred maltose to sucrose. This mutant can potentially be used for malt hydrolysis because it can efficiently hydrolyze maltose. In this work, to elucidate important factors contributing to substrate specificity of this enzyme and gain insight into how the Y227H mutation causes substrate specificity change, WT and MT homology models were constructed, and sucrose/maltose was docked into active sites of the WT and MT. AMBER14 was employed to perform three independent molecular dynamics runs for these four complexes. Based on the relative binding free energies calculated by the MM-GBSA method, sucrose is better than maltose for WT binding, while maltose is better than sucrose for MT binding. These rankings support the experimentally observed substrate specificity that WT preferred sucrose to maltose as a substrate, while MT preferred maltose to sucrose, suggesting the importance of binding affinity for substrate specificity. We also found that the Y227H mutation caused changes in the proximities between the atoms necessary for sucrose/maltose hydrolysis that may affect enzyme efficiency in the hydrolysis of sucrose/maltose. Moreover, the per-residue binding free energy decomposition results show that Y227/H227 may be a key residue for preference binding of sucrose/maltose in the WT/MT active site. Our study provides important and novel insight into the binding of sucrose/maltose in the active site of Apis mellifera HBGase III and into how the Y227H mutation leads to the substrate specificity change at the molecular level. This knowledge could be beneficial in the design of this enzyme for increased production of desired products.