Free‐energy function based on an all‐atom model for proteins

We have developed a free‐energy function based on an all‐atom model for proteins. It comprises two components, the hydration entropy (HE) and the total dehydration penalty (TDP). Upon a transition to a more compact structure, the number of accessible configurations arising from the translational displacement of water molecules in the system increases, leading to a water‐entropy gain. To fully account for this effect, the HE is calculated using a statistical‐mechanical theory applied to a molecular model for water. The TDP corresponds to the sum of the hydration energy and the protein intramolecular energy when a fully extended structure, which possesses the maximum number of hydrogen bonds with water molecules and no intramolecular hydrogen bonds, is chosen as the standard one. When a donor and an acceptor (e.g., N and O, respectively) are buried in the interior after the break of hydrogen bonds with water molecules, if they form an intramolecular hydrogen bond, no penalty is imposed. When a donor or an acceptor is buried with no intramolecular hydrogen bond formed, an energetic penalty is imposed. We examine all the donors and acceptors for backbone‐backbone, backbone‐side chain, and side chain‐side chain intramolecular hydrogen bonds and calculate the TDP. Our free‐energy function has been tested for three different decoy sets. It is better than any other physics‐based or knowledge‐based potential function in terms of the accuracy in discriminating the native fold from misfolded decoys and the achievement of high Z‐scores. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

On the satisfaction of backbone‐carbonyl lone pairs of electrons in protein structures

Protein structures are stabilized by a variety of noncovalent interactions (NCIs), including the hydrophobic effect, hydrogen bonds, electrostatic forces and van der Waals’ interactions. Our knowledge of the contributions of NCIs, and the interplay between them remains incomplete. This has implications for computational modeling of NCIs, and our ability to understand and predict protein structure, stability, and function. One consideration is the satisfaction of the full potential for NCIs made by backbone atoms. Most commonly, backbone‐carbonyl oxygen atoms located within α‐helices and β‐sheets are depicted as making a single hydrogen bond. However, there are two lone pairs of electrons to be satisfied for each of these atoms. To explore this, we used operational geometric definitions to generate an inventory of NCIs for backbone‐carbonyl oxygen atoms from a set of high‐resolution protein structures and associated molecular‐dynamics simulations in water. We included more‐recently appreciated, but weaker NCIs in our analysis, such as n→π* interactions, Cα‐H bonds and methyl‐H bonds. The data demonstrate balanced, dynamic systems for all proteins, with most backbone‐carbonyl oxygen atoms being satisfied by two NCIs most of the time. Combinations of NCIs made may correlate with secondary structure type, though in subtly different ways from traditional models of α‐ and β‐structure. In addition, we find examples of under‐ and over‐satisfied carbonyl‐oxygen atoms, and we identify both sequence‐dependent and sequence‐independent secondary‐structural motifs in which these reside. Our analysis provides a more‐detailed understanding of these contributors to protein structure and stability, which will be of use in protein modeling, engineering and design.     

Asparagine and glutamine differ in their propensities to form specific side chain-backbone hydrogen bonded motifs in proteins

Short range side chain‐backbone hydrogen bonded motifs involving Asn and Gln residues have been identified from a data set of 1370 protein crystal structures (resolution ≤ 1.5 Å). Hydrogen bonds involving residues i ? 5 to i + 5 have been considered. Out of 12,901 Asn residues, 3403 residues (26.4%) participate in such interactions, while out of 10,934 Gln residues, 1780 Gln residues (16.3%) are involved in these motifs. Hydrogen bonded ring sizes (C_n, where n is the number of atoms involved), directionality and internal torsion angles are used to classify motifs. The occurrence of the various motifs in the contexts of protein structure is illustrated. Distinct differences are established between the nature of motifs formed by Asn and Gln residues. For Asn, the most highly populated motifs are the C₁₀ (CO^δ_i …NH_{i + 2}), C₁₃ (CO^δ_i …NH_{i + 3}) and C₁₇ (N^δH_i …CO_{i ? 4}) structures. In contrast, Gln predominantly forms C₁₆ (CO^ε_i …NH_{i ? 3}), C₁₂ (N^εH_i …CO_{i ? 2}), C₁₅ (N^εH_i …CO_{i ? 3}) and C₁₈ (N^εH_i …CO_{i ? 4}) motifs, with only the C₁₈motif being analogous to the Asn C₁₇structure. Specific conformational types are established for the Asn containing motifs, which mimic backbone β‐turns and α‐turns. Histidine residues are shown to serve as a mimic for Asn residues in side chain‐backbone hydrogen bonded ring motifs. Illustrative examples from protein structures are considered. Proteins 2012; © 2011 Wiley Periodicals, Inc.     

Geometric characteristics of hydrogen bonds involving sulfur atoms in proteins

Sulfur atoms have been known to participate in hydrogen bonds (H‐bonds) and these sulfur‐containing H‐bonds (SCHBs) are suggested to play important roles in certain biological processes. This study aims to comprehensively characterize all the SCHBs in 500 high‐resolution protein structures (≤1.8 Å). We categorized SCHBs into six types according to donor/acceptor behaviors and used explicit hydrogen approach to distinguish SCHBs from those of nonhydrogen bonding interactions. It is revealed that sulfur atom is a very poor H‐bond acceptor, but a moderately good H‐bond donor. In α‐helix, considerable SCHBs were found between the sulphydryl group of cysteine residue i and the carbonyl oxygen of residue i‐4, and these SCHBs exert effects in stabilizing helices. Although for other SCHBs, they possess no specific secondary structural preference, their geometric characteristics in proteins and in free small compounds are significantly distinct, indicating the protein SCHBs are geometrically distorted. Interestingly, sulfur atom in the disulfide bond tends to form bifurcated H‐bond whereas in cysteine‐cysteine pairs prefer to form dual H‐bond. These special H‐bonds remarkably boost the interaction between H‐bond donor and acceptor. By oxidation/reduction manner, the mutual transformation between the dual H‐bonds and disulfide bonds for cysteine‐cysteine pairs can accurately adjust the structural stability and biological function of proteins in different environments. Furthermore, few loose H‐bonds were observed to form between the sulphydryl groups and aromatic rings, and in these cases the donor H is almost over against the rim rather than the center of the aromatic ring. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

Helicity of short E-R/K peptides

Understanding the secondary structure of peptides is important in protein folding, enzyme function, and peptide‐based drug design. Previous studies of synthetic Ala‐based peptides (>12 a.a.) have demonstrated the role for charged side chain interactions involving Glu/Lys or Glu/Arg spaced three (i, i + 3) or four (i, i + 4) residues apart. The secondary structure of short peptides (<9 a.a.), however, has not been investigated. In this study, the effect of repetitive Glu/Lys or Glu/Arg side chain interactions, giving rise to E‐R/K helices, on the helicity of short peptides was examined using circular dichroism. Short E‐R/K–based peptides show significant helix content. Peptides containing one or more E‐R interactions display greater helicity than those with similar E‐K interactions. Significant helicity is achieved in Arg‐based E‐R/K peptides eight, six, and five amino acids long. In these short peptides, each additional i + 3 and i + 4 salt bridge has substantial contribution to fractional helix content. The E‐R/K peptides exhibit a strongly linear melt curve indicative of noncooperative folding. The significant helicity of these short peptides with predictable dependence on number, position, and type of side chain interactions makes them an important consideration in peptide design.     

Delineation of a new structural motif involving NHN γ-turn

Macromolecules are characterized by distinctive arrangement of hydrogen bonds. Different patterns of hydrogen bonds give rise to distinct and stable structural motifs. An analysis of 4114 non-redundant protein chains reveals the existence of a three-residue, (i − 1) to (i + 1), structural motif, having two hydrogen-bonded five-membered pseudo rings (the first, an N H···OC involving the first residue, and the second being N H∙∙∙N involving the last two residues), separated by a peptide bond. There could be an additional hydrogen bond between the side-chain at (i-1) and the main-chain NH of (i + 1). The average backbone torsion angles of −76(±21)° and – 12(±17)° at i creates a tight turn in the polypeptide chain, akin to a γ-turn. Indeed, a search of three-residue fragments with restriction on the terminal C^α···C^α distance and the existence of the two pseudo rings on either side revealed the presence 14 846 cases of a variant, termed NHN γ-turn, distinct from the NHO γ-turn (2032 cases) that has traditionally been characterized by the presence of NHO hydrogen bond linking the terminal main-chain atoms. As in the latter, the newly identified γ-turns are also of two types—classical and inverse, occurring in the ratio of 1:6. The propensities of residues to occur in these turns and their secondary structural features have been enumerated. An understanding of these turns would be useful for structure prediction and loop modeling, and may serve as models to represent some of the unfolded state or disordered region in proteins.     

Water and side‐chain embedded π‐turns

Elucidating protein function from its structure is central to the understanding of cellular mechanisms. This involves deciphering the dependence of local structural motifs on sequence. These structural motifs may be stabilized by direct or water‐mediated hydrogen bonding among the constituent residues. π‐Turns, defined by interactions between (i) and (i + 5) positions, are large enough to contain a central space that can embed a water molecule (or a protein moiety) to form a stable structure. This work is an analysis of such embedded π‐turns using a nonredundant dataset of protein structures. A total of 2965 embedded π‐turns have been identified, as also 281 embedded Schellman motif, a type of π‐turn which occurs at the C‐termini of α‐helices. Embedded π‐turns and Schellman motifs have been classified on the basis of the protein atoms of the terminal turn residues that are linked by the embedded moiety, conformation, residue composition, and compared with the turns that have terminal residues connected by direct hydrogen bonds. Geometrically, the turns have been fitted to a circle and the position of the linker relative to its center analyzed. The hydroxyl group of Ser and Thr, located at (i + 3) position, is the most prominent linker for the side‐chain mediated π‐turns. Consideration of residue conservation among homologous sequences indicates the terminal and the linker positions to be the most conserved. The embedded π‐turn as a binding site (for the linker) is discussed in the context of “nest,” a concave depression that is formed in protein structures with adjacent residues having enantiomeric main‐chain conformations. © 2013 Wiley Periodicals, Inc. Biopolymers 101: 441–453, 2014.     

Green‐lighting green fluorescent protein: Faster and more efficient folding by eliminating a cis–trans peptide isomerization event

Wild‐type green fluorescent protein (GFP) folds on a time scale of minutes. The slow step in folding is a cis–trans peptide bond isomerization. The only conserved cis‐peptide bond in the native GFP structure, at P89, was remodeled by the insertion of two residues, followed by iterative energy minimization and side chain design. The engineered GFP was synthesized and found to fold faster and more efficiently than its template protein, recovering 50% more of its fluorescence upon refolding. The slow phase of folding is faster and smaller in amplitude, and hysteresis in refolding has been eliminated. The elimination of a previously reported kinetically trapped state in refolding suggests that X‐P89 is trans in the trapped state. A 2.55 Å resolution crystal structure revealed that the new variant contains only trans‐peptide bonds, as designed. This is the first instance of a computationally remodeled fluorescent protein that folds faster and more efficiently than wild type.     

Computational characterization of the chemical step in the GTP hydrolysis by Ras‐GAP for the wild‐type and G13V mutated Ras

The free energy profiles for the chemical reaction of the guanosine triphosphate hydrolysis GTP + H₂O → GDP + Pi by Ras‐GAP for the wild‐type and G13V mutated Ras were computed by using molecular dynamics protocols with the QM(ab initio)/MM potentials. The results are consistent with the recent measurements of reaction kinetics in Ras‐GAP showing about two‐order reduction of the rate constant upon G13V mutation in Ras: the computed activation barrier on the free energy profile is increased by 3 kcal/mol upon the G13V replacement. The major reason for a higher energy barrier is a shift of the “arginine finger” (R789 from GAP) from the favorable position in the active site. The results of simulations provide support for the mechanism of the reference reaction according to which the Q61 side chain directly participates in chemical transformations at the proton transfer stage. Proteins 2015; 83:1046–1053. © 2015 Wiley Periodicals, Inc.     

ω‐Turn: A novel β‐turn mimic in globular proteins stabilized by main‐chain to side‐chain CH···O interaction

Mimicry of structural motifs is a common feature in proteins. The 10‐membered hydrogen‐bonded ring involving the main‐chain C?O in a β‐turn can be formed using a side‐chain carbonyl group leading to Asx‐turn. We show that the N? H component of hydrogen bond can be replaced by a C^γ‐H group in the side chain, culminating in a nonconventional C? H···O interaction. Because of its shape this β‐turn mimic is designated as ω‐turn, which is found to occur ～three times per 100 residues. Three residues (i to i + 2) constitute the turn with the C? H···O interaction occurring between the terminal residues, constraining the torsion angles ?_{i + 1}, ψ_{i + 1}, ?_{i + 2} and χ′_{1(i + 2)} (using the interacting C^γ atom). Based on these angles there are two types of ω‐turns, each of which can be further divided into two groups. C^β‐branched side‐chains, and Met and Gln have high propensities to occur at i + 2; for the last two residues the carbonyl oxygen may participate in an additional interaction involving the S and amino group, respectively. With Cys occupying the i + 1 position, such turns are found in the metal‐binding sites. N‐linked glycosylation occurs at the consensus pattern Asn‐Xaa‐Ser/Thr; with Thr at i + 2, the sequence can adopt the secondary structure of a ω‐turn, which may be the recognition site for protein modification. Location between two β‐strands is the most common occurrence in protein tertiary structure, and being generally exposed ω‐turn may constitute the antigenic determinant site. It is a stable scaffold and may be used in protein engineering and peptide design. Proteins 2015; 83:203–214. © 2014 Wiley Periodicals, Inc.     

Structural changes induced by the deamidation and isomerization of asparagine revealed by the crystal structure of Ustilago sphaerogena ribonuclease U2B

Under physiological conditions, the deamidation and isomerization of asparagine to isoaspartate (isoAsp) proceeds nonenzymatically via succinimide. Although a large number of proteins have been reported to contain isoAsp, information concerning the three‐dimensional structure of proteins containing isoaspartate is still limited. We have crystallized isoAsp containing Ustilago sphaerogena ribonuclease U2B, and determined the crystal structure at 1.32 Å resolution. The structure revealed that the formation of isoAsp32 induces a single turn unfolding of the α‐helix from Asp29 to Asp34, and the region from Asp29 to Arg35 forms a U‐shaped loop structure. The electron density map shows that isoAsp32 retained the L‐configuration at the C_α atom. IsoAsp32 is in gauche conformation about a C_α? C_β bond, and the polypeptide chain bends by ～90° at isoAsp32. IsoAsp32 protrudes from the surface of the protein, and the abnormal β‐peptide bond in the main‐chain and α‐carboxylate in the side‐chain is fully exposed. The structure suggests that the deamidation of the Asn and the isoAsp formation in proteins could confer immunogenicity. © 2010 Wiley Periodicals, Inc. Biopolymers 93: 1003–1010, 2010.     

Characterization of a type of β‐bend ribbon spiral generated by the repeating (Xaa–Yaa–Aib–Pro) motif: The solution structure of harzianin HC IX,a 14‐residue peptaibol forming voltage‐dependent ion channels

The three‐dimensional solution structure of harzianin HC IX, a peptaibol antibiotic isolated from the fungus Trichoderma harzianum, was determined using CD, homonuclear, and heteronuclear two‐dimensional nmr spectroscopy combined with molecular modeling. This 14‐residue peptide, Ac Aib¹ Asn² Leu³ Aib⁴ Pro⁵ Ala⁶ Ile⁷ Aib⁸ Pro⁹ Iva¹⁰ Leu¹¹ Aib¹² Pro¹³ Leuol¹⁴ (Aib, α‐aminoisobutyric acid; Iva, isovaline; Leuol, leucinol), is a main representative of a short‐sequence peptaibol class characterized by an acetylated N‐terminus, a C‐terminal amino alcohol, and the presence of three Aib‐L ‐Pro motifs at positions 4–5, 8–9, and 12–13, separated by two dipeptide units. In spite of a lower number of residues, compared to the 18/20‐residue peptaibols such as alamethicin, harzianin HC IX exhibits remarkable membrane‐perturbing properties. It interacts with phospholipid bilayers, increasing their permeability and forming voltage‐gated ion channels through a mechanism slightly differing from that proposed for alamethicin. Sequence‐specific ¹H‐ and ¹³C‐nmr assignments and conformational nmr parameters (³J_NHCαH coupling constants, quantitative nuclear Overhauser enhancement data, temperature coefficients of amide and carbonyl groups, NH–ND exchange rates) were obtained in methanol solution. Sixty structures were calculated based on 98 interproton distance restraints and 6 Φ dihedral angle restraints, using high temperature restrained molecular dynamics and energy minimization. Thirty‐seven out of the sixty generated structures were consistent with the nmr data and were convergent. The peptide backbone consists in a ribbon of overlapping β‐turns twisted into a continuous spiral from Asn² to Leuol¹⁴ and forming a 26 Å long helix‐like structure. This structure is slightly amphipathic, with the three Aib–Pro motifs aligned on the less hydrophobic face of the spiral where the Asn² side chain is also present, while the more hydrophobic bulky side chains of leucines, isoleucine, isovaline, and leucinol are located on the concave side. The repetitive (Xaa–Yaa–Aib–Pro) tetrapeptide subunit, making up the peptide sequence, is characterized by four sets of (Φ,Ψ) torsional angles, with the following mean values: Φ_i = −90°, Ψ_i = −27°; Φ_i+1 = −98°, Ψ_i+1 = −17°; Φ_i+2 = −49°, Ψ_i+2 = −50°; Φ_i+3 = −78°, Ψ_i+3 = +3°. We term this particular structure, specifically occurring in the case of (Xaa–Yaa–Aib–Pro)_n sequences, the (Xaa–Yaa–Aib–Pro)‐β‐bend ribbon spiral. It is stabilized by 4 → 1 intramolecular hydrogen bonds and differs from both the canonical 3₁₀‐helix made of a succession of type III β‐turns and from the β‐bend ribbon spiral that has been described in the case of (Aib–Pro)n peptide segments. © 1999 John Wiley & Sons, Inc. Biopoly 50: 71–85, 1999     

Physical�Cchemical determinants of coil conformations in globular proteins

We present a method with the potential to generate a library of coil segments from first principles. Proteins are built from α‐helices and/or β‐strands interconnected by these coil segments. Here, we investigate the conformational determinants of short coil segments, with particular emphasis on chain turns. Toward this goal, we extracted a comprehensive set of two‐, three‐, and four‐residue turns from X‐ray–elucidated proteins and classified them by conformation. A remarkably small number of unique conformers account for most of this experimentally determined set, whereas remaining members span a large number of rare conformers, many occurring only once in the entire protein database. Factors determining conformation were identified via Metropolis Monte Carlo simulations devised to test the effectiveness of various energy terms. Simulated structures were validated by comparison to experimental counterparts. After filtering rare conformers, we found that 98% of the remaining experimentally determined turn population could be reproduced by applying a hydrogen bond energy term to an exhaustively generated ensemble of clash‐free conformers in which no backbone polar group lacks a hydrogen‐bond partner. Further, at least 90% of longer coil segments, ranging from 5‐ to 20 residues, were found to be structural composites of these shorter primitives. These results are pertinent to protein structure prediction, where approaches can be divided into either empirical or ab initio methods. Empirical methods use database‐derived information; ab initio methods rely on physical–chemical principles exclusively. Replacing the database‐derived coil library with one generated from first principles would transform any empirically based method into its corresponding ab initio homologue.     

In silico study on the effect of surface lysines and arginines on the electrostatic interactions and protein stability

Charged amino acids are mostly exposed on a protein surface, thereby forming a network of interactions with the surrounding amino acids as well as with water. In particular, positively charged arginine and lysine have different side chain geometries and provide a different number of potential electrostatic interactions. This study reports a comparative analysis of the difference in the number of two representative electrostatic interactions, such as salt-bridges and hydrogen bonds, contributed by surface arginine and lysine, as well as their effect on protein stability using molecular modeling and dynamics simulation techniques. Two in silico variants, the R variant with all arginines and the K variant with all lysines on the protein surface, were modeled by mutating all the surface lysines to arginines and the surface arginines to lysines, respectively, for each of the 10 model proteins. A structural comparison of the respective two variants showed that the majority of R variants possessed more salt-bridges and hydrogen bond interactions than the K variants, indicating that arginine provides a higher probability of electrostatic interactions than lysine owing to its side chain geometry. Molecular dynamics simulations of these variants revealed the R variants to be more stable than the K variants at room temperature but this effect was not prominent under protein denaturating conditions, such as 353 and 333 K with 8 M urea. These results suggest that the arginine residues on a protein surface contribute to the protein stability slightly more than lysine by enhancing the electrostatic interactions.     

Studying the Structural and Folding Features of Long‐Sequence Trichobrachin Peptides

In this theoretical study, the folding processes of long‐sequence trichobrachin peptides (i.e., TB IIb peptides) were investigated by molecular dynamics methods. The formation of various helical structures (i.e., 3₁₀‐, α‐, and left‐handed α‐helices) was studied with regard to the entire sequence of peptides, as well as to each amino acid. The results pointed out that TB IIb molecules showed a propensity to form helical conformations, and they could be characterized by 3₁₀‐helical structure rather than by α‐helical structure. The formation of local (i.e., i←i+3 and i←i+4) as well as of non‐local (i.e., i←i+n, where n>4; and all i→i+n) H‐bonds was also examined. The results revealed that the occurrence of local, helix‐stabilizing H‐bonds was in agreement with the appearance of helical conformations, and the non‐local H‐bonds did not produce relevant effects on the evolution of helical structures. Based on the data obtained by our structural investigation, differences were observed between the TB IIb peptides, according to the type of amino acid located in the 17th position of their sequences. In summary, the folding processes were explored for TB IIb molecules, and our theoretical study led to the conclusion that these long‐sequence peptaibols showed characteristic structural and folding features.     

Structural insights into substrate and coenzyme preference by SDR family protein Gox2253 from Gluconobater oxydans

Gox2253 from Gluconobacter oxydans belongs to the short‐chain dehydrogenases/reductases family, and catalyzes the reduction of heptanal, octanal, nonanal, and decanal with NADPH. To develop a robust working platform to engineer novel G. oxydans oxidoreductases with designed coenzyme preference, we adopted a structure based rational design strategy using computational predictions that considers the number of hydrogen bonds formed between enzyme and docked coenzyme. We report the crystal structure of Gox2253 at 2.6 Å resolution, ternary models of Gox2253 mutants in complex with NADH/short‐chain aldehydes, and propose a structural mechanism of substrate selection. Molecular dynamics simulation shows that hydrogen bonds could form between 2′‐hydroxyl group in the adenosine moiety of NADH and the side chain of Gox2253 mutant after arginine at position 42 is replaced with tyrosine or lysine. Consistent with the molecular dynamics prediction, Gox2253‐R42Y/K mutants can use both NADH and NADPH as a coenzyme. Hence, the strategies here could provide a practical platform to engineer coenzyme selectivity for any given oxidoreductase and could serve as an additional consideration to engineer substrate‐binding pockets. Proteins 2014; 82:2925–2935. © 2014 Wiley Periodicals, Inc.     

VCD Robustness of the Amide‐I and Amide‐II Vibrational Modes of Small Peptide Models

The rotational strengths and the robustness values of amide‐I and amide‐II vibrational modes of For(AA)_nNHMe (where AA is Val, Asn, Asp, or Cys, n = 1–5 for Val and Asn; n = 1 for Asp and Cys) model peptides with α‐helix and β‐sheet backbone conformations were computed by density functional methods. The robustness results verify empirical rules drawn from experiments and from computed rotational strengths linking amide‐I and amide‐II patterns in the vibrational circular dichroism (VCD) spectra of peptides with their backbone structures. For peptides with at least three residues (n ≥ 3) these characteristic patterns from coupled amide vibrational modes have robust signatures. For shorter peptide models many vibrational modes are nonrobust, and the robust modes can be dependent on the residues or on their side chain conformations in addition to backbone conformations. These robust VCD bands, however, provide information for the detailed structural analysis of these smaller systems. Chirality 27:625–634, 2015 © 2015 Wiley Periodicals, Inc.     

Constructing arabinofuranosidases for dual arabinoxylan debranching activity

Enzymatic conversion of arabinoxylan requires α‐L‐arabinofuranosidases able to remove α‐L‐arabinofuranosyl residues (α‐L‐Araf) from both mono‐ and double‐substituted D‐xylopyranosyl residues (Xylp) in xylan (i.e., AXH‐m and AXH‐d activity). Herein, SthAbf62A (a family GH62 α‐L‐arabinofuranosidase with AXH‐m activity) and BadAbf43A (a family GH43 α‐L‐arabinofuranosidase with AXH‐d3 activity), were fused to create SthAbf62A_BadAbf43A and BadAbf43A_SthAbf62A. Both fusion enzymes displayed dual AXH‐m,d and synergistic activity toward native, highly branched wheat arabinoxylan (WAX). When using a customized arabinoxylan substrate comprising mainly α‐(1 → 3)‐L‐Araf and α‐(1 → 2)‐L‐Araf substituents attached to disubstituted Xylp (d‐2,3‐WAX), the specific activity of the fusion enzymes was twice that of enzymes added as separate proteins. Moreover, the SthAbf62A_BadAbf43A fusion removed 83% of all α‐L‐Araf from WAX after a 20 hr treatment. ¹H NMR analyses further revealed differences in SthAbf62A_BadAbf43 rate of removal of specific α‐L‐Araf substituents from WAX, where 9.4 times higher activity was observed toward d‐α‐(1 → 3)‐L‐Araf compared to m‐α‐(1 → 3)‐L‐Araf positions.     

Tacticity effects on conformational structure and hydration of poly-(methacrylic acid) in aqueous solutions-a molecular dynamics simulation study

The influence of tacticity on chain dimensions, backbone and side-group conformational states and relaxation dynamics, intermolecular hydrogen bonding and its relaxation dynamics was investigated for 30 repeat unit poly(methacrylic acid) (PMA) as a function of the degree-of-neutralisation, f (i.e. charge density) [0 < f < 1 range] in explicit water with Na⁺ neutralising counter-ions, using molecular dynamics simulations. Chain expansion with increase in charge density is observed. Simulation results, where applicable, compare well with experimental results in the literature. Isotactic (i-PMA) chain exhibits the largest expansion (89% change in R_g), and syndiotactic (s-PMA-RR-0.7) chain shows 31%. For fully neutralised chain (high charge density), the probability for trans conformation at backbone bonds follows the trend i-PMA > a-PMA > s-PMA. At high charge density, a higher probability of trans state is obtained at meso dyads as compared to racemic dyads and the side group in i-PMA relative to other tacticity is conformationally more extended. RR triads impart better hydrogen bonding with water. In the COOH side group, greater rotational mobility is observed for i-PMA as compared to s-PMA. Water coordination to PMA is invariant with tacticity at low f. For f = 1, s-PMA exhibits the best level of H-bonding, due to its 3D chemical structural configuration. Binding distance between water and PMA atoms and their coordination number remain unchanged with respect to ionisation of COOH groups. Isotactic (i-PMA) exhibits the fastest relaxation of polymer-water H-bonds. While counter-ion coordination on to PMA is unaffected by tacticity, the intensity of ion condensation as well as hydrogen bond relaxation time at f = 1 varies in the order i-PMA > a-PMA > s-PMA.