Initiation of fatty acid biosynthesis in Pseudomonas putida KT2440.

Deciphering the mechanisms of bacterial fatty acid biosynthesis is crucial for both the engineering of bacterial hosts to produce fatty acid-derived molecules and the development of new antibiotics. However, gaps in our understanding of the initiation of fatty acid biosynthesis remain. Here, we demonstrate that the industrially relevant microbe Pseudomonas putida KT2440 contains three distinct pathways to initiate fatty acid biosynthesis. The first two routes employ conventional ß-ketoacyl-ACP synthase III enzymes, FabH1 and FabH2, that accept short- and medium-chain-length acyl-CoAs, respectively. The third route utilizes a malonyl-ACP decarboxylase enzyme, MadB. A combination of exhaustive in vivo alanine-scanning mutagenesis, in vitro biochemical characterization, X-ray crystallography, and computational modeling elucidate the presumptive mechanism of malonyl-ACP decarboxylation via MadB. Given that functional homologs of MadB are widespread throughout domain Bacteria, this ubiquitous alternative fatty acid initiation pathway provides new opportunities to target a range of biotechnology and biomedical applications.

Potential Pathogenicity Determinants Identified from Structural Proteomics of SARS-CoV and SARS-CoV-2.

Despite SARS-CoV and SARS-CoV-2 being equipped with highly similar protein arsenals, the corresponding zoonoses have spread among humans at extremely different rates. The specific characteristics of these viruses that led to such distinct outcomes remain unclear. Here, we apply proteome-wide comparative structural analysis aiming to identify the unique molecular elements in the SARS-CoV-2 proteome that may explain the differing consequences. By combining protein modeling and molecular dynamics simulations, we suggest nonconservative substitutions in functional regions of the spike glycoprotein (S), nsp1, and nsp3 that are contributing to differences in virulence. Particularly, we explain why the substitutions at the receptor-binding domain of S affect the structure-dynamics behavior in complexes with putative host receptors. Conservation of functional protein regions within the two taxa is also noteworthy. We suggest that the highly conserved main protease, nsp5, of SARS-CoV and SARS-CoV-2 is part of their mechanism of circumventing the host interferon antiviral response. Overall, most substitutions occur on the protein surfaces and may be modulating their antigenic properties and interactions with other macromolecules. Our results imply that the striking difference in the pervasiveness of SARS-CoV-2 and SARS-CoV among humans seems to significantly derive from molecular features that modulate the efficiency of viral particles in entering the host cells and blocking the host immune response.

Structural insights into ß-1,3-glucan cleavage by a glycoside hydrolase family.

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.

Publisher Correction: Structural insights into ß-1,3-glucan cleavage by a glycoside hydrolase family.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Characterization of aromatic acid/proton symporters in Pseudomonas putida KT2440 toward efficient microbial conversion of lignin-related aromatics.

Pseudomonas putida KT2440 (hereafter KT2440) is a well-studied platform bacterium for the production of industrially valuable chemicals from heterogeneous mixtures of aromatic compounds obtained from lignin depolymerization. KT2440 can grow on lignin-related monomers, such as ferulate (FA), 4-coumarate (4CA), vanillate (VA), 4-hydroxybenzoate (4HBA), and protocatechuate (PCA). Genes associated with their catabolism are known, but knowledge about the uptake systems remains limited. In this work, we studied the KT2440 transporters of lignin-related monomers and their substrate selectivity. Based on the inhibition by protonophores, we focused on five genes encoding aromatic acid/H+ symporter family transporters categorized into major facilitator superfamily that uses the proton motive force. The mutants of PP_1376 (pcaK) and PP_3349 (hcnK) exhibited significantly reduced growth on PCA/4HBA and FA/4CA, respectively, while no change was observed on VA for any of the five gene mutants. At pH 9.0, the conversion of these compounds by hcnK mutant (FA/4CA) and vanK mutant (VA) was dramatically reduced, revealing that these transporters are crucial for the uptake of the anionic substrates at high pH. Uptake assays using 14C-labeled substrates in Escherichia coli and biosensor-based assays confirmed that PcaK, HcnK, and VanK have ability to take up PCA, FA/4CA, and VA/PCA, respectively. Additionally, analyses of the predicted protein structures suggest that the size and hydropathic properties of the substrate-binding sites of these transporters determine their substrate preferences. Overall, this study reveals that at physiological pH, PcaK and HcnK have a major role in the uptake of PCA/4HBA and FA/4CA, respectively, and VanK is a VA/PCA transporter. This information can contribute to the engineering of strains for the efficient conversion of lignin-related monomers to value-added chemicals.

Molecular basis of substrate recognition and specificity revealed in family 12 glycoside hydrolases.

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.

Aspergillus niger ß-glucosidase has a cellulase-like tadpole molecular shape: insights into glycoside hydrolase family 3 (GH3) ß-glucosidase structure and function.

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.

Structural variants identified using non-Mendelian inheritance patterns advance the mechanistic understanding of autism spectrum disorder.

The heritability of autism spectrum disorder (ASD), based on 680,000 families and five countries, is estimated to be nearly 80%, yet heritability reported from SNP-based studies are consistently lower, and few significant loci have been identified with genome-wide association studies. This gap in genomic information may reside in rare variants, interaction among variants (epistasis), or cryptic structural variation (SV) and may provide mechanisms that underlie ASD. Here we use a method to identify potential SVs based on non-Mendelian inheritance patterns in pedigrees using parent-child genotypes from ASD families and demonstrate that they are enriched in ASD-risk genes. Most are in non-coding genic space and are over-represented in expression quantitative trait loci, suggesting that they affect gene regulation, which we confirm with their overlap of differentially expressed genes in postmortem brain tissue of ASD individuals. We then identify an SV in the GRIK2 gene that alters RNA splicing and a regulatory region of the ACMSD gene in the kynurenine pathway as significantly associated with a non-verbal ASD phenotype, supporting our hypothesis that these currently excluded loci can provide a clearer mechanistic understanding of ASD. Finally, we use an explainable artificial intelligence approach to define subgroups demonstrating their use in the context of precision medicine.

Exploring the role of plant lysin motif receptor-like kinases in regulating plant-microbe interactions in the bioenergy crop Populus.

For plants, distinguishing between mutualistic and pathogenic microbes is a matter of survival. All microbes contain microbe-associated molecular patterns (MAMPs) that are perceived by plant pattern recognition receptors (PRRs). Lysin motif receptor-like kinases (LysM-RLKs) are PRRs attuned for binding and triggering a response to specific MAMPs, including chitin oligomers (COs) in fungi, lipo-chitooligosaccharides (LCOs), which are produced by mycorrhizal fungi and nitrogen-fixing rhizobial bacteria, and peptidoglycan in bacteria. The identification and characterization of LysM-RLKs in candidate bioenergy crops including Populus are limited compared to other model plant species, thus inhibiting our ability to both understand and engineer microbe-mediated gains in plant productivity. As such, we performed a sequence analysis of LysM-RLKs in the Populus genome and predicted their function based on phylogenetic analysis with known LysM-RLKs. Then, using predictive models, molecular dynamics simulations, and comparative structural analysis with previously characterized CO and LCO plant receptors, we identified probable ligand-binding sites in Populus LysM-RLKs. Using several machine learning models, we predicted remarkably consistent binding affinity rankings of Populus proteins to CO. In addition, we used a modified Random Walk with Restart network-topology based approach to identify a subset of Populus LysM-RLKs that are functionally related and propose a corresponding signal transduction cascade. Our findings provide the first look into the role of LysM-RLKs in Populus-microbe interactions and establish a crucial jumping-off point for future research efforts to understand specificity and redundancy in microbial perception mechanisms.

Antiviral Strategies Against SARS-CoV-2: A Systems Biology Approach.

The unprecedented scientific achievements in combating the COVID-19 pandemic reflect a global response informed by unprecedented access to data. We now have the ability to rapidly generate a diversity of information on an emerging pathogen and, by using high-performance computing and a systems biology approach, we can mine this wealth of information to understand the complexities of viral pathogenesis and contagion like never before. These efforts will aid in the development of vaccines, antiviral medications, and inform policymakers and clinicians. Here we detail computational protocols developed as SARS-CoV-2 began to spread across the globe. They include pathogen detection, comparative structural proteomics, evolutionary adaptation analysis via network and artificial intelligence methodologies, and multiomic integration. These protocols constitute a core framework on which to build a systems-level infrastructure that can be quickly brought to bear on future pathogens before they evolve into pandemic proportions.

Structural and functional characterization of NEMO cleavage by SARS-CoV-2 3CLpro.

In addition to its essential role in viral polyprotein processing, the SARS-CoV-2 3C-like (3CLpro) protease can cleave human immune signaling proteins, like NF-κB Essential Modulator (NEMO) and deregulate the host immune response. Here, in vitro assays show that SARS-CoV-2 3CLpro cleaves NEMO with fine-tuned efficiency. Analysis of the 2.14 Å resolution crystal structure of 3CLpro C145S bound to NEMO 226-235 reveals subsites that tolerate a range of viral and host substrates through main chain hydrogen bonds while also enforcing specificity using side chain hydrogen bonds and hydrophobic contacts. Machine learning- and physics-based computational methods predict that variation in key binding residues of 3CLpro- NEMO helps explain the high fitness of SARS-CoV-2 in humans. We posit that cleavage of NEMO is an important piece of information to be accounted for in the pathology of COVID-19.

A mechanistic model and therapeutic interventions for COVID-19 involving a RAS-mediated bradykinin storm.

Neither the disease mechanism nor treatments for COVID-19 are currently known. Here, we present a novel molecular mechanism for COVID-19 that provides therapeutic intervention points that can be addressed with existing FDA-approved pharmaceuticals. The entry point for the virus is ACE2, which is a component of the counteracting hypotensive axis of RAS. Bradykinin is a potent part of the vasopressor system that induces hypotension and vasodilation and is degraded by ACE and enhanced by the angiotensin1-9 produced by ACE2. Here, we perform a new analysis on gene expression data from cells in bronchoalveolar lavage fluid (BALF) from COVID-19 patients that were used to sequence the virus. Comparison with BALF from controls identifies a critical imbalance in RAS represented by decreased expression of ACE in combination with increases in ACE2, renin, angiotensin, key RAS receptors, kinogen and many kallikrein enzymes that activate it, and both bradykinin receptors. This very atypical pattern of the RAS is predicted to elevate bradykinin levels in multiple tissues and systems that will likely cause increases in vascular dilation, vascular permeability and hypotension. These bradykinin-driven outcomes explain many of the symptoms being observed in COVID-19.

In late 2019, a new virus named SARS-CoV-2, which causes a disease in humans called COVID-19, emerged in China and quickly spread around the world. Many individuals infected with the virus develop only mild, symptoms including a cough, high temperature and loss of sense of smell; while others may develop no symptoms at all. However, some individuals develop much more severe, life-threatening symptoms affecting the lungs and other parts of the body including the heart and brain. SARS-CoV-2 uses a human enzyme called ACE2 like a 'Trojan Horse' to sneak into the cells of its host. ACE2 lowers blood pressure in the human body and works against another enzyme known as ACE (which has the opposite effect). Therefore, the body has to balance the levels of ACE and ACE2 to maintain a normal blood pressure. It remains unclear whether SARS-CoV-2 affects how ACE2 and ACE work. When COVID-19 first emerged, a team of researchers in China studied fluid and cells collected from the lungs of patients to help them identify the SARS-CoV-2 virus. Here, Garvin et al. analyzed the data collected in the previous work to investigate whether changes in how the body regulates blood pressure may contribute to the life-threatening symptoms of COVID-19. The analyses found that SARS-CoV-2 caused the levels of ACE in the lung cells to decrease, while the levels of ACE2 increased. This in turn increased the levels of a molecule known as bradykinin in the cells (referred to as a 'Bradykinin Storm'). . Previous studies have shown that bradykinin induces pain and causes blood vessels to expand and become leaky which will lead to swelling and inflammation of the surrounding tissue. In addition, the analyses found that production of a substance called hyaluronic acid was increased and the enzymes that could degrade it greatly decreased. Hyaluronic acid can absorb more than 1,000 times its own weight in water to form a hydrogel. The Bradykinin-Storm-induced leakage of fluid into the lungs combined with the excess hyaluronic acid would likely result in a Jello-like substance that is preventing oxygen uptake and carbon dioxide release in the lungs of severely affected COVID-19 patients. Therefore, the findings of Garvin et al. suggest that the Bradykinin Storm may be responsible for the more severe symptoms of COVID-19. Further experiments identified several existing medicinal drugs that have the potential to be re-purposed to treat the Bradykinin Storm. A possible next step would be to carry out clinical trials to assess how effective these drugs are in treating patients with COVID-19. In addition, understanding how SARS-Cov-2 affects the body will help researchers and clinicians identify individuals who are most at risk of developing life-threatening symptoms.

Crystallographic structure and molecular dynamics simulations of the major endoglucanase from Xanthomonas campestris pv. campestris shed light on its oligosaccharide products release pattern.

Cellulases are essential enzymatic components for the transformation of plant biomass into fuels, renewable materials 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.

The impact of O-glycan chemistry on the stability of intrinsically disordered proteins.

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.

Molecular characterization of a family 5 glycoside hydrolase suggests an induced-fit enzymatic mechanism.

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.

X-ray structure and molecular dynamics simulations of endoglucanase 3 from Trichoderma harzianum: structural organization and substrate recognition by endoglucanases that lack cellulose binding module.

Plant biomass holds a promise for the production of second-generation ethanol via enzymatic hydrolysis, but its utilization as a biofuel resource is currently limited to a large extent by the cost and low efficiency of the cellulolytic enzymes. Considerable efforts have been dedicated to elucidate the mechanisms of the enzymatic process. It is well known that most cellulases possess a catalytic core domain and a carbohydrate binding module (CBM), without which the enzymatic activity can be drastically reduced. However, Cel12A members of the glycosyl hydrolases family 12 (GHF12) do not bear a CBM and yet are able to hydrolyze amorphous cellulose quite efficiently. Here, we use X-ray crystallography and molecular dynamics simulations to unravel the molecular basis underlying the catalytic capability of endoglucanase 3 from Trichoderma harzianum (ThEG3), a member of the GHF12 enzymes that lacks a CBM. A comparative analysis with the Cellulomonas fimi CBM identifies important residues mediating interactions of EG3s with amorphous regions of the cellulose. For instance, three aromatic residues constitute a harboring wall of hydrophobic contacts with the substrate in both ThEG3 and CfCBM structures. Moreover, residues at the entrance of the active site cleft of ThEG3 are identified, which might hydrogen bond to the substrate. We advocate that the ThEG3 residues Asn152 and Glu201 interact with the substrate similarly to the corresponding CfCBM residues Asn81 and Arg75. Altogether, these results show that CBM motifs are incorporated within the ThEG3 catalytic domain and suggest that the enzymatic efficiency is associated with the length and position of the substrate chain, being higher when the substrate interact with the aromatic residues at the entrance of the cleft and the catalytic triad. Our results provide guidelines for rational protein engineering aiming to improve interactions of GHF12 enzymes with cellulosic substrates.

Molecular basis of the thermostability and thermophilicity of laminarinases: X-ray structure of the hyperthermostable laminarinase from Rhodothermus marinus and molecular dynamics simulations.

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.