RESUMO
The carbohydrate binding module 21 (CBM21) from Rhizopus oryzae is a dual-site CBM proposed to disrupt polysaccharide structures. Additionally, it serves as a purification tag in industry. CBM21 crystal structure features a Glc residue in an unusual 1S3 conformation, whose relevance for the CBM mechanism of action is unclear. In this context, we seek to contribute for the understanding of CBM21 mechanism of action by: i) investigating the role of the 1S3 conformation on carbohydrate recognition, and ii) characterize the protein-carbohydrate binding dynamics using molecular dynamics and metadynamics simulations at MM and QM/MM levels. Results indicate the 1S3 Glc conformation is unlikely to occur under biological conditions, being originated from the crystallographic environment. CBM21 binding to small ligands appears transient and unstable, while protein dimerization and polysaccharide chain size influence complex stability. In interactions with amylose, CBM21 exhibits a repeated unbinding followed by re-binding, while simultaneously alternating between binding sites I and II. These results suggest that CBM21 acts through transient interactions, directing carbohydrates to the catalytic center rather than forming strong and long-lasting bonds with carbohydrates. Accordingly, we expect such atomistic depiction of CBM21 mechanism could aid in CBM design targeting biotechnological applications.
Assuntos
Amilose , Módulos de Ligação de Carboidratos , Carboidratos/química , Polissacarídeos/química , Sítios de Ligação , Ligação ProteicaRESUMO
Human leishmaniasis is a neglected tropical disease which affects nearly 1.5 million people every year, with Mexico being an important endemic region. One of the major defense mechanisms of these parasites is based in the polyamine metabolic pathway, as it provides the necessary compounds for its survival. Among the enzymes in this route, trypanothione reductase (TryR), an oxidoreductase enzyme, is crucial for the Leishmania genus' survival against oxidative stress. Thus, it poses as an attractive drug target, yet due to the size and features of its catalytic pocket, modeling techniques such as molecular docking focusing on that region is not convenient. Herein, we present a computational study using several structure-based approaches to assess the druggability of TryR from L. mexicana, the predominant Leishmania species in Mexico, beyond its catalytic site. Using this consensus methodology, three relevant pockets were found, of which the one we call σ-site promises to be the most favorable one. These findings may help the design of new drugs of trypanothione-related diseases.
Assuntos
Antiprotozoários , Leishmania , Leishmaniose , Humanos , Simulação de Acoplamento Molecular , Leishmania/metabolismo , NADH NADPH Oxirredutases/metabolismo , Leishmaniose/parasitologia , Antiprotozoários/uso terapêuticoRESUMO
Phytases [myo-inositol(1,2,3,4,5,6) hexakisphosphate phosphohydrolases] are phytate-specific phosphatases not present in monogastric animals. Nevertheless, they are an essential supplement to feeding such animals and for human special diets. It is crucial, hence, the biotechnological use of phytases with intrinsic stability and activity at the acid pHs from gastric environments. Here we use Metadynamics (METADY) simulations to probe the conformational space of the Aspergillus nidulans phytase and the differential effects of pH and glycosylation in this same space. The results suggest that strategic combinations of pH and glycosylation affect the stability of native-like conformations and alternate these structures from a metastable to a stable profile. Furthermore, the protein segments previously reported as more thermosensitive in phytases from this family present a pivotal role in the conformational changes at different conditions, especially H2, H5-7, L8, L10, L12, and L17. Also, the glycosylations and the pH-dependent charge balance modulate the mobility and interactions at these same regions, with consequences for the surface solvation and active site exposition. Finally, although the glycosylations have stabilized the native structure and improved the substrate docking at all the studied pHs, the data suggest a higher phytate receptivity at catalytic poses for the unglycosylated structure at pH 6.5 and the glycosylated one at pH 4.5. This behavior agrees with the exact change in optimum pH reported for this enzyme, expressed on low or high glycosylating systems. We hope the results and insights presented here will be helpful in future approaches for rational engineering of technologically promising phytases and intelligent planning of their heterologous expression systems and conditions for use.Communicated by Ramaswamy H. Sarma.
RESUMO
The enzymatic breakdown of carbohydrates plays a critical role in several biological events and enables the development of sustainable processes to obtain bioproducts and biofuels. In this scenario, the design of efficient inhibitors for glycosidases that can act as drug targets and the engineering of carbohydrate-active enzymes with tailored catalytic properties is of remarkable importance. To guide rational approaches, it is necessary to elucidate enzyme molecular mechanisms, in particular understanding how the microenvironment modulates the conformational space explored by the substrate. Computer simulations, especially those based on ab initio methods, have provided a suitable atomic description of carbohydrate conformations and catalytic reactions in several glycosidase families. In this review, we will focus on how the active-site topology (pocket or cleft) and mode of cleavage (endo or exo) can affect the catalytic mechanisms adopted by glycosidases, in particular the substrate conformations along the reaction coordinate.
Assuntos
Carboidratos , Glicosídeo Hidrolases , Humanos , Glicosídeo Hidrolases/metabolismo , Configuração de Carboidratos , Domínio Catalítico , AçúcaresRESUMO
ß-glucosidases play a pivotal role in second-generation biofuel (2G-biofuel) production. For this application, thermostable enzymes are essential due to the denaturing conditions on the bioreactors. Random amino acid substitutions have originated new thermostable ß-glucosidases, but without a clear understanding of their molecular mechanisms. Here, we probe by different molecular dynamics simulation approaches with distinct force fields and submitting the results to various computational analyses, the molecular bases of the thermostabilization of the Paenibacillus polymyxa GH1 ß-glucosidase by two-point mutations E96K (TR1) and M416I (TR2). Equilibrium molecular dynamic simulations (eMD) at different temperatures, principal component analysis (PCA), virtual docking, metadynamics (MetaDy), accelerated molecular dynamics (aMD), Poisson-Boltzmann surface analysis, grid inhomogeneous solvation theory and colony method estimation of conformational entropy allow to converge to the idea that the stabilization carried by both substitutions depend on different contributions of three classic mechanisms: (i) electrostatic surface stabilization; (ii) efficient isolation of the hydrophobic core from the solvent, with energetic advantages at the solvation cap; (iii) higher distribution of the protein dynamics at the mobile active site loops than at the protein core, with functional and entropic advantages. Mechanisms i and ii predominate for TR1, while in TR2, mechanism iii is dominant. Loop A integrity and loops A, C, D, and E dynamics play critical roles in such mechanisms. Comparison of the dynamic and topological changes observed between the thermostable mutants and the wildtype protein with amino acid co-evolutive networks and thermostabilizing hotspots from the literature allow inferring that the mechanisms here recovered can be related to the thermostability obtained by different substitutions along the whole family GH1. We hope the results and insights discussed here can be helpful for future rational approaches to the engineering of optimized ß-glucosidases for 2G-biofuel production for industry, biotechnology, and science.
Assuntos
Biocombustíveis , beta-Glucosidase , beta-Glucosidase/genética , beta-Glucosidase/química , beta-Glucosidase/metabolismo , Substituição de Aminoácidos , Simulação de Dinâmica Molecular , Domínio CatalíticoRESUMO
Angiotensin-(1-7) re-balance the Renin-Angiotensin system affected during several pathologies, including the new COVID-19; cardiovascular diseases; and cancer. However, one of the limiting factors for its therapeutic use is its short half-life, which might be overcome with the use of dendrimers as nanoprotectors. In this work, we addressed the following issues: (1) the capacity of our computational protocol to reproduce the experimental structural features of the (hydroxyl/amino)-terminated PAMAM dendrimers as well as the Angiotensin-(1-7) peptide; (2) the coupling of Angiotensin-(1-7) to (hydroxyl/amino)-terminated PAMAM dendrimers in order to gain insight into the structural basis of its molecular binding; (3) the capacity of the dendrimers to protect Angiotensin-(1-7); and (4) the effect of pH changes on the peptide binding and covering. Our Molecular-Dynamics/Metadynamics-based computational protocol well modeled the structural experimental features reported in the literature and our double-docking approach was able to provide reasonable initial structures for stable complexes. At neutral pH, PAMAM dendrimers with both terminal types were able to interact stably with 3 Angiotensin-(1-7) peptides through ASP1, TYR4 and PRO7 key amino acids. In general, they bind on the surface in the case of the hydroxyl-terminated compact dendrimer and in the internal zone in the case of the amino-terminated open dendrimer. At acidic pH, PAMAM dendrimers with both terminal groups are still able to interact with peptides either internalized or in its periphery, however, the number of contacts, the percentage of coverage and the number of hydrogen bonds are lesser than at neutral pH, suggesting a state for peptide release. In summary, amino-terminated PAMAM dendrimer showed slightly better features to bind, load and protect Angiotensin-(1-7) peptides.
Assuntos
COVID-19 , Dendrímeros , Aminoácidos , Angiotensina I , Dendrímeros/química , Humanos , Simulação de Dinâmica Molecular , Fragmentos de Peptídeos , PeptídeosRESUMO
Crystallin aggregation in the eye lens is one of the leading causes of cataract formation. The increase in the human γD-crystallin (HγDC) aggregation propensity has been associated with the oligomerization of its partially folded and fully unfolded structure. A recent study demonstrated that the binding of flavonoid morin (MOR) to HγDC inhibits the fibrillation of this protein. In this work, we carry out an exhaustive search for the possible binding site of MOR on HγDC by combining an ensemble docking approach with the Wrap 'N' Shake protocol. In agreement with previous results, we found a potential MOR-binding site in the cleft formed between the N-terminal and C-terminal domains of HγDC. MOR preference for the cleft residues was observed even with the interface-opened intermediate conformers of HγDC. In addition, metadynamics simulations were carried out to corroborate the stabilizing activity of MOR on HγDC structure and to identify the structural regions implicated during the unfolding inhibition. Overall, this study provides relevant insights into the identification of new HγDC aggregation inhibitors.
Assuntos
Flavonoides , Sítios de Ligação , HumanosRESUMO
The main-protease (Mpro) catalyzes a crucial step for the SARS-CoV-2 life cycle. The recent SARS-CoV-2 presents the main protease (MCoV2pro) with 12 mutations compared to SARS-CoV (MCoV1pro). Recent studies point out that these subtle differences lead to mobility variances at the active site loops with functional implications. We use metadynamics simulations and a sort of computational analysis to probe the dynamic, pharmacophoric and catalytic environment differences between the monomers of both enzymes. So, we verify how much intrinsic distinctions are preserved in the functional dimer of MCoV2pro, as well as its implications for ligand accessibility and optimized drug screening. We find a significantly higher accessibility to open binding conformers in the MCoV2pro monomer compared to MCoV1pro. A higher hydration propensity for the MCoV2pro S2 loop with the A46S substitution seems to exercise a key role. Quantum calculations suggest that the wider conformations for MCoV2pro are less catalytically active in the monomer. However, the statistics for contacts involving the N-finger suggest higher maintenance of this activity at the dimer. Docking analyses suggest that the ability to vary the active site width can be important to improve the access of the ligand to the active site in different ways. So, we carry out a multiconformational virtual screening with different ligand bases. The results point to the importance of taking into account the protein conformational multiplicity for new promissors anti MCoV2pro ligands. We hope these results will be useful in prospecting, repurposing and/or designing new anti SARS-CoV-2 drugs.Communicated by Ramaswamy H. Sarma.
Assuntos
COVID-19 , SARS-CoV-2 , Humanos , SARS-CoV-2/metabolismo , Domínio Catalítico , Ligantes , Inibidores de Proteases/farmacologia , Inibidores de Proteases/química , Proteínas não Estruturais Virais/química , Antivirais/farmacologia , Antivirais/química , Simulação de Acoplamento Molecular , Simulação de Dinâmica Molecular , Cisteína Endopeptidases/químicaRESUMO
Cyclin-Dependent Kinase 2 (CDK2) and Vascular-Endothelial Growth Factor Receptor 2 (VEGFR2) are promising targets for the design of novel inhibitors in anticancer therapeutics. In a recent work, our group designed a set of potential dual inhibitors predicted to occupy an allosteric back pocket near the active site of both enzymes, but their dynamic and unbinding behavior was unclear. Here, we used molecular dynamics (MD) and metadynamics (meta-D) simulations to study two of these virtual candidates (herein called IQ2 and IQ3). Their binding mode was predicted to be similar to that observed in LQ5 and BAX, well-known back-pocket binders of CDK2 and VEGFR2, respectively, including H-bonding with critical residues such as Leu83/Cys113 and Asp145/Asp190 (but excepting H-bonding with Glu51/Glu111) in CDK2/VEGFR2, correspondingly. Likewise, while LQ5 and BAX unbound through the allosteric channel as expected for type-IIA inhibitors, IQ2 and IQ3 unbound via the ATP channel (except for CDK2-IQ2) as expected for type-I½A inhibitors. Interestingly, a C-C single/double bond difference between IQ2/IQ3, respectively, resulted associated with differences in the AS/T loop flexibility observed for CDK2. These insights will help developing scaffold modifications during an optimization stage, serving as a starting point to develop dual kinase inhibitors in challenging biological targets with a promising anticancer potential.Communicated by Ramaswamy H. Sarma.
Assuntos
Simulação de Dinâmica Molecular , Quinase 2 Dependente de Ciclina/química , Ligação Proteica , Sítios de LigaçãoRESUMO
Since the beginning of oil exploration, whole ecosystems have been affected by accidents and bad practices involving petroleum compounds. In this sense, bioremediation stands out as the cheapest and most eco-friendly alternatives to reverse the damage done in oil-impacted areas. However, more efforts must be made to engineer enzymes that could be used in the bioremediation process. Interestingly, a recent work described that α-amylase, one of the most evolutionary conserved enzymes, was able to promiscuously degrade n-alkanes, a class of molecules abundant in the petroleum admixture. Considering that α-amylase is expressed in almost all known organisms, and employed in numerous biotechnological processes, using it can be a great leap toward more efficient applications of enzyme or microorganism-consortia bioremediation approaches. In this work, we employed a strict computational approach to design new α-amylase mutants with potentially enhanced catalytic efficiency toward n-alkanes. Using in silico techniques, such as molecular docking, molecular dynamics, metadynamics, and residue-residue interaction networks, we generated mutants potentially more efficient for degrading n-alkanes, L183Y, and N314A. Our results indicate that the new mutants have an increased binding rate for tetradecane, the longest n-alkane previously tested, which can reside in the catalytic center for more extended periods. Additionally, molecular dynamics and network analysis showed that the new mutations have no negative impact on protein structure than the WT. Our results aid in solidifying this enzyme as one more tool in the petroleum bioremediation toolbox.
Assuntos
Alcanos/metabolismo , Simulação de Acoplamento Molecular , alfa-Amilases/metabolismo , Alcanos/química , Bacillus subtilis/enzimologia , Biocatálise , Biodegradação Ambiental , alfa-Amilases/química , alfa-Amilases/genéticaRESUMO
Flexibility is involved in a wide range of biological processes, such as protein assembly and binding recognition. EPSP synthase is an enzyme that must undergo a large conformational change to accommodate its ligands into its binding cavity. However, although the structure of EPSP synthase has been determined, its plasticity has not been explored in depth. Therefore, in this work, we extensively examined the influence of the flexibility of Mycobacterium tuberculosis EPSP (MtEPSP) synthase on the function of this protein using classical and replica-exchange metadynamics simulations. We were able to identify five well-populated conformational clusters for MtEPSP synthase: two corresponding to open, one to ajar, and two to closed conformations. We also pinpointed three hydrophobic regions that are responsible for guiding transitions among these states. Taken together, the new findings presented here indicate how the hydrophobic regions modulate the flexibility of MtEPSP synthase, and they highlight the importance of considering these dynamic features in drug design projects employing this enzyme as a target. Graphical abstract The flexibility of EPSP synthase as a function of the pincer angles.
Assuntos
3-Fosfoshikimato 1-Carboxiviniltransferase/química , Proteínas de Bactérias/química , Mycobacterium tuberculosis/enzimologia , Domínios Proteicos , Relação Estrutura-AtividadeRESUMO
New strategies to control Leishmania disease demand an extensive knowledge about several aspects of infection including the understanding of its molecular events. In murine models, cysteine proteinase B from Leishmania amazonensis promotes regulation of immune response, and fragments from its C-terminus extension (cyspep) can play a decisive role in the host-parasite interaction. The interaction between cyspep-derived peptides and major histocompatibility complex (MHC) proteins is a crucial factor in Leishmania infections. Seven cyspep-derived peptides, previously identified as capable of interacting with H-2 (murine) MHC class I proteins, were studied in this work. We established a protocol to simulate the unbinding of these peptides from the cleft of H-2 receptors. From the simulations, we estimated the corresponding free energy of dissociation (ΔGd ) and described the molecular events that occur during the exit of peptides from the cleft. To test the reliability of this method, we first applied it to a calibration set of four crystallographic MHC/peptide complexes. Next, we explored the unbinding of the seven complexes mentioned above. Results were consistent with ΔGd values obtained from surface plasmon resonance (SPR) experiments. We also identified some of the primary interactions between peptides and H-2 receptors, and we detected three regions of influence for the interaction. This pattern was systematically observed for the peptides and helped determine a minimum distance for the real interaction between peptides and H-2 proteins occurring at â¼ 25 Å.