ABSTRACT
This study characterized the binding mechanisms of the lectin cMoL (from Moringa oleifera seeds) to carbohydrates using spectroscopy and molecular dynamics (MD). The interaction with carbohydrates was studied by evaluating lectin fluorescence emission after titration with glucose or galactose (2.0-11 mM). The Stern-Volmer constant (Ksv), binding constant (Ka), Gibbs free energy (∆G), and Hill coefficient were calculated. After the urea-induced denaturation of cMoL, evaluations were performed using fluorescence spectroscopy, circular dichroism (CD), and hemagglutinating activity (HA) evaluations. The MD simulations were performed using the Amber 20 package. The decrease in Ksv revealed that cMoL interacts with carbohydrates via a static mechanism. The cMoL bound carbohydrates spontaneously (ΔG < 0) and presented a Ka on the order of 102, with high selectivity for glucose. Protein-ligand complexes were stabilized by hydrogen bonds and hydrophobic interactions. The Hill parameter (h~2) indicated that the binding occurs through the cMoL dimer. The loss of HA at urea concentrations at which the fluorescence and CD spectra indicated protein monomerization confirmed these results. The MD simulations revealed that glucose bound to the large cavity formed between the monomers. In conclusion, the biotechnological application of cMoL lectin requires specific methods or media to improve its dimeric protein structure.
Subject(s)
Molecular Dynamics Simulation , Moringa oleifera , Protein Binding , Seeds , Moringa oleifera/chemistry , Seeds/chemistry , Plant Lectins/chemistry , Protein Multimerization , Carbohydrates/chemistry , Circular Dichroism , Lectins/chemistry , Lectins/metabolism , Spectrometry, Fluorescence , Protein Conformation , Thermodynamics , Hydrogen BondingABSTRACT
Fold-switching enables metamorphic proteins to reversibly interconvert between two highly dissimilar native states to regulate their protein functions. While about 100 proteins have been identified to undergo fold-switching, unveiling the key residues behind this mechanism for each protein remains challenging. Reasoning that fold-switching in proteins is driven by dynamic changes in local energetic frustration, we combined fold-switching simulations generated using simplified structure-based models with frustration analysis to identify key residues involved in this process based on the change in the density of minimally frustrated contacts during refolding. Using this approach to analyze the fold-switch of the bacterial transcription factor RfaH, we identified 20 residues that significantly change their frustration during its fold-switch, some of which have been experimentally and computationally reported in previous works. Our approach, which we developed as an additional module for the FrustratometeR package, highlights the role of local frustration dynamics in protein fold-switching and offers a robust tool to enhance our understanding of other proteins with significant conformational shifts.
Subject(s)
Escherichia coli Proteins , Protein Folding , Escherichia coli Proteins/chemistry , Escherichia coli Proteins/metabolism , Trans-Activators/chemistry , Trans-Activators/metabolism , Trans-Activators/genetics , Molecular Dynamics Simulation , Peptide Elongation Factors/chemistry , Peptide Elongation Factors/metabolism , Models, Molecular , Protein Conformation , ThermodynamicsABSTRACT
Deep learning methods, trained on the increasing set of available protein 3D structures and sequences, have substantially impacted the protein modeling and design field. These advancements have facilitated the creation of novel proteins, or the optimization of existing ones designed for specific functions, such as binding a target protein. Despite the demonstrated potential of such approaches in designing general protein binders, their application in designing immunotherapeutics remains relatively underexplored. A relevant application is the design of T cell receptors (TCRs). Given the crucial role of T cells in mediating immune responses, redirecting these cells to tumor or infected target cells through the engineering of TCRs has shown promising results in treating diseases, especially cancer. However, the computational design of TCR interactions presents challenges for current physics-based methods, particularly due to the unique natural characteristics of these interfaces, such as low affinity and cross-reactivity. For this reason, in this study, we explored the potential of two structure-based deep learning protein design methods, ProteinMPNN and ESM-IF1, in designing fixed-backbone TCRs for binding target antigenic peptides presented by the MHC through different design scenarios. To evaluate TCR designs, we employed a comprehensive set of sequence- and structure-based metrics, highlighting the benefits of these methods in comparison to classical physics-based design methods and identifying deficiencies for improvement.
Subject(s)
Computational Biology , Deep Learning , Receptors, Antigen, T-Cell , Receptors, Antigen, T-Cell/immunology , Receptors, Antigen, T-Cell/chemistry , Receptors, Antigen, T-Cell/metabolism , Computational Biology/methods , Humans , Protein Engineering/methods , Models, Molecular , Protein Conformation , Protein BindingABSTRACT
Molecular dynamics (MD) simulations produce a substantial volume of high-dimensional data, and traditional methods for analyzing these data pose significant computational demands. Advances in MD simulation analysis combined with deep learning-based approaches have led to the understanding of specific structural changes observed in MD trajectories, including those induced by mutations. In this study, we model the trajectories resulting from MD simulations of the SARS-CoV-2 spike protein-ACE2, specifically the receptor-binding domain (RBD), as interresidue distance maps, and use deep convolutional neural networks to predict the functional impact of point mutations, related to the virus's infectivity and immunogenicity. Our model was successful in predicting mutant types that increase the affinity of the S protein for human receptors and reduce its immunogenicity, both based on MD trajectories (precision = 0.718; recall = 0.800; [Formula: see text] = 0.757; MCC = 0.488; AUC = 0.800) and their centroids. In an additional analysis, we also obtained a strong positive Pearson's correlation coefficient equal to 0.776, indicating a significant relationship between the average sigmoid probability for the MD trajectories and binding free energy (BFE) changes. Furthermore, we obtained a coefficient of determination of 0.602. Our 2D-RMSD analysis also corroborated predictions for more infectious and immune-evading mutants and revealed fluctuating regions within the receptor-binding motif (RBM), especially in the [Formula: see text] loop. This region presented a significant standard deviation for mutations that enable SARS-CoV-2 to evade the immune response, with RMSD values of 5Å in the simulation. This methodology offers an efficient alternative to identify potential strains of SARS-CoV-2, which may be potentially linked to more infectious and immune-evading mutations. Using clustering and deep learning techniques, our approach leverages information from the ensemble of MD trajectories to recognize a broad spectrum of multiple conformational patterns characteristic of mutant types. This represents a strategic advantage in identifying emerging variants, bypassing the need for long MD simulations. Furthermore, the present work tends to contribute substantially to the field of computational biology and virology, particularly to accelerate the design and optimization of new therapeutic agents and vaccines, offering a proactive stance against the constantly evolving threat of COVID-19 and potential future pandemics.
Subject(s)
Angiotensin-Converting Enzyme 2 , Deep Learning , Molecular Dynamics Simulation , SARS-CoV-2 , Spike Glycoprotein, Coronavirus , Spike Glycoprotein, Coronavirus/chemistry , Spike Glycoprotein, Coronavirus/genetics , Spike Glycoprotein, Coronavirus/metabolism , Humans , SARS-CoV-2/genetics , SARS-CoV-2/chemistry , SARS-CoV-2/metabolism , Angiotensin-Converting Enzyme 2/chemistry , Angiotensin-Converting Enzyme 2/metabolism , COVID-19/virology , Protein Binding , Protein Conformation , Mutation , Binding Sites , Protein DomainsABSTRACT
Plasmodium falciparum apical membrane antigen 1 (PfAMA1) is a surface protein found in two stages of the malaria life cycle. This is a protein involved in a reorientation movement of the parasite so that cell invasion occurs in the so-called "moving junction", relevant when the membranes of the parasite and the host are in contact. The structure of a conformational epitope of domain III of PfAMA1 in complex with the monoclonal antibody Fab F8.12.19 is experimentally known. Here, we used molecular dynamics with enhanced sampling by Hamiltonian replica exchange molecular dynamics (HREMD) to understand the effect of intermolecular interactions, conformational variability, and intrinsically disordered regions on the mechanism of antigen-antibody interaction. Clustering methods and the analysis of conformational variability were used in order to understand the influence of the presence of the partner protein in the complex. The free-state epitope accesses a broader conformational pool, including disordered conformations not seen in the bound state. The simulations suggest an extended conformational selection mechanism in which the antibody stabilizes a conformational set of the epitope existing in the free state. The stabilization of the active conformation occurs mainly through hydrogen bonds: Tyr(H33)-Asp493, His(L94)-Val510, Ser(L93)-Glu511, Tyr(H56)-Asp485, and Tyr(H35)-Asp493. The antibody has a structure with few flexible regions, and only the complementarity determining region (CDR) H3 shows greater plasticity in the presence of the epitope.
Subject(s)
Antigens, Protozoan , Membrane Proteins , Molecular Dynamics Simulation , Plasmodium falciparum , Protozoan Proteins , Protozoan Proteins/chemistry , Protozoan Proteins/immunology , Antigens, Protozoan/chemistry , Antigens, Protozoan/immunology , Membrane Proteins/chemistry , Membrane Proteins/immunology , Plasmodium falciparum/immunology , Plasmodium falciparum/chemistry , Protein Conformation , Antibodies, Monoclonal/chemistry , Antibodies, Monoclonal/immunologyABSTRACT
Dienelactone hydrolase (DLH) is one of numerous hydrolytic enzymes with an α/ß-hydrolase fold, which catalyze the hydrolysis of dienelactone to maleylacetate. The DLHs share remarkably similar tertiary structures and a conserved arrangement of catalytic residues. This study presents the crystal structure and comprehensive functional characterization of a novel thermostable DLH from the bacterium Hydrogenobacter thermophilus (HtDLH). The crystal structure of the HtDLH, solved at a resolution of about 1.67â¯Å, exhibits a canonical α/ß-hydrolase fold formed by eight ß-sheet strands in the core, with one buried α-helix and six others exposed to the solvent. The structure also confirmed the conserved catalytic triad of DHLs formed by Cys121, Asp170, and His202 residues. The HtDLH forms stable homodimers in solution. Functional studies showed that HtDLH has the expected esterase activity over esters with short carbon chains, such as p-nitrophenyl acetate, reaching optimal activity at pH 7.5 and 70⯰C. Furthermore, HtDLH maintains more than 50â¯% of its activity even after incubation at 90⯰C for 16â¯h. Interestingly, HtDLH exhibits catalytic activity towards polyethylene terephthalate (PET) monomers, including bis-1,2-hydroxyethyl terephthalate (BHET) and 1-(2-hydroxyethyl) 4-methyl terephthalate, as well as other aliphatic and aromatic esters. These findings associated with the lack of activity on amorphous PET indicate that HtDLH has characteristic of a BHET-degrading enzyme. This work expands our understanding of enzyme families involved in PET degradation, providing novel insights for plastic biorecycling through protein engineering, which could lead to eco-friendly solutions to reduce the accumulation of plastic in landfills and natural environments.
Subject(s)
Carboxylic Ester Hydrolases , Enzyme Stability , Substrate Specificity , Crystallography, X-Ray , Carboxylic Ester Hydrolases/metabolism , Carboxylic Ester Hydrolases/chemistry , Carboxylic Ester Hydrolases/genetics , Phthalic Acids/metabolism , Phthalic Acids/chemistry , Esters/metabolism , Esters/chemistry , Models, Molecular , Bacterial Proteins/metabolism , Bacterial Proteins/chemistry , Bacterial Proteins/genetics , Protein Conformation , Hydrogen-Ion Concentration , Kinetics , Hydrolysis , Catalytic Domain , TemperatureABSTRACT
To get a deeper understanding of the structural bases of the allosteric transition between T and R states of plant and bacterial phosphoenolpyruvate carboxylases (PEPCs), we obtained the first T-state crystal structures of the maize photosynthetic PEPC (ZmPEPC-C4) and exhaustively compared them with the previously reported R-state ZmPEPC-C4 and other T-state structures. We identified previously unrecognized significant conformational changes in the T state: that of the α8-α9 loop, which connects the two kinds of activator allosteric sites with the active site, the conversion of the α30 helix into a 310 helix, leading to the disorganization of the active site lid and activators allosteric sites, and the closure of the inhibitor allosteric-site lid. Additionally, we identified previously overlooked, highly conserved residues of potential interest in the allosteric transition, including two histidines whose protonation might stabilize the T state. The crystal structures reported here also suggest similar tetrameric quaternary arrangements of PEPC enzymes in the R and T states, and the location of the bicarbonate binding site, as well as the conformational changes required for the carboxylation step. Our findings and working hypothesis advance the understanding of the structural features of the allosteric PEPC enzymes and provide a foundation for future experiments.
Subject(s)
Models, Molecular , Phosphoenolpyruvate Carboxylase , Zea mays , Zea mays/enzymology , Zea mays/chemistry , Allosteric Regulation , Crystallography, X-Ray , Phosphoenolpyruvate Carboxylase/chemistry , Phosphoenolpyruvate Carboxylase/metabolism , Allosteric Site , Catalytic Domain , Protein Conformation , Amino Acid SequenceABSTRACT
Faced with the emergence of multiresistant microorganisms that affect human health, microbial agents have become a serious global threat, affecting human health and plant crops. Antimicrobial peptides have attracted significant attention in research for the development of new microbial control agents. This work's goal was the structural characterization and analysis of antifungal activity of chitin-binding peptides from Capsicum baccatum and Capsicum frutescens seeds on the growth of Candida and Fusarium species. Proteins were initially submitted to extraction in phosphate buffer pH 5.4 and subjected to chitin column chromatography. Posteriorly, two fractions were obtained for each species, Cb-F1 and Cf-F1 and Cb-F2 and Cf-F2, respectively. The Cb-F1 (C. baccatum) and Cf-F1 (C. frutescens) fractions did not bind to the chitin column. The electrophoresis results obtained after chromatography showed two major protein bands between 3.4 and 14.2 kDa for Cb-F2. For Cf-F2, three major bands were identified between 6.5 and 14.2 kDa. One band from each species was subjected to mass spectrometry, and both bands showed similarity to nonspecific lipid transfer protein. Candida albicans and Candida tropicalis had their growth inhibited by Cb-F2. Cf-F2 inhibited the development of C. albicans but did not inhibit the growth of C. tropicalis. Both fractions were unable to inhibit the growth of Fusarium species. The toxicity of the fractions was tested in vivo on Galleria mellonella larvae, and both showed a low toxicity rate at high concentrations. As a result, the fractions have enormous promise for the creation of novel antifungal compounds.
Subject(s)
Antifungal Agents , Candida , Chitin , Fusarium , Molecular Docking Simulation , Antifungal Agents/pharmacology , Antifungal Agents/chemistry , Antifungal Agents/metabolism , Chitin/chemistry , Chitin/metabolism , Fusarium/drug effects , Candida/drug effects , Carrier Proteins/chemistry , Carrier Proteins/metabolism , Animals , Capsicum/chemistry , Plant Proteins/chemistry , Plant Proteins/metabolism , Plant Proteins/pharmacology , Microbial Sensitivity Tests , Protein Binding , Protein ConformationABSTRACT
BACKGROUND: Humoral immune response against the pre-fusion (pre-F) conformation of respiratory syncytial virus (RSV) F protein has been proposed to play a protective role against infection. An RSV pre-F maternal vaccine has been recently approved in several countries to protect young infants against RSV. We aimed to assess serum IgG titers against the pre-F and post-F conformations of RSV F protein and their association with life-threatening RSV disease (LTD) in previously healthy infants. METHODS: A prospective cohort study including hospitalized infants <12 months with a first RSV infection was conducted during 2017-2019. Patients with LTD required intensive care and mechanical respiratory assistance. RSV pre-F exclusive and post-F antibody responses were determined by post-F competition and non-competition immunoassays, respectively, and neutralizing activity was measured by plaque reduction neutralization test. RESULTS: Fifty-eight patients were included; the median age was 3.5 months and 41 % were females. Fifteen patients developed LTD. RSV F-specific antibody titers positively correlated with neutralizing antibody titers in acute and convalescent phases but, importantly, they did not associate with LTD. Acute RSV pre-F exclusive and post-F IgG titers negatively correlated with patient age (P = 0.0007 and P < 0.0001), while a positive correlation was observed between the fold changes in RSV F-specific antibody titers between convalescent and acute phase and patient age (P = 0.0014 and P < 0.0001). Infants ≤2 months exhibited significantly lower fold-changes in RSV F-specific and neutralizing antibody titers between convalescence and acute phase than older infants. Additionally, acute RSV antibody titers showed no correlation with nasal RSV load and, furthermore, nasal viral load was not associated with the development of LTD. CONCLUSIONS: This study highlights that protection against life-threatening RSV disease is not necessarily antibody-dependent. Further characterization of the immune response against RSV and its role in protection against severe disease is important for the development of the safest possible preventive strategies.
Subject(s)
Antibodies, Neutralizing , Antibodies, Viral , Immunoglobulin G , Respiratory Syncytial Virus Infections , Respiratory Syncytial Virus, Human , Viral Fusion Proteins , Humans , Respiratory Syncytial Virus Infections/immunology , Respiratory Syncytial Virus Infections/prevention & control , Female , Infant , Antibodies, Viral/blood , Antibodies, Viral/immunology , Viral Fusion Proteins/immunology , Prospective Studies , Respiratory Syncytial Virus, Human/immunology , Male , Antibodies, Neutralizing/immunology , Antibodies, Neutralizing/blood , Immunoglobulin G/blood , Immunoglobulin G/immunology , Protein Conformation , Respiratory Syncytial Virus Vaccines/immunology , Infant, NewbornABSTRACT
Septins are filamentous nucleotide-binding proteins which can associate with membranes in a curvature-dependent manner leading to structural remodelling and barrier formation. Ciona intestinalis, a model for exploring the development and evolution of the chordate lineage, has only four septin-coding genes within its genome. These represent orthologues of the four classical mammalian subgroups, making it a minimalist non-redundant model for studying the modular assembly of septins into linear oligomers and thereby filamentous polymers. Here, we show that C. intestinalis septins present a similar biochemistry to their human orthologues and also provide the cryo-EM structures of an octamer, a hexamer and a tetrameric sub-complex. The octamer, which has the canonical arrangement (2-6-7-9-9-7-6-2) clearly shows an exposed NC-interface at its termini enabling copolymerization with hexamers into mixed filaments. Indeed, only combinations of septins which had CiSEPT2 occupying the terminal position were able to assemble into filaments via NC-interface association. The CiSEPT7-CiSEPT9 tetramer is the smallest septin particle to be solved by Cryo-EM to date and its good resolution (2.7 Å) provides a well-defined view of the central NC-interface. On the other hand, the CiSEPT7-CiSEPT9 G-interface shows signs of fragility permitting toggling between hexamers and octamers, similar to that seen in human septins but not in yeast. The new structures provide insights concerning the molecular mechanism for cross-talk between adjacent interfaces. This indicates that C. intestinalis may represent a valuable tool for future studies, fulfilling the requirements of a complete but simpler system to understand the mechanisms behind the assembly and dynamics of septin filaments.
Subject(s)
Ciona intestinalis , Cryoelectron Microscopy , Models, Molecular , Protein Multimerization , Septins , Ciona intestinalis/metabolism , Ciona intestinalis/chemistry , Ciona intestinalis/genetics , Septins/metabolism , Septins/chemistry , Septins/genetics , Animals , Humans , Nucleotides/metabolism , Nucleotides/chemistry , Protein Conformation , Protein BindingABSTRACT
CONTEXT: Geometrical knots are rare structural arrangements in proteins in which the polypeptide chain ties itself into a knot, which is very intriguing due to the uncertainty of their impact on the protein properties. Presently, classical molecular dynamics is the most employed technique in the few studies found on this topic, so any information on how the presence of knots affects the reactivity and electronic properties of proteins is even scarcer. Using the electronic structure methods and quantum chemical descriptors analysis, we found that the same amino-acid residues in the knot core have statistically larger values for the unknotted protein, for both hard-hard and soft-soft interaction descriptors. In addition, we present a computationally feasible protocol, where we show it is possible to separate the contribution of the geometrical knot to the reactivity and other electronic structure properties. METHODS: In order to investigate these systems, we used PRIMoRDiA, a new software developed by our research group, to explore the electronic structure of biological macromolecules. We evaluated several local quantum chemical descriptors to unveil relevant patterns potentially originating from the presence of the geometrical knot in two proteins, belonging to the ornithine transcarbamylase family. We compared several sampled structures from these two enzymes that are highly similar in both tertiary structure and function, but one of them has a knot whereas the other does not. The sampling was carried out through molecular dynamics simulations using ff14SB force field along 50 ns, and the semiempirical convergence was performed with PM7 Hamiltonian.
Subject(s)
Molecular Dynamics Simulation , Ornithine Carbamoyltransferase , Ornithine Carbamoyltransferase/chemistry , Ornithine Carbamoyltransferase/metabolism , Protein Conformation , Models, MolecularABSTRACT
Hemoglobin (Hb) is a hemeprotein found inside erythrocytes and is crucial in transporting oxygen and carbon dioxide in our bodies. In erythrocytes (Ery), the main energy source is glucose metabolized through glycolysis. However, a fraction of Hb can undergo glycation, in which a free amine group from the protein spontaneously binds to the carbonyl of glucose in the bloodstream, resulting in the formation of glycated hemoglobin (HbA1c), widely used as a marker for diabetes. Glycation leads to structural and conformational changes, compromising the function of proteins, and is intensified in the event of hyperglycemia. The main changes in Hb include structural alterations to the heme group, compromising its main function (oxygen transport). In addition, amyloid aggregates can form, which are strongly related to diabetic complications and neurodegenerative diseases. Therefore, this chapter discusses in vitro protocols for producing glycated Hb, as well as the main techniques and biophysical assays used to assess changes in the protein's structure before and after the glycation process. This more complete understanding of the effects of glycation on Hb is fundamental for understanding the complications associated with hyperglycemia and for developing more effective prevention and treatment strategies.
Subject(s)
Hemoglobins , Humans , Glycosylation , Hemoglobins/metabolism , Hemoglobins/chemistry , Glycated Hemoglobin/metabolism , Protein Conformation , AnimalsABSTRACT
This work focuses on the δ receptor (DOR), a G protein-coupled receptor (GPCR) belonging to the opioid receptor group. DOR is expressed in numerous tissues, particularly within the nervous system. Our study explores computationally the receptor's interactions with various ligands, including opiates and opioid peptides. It elucidates how these interactions influence the δ receptor response, relevant in a wide range of health and pathological processes. Thus, our investigation aims to explore the significance of DOR as an incoming drug target for pain relief and neurodegenerative diseases and as a source for novel opioid non-narcotic analgesic alternatives. We analyze the receptor's structural properties and interactions using Molecular Dynamics (MD) simulations and Gaussian-accelerated MD across different functional states. To thoroughly assess the primary differences in the structural and conformational ensembles across our different simulated systems, we initiated our study with 1 µs of conventional Molecular Dynamics. The strategy was chosen to encompass the full activation cycle of GPCRs, as activation processes typically occur within this microsecond range. Following the cMD, we extended our study with an additional 100 ns of Gaussian accelerated Molecular Dynamics (GaMD) to enhance the sampling of conformational states. This simulation approach allowed us to capture a comprehensive range of dynamic interactions and conformational changes that are crucial for GPCR activation as influenced by different ligands. Our study includes comparing agonist and antagonist complexes to uncover the collective patterns of their functional states, regarding activation, blocking, and inactivation of DOR, starting from experimental data. In addition, we also explored interactions between agonist and antagonist molecules from opiate and opioid classifications to establish robust structure-activity relationships. These interactions have been systematically quantified using a Quantitative Structure-Activity Relationships (QSAR) model. This research significantly contributes to our understanding of this significant pharmacological target, which is emerging as an attractive subject for drug development.
Subject(s)
Molecular Dynamics Simulation , Receptors, Opioid, delta , Receptors, Opioid, delta/metabolism , Receptors, Opioid, delta/chemistry , Humans , Ligands , Analgesics, Opioid/pharmacology , Analgesics, Opioid/chemistry , Protein Binding , Protein ConformationABSTRACT
The MC1R protein is a receptor found in melanocytes that plays a role in melanin synthesis. Mutations in this protein can impact hair color, skin tone, tanning ability, and increase the risk of skin cancer. The MC1R protein is activated by the alpha-melanocyte-stimulating hormone (α-MSH). Previous studies have shown that mutations affect the interaction between MC1R and α-MSH; however, the mechanism behind this process is poorly understood. Our study aims to shed light on this mechanism using molecular dynamics (MD) simulations to analyze the Asp84Glu and Asp294His variants. We simulated both the wild-type (WT) protein and the mutants with and without ligand. Our results reveal that mutations induce unique conformations during state transitions, hindering the switch between active and inactive states and decreasing cellular levels of cAMP. Interestingly, Asp294His showed increased ligand affinity but decreased protein activity, highlighting that tighter binding does not always lead to increased activation. Our study provides insights into the molecular mechanisms underlying the impact of MC1R mutations on protein activity.
Subject(s)
Cyclic AMP , Mutation , Receptor, Melanocortin, Type 1 , alpha-MSH , Humans , alpha-MSH/chemistry , alpha-MSH/metabolism , alpha-MSH/genetics , Binding Sites , Cyclic AMP/metabolism , Cyclic AMP/chemistry , Ligands , Molecular Dynamics Simulation , Protein Binding , Protein Conformation , Receptor, Melanocortin, Type 1/genetics , Receptor, Melanocortin, Type 1/chemistry , Receptor, Melanocortin, Type 1/metabolismABSTRACT
Glycoside hydrolase family 5 (GH5) encompasses enzymes with several different activities, including endo-1,4-ß-mannosidases. These enzymes are involved in mannan degradation, and have a number of biotechnological applications, such as mannooligosaccharide prebiotics production, stain removal and dyes decolorization, to name a few. Despite the importance of GH5 enzymes, only a few members of subfamily 7 were structurally characterized. In the present work, biochemical and structural characterization of Bacillus licheniformis GH5 mannanase, BlMan5_7 were performed and the enzyme cleavage pattern was analyzed, showing that BlMan5_7 requires at least 5 occupied subsites to perform efficient hydrolysis. Additionally, crystallographic structure at 1.3 Å resolution was determined and mannoheptaose (M7) was docked into the active site to investigate the interactions between substrate and enzyme through molecular dynamic (MD) simulations, revealing the existence of a - 4 subsite, which might explain the generation of mannotetraose (M4) as an enzyme product. Biotechnological application of the enzyme in stain removal was investigated, demonstrating that BlMan5_7 addition to washing solution greatly improves mannan-based stain elimination.
Subject(s)
Bacillus licheniformis , Catalytic Domain , Mutagenesis, Site-Directed , Bacillus licheniformis/enzymology , Bacillus licheniformis/genetics , Crystallography, X-Ray , Molecular Dynamics Simulation , Mannosidases/chemistry , Mannosidases/genetics , Mannosidases/metabolism , Substrate Specificity , Hydrolysis , Tetroses/chemistry , Tetroses/metabolism , Bacterial Proteins/chemistry , Bacterial Proteins/genetics , Bacterial Proteins/metabolism , Protein Conformation , Mannans/chemistry , Mannans/metabolism , beta-Mannosidase/chemistry , beta-Mannosidase/genetics , beta-Mannosidase/metabolism , Models, Molecular , Molecular Docking Simulation , OligosaccharidesABSTRACT
We analyzed the thermal stability of the BstHPr protein through the site-directed point mutation Lys62 replaced by Ala residue using molecular dynamics simulations at five different temperatures: 298, 333, 362, 400, and 450 K, for periods of 1 µs and in triplicate. The results from the mutant thermophilic BstHPrm protein were compared with those of the wild-type thermophilic BstHPr protein and the mesophilic BsHPr protein. Structural and molecular interaction analyses show that proteins lose stability as temperature increases. Mutant and wild-type proteins behave similarly up to 362 K. However, at 400 K the mutant protein shows greater structural instability, losing more buried hydrogen bonds and exposing more of its non-polar residues to the solvent. Therefore, in this study, we confirmed that the salt bridge network of the Glu3-Lys62-Glu36 triad, made up of the Glu3-Lys62 and Glu36-Lys62 ion pairs, provides thermal stability to the thermophilic BstHPr protein.
Subject(s)
Molecular Dynamics Simulation , Protein Stability , Hydrogen Bonding , Temperature , Mutation , Bacterial Proteins/chemistry , Bacterial Proteins/genetics , Bacterial Proteins/metabolism , Amino Acid Substitution , Protein Conformation , Mutagenesis, Site-DirectedABSTRACT
The accurate experimental estimation of protein-ligand systems' residence time (τ) has become very relevant in drug design projects due to its importance in the last stages of refinement of the drug's pharmacodynamics and pharmacokinetics. It is now well-known that it is not sufficient to estimate the affinity of a protein-drug complex in the thermodynamic equilibrium process in in vitro experiments (closed systems), where the concentrations of the drug and protein remain constant. On the contrary, it is mandatory to consider the conformational dynamics of the system in terms of the binding and unbinding processes between protein and drugs in in vivo experiments (open systems), where their concentrations are in constant flux. This last model has been proven to dictate much of several drugs' pharmacological activities in vivo. At the atomistic level, molecular dynamics simulations can explain why some drugs are more effective than others or unveil the molecular aspects that make some drugs work better in one molecular target. Here, the protein kinases Aurora A/B, complexed with its inhibitor Danusertib, were studied using conventional and enhanced molecular dynamics (MD) simulations to estimate the dissociation paths and, therefore, the computational τ values and their comparison with experimental ones. Using classical molecular dynamics (cMD), three differential residues within the Aurora A/B active site, which seems to play an essential role in the observed experimental Danusertib's residence time against these kinases, were characterized. Then, using WT-MetaD, the relative Danusertib's residence times against Aurora A/B kinases were measured in a nanosecond time scale and were compared to those τ values observed experimentally. In addition, the potential dissociation paths of Danusertib in Aurora A and B were characterized, and differences that might be explained by the differential residues in the enzyme's active sites were found. In perspective, it is expected that this computational protocol can be applied to other protein-ligand complexes to understand, at the molecular level, the differences in residence times and amino acids that may contribute to it.
Subject(s)
Aurora Kinase A , Aurora Kinase B , Molecular Dynamics Simulation , Aurora Kinase B/metabolism , Aurora Kinase B/chemistry , Aurora Kinase B/antagonists & inhibitors , Aurora Kinase A/metabolism , Aurora Kinase A/chemistry , Aurora Kinase A/antagonists & inhibitors , Pyrazoles/chemistry , Pyrazoles/metabolism , Protein Conformation , Protein Kinase Inhibitors/chemistry , Protein Kinase Inhibitors/pharmacology , Protein Kinase Inhibitors/metabolism , Protein Binding , Humans , Benzamides/chemistry , Benzamides/metabolism , Benzamides/pharmacology , ThermodynamicsABSTRACT
Assessing membrane protein stability is among the major challenges in protein science due to their inherent complexity, which complicates the application of conventional biophysical tools. In this work, sodium dodecyl sulfate-induced denaturation of AfCopA, a Cu(I)-transport ATPase from Archaeoglobus fulgidus, was explored using a combined model-free spectral phasor analysis and a model-dependent thermodynamic analysis. Decrease in tryptophan and 1-anilino-naphthalene-8-sulfonate fluorescence intensity, displacements in the spectral phasor space, and the loss of ATPase activity were reversibly induced by this detergent. Refolding from the SDS-induced denatured state yields an active enzyme that is functionally and spectroscopically indistinguishable from the native state of the protein. Phasor analysis of Trp spectra allowed us to identify two intermediate states in the SDS-induced denaturation of AfCopA, a result further supported by principal component analysis. In contrast, traditional thermodynamic analysis detected only one intermediate state, and including the second one led to overparameterization. Additionally, ANS fluorescence spectral analysis detected one more intermediate and a gradual change at the level of the hydrophobic transmembrane surface of the protein. Based on this evidence, a model for acquiring the native structure of AfCopA in a membrane-like environment is proposed.
Subject(s)
Archaeoglobus fulgidus , Membrane Proteins , Protein Denaturation , Sodium Dodecyl Sulfate , Thermodynamics , Sodium Dodecyl Sulfate/chemistry , Sodium Dodecyl Sulfate/pharmacology , Archaeoglobus fulgidus/enzymology , Membrane Proteins/chemistry , Membrane Proteins/metabolism , Adenosine Triphosphatases/chemistry , Adenosine Triphosphatases/metabolism , Spectrometry, Fluorescence , Protein Stability , Archaeal Proteins/chemistry , Archaeal Proteins/metabolism , Anilino Naphthalenesulfonates/chemistry , Anilino Naphthalenesulfonates/metabolism , Tryptophan/chemistry , Tryptophan/metabolism , Copper/chemistry , Copper/metabolism , Protein Folding , Protein ConformationABSTRACT
DapE is a Zn2+-metallohydrolase recognized as a drug target for bacterial control. It is a homodimer that requires the exchange of interface strands by an induced fit essential for catalysis. Identifying novel anti-DapE agents requires greater structural details. Most of the characterized DapEs are from the Gram-negative group. Here, two high-resolution DapE crystal structures from Enterococcus faecium are presented for the first time with novel aspects. A loosened enzyme intermediate between the open and closed conformations is observed. Substrates may bind to loose state, subsequently it closes, where hydrolysis occurs, and finally, the change to the open state leads to the release of the products. Mutation of His352 suggests a role, along with His194, in the oxyanion stabilization in the mono-metalated Zn2+ isoform, while in the di-metalated isoform, the metal center 2 complements it function. An aromatic-π box potentially involved in the interaction of DapE with other proteins, and a peptide flip could determine the specificity in the Gram-positive ArgE/DapE group. Finally, details of two extra-catalytic cavities whose geometry changes depending on the conformational state of the enzyme are presented. These cavities could be a target for developing non-competitive agents that trap the enzyme in an inactive state.
Subject(s)
Bacterial Proteins , Enterococcus faecium , Amidohydrolases/chemistry , Amidohydrolases/metabolism , Amidohydrolases/genetics , Amino Acid Sequence , Bacterial Proteins/chemistry , Bacterial Proteins/genetics , Bacterial Proteins/metabolism , Catalytic Domain , Crystallography, X-Ray , Enterococcus faecium/enzymology , Ligands , Models, Molecular , Protein Binding , Protein Conformation , Substrate Specificity , Zinc/chemistry , Zinc/metabolismABSTRACT
Ubiquitin-specific protease 7 (USP7) is a deubiquitinase enzyme that plays a critical role in regulating various cellular processes by cleaving ubiquitin molecules from target proteins. The C-terminal loop (CTL) motif is a specific region at the C-terminal end of the USP7 enzyme. Recent experiments suggest that the CTL motif plays a role in USP7's catalytic activity by contributing to the enzyme's structural stability, substrate recognition, and catalytic efficiency. The objective of this work is to elucidate these roles through the utilization of computational methods for molecular simulations. For this, we conducted extensive molecular dynamics (MD) simulations to investigate the conformational dynamics and protein-protein interactions within the USP7 enzyme-substrate complex with the substrate consisting of the ubiquitin tagged with the fluorescent label rhodamine 110-gly (Ub-Rho). Our results demonstrate that the CTL motif plays a crucial role in stabilizing the Ubl domains' conformation and augmenting the stability of active conformations within the enzyme-substrate complex. Conversely, the absence of the CTL motif results in increased flexibility and variability in Ubl domains' motion, leading to a reduced percentage of active conformations. Furthermore, our analysis of protein-protein interactions highlights the significance of the CTL motif in anchoring the Ubl45 domains to the catalytic domain (CD), thereby facilitating stable interactions with the substrate. Overall, our findings provide valuable insights into the conformational dynamics and protein-protein interactions inherent in the USP7 enzyme-substrate complex. These insights shed light on some mechanistic details of USP7 concerning the substrate's recognition before its catalytic action.