Search

Water is one of the principal constituents by mass of living plant cell walls. However, its role and interactions with secondary cell wall polysaccharides and the impact of dehydration and subsequent rehydration on the molecular architecture are still to be elucidated. This work combines multidimensional solid-state 13C magic-angle-spinning (MAS) nuclear magnetic resonance (NMR) with molecular dynamics modeling to decipher the role of water in the molecular architecture of softwood secondary cell walls. The proximities between all main polymers, their molecular conformations, and interaction energies are compared in never-dried, oven-dried, and rehydrated states. Water is shown to play a critical role at the hemicellulose–cellulose interface. After significant molecular shrinkage caused by dehydration, the original molecular conformation is not fully recovered after rehydration. The changes include xylan becoming more closely and irreversibly associated with cellulose and some mannan becoming more mobile and changing conformation. These irreversible nanostructural changes provide a basis for explaining and improving the properties of wood-based materials.

Bones are responsible for body support, structure, motion, and several other functions that enable and facilitate life for many different animal species. They exhibit a complex network of distinct physical structures and mechanical properties, which ultimately depend on the fraction of their primary constituents at the molecular scale. However, the relationship between structure and mechanical properties in bones are still not fully understood. Here, we investigate structural and mechanical properties of all-atom bone molecular models composed of type-I collagen, hydroxyapatite (HA), and water by means of fully atomistic molecular dynamics simulations. Our models encompass an extrafibrillar volume (EFV) and consider mineral content in both the EFV and intrafibrillar volume (IFV), consistent with experimental observations. We investigate solvation structures and elastic properties of bone microfibril models with different degrees of mineralization, ranging from highly mineralized to weakly mineralized and nonmineralized models. We find that the local tetrahedral order of water is lost in similar ways in the EFV and IFV regions for all HA containing models, as calcium and phosphate ions are strongly coordinated with water molecules. We also subject our models to tensile loads and analyze the spatial stress distribution over the nanostructure of the material. Our results show that both mineral and water contents accumulate significantly higher stress levels, most notably in the EFV, thus revealing that this region, which has been only recently incorporated in all-atom molecular models, is fundamental for studying the mechanical properties of bones at the nanoscale. Furthermore, our results corroborate the well-established finding that high mineral content makes bone stiffer.

Glycoside hydrolases (GH) cleave carbohydrate glycosidic bonds and play pivotal roles in living organisms and in many industrial processes. Unlike acid-catalyzed hydrolysis of carbohydrates in solution, which can occur either via cyclic or acyclic oxocarbenium-like transition states, it is widely accepted that GH-catalyzed hydrolysis proceeds via a general acid mechanism involving a cyclic oxocarbenium-like transition state with protonation of the glycosidic oxygen. The GH45 subfamily C inverting endoglucanase from Phanerochaete chrysosporium (PcCel45A) defies the classical inverting mechanism as its crystal structure conspicuously lacks a general Asp or Glu base residue. Instead, PcCel45A has an Asn residue, a notoriously weak base in solution, as one of its catalytic residues at position 92. Moreover, unlike other inverting GHs, the relative position of the catalytic residues in PcCel45A impairs the proton abstraction from the nucleophilic water that attacks the anomeric carbon, a key step in the classical mechanism. Here, we investigate the viability of an endocyclic mechanism for PcCel45A using hybrid quantum mechanics/molecular mechanics (QM/MM) simulations, with the QM region treated with the self-consistent-charge density-functional tight-binding level of theory. In this mechanism, an acyclic oxocarbenium-like transition state is stabilized leading to the opening of the glucopyranose ring and formation of an unstable acyclic hemiacetal that can be readily decomposed into hydrolysis product. In silico characterization of the Michaelis complex shows that PcCel45A significantly restrains the sugar ring to the 4C1 chair conformation at the −1 subsite of the substrate binding cleft, in contrast to the classical exocyclic mechanism in which ring puckering is critical. We also show that PcCel45A provides an environment where the catalytic Asn92 residue in its standard amide form participates in a cooperative hydrogen bond network resulting in its increased nucleophilicity due to an increased negative charge on the oxygen atom. Our results for PcCel45A suggest that carbohydrate hydrolysis catalyzed by GHs may take an alternative route from the classical mechanism.

Serine hydrolases cleave peptide and ester bonds and are ubiquitous in nature, with applications in biotechnology, in materials, and as drug targets. The serine hydrolase two-step mechanism employs a serine-histidine-aspartate/glutamate catalytic triad, where the histidine residue acts as a base to activate poor nucleophiles (a serine residue or a water molecule) and as an acid to allow the dissociation of poor leaving groups. This mechanism has been the subject of debate regarding how histidine shuttles the proton from the nucleophile to the leaving group. To elucidate the reaction mechanism of serine hydrolases, we employ quantum mechanics/molecular mechanics-based transition path sampling to obtain the reaction coordinate using the Aspergillus niger feruloyl esterase A (AnFaeA) as a model enzyme. The optimal reaction coordinates include terms involving nucleophilic attack on the carbonyl carbon and proton transfer to, and dissociation of, the leaving group. During the reaction, the histidine residue undergoes a reorientation on the time scale of hundreds of femtoseconds that supports the "moving histidine" mechanism, thus calling into question the "ring flip" mechanism. We find a concerted mechanism, where the transition state coincides with the tetrahedral intermediate with the histidine residue pointed between the nucleophile and the leaving group. Moreover, motions of the catalytic aspartate toward the histidine occur concertedly with proton abstraction by the catalytic histidine and help stabilize the transition state, thus partially explaining how serine hydrolases enable poor nucleophiles to attack the substrate carbonyl carbon. Rate calculations indicate that the second step (deacylation) is rate-determining, with a calculated rate constant of 66 s-1. Overall, these results reveal the pivotal role of active-site dynamics in the catalytic mechanism of AnFaeA, which is likely similar in other serine hydrolases.

SARS-CoV-2, the pathogenic agent of COVID-19, employs angiotensin converting enzyme-2 (ACE2) as its cell entry receptor. Clinical data reveal that in severe COVID-19, SARS-CoV-2 infects the lung, leading to a frequently lethal triad of respiratory insufficiency, acute cardiovascular failure, and coagulopathy. Physiologically, ACE2 plays a role in the regulation of three systems that could potentially be involved in the pathogenesis of severe COVID-19: the kinin–kallikrein system, resulting in acute lung inflammatory edema; the renin–angiotensin system, promoting cardiovascular instability; and the coagulation system, leading to thromboembolism. Here we assembled a healthy human lung cell atlas meta-analysis with ~ 130,000 public single-cell transcriptomes and show that key elements of the bradykinin, angiotensin and coagulation systems are co-expressed with ACE2 in alveolar cells and associated with their differentiation dynamics, which could explain how changes in ACE2 promoted by SARS-CoV-2 cell entry result in the development of the three most severe clinical components of COVID-19.

Aquaporins are membrane proteins responsible for permeating water, ions, dissolved gases, and other small molecular weight compounds through the protective cell membranes of living organisms. These proteins have been gaining increased importance as targets for treating a variety of parasitic diseases, since they control key physiological processes in the life cycle of parasitic protozoans, such as the uptake of nutrients, release of metabolites, and alleviation of osmotic stress. In this work, we use homology modeling to build three-dimensional structures for the four main aquaporins encoded and expressed by Leishmania major, a protozoan that causes leishmaniasis and affects millions of people worldwide. Physico-chemical properties of the proposed models for LmAQP1, LmAQPα, LmAQPβ, and LmAQPγ are then investigated using molecular dynamics simulations and the reference interaction site model (RISM) molecular theory of solvation. Pore characteristics, water permeation, and potential of mean force across the AQP channels for water, methanol, urea, ammonia, and carbon dioxide are examined and compared with results obtained for a protozoan (Plasmodium falciparum) aquaporin for which a crystal structure is available.

The fundamental and assorted roles of β-1,3-glucans in nature are underpinned on diverse chemistry and molecular structures, demanding sophisticated and intricate enzymatic systems for their processing. In this work, the selectivity and modes of action of a glycoside hydrolase family active on β-1,3-glucans were systematically investigated combining sequence similarity network, phylogeny, X-ray crystallography, enzyme kinetics, mutagenesis and molecular dynamics. This family exhibits a minimalist and versatile (α/β)-barrel scaffold, which can harbor distinguishing exo or endo modes of action, including an ancillary-binding site for the anchoring of triple-helical β-1,3-glucans. The substrate binding occurs via a hydrophobic knuckle complementary to the canonical curved conformation of β-1,3-glucans or through a substrate conformational change imposed by the active-site topology of some fungal enzymes. Together, these findings expand our understanding of the enzymatic arsenal of bacteria and fungi for the breakdown and modification of β-1,3-glucans, which can be exploited for biotechnological applications.

Metal–organic frameworks (MOFs) comprise a class of highly porous nanomaterials formed by the assembly of organic molecular templates connected by metal ions. These materials exhibit a large diversity of pore size and geometry, topology, surface area, and chemical functionality. MOFs are particularly promising materials for developing new technologies for capture and storage of greenhouse gases such as methane and carbon dioxide. Here we apply the three dimensional reference interaction site model (3D-RISM) molecular theory of solvation to study the interactions of CO2 and CH4 with the metal–organic material MIL-47. The 3D-RISM integral equations were solved to determine the three dimensional density correlation functions of the gas (solvent) relative to the atomic positions of the MIL-47 framework, treated as static solute sites. The computed solvent spatial distributions inside MIL-47 pores were used to identify whether or not there exist preferable binding sites and the binding free energy landscape for the gas of interest at low computational costs compared with other molecular modeling techniques, such as grand canonical Monte Carlo and molecular dynamics simulations. The 3D-RISM formalism was applied to pure CO2, pure CH4, and binary mixtures of these gases of various compositions under different pressure conditions. The results indicate that both gases bind very weakly to MIL-47 and that this material exhibits nearly vanishing CO2/CH4 selectivity. The 3D-RISM computations presented here can be extended to investigate the physical adsorption of gases on other MOFs and nanoporous materials, providing an alternative low-cost computational approach to study gas capture and storage in nanoporous materials in general and, in particular, to determine the binding free-energy landscape in these systems.

Cellulose possesses considerable potential for a wide range of sustainable applications. Nanocellulose-based material properties are primarily dependent on the structural surface characteristics of its crystalline planes. Experimental measurements of the affinity of crystalline nanocellulose surfaces with water are scarce and challenging to obtain. Therefore, the relative hydrophilicity of different cellulose allomorphs crystalline planes is often inferred from qualitative assessments of their surface and the exposition of polar groups to the solvent. This work investigates the relative hydrophilicity of cellulose surfaces using molecular dynamics simulations. The behavior of a water droplet laid on different crystal planes was used to determine their relative hydrophilicity. The water molecules fully spread onto highly hydrophilic surfaces. However, a water droplet placed on less hydrophilic surfaces equilibrates as an oblate spheroidal cap allowing the measurement of a contact angle. The results indicate that the Iα (010), Iα (11̅0), Iβ (010), and Iβ (110) faces, as well as the faces of human-made celluloses II and III_I (100), (11̅0), (010), and (110) are all highly hydrophilic. They all have a contact angle value inferior to 11°. Not unexpectedly, the Iα (001) and Iβ (100) surfaces are less hydrophilic with contact angles of 48 and 34°, respectively. However, the Iβ (11̅0) plane, often referred to as a hydrophilic surface, forms a contact angle of about 32°. The results are rationalized in terms of structure, exposure of hydroxyl groups to the solvent, and degree of cellulose-cellulose versus cellulose-water hydrogen bonds on each face. The simulations also show that the surface oxidation degree tunes the surface hydrophilicity in a nonlinear manner due to cooperative effects involving water-cellulose interactions. Our study helps us to understand how the degree of hydrophilicity of cellulose emerges from specific structural features of each crystalline surface.

Glucuronoyl esterases (GEs) are α/β serine hydrolases and a relatively new addition in the toolbox to reduce the recalcitrance of lignocellulose, the biggest obstacle in cost-effective utilization of this important renewable resource. While biochemical and structural characterization of GEs have progressed greatly recently, there have yet been no mechanistic studies shedding light onto the rate-limiting steps relevant for biomass conversion. The bacterial GE OtCE15A possesses a classical yet distinctive catalytic machinery, with easily identifiable catalytic Ser/His completed by two acidic residues (Glu and Asp) rather than one as in the classical triad, and an Arg side chain participating in the oxyanion hole. By QM/MM calculations, we identified deacylation as the decisive step in catalysis, and quantified the role of Asp, Glu and Arg, showing the latter to be particularly important. The results agree well with experimental and structural data. We further calculated the free-energy barrier of post-catalysis dissociation from a complex natural substrate, suggesting that in industrial settings non-catalytic processes may constitute the rate-limiting step, and pointing to future directions for enzyme engineering in biomass utilization.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the agent of a major global outbreak of respiratory tract disease known as coronavirus disease-2019 (COVID-19). SARS-CoV-2 infects mainly lungs and may cause several immune-related complications, such as lymphocytopenia and cytokine storm, which are associated with the severity of the disease and predict mortality1,2. The mechanism by which SARS-CoV-2 infection may result in immune system dysfunction is still not fully understood. Here we show that SARS-CoV-2 infects human CD4+T helper cells, but not CD8+T cells, and is present in blood and bronchoalveolar lavage T helper cells of severe COVID-19 patients. We demonstrated that SARS-CoV-2 spike glycoprotein (S) directly binds to the CD4 molecule, which in turn mediates the entry of SARS-CoV-2 in T helper cells. This leads to impaired CD4 T cell function and may cause cell death. SARS-CoV-2-infected T helper cells express higher levels of IL-10, which is associated with viral persistence and disease severity. Thus, CD4-mediated SARS-CoV-2 infection of T helper cells may contribute to a poor immune response in COVID-19 patients.

Metal-organic frameworks (MOFs) comprise a class of highly porous nanomaterials formed by the assembly of organic molecular templates connected by metal ions. These materials exhibit a large diversity of pore size and geometry, topology, surface area, and chemical functionality. MOFs are particularly promising materials for developing new technologies for capture and storage of greenhouse gases such as methane and carbon dioxide. Here we apply the three dimensional reference interaction site model (3D-RISM) molecular theory of solvation to study the interactions of CO

SARS-CoV-2, the pathogenic agent of COVID-19, employs angiotensin converting enzyme-2 (ACE2) as its cell entry receptor. Clinical data reveal that in severe COVID-19, SARS-CoV-2 infects the lung, leading to a frequently lethal triad of respiratory insufficiency, acute cardiovascular failure, and coagulopathy. Physiologically, ACE2 plays a role in the regulation of three systems that could potentially be involved in the pathogenesis of severe COVID-19: the kinin-kallikrein system, resulting in acute lung inflammatory edema; the renin-angiotensin system, promoting cardiovascular instability; and the coagulation system, leading to thromboembolism. Here we analyzed ~130,000 human lung single-cell transcriptomes and show that key elements of the kinin-kallikrein, renin-angiotensin and coagulation systems are co-expressed with ACE2 in alveolar cells, which could explain how changes in ACE2 promoted by SARS- CoV-2 cell entry result in the development of the three most severe clinical components of COVID-19.

The linker of the endoglucanase from Xanthomonas campestris pv. campestris ((PT)12) has a specific sequence, a repeating proline-threonine motif. In order to understand its role, it has been compared to a regular sequence linker, in this work-the cellobiohydrolase 2 from Trichoderma reesei (CBH2). Elastic properties of the two linkers have been estimated by calculating free energy profile along the linker length from an enhanced sampling molecular dynamics simulation. The (PT)12 exhibits more pronounced elastic behaviour than CBH2. The PT repeating motif results in a two-mode energy profile which could be very useful in the enzyme motions along the substrate during hydrolytic catalysis.

The interaction between mannan polysaccharides and cellulose microfibrils contributes to cell wall properties in some vascular plants, but the molecular arrangement of mannan in the cell wall and the nature of the molecular bonding between mannan and cellulose remain unknown. Previous studies have shown that mannan is important in maintaining Arabidopsis (Arabidopsis thaliana) seed mucilage architecture, and that Cellulose Synthase-Like A2 (CSLA2) synthesizes a glucomannan backbone, which Mannan α-Galactosyl Transferase1 (MAGT1/GlycosylTransferase-Like6/Mucilage Related10) might decorate with single α-Gal branches. Here, we investigated the ratio and sequence of Man and Glc and the arrangement of Gal residues in Arabidopsis mucilage mannan using enzyme sequential digestion, carbohydrate gel electrophoresis, and mass spectrometry. We found that seed mucilage galactoglucomannan has a backbone consisting of the repeating disaccharide [4)-β-Glc-(1,4)-β-Man-(1,], and most of the Man residues in the backbone are substituted by single α-1,6-Gal. CSLA2 is responsible for the synthesis of this patterned glucomannan backbone and MAGT1 catalyses the addition of α-Gal. In vitro activity assays revealed that MAGT1 transferred α-Gal from UDP-Gal only to Man residues within the CSLA2 patterned glucomannan backbone acceptor. These results indicate that CSLAs and galactosyltransferases are able to make precisely defined galactoglucomannan structures. Molecular dynamics simulations suggested this patterned galactoglucomannan is able to bind stably to some hydrophilic faces and to hydrophobic faces of cellulose microfibrils. A specialization of the biosynthetic machinery to make galactoglucomannan with a patterned structure may therefore regulate the mode of binding of this hemicellulose to cellulose fibrils.

Since the commercial introduction of Ion Mobility coupled with Mass Spectrometry (IM-MS) devices in 2003, a large number of research laboratories have embraced the technique. IM-MS is a fairly rapid experiment used as a molecular separation tool and to obtain structural information. The interpretation of IM-MS data is still challenging and relies heavily on theoretical calculations of the molecule's collision cross section (CCS) against a buffer gas. Here, a new software (HPCCS) is presented, which performs CCS calculations using high perfomance computing techniques. Based on the trajectory method, HPCCS can accurately calculate CCS for a great variety of molecules, ranging from small organic molecules to large protein complexes, using helium or nitrogen as buffer gas with considerable gains in computer time compared to publicly available codes under the same level of theory. HPCCS is available as free software under the Academic Use License at https://github.com/cepid-cces/hpccs. © 2018 Wiley Periodicals, Inc.

Human hemoglobin (Hb) Coimbra (βAsp99Glu) is one of the seven βAsp99 Hb variants described to date. All βAsp99 substitutions result in increased affinity for O2 and decreased heme-heme cooperativity and their carriers are clinically characterized by erythrocytocis, caused by tissue hypoxia. Since βAsp99 plays an important role in the allosteric α1β2 interface and the mutation in Hb Coimbra only represents the insertion of a CH2 group in this interface, the present study of Hb Coimbra is important for a better understanding of the global impact of small modifications in this allosteric interface. We carried out functional, kinetic and dynamic characterization of this hemoglobin, focusing on the interpretation of these results in the context of a growth of the position 99 side chain length in the α1β2 interface. Oxygen affinity was evaluated by measuring p50 values in distinct pHs (Bohr effect), and the heme-heme cooperativity was analyzed by determining the Hill coefficient (n), in addition to the effect of the allosteric effectors inositol hexaphosphate (IHP) and 2,3-bisphosphoglyceric acid (2,3-BPG). Computer simulations revealed a stabilization of the R state in the Coimbra variant with respect to the wild type, and consistently, the T-to-R quaternary transition was observed on the nanosecond time scale of classical molecular dynamics simulations. Fil: Jorge, S. E.. Universidade Estadual de Campinas; Brasil Fil: Bringas, Mauro. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Química, Física de los Materiales, Medioambiente y Energía. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Química, Física de los Materiales, Medioambiente y Energía; Argentina Fil: Petruk, Ariel Alcides. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Química, Física de los Materiales, Medioambiente y Energía. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Química, Física de los Materiales, Medioambiente y Energía; Argentina Fil: Arrar, Mehrnoosh. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Química, Física de los Materiales, Medioambiente y Energía. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Química, Física de los Materiales, Medioambiente y Energía; Argentina Fil: Marti, Marcelo Adrian. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales; Argentina Fil: Skaf, M. S.. Universidade Estadual de Campinas; Brasil Fil: Costa, F. F.. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Química, Física de los Materiales, Medioambiente y Energía. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Química, Física de los Materiales, Medioambiente y Energía; Argentina Fil: Capece, Luciana. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Química, Física de los Materiales, Medioambiente y Energía. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Química, Física de los Materiales, Medioambiente y Energía; Argentina Fil: Sonati, M. F.. Universidade Estadual de Campinas; Brasil Fil: Estrin, Dario Ariel. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Química, Física de los Materiales, Medioambiente y Energía. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Química, Física de los Materiales, Medioambiente y Energía; Argentina

Cellobiohydrolases (CBHs) are key enzymes for the saccharification of cellulose and play major roles in industrial settings for biofuel production. The catalytic core domain of these enzymes exhibits a long and narrow binding tunnel capable of binding glucan chains from crystalline cellulose and processively hydrolyze them. The binding cleft is topped by a set of loops, which are believed to play key roles in substrate binding and cleavage processivity. Here, we present an analysis of the loop motions of the Trichoderma reesei Cel7A catalytic core domain (TrCel7A) using conventional and accelerated molecular dynamics simulations. We observe that the loops exhibit highly coupled fluctuations and cannot move independently of each other. In the absence of a substrate, the characteristic large amplitude dynamics of TrCel7A consists of breathing motions, where the loops undergo open-and-close fluctuations. Upon substrate binding, the open-close fluctuations of the loops are quenched and one of the loops moves parallel to the binding site, possibly to allow processive motion along the glucan chain. Using microsecond accelerated molecular dynamics, we observe large-scale fluctuations of the loops (up to 37 Å) and the entire exposure of the TrCel7A binding site in the absence of the substrate, resembling an endoglucanase. These results suggest that the initial CBH-substrate contact and substrate recognition by the enzyme are similar to that of endoglucanases and, once bound to the substrate, the loops remain closed for proper enzymatic activity.

Protein glycosylation is a diverse post-translational modification that serves myriad biological functions. O-linked glycans in particular vary widely in extent and chemistry in eukaryotes, with secreted proteins from fungi and yeast commonly exhibiting O-mannosylation in intrinsically disordered regions of proteins, likely for proteolysis protection, among other functions. However, it is not well understood why mannose is often the preferred glycan, and more generally, if the neighboring protein sequence and glycan have coevolved to protect against proteolysis in glycosylated intrinsically disordered proteins (IDPs). Here, we synthesized variants of a model IDP, specifically a natively O-mannosylated linker from a fungal enzyme, with α-O-linked mannose, glucose, and galactose moieties, along with a non-glycosylated linker. Upon exposure to thermolysin, O-mannosylation, by far, provides the highest extent of proteolysis protection. To explain this observation, extensive molecular dynamics simulations were conducted, revealing that the axial configuration of the C2-hydroxyl group (2-OH) of α-mannose adjacent to the glycan–peptide bond strongly influences the conformational features of the linker. Specifically, α-mannose restricts the torsions of the IDP main chain more than other glycans whose equatorial 2-OH groups exhibit interactions that favor perpendicular glycan–protein backbone orientation. We suggest that IDP stiffening due to O-mannosylation impairs protease action, with contributions from protein–glycan interactions, protein flexibility, and protein stability. Our results further imply that resistance to proteolysis is an important driving force for evolutionary selection of α-mannose in eukaryotic IDPs, and more broadly, that glycan motifs for proteolysis protection likely coevolve with the protein sequence to which they attach.

Processive cellulases are highly efficient molecular engines involved in the cellulose breakdown process. However, the mechanism that processive bacterial enzymes utilize to recruit and retain cellulose strands in the catalytic site remains poorly understood. Here, integrated enzymatic assays, protein crystallography and computational approaches were combined to study the enzymatic properties of the processive BlCel48B cellulase from Bacillus licheniformis. Hydrolytic efficiency, substrate binding affinity, cleavage patterns, and the apparent processivity of bacterial BlCel48B are significantly impacted by the cellulose size and its surface morphology. BlCel48B crystallographic structure was solved with ligands spanning -5 to -2 and +1 to +2 subsites. Statistical coupling analysis and molecular dynamics show that co-evolved residues on active site are critical for stabilizing ligands in the catalytic tunnel. Our results provide mechanistic insights into BlCel48B molecular-level determinants of activity, substrate binding, and processivity on insoluble cellulose, thus shedding light on structure-activity correlations of GH48 family members in general.

A technical overview of the High Performance Collision Cross Section (HPCCS) software for accurate and efficient calculations of collision cross sections for molecular ions ranging from small organic molecules to large protein complexes is presented. The program uses helium or nitrogen as buffer gas with considerable gains in computer time compared to publicly available codes under the Trajectory Method approximation. HPCCS is freely available under the Academic Use License at https://github.com/cepid-cces/hpccs .

Given the abundance of lignin in nature, multiple enzyme systems have been discovered to cleave the β-O-4 bonds, the most prevalent intermonomer linkage. In particular, stereospecific cleavage of lignin oligomers by glutathione S-transferases (GSTs) has been reported in several sphingomonads. Here, we apply quantum mechanics/molecular mechanics simulations to study the mechanism of two glutathione-dependent enzymes in the β-aryl ether catabolic pathway of Sphingomonas sp. SYK-6, namely, LigF, a β-etherase, and LigG, a lyase. For LigF, the free-energy landscape supports a SN2 reaction mechanism, with the monoaromatic leaving group being promptly neutralized upon release. Specific interactions with conserved residues are responsible for stereoselectivity and for activation of the cofactor as a nucleophile. A glutathione conjugate is also released by LigF and serves the substrate of LigG, undergoing a SN2-like reaction, in which Cys15 acts as the nucleophile, to yield the second monoaromatic product. The simulations suggest that the electron-donating substituent at the para-position found in lignin-derived aromatics and the interaction with Tyr217 are essential for reactivity in LigG. Overall, this work deepens the understanding of the stereospecific enzymatic mechanisms in the β-aryl ether cleavage pathway and reveals key structural features underpinning the ligninolytic activity detected in several sphingomonad GSTs.

The fundamental and assorted roles of β-1,3-glucans in nature are underpinned on diverse chemistry and molecular structures, demanding sophisticated and intricate enzymatic systems for their processing. In this work, the selectivity and modes of action of a glycoside hydrolase family active on β-1,3-glucans were systematically investigated combining sequence similarity network, phylogeny, X-ray crystallography, enzyme kinetics, mutagenesis and molecular dynamics. This family exhibits a minimalist and versatile (α/β)-barrel scaffold, which can harbor distinguishing exo or endo modes of action, including an ancillary-binding site for the anchoring of triple-helical β-1,3-glucans. The substrate binding occurs via a hydrophobic knuckle complementary to the canonical curved conformation of β-1,3-glucans or through a substrate conformational change imposed by the active-site topology of some fungal enzymes. Together, these findings expand our understanding of the enzymatic arsenal of bacteria and fungi for the breakdown and modification of β-1,3-glucans, which can be exploited for biotechnological applications.

Made available in DSpace on 2020-12-12T02:23:44Z (GMT). No. of bitstreams: 0 Previous issue date: 2019-11-14 Background: β-Xylosidases are glycoside hydrolases (GHs) that cleave xylooligosaccharides and/or xylobiose into shorter oligosaccharides and xylose. Aspergillus nidulans is an established genetic model and good source of carbohydrate-active enzymes (CAZymes). Most fungal enzymes are N-glycosylated, which influences their secretion, stability, activity, signalization, and protease protection. A greater understanding of the N-glycosylation process would contribute to better address the current bottlenecks in obtaining high secretion yields of fungal proteins for industrial applications. Results: In this study, BxlB-a highly secreted GH3 β-xylosidase from A. nidulans, presenting high activity and several N-glycosylation sites-was selected for N-glycosylation engineering. Several glycomutants were designed to investigate the influence of N-glycans on BxlB secretion and function. The non-glycosylated mutant (BxlBnon-glyc) showed similar levels of enzyme secretion and activity compared to the wild-type (BxlBwt), while a partially glycosylated mutant (BxlBN1;5;7) exhibited increased activity. Additionally, there was no enzyme secretion in the mutant in which the N-glycosylation context was changed by the introduction of four new N-glycosylation sites (BxlBCC), despite the high transcript levels. BxlBwt, BxlBnon-glyc, and BxlBN1;5;7 formed similar secondary structures, though the mutants had lower melting temperatures compared to the wild type. Six additional glycomutants were designed based on BxlBN1;5;7, to better understand its increased activity. Among them, the two glycomutants which maintained only two N-glycosylation sites each (BxlBN1;5 and BxlBN5;7) showed improved catalytic efficiency, whereas the other four mutants' catalytic efficiencies were reduced. The N-glycosylation site N5 is important for improved BxlB catalytic efficiency, but needs to be complemented by N1 and/or N7. Molecular dynamics simulations of BxlBnon-glyc and BxlBN1;5 reveals that the mobility pattern of structural elements in the vicinity of the catalytic pocket changes upon N1 and N5 N-glycosylation sites, enhancing substrate binding properties which may underlie the observed differences in catalytic efficiency between BxlBnon-glyc and BxlBN1;5. Conclusions: This study demonstrates the influence of N-glycosylation on A. nidulans BxlB production and function, reinforcing that protein glycoengineering is a promising tool for enhancing thermal stability, secretion, and enzymatic activity. Our report may also support biotechnological applications for N-glycosylation modification of other CAZymes. Department of Biochemistry and Tissue Biology Institute of Biology University of Campinas (UNICAMP), Rua Monteiro Lobato, 255 Cidade Universitária Zeferino Vaz Department of Physics Institute of Biosciences Humanities and Exact Sciences São Paulo State University (UNESP) Department of Medical Science Faculty of Medicine University of Campinas (UNICAMP) Department of Biochemistry and Immunology Ribeirão Preto Medical School University of São Paulo (USP) Institute of Chemistry and Center for Computing in Engineering and Sciences University of Campinas (UNICAMP) Department of Physics Institute of Biosciences Humanities and Exact Sciences São Paulo State University (UNESP)

Materials for vehicular methane storage have been extensively studied, although no suitable material has been found. In this work, we use molecular simulation to investigate three types of carbon-b...

Molecular dynamic simulations were used to investigate the wetting behavior of nanocrystalline cellulose allomorphs. Three 40ns independent MD simulations using CHARMM force field gave droplet profiles, whose contact angles were obtained with LBADSA plugin for ImageJ. Surface hydrophilicity was related to hydroxyl group availability for hydrogen bonding. Oxidation of C6 hydroxyl groups significantly boosts surface hydrophilicity, and even low levels of modification (~7%) leads to contact angles close to zero degrees. These data can be useful for biotechnological applications.

Cellulases are essential enzymatic components for the transformation of plant biomass into fuels, renewable ma- terials and green chemicals. Here, we determined the crystal structure, pattern of hydrolysis products release, and conducted molecular dynamics simulations of the major endoglucanase from the Xanthomonas campestris pv. campestris (XccCel5A). XccCel5A has a TIM barrel fold with the catalytic site centrally placed in a binding groove surrounded by aromatic side chains. Molecular dynamics simulations show that productive position of the substrate is secured by a network of hydrogen bonds in the four main subsites, which differ in details from homologous structures. Capillary zone electrophoresis and computational studies reveal XccCel5A can act both as endoglucanase and licheninase, but there are preferable arrangements of substrate regarding β-1,3 and β- 1,4 bonds within the binding cleft which are related to the enzymatic efficiency.

Glycoside hydrolases (GHs) are involved in the degradation of a wide diversity of carbohydrates and present several biotechnological applications. Many GH families are composed of enzymes with a single well-defined specificity. In contrast, enzymes from the GH16 family can act on a range of different polysaccharides, including β-glucans and galactans. SCLam, a GH16 member derived from a soil metagenome, an endo-β-1,3(4)-glucanase (EC 3.2.1.6), can cleave both β-1,3 and β-1,4 glycosidic bonds in glucans, such as laminarin, barley β-glucan, and cello-oligosaccharides. A similar cleavage pattern was previously reported for other GH16 family members. However, the molecular mechanisms for this dual cleavage activity on (1,3)- and (1,4)-β-D-glycosidic bonds by laminarinases have not been elucidated. In this sense, we determined the X-ray structure of a presumably inactive form of SCLam cocrystallized with different oligosaccharides. The solved structures revealed general bound products that are formed owing to residual activities of hydrolysis and transglycosylation. Biochemical and biophysical analyses and molecular dynamics simulations help to rationalize differences in activity toward different substrates. Our results depicted a bulky aromatic residue near the catalytic site critical to select the preferable configuration of glycosidic bonds in the binding cleft. Altogether, these data contribute to understanding the structural basis of recognition and hydrolysis of β-1,3 and β-1,4 glycosidic linkages of the laminarinase enzyme class, which is valuable for future studies on the GH16 family members and applications related to biomass conversion into feedstocks and bioproducts.

Cellulose, the most abundant biopolymer on Earth, represents a resource for sustainable production of biofuels. Thermochemical treatments make lignocellulosic biomaterials more amenable to depolymerization by exposing cellulose microfibrils to enzymatic or chemical attacks. In such treatments, the solvent plays fundamental roles in biomass modification, but the molecular events underlying these changes are still poorly understood. Here, the 3D-RISM-KH molecular theory of solvation has been employed to analyze the role of water in cellulose aggregation under different thermodynamic conditions. The results show that, under ambient conditions, highly structured hydration shells around cellulose create repulsive forces that protect cellulose microfibrils from aggregating. Under hydrothermal pretreatment conditions, however, the hydration shells lose structure, and cellulose aggregation is favored. These effects are largely due to a decrease in cellulose-water interactions relative to those at ambient conditions, so that cellulose-cellulose attractive interactions become prevalent. Our results provide an explanation to the observed increase in the lateral size of cellulose crystallites when biomass is subject to pretreatments and deepen the current understanding of the mechanisms of biomass modification.

Fungal GH12 enzymes are classified as xyloglucanases when they specifically target xyloglucans, or promiscuous endoglucanases when they exhibit catalytic activity against xyloglucan and β-glucan chains. Several structural and functional studies involving GH12 enzymes tried to explain the main patterns of xyloglucan activity, but what really determines xyloglucanase specificity remains elusive. Here, three fungal GH12 enzymes from Aspergillus clavatus (AclaXegA), A. zonatus (AspzoGH12), and A. terreus (AtEglD) were studied to unveil the molecular basis for substrate specificity. Using functional assays, site-directed mutagenesis, and molecular dynamics simulations, we demonstrated that three main regions are responsible for substrate selectivity: (i) the YSG group in loop 1; (ii) the SST group in loop 2; and (iii) loop A3-B3 and neighboring residues. Functional assays and sequence alignment showed that while AclaXegA is specific to xyloglucan, AtEglD cleaves β-glucan, and xyloglucan. However, AspzoGH12 was also shown to be promiscuous contrarily to a sequence alignment-based prediction. We find that residues Y111 and R93 in AtEglD harbor the substrate in an adequate orientation for hydrolysis in the catalytic cleft entrance and that residues Y19 in AclaXegA and Y30 in AspzoGH12 partially compensate the absence of the YSG segment, typically found in promiscuous enzymes. The results point out the multiple structural factors underlying the substrate specificity of GH12 enzymes. Biotechnol. Bioeng. 2016;113: 2577-2586. © 2016 Wiley Periodicals, Inc.

Carbon nanoscrolls are similar to multi-walled carbon nanotubes but constructed from rolled graphene sheets into papyrus-like structures. In this work, molecular simulations are used to evaluate the post-combustion CO2 capture properties of nanoscrolls made of graphene, α-, β-, and γ-graphyne, boron nitride, and three types of carbon nitride. The CO2 uptake capacity, CO2/N2 selectivity and CO2 working capacity were computed with grand canonical Monte Carlo simulations at conditions relevant to post-combustion CO2 capture. The interlayer spacing of the nanoscrolls was optimized for each property and sheet material. For graphene nanoscrolls, the optimal interlayer spacing of 7.3 A was identified for both the CO2 uptake and selectivity, while for working capacity the optimal interlayer spacing was determined to be 8.6 A. It was found that the CO2 uptake capacity of the materials correlated to the density of the sheets from which they were formed. Nanoscrolls made from graphene and boron nitride, which have the highest number of atoms per unit area, also showed the highest CO2 uptakes. At 0.15 bar CO2, 313 K, graphene and boron nitride nanoscrolls exhibited exceptional CO2 uptake capacities of 7.7 and 8.2 mmol/g, respectively, while also exhibiting high CO2/N2 selectivities of 135 and 153, respectively. Molecular dynamics simulations were used to examine the adsorption kinetics. The simulations showed that an empty graphene nanoscroll with a roll length of 200 A could adsorb CO2 into the center of the roll within 10 ns. Materials with pores that can allow CO2 to pass through, such as graphynes, showed much faster adsorption times.

Expansins are disruptive proteins that loosen growing plant cell walls and can enhance the enzymatic hydrolysis of cellulose. The canonical expansin structure consists of one domain responsible for substrate binding (D2) and another domain (D1) of unknown function, but essential for activity. Although the effects of expansins on cell walls and cellulose fibrils are known, the molecular mechanism underlying their biophysical function is poorly understood. Here, we use molecular dynamics simulations to gain insights into the mechanism of action of the Bacillus subtilis expansin BsEXLX1. We show that BsEXLX1 can slide on the hydrophobic surface of crystalline cellulose via the flat aromatic surface of its binding domain D2, comprised mainly of residues Trp125 and Trp126. Also, we observe that BsEXLX1 can hydrogen bond a free glucan chain in a twisted conformation and that the twisting is chiefly induced by means of residue Asp82 located on D1, which has been shown to be essential for expansin activity. These results suggest that BsEXLX1 could move on the surface of cellulose and disrupt hydrogen bonds by twisting glucan chains. Simulations of the inactive BsEXLX1 mutants Asp82Asn and Tyr73Ala indicate structural alterations around the twisting center in the domain D1, which suggest a molecular basis for the lack of activity of these mutants and corroborate the idea that BsEXLX1 works by inducing twists on glucan chains. Moreover, simulations of the double mutant Asp82Asn/Tyr73Leu predict the recovery of the lost activity of BsEXLX1-Asp82Asn. Our results provide a dynamical view of the expansin-substrate interactions at the molecular scale and help shed light on the expansin mechanism.

Microbial aromatic catabolism offers a promising approach to convert lignin, a vast source of renewable carbon, into useful products. Aryl-O-demethylation is an essential biochemical reaction to ultimately catabolize coniferyl and sinapyl lignin-derived aromatic compounds, and is often a key bottleneck for both native and engineered bioconversion pathways. Here, we report the comprehensive characterization of a promiscuous P450 aryl-O-demethylase, consisting of a cytochrome P450 protein from the family CYP255A (GcoA) and a three-domain reductase (GcoB) that together represent a new two-component P450 class. Though originally described as converting guaiacol to catechol, we show that this system efficiently demethylates both guaiacol and an unexpectedly wide variety of lignin-relevant monomers. Structural, biochemical, and computational studies of this novel two-component system elucidate the mechanism of its broad substrate specificity, presenting it as a new tool for a critical step in biological lignin conversion., Catabolizing lignin-derived aromatic compounds requires an aryl-O-demethylation step. Here the authors present the structures of GcoA and GcoB, a cytochrome P450-reductase pair that catalyzes aryl-O-demethylations and show that GcoA displays broad substrate specificity, which is of interest for biotechnology applications.

Significance Synthetic polymers are ubiquitous in the modern world but pose a global environmental problem. While plastics such as poly(ethylene terephthalate) (PET) are highly versatile, their resistance to natural degradation presents a serious, growing risk to fauna and flora, particularly in marine environments. Here, we have characterized the 3D structure of a newly discovered enzyme that can digest highly crystalline PET, the primary material used in the manufacture of single-use plastic beverage bottles, in some clothing, and in carpets. We engineer this enzyme for improved PET degradation capacity and further demonstrate that it can also degrade an important PET replacement, polyethylene-2,5-furandicarboxylate, providing new opportunities for biobased plastics recycling., Poly(ethylene terephthalate) (PET) is one of the most abundantly produced synthetic polymers and is accumulating in the environment at a staggering rate as discarded packaging and textiles. The properties that make PET so useful also endow it with an alarming resistance to biodegradation, likely lasting centuries in the environment. Our collective reliance on PET and other plastics means that this buildup will continue unless solutions are found. Recently, a newly discovered bacterium, Ideonella sakaiensis 201-F6, was shown to exhibit the rare ability to grow on PET as a major carbon and energy source. Central to its PET biodegradation capability is a secreted PETase (PET-digesting enzyme). Here, we present a 0.92 Å resolution X-ray crystal structure of PETase, which reveals features common to both cutinases and lipases. PETase retains the ancestral α/β-hydrolase fold but exhibits a more open active-site cleft than homologous cutinases. By narrowing the binding cleft via mutation of two active-site residues to conserved amino acids in cutinases, we surprisingly observe improved PET degradation, suggesting that PETase is not fully optimized for crystalline PET degradation, despite presumably evolving in a PET-rich environment. Additionally, we show that PETase degrades another semiaromatic polyester, polyethylene-2,5-furandicarboxylate (PEF), which is an emerging, bioderived PET replacement with improved barrier properties. In contrast, PETase does not degrade aliphatic polyesters, suggesting that it is generally an aromatic polyesterase. These findings suggest that additional protein engineering to increase PETase performance is realistic and highlight the need for further developments of structure/activity relationships for biodegradation of synthetic polyesters.