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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

The glycoside hydrolase family 45 (GH45) of carbohydrate modifying enzymes is mostly comprised of β-1,4-endoglucanases. Significant diversity between the GH45 members has prompted the division of this family into three subfamilies: A, B and C, which may differ in terms of the mechanism, general architecture, substrate binding and cleavage. Here, we use a combination of X-ray crystallography, bioinformatics, enzymatic assays, molecular dynamics simulations and site-directed mutagenesis experiments to characterize the structure, substrate binding and enzymatic specificity of the GH45 subfamily C endoglucanase from Phanerochaete chrysosporium (PcCel45A). We investigated the role played by different residues in the binding of the enzyme to cellulose oligomers of different lengths and examined the structural characteristics and dynamics of PcCel45A that make subfamily C so dissimilar to other members of the GH45 family. Due to the structural similarity shared between PcCel45A and domain I of expansins, comparative analysis of their substrate binding was also carried out. Our bioinformatics sequence analyses revealed that the hydrolysis mechanisms in GH45 subfamily C is not restricted to use of the imidic asparagine as a general base in the “Newton’s cradle” catalytic mechanism recently proposed for this subfamily.

Carbohydrate-Active Enzymes are key enzymes for biomass-to-bioproducts conversion. α- l -Arabinofuranosidases that belong to the Glycoside Hydrolase family 62 (GH62) have important applications in biofuel production from plant biomass by hydrolyzing arabinoxylans, found in both the primary and secondary cell walls of plants. In this work, we identified a GH62 α- l -arabinofuranosidase (AnAbf62Awt) that was highly secreted when Aspergillus nidulans was cultivated on sugarcane bagasse. The gene AN7908 was cloned and transformed in A. nidulans for homologous production of AnAbf62Awt, and we confirmed that the enzyme is N-glycosylated at asparagine 83 by mass spectrometry analysis. The enzyme was also expressed in Escherichia coli and the studies of circular dichroism showed that the melting temperature and structural profile of AnAbf62Awt and the non-glycosylated enzyme from E. coli (AnAbf62Adeglyc) were highly similar. In addition, the designed glycomutant AnAbf62AN83Q presented similar patterns of secretion and activity to the AnAbf62Awt, indicating that the N-glycan does not influence the properties of this enzyme. The crystallographic structure of AnAbf62Adeglyc was obtained and the 1.7 A resolution model showed a five-bladed β-propeller fold, which is conserved in family GH62. Mutants AnAbf62AY312F and AnAbf62AY312S showed that Y312 was an important substrate-binding residue. Molecular dynamics simulations indicated that the loop containing Y312 could access different conformations separated by moderately low energy barriers. One of these conformations, comprising a local minimum, is responsible for placing Y312 in the vicinity of the arabinose glycosidic bond, and thus, may be important for catalytic efficiency.

Lignocellulosic biomass is mainly constituted by cellulose, hemicellulose, and lignin and represents an important resource for the sustainable production of biofuels and green chemistry materials. Xylans, a common hemicellulose, interact with cellulose and often exhibit various side chain substitutions including acetate, (4-O-methyl) glucuronic acid, and arabinose. Recent studies have shown that the distribution of xylan substitutions is not random, but follows patterns that are dependent on the plant taxonomic family and cell wall type. Here, we use molecular dynamics simulations to investigate the role of substitutions on xylan interactions with the hydrophilic cellulose face, using the recently discovered xylan decoration pattern of the conifer gymnosperms as a model. The results show that α-1,2-linked substitutions stabilize the binding of single xylan chains independently of the nature of the substitution and that Ca2+ ions can mediate cross-links between glucuronic acid substitutions of two neighbor...

Enzymatic conversion of lignocellulosic biomass into biofuels and chemicals constitutes a potential route for sustainable development. Cellobiohydrolases are key enzymes used in industrial cocktails for depolymerization of crystalline cellulose, and their mechanism of action has been intensely studied in the past several years. Provided with a tunnel-like substrate-binding cavity, cellobiohydrolases possess the ability to processively hydrolyze glycosidic bonds of crystalline cellulose, yielding one molecule of cellobiose per catalytic cycle. As such, cellobiose expulsion from the product binding site is a necessary step in order to allow for the processive hydrolysis mechanism. However, the high-affinity binding of cellobiose to the enzyme impairs the process and causes activity inhibition due to reaction products. Here, we use molecular dynamics simulations to study the binding of cellobiose to the Trichoderma reesei Cel7A (TrCel7A) cellobiohydrolase and the effects of mutations that reduce cellobiose binding, without affecting the structural and dynamical integrities of the enzyme. We observe that the product binding site exhibits an intrinsic flexibility that can sterically hinder cellobiose release. Several point mutations in the product binding site reduce cellobiose-enzyme interactions, but not all modifications are able to maintain the structural integrity of the enzyme. In particular, mutation of charged residues in the TrCel7A product binding site causes perturbations that affect the structure of the loops that form the substrate-binding tunnel of the enzyme and, hence, may affect TrCel7A function in other steps of the hydrolysis mechanism. Our results suggest there is a trade-off between product inhibition and catalytic efficiency, and they provide directions for cellulases engineering.

Graphene oxide (GO) is a very promising material because it is easy to process, water-soluble, and chemically versatile due to the presence of oxygenated groups on its surface. GO has been used in different areas such as electronics, biosensing, and environmental remediation. To design efficient materials, especially for biosensing and for remediating pollutants, the knowledge of surface noncovalent interaction and functionalization is crucial. Recently, it has been suggested revisions on the structural models of GO because the presence of highly oxidized polyaromatic carboxylated fragments (oxidative debris) on the GO surfaces. These debris are produced during acid treatments commonly employed in GO synthesis and purification. Here we applied chemical analysis, bioassays, and atomistic simulations to study the influence of oxidative debris on the noncovalent interaction of GO sheets with an important organic pollutant (e.g., 1-nitropyrene). GO samples without oxidative debris were found to be 75% more ef...

Aspergillus niger is known to secrete large amounts of β-glucosidases, which have a variety of biotechnological and industrial applications. Here, we purified an A. niger β-glucosidase (AnBgl1) and conducted its biochemical and biophysical analyses. Purified enzyme with an apparent molecular mass of 116 kDa forms monomers in solution as judged by native gel electrophoresis and has a pI value of 4.55, as found for most of the fungi of β-glucosidases. Surprisingly, the small angle x-ray experiments reveal that AnBgl1 has a tadpole-like structure, with the N-terminal catalytic domain and C-terminal fibronectin III-like domain (FnIII) connected by the long linker peptide (∼100 amino acid residues) in an extended conformation. This molecular organization resembles the one adopted by other cellulases (such as cellobiohydrolases, for example) that frequently contain a catalytic domain linked to the cellulose-binding module that mediates their binding to insoluble and polymeric cellulose. The reasons why AnBgl1, which acts on the small soluble substrates, has a tadpole molecular shape are not entirely clear. However, our enzyme pulldown assays with different polymeric substrates suggest that AnBgl1 has little or no capacity to bind to and to adsorb cellulose, xylan, and starch, but it has high affinity to lignin. Molecular dynamics simulations suggested that clusters of residues located in the C-terminal FnIII domain interact strongly with lignin fragments. The simulations showed that numerous arginine residues scattered throughout the FnIII surface play an important role in the interaction with lignin by means of cation-π stacking with the lignin aromatic rings. These results indicate that the C-terminal FnIII domain could be operational for immobilization of the enzyme on the cell wall and for the prevention of unproductive binding of cellulase to the biomass lignin.

Peroxisome proliferator-activated receptors (PPARs) are members of a superfamily of nuclear transcription factors. They are involved in mediating numerous physiological effects in humans, including glucose and lipid metabolism. PPAR alpha ligands effectively treat dyslipidemia and have significant antiinflammatory and anti-atherosclerotic activities. These effects and their ligand-dependent activity make nuclear receptors obvious targets for drug design. Here, we present the structure of the human PPAR alpha in complex with WY14643, a member of fibrate class of drug, and a widely used PPAR activator. The crystal structure of this complex suggests that WY14643 induces activation of PPAR alpha in an unusual bipartite mechanism involving conventional direct helix 12 stabilization and an alternative mode that involves a second ligand in the pocket. We present structural observations, molecular dynamics and activity assays that support the importance of the second site in WY14643 action. The unique binding mode of WY14643 reveals a new pattern of nuclear receptor ligand recognition and suggests a novel basis for ligand design, offering clues for improving the binding affinity and selectivity of ligand. We show that binding of WY14643 to PPAR alpha was associated with antiinflammatory disease in a human corneal cell model, suggesting possible applications for PPAR alpha ligands. (C) 2013 Elsevier Ltd. All rights reserved.

The peroxisome proliferator-activated receptor γ (PPARγ) ligands are important therapeutic drugs for the treatment of type 2 diabetes, obesity and cardiovascular diseases. In particular, partial agonists and non-agonists are interesting targets to reduce glucose levels, presenting few side effects in comparison to full agonists. In this work, we present a set of CHARMM-based parameters of a molecular mechanics force field for two PPARγ ligands, GQ16 and SR1664. GQ16 belongs to the thiazolidinedione class of drugs and it is a PPARγ partial agonist that has been shown to promote the "browning" of white adipose tissue. SR1664 is the precursor of the PPARγ non-agonist class of ligands that activates PPARγ in a non-classical manner. Here, we use quantum chemical calculations consistent with the CHARMM protocol to obtain bonded and non-bonded parameters, including partial atomic charges and effective torsion potentials for both molecules. The newly parameterized models were evaluated by examining the behavior of GQ16 and SR1664 free in water and bound to the ligand binding pocket of PPARγ using molecular dynamics simulations. The potential parameters derived here are readily transferable to a variety of pharmaceutical compounds and similar PPARγ ligands.

Cellulose is an abundant, cheap, renewable, yet recalcitrant, material, which, if dissolved, may be formed into a wide range of materials, composites, and mixtures. Much attention has recently been focused on the use of mixtures of ionic liquids and some solvents (so-called organic electrolyte solutions, OESs) as efficient cellulose dissolution solvents, but many of the cosolvents used lack green credentials - a perennial problem where dipolar aprotic solvents are the solvents of choice. We present a rational approach, based on definition of ranges of solvent parameters gathered together in recently published databases, to find "greener" cosolvents for OES formation. Thus, γ-butyrolactone is identified as a suitable OES former for dissolution of microcrystalline cellulose and biobased γ-valerolactone as a marginally less efficient, but significantly safer, alternative. Comparison of cosolvent efficiency reveals that previous use of measures of mass, or concentration, of cellulose dissolved may have masked the similarities between 1-methylimidazole, dimethyl sulfoxide (DMSO), N,N-dimethylformamide, N-N′-dimethylimidazolidinone, N,N-dimethylacetamide, N-methylpyrrolidinone, and sulfolane (seldom considered), while comparison on a molar basis reveals that the molar volume of the solvent is an important factor. Reference-interaction site model (RISM) calculations for the DMSO/1-ethyl-3-methylimidazolium acetate OES suggest competition between DMSO and the acetate anion and preferential solvation of cellulose by the ionic liquid.

Made available in DSpace on 2018-11-26T15:29:11Z (GMT). No. of bitstreams: 0 Previous issue date: 2016-04-01 Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) Glycoside hydrolases (GHs) play fundamental roles in the decomposition of lignocellulosic biomaterials. Here, we report the full-length structure of a cellulase from Bacillus licheniformis (BlCel5B), a member of the GH5 subfamily 4 that is entirely dependent on its two ancillary modules (Ig-like module and CBM46) for catalytic activity. Using X-ray crystallography, small-angle X-ray scattering and molecular dynamics simulations, we propose that the C-terminal CBM46 caps the distal N-terminal catalytic domain (CD) to establish a fully functional active site via a combination of large-scale multidomain conformational selection and induced-fit mechanisms. The Ig-like module is pivoting the packing and unpacking motions of CBM46 relative to CD in the assembly of the binding subsite. This is the first example of a multidomain GH relying on large amplitude motions of the CBM46 for assembly of the catalytically competent form of the enzyme. Univ Sao Paulo, Sao Carlos Inst Phys, BR-13566590 Sao Paulo, Brazil Univ Estadual Campinas, Inst Chem, BR-13084862 Sao Paulo, Brazil State Univ Sao Paulo, Inst Biosci, BR-18618970 Sao Paulo, Brazil European Synchrotron Radiat Facil, CS40220, Grenoble, France State Univ Sao Paulo, Inst Biosci, BR-18618970 Sao Paulo, Brazil FAPESP: 2008/56255-9 FAPESP: 2009/52840-7 FAPESP: 2010/18773-8 FAPESP: 2013/08293-7 FAPESP: 2013/15582-5 FAPESP: 2014/10448-1 CNPq: 490022/2009-0 CNPq: 301981/2011-6 CNPq: 500091/2014-5 CNPq: 310177/2011-1

Supercritical carbon dioxide is a promising green-chemistry solvent for many enzyme-catalyzed chemical reactions, yet the striking stability of some enzymes in such unconventional environments is not well understood. Here, we investigate the stabilization of the Candida antarctica Lipase B (CALB) in supercritical carbon dioxide-water biphasic systems using molecular dynamics simulations. The preservation of the enzyme structure and optimal activity depend on the presence of small amounts of water in the supercritical dispersing medium. When the protein is at least partially hydrated, water molecules bind to specific sites on the enzyme surface and prevent carbon dioxide from penetrating its catalytic core. Strikingly, water and supercritical carbon dioxide cover the protein surface quite heterogeneously. In the first solvation layer, the hydrophilic residues at the surface of the protein are able to pin down patches of water, whereas carbon dioxide solvates preferentially hydrophobic surface residues. In the outer solvation shells, water molecules tend to cluster predominantly on top of the larger water patches of the first solvation layer instead of spreading evenly around the remainder of the protein surface. For CALB, this exposes the substrate-binding region of the enzyme to carbon dioxide, possibly facilitating diffusion of nonpolar substrates into the catalytic funnel. Therefore, by means of microheterogeneous solvation, enhanced accessibility of hydrophobic substrates to the active site can be achieved, while preserving the functional structure of the enzyme. Our results provide a molecular picture on the nature of the stability of proteins in nonaqueous media.

Cellulose-builder is a user-friendly program that builds crystalline structures of cellulose of different sizes and geometries. The program generates Cartesian coordinates for all atoms of the specified structure in the Protein Data Bank format, suitable for using as starting configurations in molecular dynamics simulations and other calculations. Crystalline structures of cellulose polymorphs Iα, Iβ, II, and III(I) of practically any size are readily constructed which includes parallelepipeds, plant cell wall cellulose elementary fibrils of any length, and monolayers. Periodic boundary conditions along the crystallographic directions are easily imposed. The program also generates atom connectivity file in PSF format, required by well-known simulation packages such as NAMD, CHARMM, and others. Cellulose-builder is based on the Bash programming language and should run on practically any Unix-like platform, demands very modest hardware, and is freely available for download from ftp://ftp.iqm.unicamp.br/pub/cellulose-builder.

Hb S-São Paulo (SP) [HBB:c.20A > T p.Glu6Val; c.196A > G p.Lys65Glu] is a new double-mutant hemoglobin that was found in heterozygosis in an 18-month-old Brazilian male with moderate anemia. It behaves like Hb S in acid electrophoresis, isoelectric focusing and solubility testing but shows different behavior in alkaline electrophoresis, cation-exchange HPLC and RP-HPLC. The variant is slightly unstable, showed reduced oxygen affinity and also appeared to form polymers more stable than the Hb S. Molecular dynamics simulation suggests that the polymerization is favored by interfacial electrostatic interactions. This provides a plausible explanation for some of the reported experimental observations Fil: Jorge, Susan E. D. C.. Universidade Estadual de Campinas; Brasil Fil: Petruk, Ariel Alcides. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Tucumán. Instituto Superior de Investigaciones Biológicas. Universidad Nacional de Tucumán. Instituto Superior de Investigaciones Biológicas; Argentina Fil: Kimura, Elza M.. Universidade Estadual de Campinas; Brasil Fil: Oliveira, Denise M.. Universidade Estadual de Campinas; Brasil Fil: Caire, Lucas. Universidade Estadual de Campinas; Brasil Fil: Suemasu, Cintia N.. Universidade Estadual de Campinas; Brasil Fil: Silveira, Paulo A. A.. Albert Einstein Institute; Brasil Fil: Albuquerque, Dulcineia M.. Universidade Estadual de Campinas; Brasil Fil: Costa, Fernando F.. Universidade Estadual de Campinas; Brasil Fil: Skaf, Munir S.. Universidade Estadual de Campinas; Brasil Fil: Martínez, Leandro. Universidade de Sao Paulo; Brasil Fil: De Fatima Sonati, Maria. Universidade Estadual de Campinas; Brasil

Glycosyl hydrolases are enzymes capable of breaking the glycosidic linkage of polysaccharides and have considerable industrial and biotechnological applications. Driven by the later applications, it is frequently desirable that glycosyl hydrolases display stability and activity under extreme environment conditions, such as high temperatures and extreme pHs. Here, we present X-ray structure of the hyperthermophilic laminarinase from Rhodothermus marinus (RmLamR) determined at 1.95 Å resolution and molecular dynamics simulation studies aimed to comprehend the molecular basis for the thermal stability of this class of enzymes. As most thermostable proteins, RmLamR contains a relatively large number of salt bridges, which are not randomly distributed on the structure. On the contrary, they form clusters interconnecting β-sheets of the catalytic domain. Not all salt bridges, however, are beneficial for the protein thermostability: the existence of charge-charge interactions permeating the hydrophobic core of the enzymes actually contributes to destabilize the structure by facilitating water penetration into hydrophobic cavities, as can be seen in the case of mesophilic enzymes. Furthermore, we demonstrate that the mobility of the side-chains is perturbed differently in each class of enzymes. The side-chains of loop residues surrounding the catalytic cleft in the mesophilic laminarinase gain mobility and obstruct the active site at high temperature. By contrast, thermophilic laminarinases preserve their active site flexibility, and the active-site cleft remains accessible for recognition of polysaccharide substrates even at high temperatures. The present results provide structural insights into the role played by salt-bridges and active site flexibility on protein thermal stability and may be relevant for other classes of proteins, particularly glycosyl hydrolases.

We develop a CHARMM-based interaction potential for rosiglitazone, a well-known selective ligand to the c isoform of the peroxisome proliferator-activated receptor (PPARc) and widely marketed antidiabetic drug of the thiazolidinedione (TZD) class. We derive partial atomic charges and dihedral torsion potentials for seven rotations in the molecule, for which there are no analogs available in CHARMM. The potential model is validated by performing a series of molecular dynamics simulations of rosiglitazone in neat water and of a fully solvated rosiglitazone-PPARc complex. The structural and dynamical behavior of the complex is analyzed in comparison with available experimental data. The potential parameters derived here are readily transferable to a variety of pharmaceutically important TZD compounds. V C 2010 Wiley Periodicals, Inc. Int J Quantum Chem 111: 1346-1354, 2011

We present a set of CHARMM-based parameters of molecularmechanics force field for neutral articaine (ATC), a potent and widely used localanesthetic in medical and dental applications. Partial atomic charges and effectivetorsion potentials around four chemical bonds, not available in CHARMM for thismolecule, were obtained from quantum chemical calculations. The newly parameterizedmodel was evaluated by examining the behavior of ATC in water and in aphospholipid bilayer using molecular dynamics simulations. V C 2010 Wiley Periodicals, Inc.Int J Quantum Chem 111: 1339–1345, 2011 Key words: articaine; local anesthetics; CHARMM parameters; molecular dynamics 1. Introduction L ocal anesthetics (LA) are compounds thatcan relieve and prevent pain by reversiblyinterrupting nerve conduction. As the nameimplies, these drugs act locally, anywhere in thenervous system, and in all types of nerve fibers.Typically, the structure of LA comprises a hydro-philic (an amine) and a hydrophobic (an aromaticring, usually benzene) region, which are sepa-rated by an amide or ester linkage. Lidocaine andbupivacaine, for instance, which are among themost commonly used LA [1], are classified asamino amides, since the hydrophilic group is atertiary amine and the linkage with the hydropho-bic group (an aromatic ring) is an amide. Pro-caine, on the other hand, has an ester as interme-diate group and is, therefore, classified as anamino ester LA [2, 3].The anesthetic effect occurs because the LAmolecules diffuse into the phospholipid mem-branes and cause conformational changes in thevoltage-dependent sodium channel, leading to thetemporary inactivation of this transmembrane

Understanding the molecular basis of the binding modes of natural and synthetic ligands to nuclear receptors is fundamental to our comprehension of the activation mechanism of this important class of hormone regulated transcription factors and to the development of new ligands. Thyroid hormone receptors (TR) are particularly important targets for pharmaceuticals development because TRs are associated with the regulation of metabolic rates, body weight, and circulating levels of cholesterol and triglycerides in humans. While several high-affinity ligands are known, structural information is only partially available. In this work we obtain structural models of several TR-ligand complexes with unknown structure by docking high affinity ligands to the receptors' ligand binding domain with subsequent relaxation by molecular dynamics simulations. The binding modes of these ligands are discussed providing novel insights into the development of TR ligands. The experimental binding free energies are reasonably well-reproduced from the proposed models using a simple linear interaction energy free-energy calculation scheme.

Molecular dynamics simulations have been carried out to investigate the dynamics of fructose aqueous solutions up to 70 wt % concentration. We find that the hydrogen (H)-bonded network of fructose molecules extends with increasing sugar content and forms a structurally heterogeneous system around and above 45 wt % concentration, characterized as a percolated-like solute domain permeated by patchy regions of solvent. The presence of such aggregates in concentrated solutions promotes the slowing down of water translational, reorientational, and H-bonding dynamics, typical of many biomolecular environments. Analysis of the effects of the topological and energetic disorder of the sugar aggregates on vicinal water dynamics, similar to that recently carried out for the hydration layer of proteins by Pizzitutti et al. (J. Phys. Chem. B 2007, 111, 7584), reveals many similarities between the dynamical anomaly of the hydration layers of both systems. Like a protein surface, topological and energetic disorders of the sugar aggregates both contribute to the translational diffusion anomaly. However, unlike in the vicinity of a protein surface, the rotational relaxation is also hindered by the topological disorder created by the intertwined, percolating sugar clusters in concentrated solutions.

The solvation properties of 5-fluorouracil in supercritical carbon dioxide (SC-CO2) is investigated by means of molecular dynamics computer simulations using a force field that treats all atoms explicitly. Solvation shell distribution functions and solvent density maps around the solute were computed at different supercritical temperatures in the range 308–328 K and densities varying from 0.117 to 0.94 g/cm3. At the lowest density, we detect a strongly inhomogeneous local density augmentation effect manifested as the clustering of solvent molecules above and below the solute molecular plane and around its carbonyl and NH groups. The free-energy benefit for a CO2 molecule to be in the first solvation shell relative to the bulk is roughly twice the thermal energy at such conditions. The distribution of the relative orientation of CO2 with respect to the solute suggests that local bond multipolar solute–solvent interactions are responsible for the solvent density enhancement in these systems. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008

Understanding the nature of allostery in DNA-nuclear receptor (NR) complexes is of fundamental importance for drug development since NRs regulate the transcription of a myriad of genes in humans and other metazoans. Here, we investigate allostery in the peroxisome proliferator-activated/retinoid X receptor heterodimer. This important NR complex is a target for antidiabetic drugs since it binds to DNA and functions as a transcription factor essential for insulin sensitization and lipid metabolism. We find evidence of interdependent motions of Ω-loops and PPARγ-DNA binding domain with contacts susceptible to conformational changes and mutations, critical for regulating transcriptional functions in response to sequence-dependent DNA dynamics. Statistical network analysis of the correlated motions, observed in molecular dynamics simulations, shows preferential allosteric pathways with convergence centers comprised of polar amino acid residues. These findings are particularly relevant for the design of allosteric modulators of ligand-dependent transcription factors.

Molecular dynamics simulations have been performed to investigate the interactions between chlorpromazine (CPZ) and Langmuir monolayers of the zwitterionic dipalmitoylphosphatidylcholine (DPPC) and the anionic dipalmitoylphosphatidylglycerol (DPPG). Simulations for a fixed surface density and different charge states - neutral and protonated CPZ - were able to capture important features of the CPZ-phospholipid monolayer interaction. Neutral CPZ is predominantly found in the hydrophobic tail region, whereas protonated CPZ is located at the lipid-water interface. Specific interactions (hydrogen bonds) between protonated CPZ and the lipid head groups were found for both zwitterionic and anionic monolayers. We computed lipid tail order parameters and investigated the effects of the drug upon tail ordering. We also computed electrostatic surface potentials and found qualitative good agreement with experimental results.