The investigation of type-specific features of the copper coordinating AA9 proteins and their effect on the interaction with crystalline cellulose using molecular dynamics studies

AA9 proteins are metallo-enzymes which are crucial for the early stages of cellulose degradation. AA9 proteins have been suggested to cleave glycosidic bonds linking cellulose through the use of their Cu2+ coordinating active site. AA9 proteins possess different regioselectivities depending on the resulting cleavage they form and as result, are grouped accordingly. Type 1 AA9 proteins cleave the C1 carbon of cellulose while Type 2 AA9 proteins cleave the C4 carbon and Type 3 AA9 proteins cleave either C1 or C4 carbons. The steric congestion of the AA9 active site has been proposed to be a contributor to the observed regioselectivity. As such, a bioinformatics characterisation of type-specific sequence and structural features was performed. Initially AA9 protein sequences were obtained from the Pfam database and multiple sequence alignment was performed. The sequences were phylogenetically characterised and sequences were grouped into their respective types and sub-groups were identified. A selection analysis was performed on AA9 LPMO types to determine the selective pressure acting on AA9 protein residues. Motif discovery was then performed to identify conserved sequence motifs in AA9 proteins. Once type-specific sequence features were identified structural mapping was performed to assess possible effects on substrate interaction. Physicochemical property analysis was also performed to assess biochemical differences between AA9 LPMO types. Molecular dynamics (MD) simulations were then employed to dynamically assess the consequences of the discovered type-specific features on AA9-cellulose interaction. Due to the absence of AA9 specific force field parameters MD simulations were not readily applicable. As a result, Potential Energy Surface (PES) scans were performed to evaluate the force field parameters for the AA9 active site using the PM6 semi empirical approach and least squares fitting. A Type 1 AA9 active site was constructed from the crystal structure 4B5Q, encompassing only the Cu2+ coordinating residues, the Cu2+ ion and two water residues. Due to the similarity in AA9 active sites, the Type force field parameters were validated on all three AA9 LPMO types. Two MD simulations for each AA9 LPMO types were conducted using two separate Lennard-Jones parameter sets. Once completed, the MD trajectories were analysed for various features including the RMSD, RMSF, radius of gyration, coordination during simulation, hydrogen bonding, secondary structure conservation and overall protein movement. Force field parameters were successfully evaluated and validated for AA9 proteins. MD simulations of AA9 proteins were able to reveal the presence of unique type-specific binding modes of AA9 active sites to cellulose. These binding modes were characterised by the presence of unique type-specific loops which were present in Type 2 and 3 AA9 proteins but not in Type 1 AA9 proteins. The loops were found to result in steric congestion that affects how the Cu2+ ion interacts with cellulose. As a result, Cu2+ binding to cellulose was observed for Type 1 and not Type 2 and 3 AA9 proteins. In this study force field parameters have been evaluated for the Type 1 active site of AA9 proteins and this parameters were evaluated on all three types and binding. Future work will focus on identifying the nature of the reactive oxygen species and performing QM/MM calculations to elucidate the reactive mechanism of all three AA9 LPMO types

The large scale bioinformatics analysis of auxiliary activity family 9 enzymes

Biofuels have been proposed to be a suitable replacement to the already depleting fossil fuels. The complex structures of plant biomasses present a challenge the production of biofuels due to recalcitrance. The complex cellulose structure and hydrogen bonding between repeat units of cellulose is believed to be a major contributor to the recalcitrance of cellulose. Fungal organisms come equipped with various oxidative enzymes involved in degradation of plant biomass. The exact mechanism of cellulose degradation remains elusive. The GH61 is a group of proteins which are PMOs. GH61 sequences where previously described as endoglucanases due to weak endoglucanase activity. These enzymes were later found not possess any enzyme activity of their own however they could enhance the activity of other cellulose degrading enzymes. As a result reclassification of these enzymes as AA9 has been implemented. AA9 proteins have been reported to share structural homology with the bacterial AA10 group of enzymes. Based on cleavage products that are produced when AA9 proteins interact with cellulose, AA9 proteins have been grouped into three types. To date the exact mechanism and the sequence and structural basis for differentiating between the various AA9 types remains unknown. Using various bionformatic techniques sequence and structural elements were identified for distinguishing between the AA9 types. A large dataset of sequences was obtained from the Pfam database from UNIPROT entries. Due to high divergence of AA9 sequences, a smaller dataset with the more divergent sequences removed was created. The inclusion of the reference sequences to the data set was done to observe which sequences belong to a certain type. Phylogenetic analysis was able to group AA9 proteins into three distinct groups. MSA and motif analysis revealed that the N-Terminus of these proteins is mostly responsible for type specificity. Structural analysis of AA9 PDB structures and homology models allowed the effect of physicochemical properties to be gauged structurally. The presence of 310 helices and aromatic residues the surface of AA9 sequences is an observation which still warrants further investigation

AMBER force field parameters for the Zn (II) ions of the tunneling-fold enzymes GTP cyclohydrolase I and 6-pyruvoyl tetrahydropterin synthase

The folate biosynthesis pathway is an essential pathway for cell growth and survival. Folate derivatives serve as a source of the one-carbon units in several intracellular metabolic reactions. Rapidly dividing cells rely heavily on the availability of folate derivatives for their proliferation. As a result, drugs targeting this pathway have shown to be effective against tumor cells and pathogens, but drug resistance against the available antifolate drugs emerged quickly. Therefore, there is a need to develop new treatment strategies and identify alternative metabolic targets. The two de novo folate biosynthesis pathway enzymes, GTP cyclohydrolase I (GCH1) and 6-pyruvoyl tetrahydropterin synthase (PTPS), can provide an alternative strategy to overcome the drug resistance that emerged in the two primary targeted enzymes dihydrofolate reductase and dihydropteroate synthase. Both GCH1 and PTPS enzymes contain Zn2+ ions in their active sites, and to accurately study their dynamic behaviors using all-atom molecular dynamics (MD) simulations, appropriate parameters that can describe their metal sites should be developed and validated. In this study, force field parameters of the GCH1 and PTPS metal centers were generated using quantum mechanics (QM) calculations and then validated through MD simulations to ensure their accuracy in describing and maintaining the Zn2+ ion coordination environment. The derived force field parameters will provide accurate and reliable MD simulations involving these two enzymes for future in-silico identification of drug candidates against the GCH1 and PTPS enzymes

Analysis of the galactomannan binding ability of β-mannosidases, BtMan2A and CmMan5A, regarding their activity and synergism with a β-mannanase

Both β-mannanases and β-mannosidases are required for mannan-backbone degradation into mannose. In this study, two β-mannosidases of glycoside hydrolase (GH) families 2 (BtMan2A) and 5 (CmMan5A) were evaluated for their substrate specificities and galactomannan binding ability. BtMan2A preferred short manno-oligomers, while CmMan5A preferred longer ones; DP >2, and galactomannans. BtMan2A displayed irreversible galactomannan binding, which was pH-dependent, with higher binding observed at low pH, while CmMan5A had limited binding. Docking and molecular dynamics (MD) simulations showed that BtMan2A galactomannan binding was stronger under acidic conditions (-8.4 kcal/mol) than in a neutral environment (-7.6 kcal/mol), and the galactomannan ligand was more unstable under neutral conditions than acidic conditions. Qualitative surface plasmon resonance (SPR) experimentally confirmed the reduced binding capacity of BtMan2A at pH 7. Finally, synergistic β-mannanase to β-mannosidase (BtMan2A or CmMan5A) ratios required for maximal galactomannan hydrolysis were determined. All CcManA to CmMan5A combinations were synergistic (≈1.2-fold), while combinations of CcManA with BtMan2A (≈1.0-fold) yielded no hydrolysis improvement. In conclusion, the low specific activity of BtMan2A towards long and galactose-containing oligomers and its non-catalytic galactomannan binding ability led to no synergy with the mannanase, making GH2 mannosidases ineffective for use in cocktails for mannan degradation.The National Research Foundation (NRF) of South Africa, Rhodes University and the University of Pretoria.http://www.elsevier.com/locate/csbjhj2023BiochemistryGeneticsMicrobiology and Plant Patholog

AMBER force field parameters for the Zn (II) ions of the tunneling-fold enzymes GTP cyclohydrolase I and 6-pyruvoyl tetrahydropterin synthase:

The folate biosynthesis pathway is an essential pathway for cell growth and survival. Folate derivatives serve as a source of the one-carbon units in several intracellular metabolic reactions. Rapidly dividing cells rely heavily on the availability of folate derivatives for their proliferation. As a result, drugs targeting this pathway have shown to be effective against tumor cells and pathogens, but drug resistance against the available antifolate drugs emerged quickly. Therefore, there is a need to develop new treatment strategies and identify alternative metabolic targets. The two de novo folate biosynthesis pathway enzymes, GTP cyclohydrolase I (GCH1) and 6-pyruvoyl tetrahydropterin synthase (PTPS), can provide an alternative strategy to overcome the drug resistance that emerged in the two primary targeted enzymes dihydrofolate reductase and dihydropteroate synthase

The determination of CHARMM force field parameters for the Mg2+ containing HIV-1 integrase:

The HIV integrase enzyme is a validated drug target. However, its potential has remained largely unexploited until recently due to lack of structural and mechanistic information. Its catalytic core domain (CCD) is crucial for the viral-human DNA integration making integrase an ideal target for inhibitor design. However, in order to do so, force field parameters for the integrase magnesium ion need to be established. Quantum mechanical calculations were used to derive force field parameters which were validated through molecular dynamics studies. Our results show that the parameters determined accurately maintain the integrity of the metal pocket of the integrase CCD

Commento art. 29 DPR 448/1988

Additional file 3. Multiple sequence alignment of the Type 1 PMO sequences used in this study

The evaluation and validation of copper (II) force field parameters of the Auxiliary Activity family 9 enzymes:

The Auxiliary Activity family 9 (AA9) proteins are Cu2+ coordinating enzymes which are crucial for the early stages of cellulose degradation. In this study, the force field parameters for copper-containing bonds in the Type 1 AA9 protein active site were established and used in a molecular dynamics simulation on a solvated, neutralized system containing an AA9 protein, Cu2+ and a β-cellulose surface. The copper to cellulose interaction was evident during the dynamics, which could also be accelerated by the use of high Cu O van der Waals parameters. The interaction of AA9, Cu2+ and cellulose is described in detail