Pathways and mechanisms for product release in the engineered haloalkane dehalogenases explored using classical and random acceleration molecular dynamics simulations

Reinforcing Tunnel Network Exploration in Proteins Using Gaussian Accelerated Molecular Dynamics

Tunnels are structural conduits in biomolecules responsible for transporting chemical compounds and solvent molecules from the active site. They have been shown to be present in a wide variety of enzymes across all functional and structural classes. However, the study of such pathways is experimentally challenging, because they are typically transient. Computational methods, such as molecular dynamics (MD) simulations, have been successfully proposed to explore tunnels. Conventional MD (cMD) provides structural details to characterize tunnels but suffers from sampling limitations to capture rare tunnel openings on longer time scales. Therefore, in this study, we explored the potential of Gaussian accelerated MD (GaMD) simulations to improve the exploration of complex tunnel networks in enzymes. We used the haloalkane dehalogenase LinB and its two variants with engineered transport pathways, which are not only well-known for their application potential but have also been extensively studied experimentally and computationally regarding their tunnel networks and their importance in multistep catalytic reactions. Our study demonstrates that GaMD efficiently improves tunnel sampling and allows the identification of all known tunnels for LinB and its two mutants. Furthermore, the improved sampling provided insight into a previously unknown transient side tunnel (ST). The extensive conformational landscape explored by GaMD simulations allowed us to investigate in detail the mechanism of ST opening. We determined variant-specific dynamic properties of ST opening, which were previously inaccessible due to limited sampling of cMD. Our comprehensive analysis supports multiple indicators of the functional relevance of the ST, emphasizing its potential significance beyond structural considerations. In conclusion, our research proves that the GaMD method can overcome the sampling limitations of cMD for the effective study of tunnels in enzymes, providing further means for identifying rare tunnels in enzymes with the potential for drug development, precision medicine, and rational protein engineering.

Water will Find Its Way: Transport through Narrow Tunnels in Hydrolases

An aqueous environment is vital for life as we know it, and water is essential for nearly all biochemical processes at the molecular level. Proteins utilize water molecules in various ways. Consequently, proteins must transport water molecules across their internal network of tunnels to reach the desired action sites, either within them or by functioning as molecular pipes to control cellular osmotic pressure. Despite water playing a crucial role in enzymatic activity and stability, its transport has been largely overlooked, with studies primarily focusing on water transport across membrane proteins. The transport of molecules through a protein's tunnel network is challenging to study experimentally, making molecular dynamics simulations the most popular approach for investigating such events. In this study, we focused on the transport of water molecules across three different α/β-hydrolases: haloalkane dehalogenase, epoxide hydrolase, and lipase. Using a 5 μs adaptive simulation per system, we observed that only a few tunnels were responsible for the majority of water transport in dehalogenase, in contrast to a higher diversity of tunnels in other enzymes. Interestingly, water molecules could traverse narrow tunnels with subangstrom bottlenecks, which is surprising given the commonly accepted water molecule radius of 1.4 Å. Our analysis of the transport events in such narrow tunnels revealed a markedly increased number of hydrogen bonds formed between the water molecules and protein, likely compensating for the steric penalty of the process. Overall, these commonly disregarded narrow tunnels accounted for ∼20% of the total water transport observed, emphasizing the need to surpass the standard geometrical limits on the functional tunnels to properly account for the relevant transport processes. Finally, we demonstrated how the obtained insights could be applied to explain the differences in a mutant of the human soluble epoxide hydrolase associated with a higher incidence of ischemic stroke.

Decoding Inhibitor Egression from Wild-Type and G2019S Mutant LRRK2 Kinase: Insights into Unbinding Mechanisms for Precision Drug Design in Parkinson's Disease

Leucine-rich repeat kinase 2 (LRRK2) remains a viable target for drug development since the discovery of the association of its mutations with Parkinson's disease (PD). G2019S (in the kinase domain) is the most common mutation for LRRK2-based PD. Though various types of inhibitors have been developed for the kinase domain to reduce the effect of the mutation, understanding the working of these inhibitors at the molecular level is still ongoing. This study focused on the exploration of the dissociation mechanism (pathways) of inhibitors from (WT and G2019S) LRRK2 kinase (using homology model CHK1 kinase), which is one of the crucial aspects in drug discovery. Here, two ATP-competitive type I inhibitors, PF-06447475 and MLi-2 (Comp1 and Comp2 ), and one non-ATP-competitive type II inhibitor, rebastinib (Comp3), were considered for this investigation. To study the unbinding process, random accelerated molecular dynamics simulations were performed. The binding free energies of the three inhibitors for different egression paths were determined using umbrella sampling. This work found four major egression pathways that were adopted by the inhibitors Comp1 (path1, path2, and path3), Comp2 (path1, path2 and path3), and Comp3 (path3 and path4). Also, the mechanism of unbinding for each path and key residues involved in unbinding were explored. Mutation was not observed to impact the preference of the particular egression pathways for both LRRK2-Comp1 and -Comp2 systems. However, the findings suggested that the size of the inhibitor molecules might have an effect on the preference of the egression pathways. The binding energy and residence time of the inhibitors followed a similar trend to experimental observations. The findings of this work might provide insight into designing more potent inhibitors for the G2019S LRRK2 kinase.

Structure-Guided Steric Hindrance Engineering of Devosia Strain A6-243 Quinone-Dependent Dehydrogenase to Enhance Its Catalytic Efficiency

Deoxynivalenol (DON), the most widely distributed mycotoxin worldwide, causes severe health risks for humans and animals. Quinone-dependent dehydrogenase derived from Devosia strain A6-243 (DADH) can degrade DON into less toxic 3-keto-DON and then aldo-keto reductase AKR13B3 can reduce 3-keto-DON into relatively nontoxic 3-epi-DON. However, the poor catalytic efficiency of DADH made it unsuitable for practical applications, and it has become the rate-limiting step of the two-step enzymatic cascade catalysis. Here, structure-guided steric hindrance engineering was employed to enhance the catalytic efficiency of DADH. After the steric hindrance engineering, the best mutant, V429G/N431V/T432V/L434V/F537A (M5-1), showed an 18.17-fold increase in specific activity and an 11.04-fold increase in catalytic efficiency (kcat/Km) compared with that of wild-type DADH. Structure-based computational analysis provided information on the increased catalytic efficiency in the directions that attenuated steric hindrance, which was attributed to the reshaped substrate-binding pocket with an expanded catalytic binding cavity and a favorable attack distance. Tunnel analysis suggested that reshaping the active cavity by mutation might alter the shape and size of the enzyme tunnels or form one new enzyme tunnel, which might contribute to the improved catalytic efficiency of M5-1. These findings provide a promising strategy to enhance the catalytic efficiency by steric hindrance engineering.

Structural analysis of a hormone-bound Striga strigolactone receptor

Strigolactones (SLs) regulate many aspects of plant development, but ambiguities remain about how this hormone is perceived because SL-complexed receptor structures do not exist. We find that when SL binds the Striga receptor, ShHTL5, a series of conformational changes relative to the unbound state occur, but these events are not sufficient for signalling. Ligand-complexed receptors, however, form internal tunnels that posit an explanation for how SL exits its receptor after hydrolysis.

Prediction of adsorption capacity and biodegradability of polybrominated diphenyl ethers in soil

Polybrominated diphenyl ethers (PBDEs) are widely used brominated flame retardants with strong toxicity concerns. Understanding the behaviors of PBDEs in soil is essential to evaluate their environmental impact. However, the limited, incoherent, and inaccurate data has challenged predicting the adsorption capacity and biodegradability of all 209 PBDE congeners in the soil. Moreover, there are minimal studies regarding the interactions between adsorption and biodegradation behaviors of PBDEs in the soil. Herein, in this study, we adopted quantitative structure-property relationship (QSAR) modeling to predict the adsorption behavior of 209 PBDE congeners by estimating their organic carbon-water partition coefficient (K_OC) values. In addition, the biodegradability of commonly occurring PBDE congeners was evaluated by analyzing their affinity to extracellular enzymes responsible for biodegradation using molecular docking. The results highlight that the degree of bromination plays a significant role in both the absorption and biodegradation of PBDEs in the soil due to compound stability and molecular geometry. Our findings help to advance the knowledge on PBDE behaviors in the soil and facilitate PBDE remediation associated with a soil environment.

Mechanism-Based Strategy for Optimizing HaloTag Protein Labeling

HaloTag labeling technology has introduced unrivaled potential in protein chemistry and molecular and cellular biology. A wide variety of ligands have been developed to meet the specific needs of diverse applications, but only a single protein tag, DhaAHT, is routinely used for their incorporation. Following a systematic kinetic and computational analysis of different reporters, a tetramethylrhodamine- and three 4-stilbazolium-based fluorescent ligands, we showed that the mechanism of incorporating different ligands depends both on the binding step and the efficiency of the chemical reaction. By studying the different haloalkane dehalogenases DhaA, LinB, and DmmA, we found that the architecture of the access tunnels is critical for the kinetics of both steps and the ligand specificity. We showed that highly efficient labeling with specific ligands is achievable with natural dehalogenases. We propose a simple protocol for selecting the optimal protein tag for a specific ligand from the wide pool of available enzymes with diverse access tunnel architectures. The application of this protocol eliminates the need for expensive and laborious protein engineering.

Structure-guided steric hindrance engineering of Bacillus badius phenylalanine dehydrogenase for efficient L-homophenylalanine synthesis

BACKGROUND

Direct reductive amination of prochiral 2-oxo-4-phenylbutyric acid (2-OPBA) catalyzed by phenylalanine dehydrogenase (PheDH) is highly attractive in the synthesis of the pharmaceutical chiral building block L-homophenylalanine (L-HPA) given that its sole expense is ammonia and that water is the only byproduct. Current issues in this field include a poor catalytic efficiency and a low substrate loading.

RESULTS

In this study, we report a structure-guided steric hindrance engineering of PheDH from Bacillus badius to create an enhanced biocatalyst for efficient L-HPA synthesis. Mutagenesis libraries based on molecular docking, double-proximity filtering, and a degenerate codon significantly increased catalytic efficiency. Seven superior mutants were acquired, and the optimal triple-site mutant, V309G/L306V/V144G, showed a 12.7-fold higher k_cat value, and accordingly a 12.9-fold higher k_cat/K_m value, than that of the wild type. A paired reaction system comprising V309G/L306V/V144G and glucose dehydrogenase converted 1.08 M 2-OPBA to L-HPA in 210 min, and the specific space-time conversion was 30.9 mmol g^-1 L^-1 h^-1. The substrate loading and specific space-time conversion are the highest values to date. Docking simulation revealed increases in substrate-binding volume and additional degrees of freedom of the substrate 2-OPBA in the pocket. Tunnel analysis suggested the formation of new enzyme tunnels and the expansion of existing ones.

CONCLUSIONS

Overall, the results show that the mutant V309G/L306V/V144G has the potential for the industrial synthesis of L-HPA. The modified steric hindrance engineering approach can be a valuable addition to the current enzyme engineering toolbox.

An updated structural model of the A domain of the Pseudomonas putida XylR regulator poses an atypical interplay with aromatic effectors

A revised model of the aromatic binding A domain of the σ⁵⁴ -dependent regulator XylR of Pseudomonas putida mt-2 was produced based on the known 3D structures of homologous regulators PoxR, MopR and DmpR. The resulting frame was instrumental for mapping a number of mutations known to alter effector specificity, which were then reinterpreted under a dependable spatial reference. Some of these changes involved the predicted aromatic binding pocket but others occurred in distant locations, including dimerization interfaces and putative zinc binding site. The effector pocket was buried within the protein structure and accessible from the outside only through a narrow tunnel. Yet, several loop regions of the A domain could provide the flexibility required for widening such a tunnel for passage of aromatic ligands. The model was experimentally validated by treating the cells in vivo and the purified protein in vitro with benzyl bromide, which reacts with accessible nucleophilic residues on the protein surface. Structural and proteomic analyses confirmed the predicted in/out distribution of residues but also supported two additional possible scenarios of interaction of the A domain with aromatic effectors: a dynamic interaction of the fully structured yet flexible protein with the aromatic partner and/or inducer-assisted folding of the A domain.

Combination of multiple computational methods revealing specific sub-sectional recognition and hydrogen-bond dependent transportation of CKII peptide fragment in O-GlcNAc transferase

•

Mechanism of CKII peptide recognition, transportation and binding in OGT is obtained.

•

Peptide delivery is strong exothermic, highly dependent on hydrogen bond network.

•

Typical ‘spread’ & ‘V’ conformation change noticed for peptide accompanies stable OGT.

•

Specific subsection of peptide has diverse performance in its recognition and delivery.

•

Multiple methods combination may be used in other bio-system with flexible substrate.

O-linked β-N-acetyl-D-glucosamine (O-GlcNAc) transferase (OGT) is an essential enzyme in many cellular physiological catalytic reactions that regulates protein O-GlcNAcylation. Aberrant O-GlcNAcylation is related to insulin resistance, diabetic complications, cancer and neurodegenerative diseases. Understanding the peptide delivery in OGT is significant in comprehending enzymatic catalytic process, target-protein recognition and pathogenic mechanism. Herein extensive molecular dynamics (MD) simulations combined with various techniques are utilized to study the recognizing and binding mechanism of peptide fragment extracted from casein kinase II by OGT from atomic level. The residues of His496, His558, Thr633, Lys634, and Pro897 are demonstrated to play a dominant role in the peptide stabilization via hydrogen bonds and σ-π interaction, whose van der Waals and non-polar solvent effects provide the main driving force. In addition, two channels are identified. The delivery mode, mechanism together with thermodynamic and dynamic characterizations for the most favorable channel are determined. The peptide is more inclined to be recognized by OGT through the cavity comprised of residues 799–812, 893–899, and 865–871, and Tyr13-terminal is prior recognized to Met26-terminal. The transportation process is accompanied with conformation changes between the “spread” and “V” shapes. The whole process is strong exothermic that is highly dependent on the variation of hydrogen bond interactions between peptide and OGT as well as the performance of different subsections of peptide. Besides that, multiple computational methods combinations may contribute meaningfully to calculation of similar bio-systems with long and flexible substrate.

Stabilizing AqdC, a Pseudomonas Quinolone Signal-Cleaving Dioxygenase from Mycobacteria, by FRESCO-Based Protein Engineering

The mycobacterial PQS dioxygenase AqdC, a cofactor-less protein with an α/β-hydrolase fold, inactivates the virulence-associated quorum-sensing signal molecule 2-heptyl-3-hydroxy-4(1H)-quinolone (PQS) produced by the opportunistic pathogen Pseudomonas aeruginosa and is therefore a potential anti-virulence tool. We have used computational library design to predict stabilizing amino acid replacements in AqdC. While 57 out of 91 tested single substitutions throughout the protein led to stabilization, as judged by increases in T app m of >2 °C, they all impaired catalytic activity. Combining substitutions, the proteins AqdC-G40K-A134L-G220D-Y238W and AqdC-G40K-G220D-Y238W showed extended half-lives and the best trade-off between stability and activity, with increases in T app m of 11.8 and 6.1 °C and relative activities of 22 and 72 %, respectively, compared to AqdC. Molecular dynamics simulations and principal component analysis suggested that stabilized proteins are less flexible than AqdC, and the loss of catalytic activity likely correlates with an inability to effectively open the entrance to the active site.

Dual Substrate Specificity of the Rutinosidase from Aspergillus niger and the Role of Its Substrate Tunnel

Rutinosidases (α-l-rhamnopyranosyl-(1-6)-β-d-glucopyranosidases, EC 3.2.1.168, CAZy GH5) are diglycosidases that cleave the glycosidic bond between the disaccharide rutinose and the respective aglycone. Similar to many retaining glycosidases, rutinosidases can also transfer the rutinosyl moiety onto acceptors with a free -OH group (so-called transglycosylation). The recombinant rutinosidase from Aspergillus niger (AnRut) is selectively produced in Pichia pastoris. It can catalyze transglycosylation reactions as an unpurified preparation directly from cultivation. This enzyme exhibits catalytic activity towards two substrates; in addition to rutinosidase activity, it also exhibits β-d-glucopyranosidase activity. As a result, new compounds are formed by β-glucosylation or rutinosylation of acceptors such as alcohols or strong inorganic nucleophiles (NaN₃). Transglycosylation products with aliphatic aglycones are resistant towards cleavage by rutinosidase, therefore, their side hydrolysis does not occur, allowing higher transglycosylation yields. Fourteen compounds were synthesized by glucosylation or rutinosylation of selected acceptors. The products were isolated and structurally characterized. Interactions between the transglycosylation products and the recombinant AnRut were analyzed by molecular modeling. We revealed the role of a substrate tunnel in the structure of AnRut, which explained the unusual catalytic properties of this glycosidase and its specific transglycosylation potential. AnRut is attractive for biosynthetic applications, especially for the use of inexpensive substrates (rutin and isoquercitrin).

Membrane-Dependent Binding and Entry Mechanism of Dopamine into Its Receptor

Synaptic neurotransmission has recently been proposed to function via either a membrane-independent or a membrane-dependent mechanism, depending on the neurotransmitter type. In the membrane-dependent mechanism, amphipathic neurotransmitters first partition to the lipid headgroup region and then diffuse along the membrane plane to their membrane-buried receptors. However, to date, this mechanism has not been demonstrated for any neurotransmitter-receptor complex. Here, we combined isothermal calorimetry measurements with a diverse set of molecular dynamics simulation methods to investigate the partitioning of an amphipathic neurotransmitter (dopamine) and the mechanism of its entry into the ligand-binding site. Our results show that the binding of dopamine to its receptor is consistent with the membrane-dependent binding and entry mechanism. Both experimental and simulation results showed that dopamine favors binding to lipid membranes especially in the headgroup region. Moreover, our simulations revealed a ligand-entry pathway from the membrane to the binding site. This pathway passes through a lateral gate between transmembrane alpha-helices 5 and 6 on the membrane-facing side of the protein. All in all, our results demonstrate that dopamine binds to its receptor by a membrane-dependent mechanism, and this is complemented by the more traditional binding mechanism directly through the aqueous phase. The results suggest that the membrane-dependent mechanism is common in other synaptic receptors, too.

Structures of hyperstable ancestral haloalkane dehalogenases show restricted conformational dynamics

Ancestral sequence reconstruction is a powerful method for inferring ancestors of modern enzymes and for studying structure-function relationships of enzymes. We have previously applied this approach to haloalkane dehalogenases (HLDs) from the subfamily HLD-II and obtained thermodynamically highly stabilized enzymes (ΔT m up to 24 °C), showing improved catalytic properties. Here we combined crystallographic structural analysis and computational molecular dynamics simulations to gain insight into the mechanisms by which ancestral HLDs became more robust enzymes with novel catalytic properties. Reconstructed ancestors exhibited similar structure topology as their descendants with the exception of a few loop deviations. Strikingly, molecular dynamics simulations revealed restricted conformational dynamics of ancestral enzymes, which prefer a single state, in contrast to modern enzymes adopting two different conformational states. The restricted dynamics can potentially be linked to their exceptional stabilization. The study provides molecular insights into protein stabilization due to ancestral sequence reconstruction, which is becoming a widely used approach for obtaining robust protein catalysts.

The acid-base-nucleophile catalytic triad in ABH-fold enzymes is coordinated by a set of structural elements

The alpha/beta-Hydrolases (ABH) are a structural class of proteins that are found widespread in nature and includes enzymes that can catalyze various reactions in different substrates. The catalytic versatility of the ABH fold enzymes, which has been a valuable property in protein engineering applications, is based on a similar acid-base-nucleophile catalytic mechanism. In our research, we are concerned with the structure that surrounds the key units of the catalytic machinery, and we have previously found conserved structural organizations that coordinate the catalytic acid, the catalytic nucleophile and the residues of the oxyanion hole. Here, we explore the architecture that surrounds the catalytic histidine at the active sites of enzymes from 40 ABH fold families, where we have identified six conserved interactions that coordinate the catalytic histidine next to the catalytic acid and the catalytic nucleophile. Specifically, the catalytic nucleophile is coordinated next to the catalytic histidine by two weak hydrogen bonds, while the catalytic acid is directly involved in the coordination of the catalytic histidine through by two weak hydrogen bonds. The imidazole ring of the catalytic histidine is coordinated by a CH-π contact and a hydrophobic interaction. Moreover, the catalytic triad residues are connected with a residue that is located at the core of the active site of ABH fold, which is suggested to be the fourth member of a “structural catalytic tetrad”. Besides their role in the stability of the catalytic mechanism, the conserved elements of the catalytic site are actively involved in ligand binding and affect other properties of the catalytic activity, such as substrate specificity, enantioselectivity, pH optimum and thermostability of ABH fold enzymes. These properties are regularly targeted in protein engineering applications, and thus, the identified conserved structural elements can serve as potential modification sites in order to develop ABH fold enzymes with altered activities.

Molecular Cavity Topological Representation for Pattern Analysis: A NLP Analogy-Based Word2Vec Method

Cavity analysis in molecular dynamics is important for understanding molecular function. However, analyzing the dynamic pattern of molecular cavities remains a difficult task. In this paper, we propose a novel method to topologically represent molecular cavities by vectorization. First, a characterization of cavities is established through Word2Vec model, based on an analogy between the cavities and natural language processing (NLP) terms. Then, we use some techniques such as dimension reduction and clustering to conduct an exploratory analysis of the vectorized molecular cavity. On a real data set, we demonstrate that our approach is applicable to maintain the topological characteristics of the cavity and can find the change patterns from a large number of cavities.

Deciphering the Structural Basis of High Thermostability of Dehalogenase from Psychrophilic Bacterium Marinobacter sp. ELB17

Haloalkane dehalogenases are enzymes with a broad application potential in biocatalysis, bioremediation, biosensing and cell imaging. The new haloalkane dehalogenase DmxA originating from the psychrophilic bacterium Marinobacter sp. ELB17 surprisingly possesses the highest thermal stability (apparent melting temperature Tm,app = 65.9 °C) of all biochemically characterized wild type haloalkane dehalogenases belonging to subfamily II. The enzyme was successfully expressed and its crystal structure was solved at 1.45 Å resolution. DmxA structure contains several features distinct from known members of haloalkane dehalogenase family: (i) a unique composition of catalytic residues; (ii) a dimeric state mediated by a disulfide bridge; and (iii) narrow tunnels connecting the enzyme active site with the surrounding solvent. The importance of narrow tunnels in such paradoxically high stability of DmxA enzyme was confirmed by computational protein design and mutagenesis experiments.

Unraveling the Unbinding Pathways of Products Formed in Catalytic Reactions Involved in SIRT1-3: A Random Acceleration Molecular Dynamics Simulation Study

Sirtuins are a family of nicotinamide adenine dinucleotide (NAD⁺)-dependent enzymes, which undergo robust deacetylase activity, resulting in the production of nicotinamide. It is well known that nicotinamide, which is one of the products, can also act as an inhibitor for further deacetylation process by forming NAD⁺ again. Hence, the removal of nicotinamide from sirtuins is a demanding process, and the mechanistic understanding of the process remains elusive. In this investigation, we have made an attempt to unravel the unbinding pathways of nicotinamide from SIRT1, SIRT2, and SIRT3 (SIRT1-3) using Random Acceleration Molecular Dynamics (RAMD) Simulations, and we have successfully identified various unbinding channels. The selectivity of the egression channel is determined by using a thorough analysis of the frequency of egression trajectories. Similarly, various inhibitors have been docked with the active sites of SIRT1-3, and their egression pathways have been investigated to understand whether they follow the same egression pathway as that of nicotinamide. The residues that are responsible for the unbinding pathways have been determined from the analysis of RAMD trajectories. From these results, it is clear that phenylalanine and histidine residues play major roles in the egression of inhibitors. Additionally, the key residues Leu, Pro, Met, Phe, Tyr, and Ile are found to control the release by acting as gateway residues. The role of these residues from different egression channels has been studied by carrying out mutations with alanine residue. This is the first report on sirtuins, which demonstrates the novel unbinding pathways for nicotinamide/inhibitors. This work provides new insights for developing more promising SIRT1-3 inhibitors.

Engineering of lysine cyclodeaminase conformational dynamics for relieving substrate and product inhibitions in the biosynthesis of l-pipecolic acid

Efficient biocatalytic process construction by relieving substrate and product inhibitions via identification and engineering of enzyme conformational dynamics.

Uridine diphosphate release mechanism in O-N-acetylglucosamine (O-GlcNAc) transferase catalysis

O-linked N-acetylglucosamine transferase (OGT) is an essential enzyme that catalyzes the covalent bonding of N-acetylglucosamine (GlcNAc) to the hydroxyl group of a serine or threonine in the target protein. It plays an important role in many important cellular physiological catalytic reactions. Here, we determine the binding mode and the binding free energy of the OGT product (uridine diphosphate, UDP) as well as the hydrogen-bond-dependent release mechanism using extensive molecular dynamic simulations. The Lys634, Asn838, Gln839, Lys842, His901, and Asp925 residues were identified to play a major role in the UDP stabilization in the active site of OGT, where hydrogen bonding and π-π interactions mainly occur. The calculations on the mutant forms support our results. Sixteen possible release channels were identified while the two most favorable channels were determined using random acceleration molecular dynamics (RAMD) simulations combined with the constant velocity pulling (PCV) method. The thermodynamic and dynamic properties as along with the corresponding mechanism were determined and discussed according to the umbrella sampling technique. For the most optimal channel, the main free energy barrier is 13 kcal/mol, which probably originates from the hydrogen bonds between UDP and the Ala896 and Asp925 residues. Moreover, the unstable hydrogen bonds and the rollback of the ligand likely cause the other two small obstacles. This work clarifies the ligand transport mechanism in the OGT enzymatic process and is a great resource for designing inhibitors based on UDP or UDP-GlcNAc.

Insights into product release dynamics through structural analyses of thymidylate kinase

Several studies on enzyme catalysis have pointed out that the product release event could be a rate limiting step. In this study, we have compared the release event of two products, Adenosine di-phosphate (ADP) and Thymidine di-phosphate (TDP) from the active-site of human and Thermus thermophilus thymidine mono-phosphate kinase (TMPK), referred to as hTMPK and ttTMPK, respectively. TMPK catalyses the conversion of Thymidine mono-phosphate (TMP) to TDP using ATP as phosphoryl donor in the presence of Mg²⁺ ion. Most of the earlier studies on this enzyme have focused on understanding substrate binding and catalysis, but the critical product release event remains elusive. Competitive binding experiments of the substrates and the products using ttTMPK apo crystals have indicated that the substrate (TMP) can replace the bound product (TDP), even in the presence of an ADP molecule. Further, the existing random accelerated molecular dynamics (RAMD) simulation program was modified to study the release of both the products simultaneously from the active site. The RAMD simulations on product-bound structures of both ttTMPK and hTMPK, revealed that while several exit patterns of the products are permissible, the sequential exit mode is the most preferred pattern for both ttTMPK and hTMPK enzymes. Additionally, the product release from the hTMPK was found to be faster and more directional as compared to ttTMPK. Structural investigation revealed that the critical changes in the residue composition in the LID-region of ttTMPK and hTMPK have an effect on the product release and can be attributed to the observed differences during product release event. Understanding of these dissimilarities is of considerable utility in designing potent inhibitors or prodrugs that can distinguish between eukaryotic and prokaryotic homologues of thymidylate kinase.

Dynamic perspectives into the mechanisms of mutation-induced p53-DNA binding loss and inactivation using active perturbation theory: Structural and molecular insights toward the design of potent reactivators in cancer therapy

The DNA-binding ability of p53 represents the crux of its tumor suppressive activities, which involves transcriptional activation of target genes responsible for apoptosis and cell-cycle arrest. Mutational occurrences within or in close proximity to the DNA-binding surface of p53 have accounted for the loss of direct DNA-binding ability and inactivation implicated in many cases of cancer. Moreover, the design of therapeutic compounds that can restore DNA-binding ability in p53 mutants has been identified as a way forward in curtailing their oncogenic activities. However, there is still the need for more insights into evaluate the perturbations that occur at the DNA-binding interface of mp53 relative to DNA-binding loss, inactivation, and design of potent reactivators, hence the purpose of this study. Therefore, we evaluated p53-structural (R175H) and contact (R273C) mutational effects using tunnel perturbation analysis and other computational tools. We identified significant perturbations in the active tunnels of p53, which resulted in altered geometry and loss, unlike in the wild-type p53. This corroborated with structural, DNA-binding, and interaction network analysis, which showed that loss of flexibility, repulsion of DNA-interactive residues, and instability occurred at the binding interface of both mutants. Also, these mutations altered bonding interactions and network topology at the DNA-binding interface, resulting in the reduction of p53-DNA binding proximity and affinity. Therefore, these findings would aid the structure-based design of novel chemical entities capable of restoring p53-DNA binding and activation.

Tracking the route of molecular oxygen in O2-tolerant membrane-bound [NiFe] hydrogenase

[NiFe] hydrogenases catalyze the reversible splitting of H₂ into protons and electrons at a deeply buried active site. The catalytic center can be accessed by gas molecules through a hydrophobic tunnel network. While most [NiFe] hydrogenases are inactivated by O₂, a small subgroup, including the membrane-bound [NiFe] hydrogenase (MBH) of Ralstonia eutropha, is able to overcome aerobic inactivation by catalytic reduction of O₂ to water. This O₂ tolerance relies on a special [4Fe3S] cluster that is capable of releasing two electrons upon O₂ attack. Here, the O₂ accessibility of the MBH gas tunnel network has been probed experimentally using a "soak-and-freeze" derivatization method, accompanied by protein X-ray crystallography and computational studies. This combined approach revealed several sites of O₂ molecules within a hydrophobic tunnel network leading, via two tunnel entrances, to the catalytic center of MBH. The corresponding site occupancies were related to the O₂ concentrations used for MBH crystal derivatization. The examination of the O₂-derivatized data furthermore uncovered two unexpected structural alterations at the [4Fe3S] cluster, which might be related to the O₂ tolerance of the enzyme.

Reactive oxygen species leave a damage trail that reveals water channels in Photosystem II

Photosystem II (PSII), a unique membrane-bound oxidoreductase, catalyzes light-driven oxidation of water to molecular oxygen. Although high-resolution structures of PSII are known, the exact path of the substrate water molecules to the catalytic Mn₄CaO₅ center within the PSII complex remains poorly understood. PSII produces reactive oxygen species (ROS), responsible for the frequent damage and turnover of this megacomplex that occur under physiological conditions. Such ROS are known to specifically modify PSII proteins. Using high-resolution tandem mass spectrometry, we identified oxidative modifications on 36 amino acid residues on the lumenal side of PSII, in the core PSII proteins D1, D2, and CP43 of the cyanobacterium Synechocystis sp. PCC 6803. Remarkably, these oxidized residues clustered into three nearly continuous formations, tracking the pathways of ROS diffusion from the manganese center all the way out to the surface of PSII. We suggest that these profiles of oxidized residues reveal the locations of water channels within PSII. Our results provide the most comprehensive experimental evidence to date of physiologically relevant oxidized residues in PSII and illuminate three possible channels for water between the catalytic Mn cluster in the PSII complex and the bulk medium around it.

Catalytic Cycle of Haloalkane Dehalogenases Toward Unnatural Substrates Explored by Computational Modeling

The anthropogenic toxic compound 1,2,3-trichloropropane is poorly degradable by natural enzymes. We have previously constructed the haloalkane dehalogenase DhaA31 by focused directed evolution ( Pavlova, M. et al. Nat. Chem. Biol. 2009 , 5 , 727 - 733 ), which is 32 times more active than the wild-type enzyme and is currently the most active variant known against that substrate. Recent evidence has shown that the structural basis responsible for the higher activity of DhaA31 was poorly understood. Here we have undertaken a comprehensive computational study of the main steps involved in the biocatalytic hydrolysis of 1,2,3-trichloropropane to decipher the structural basis for such enhancements. Using molecular dynamics and quantum mechanics approaches we have surveyed (i) the substrate binding, (ii) the formation of the reactive complex, (iii) the chemical step, and (iv) the release of the products. We showed that the binding of the substrate and its transport through the molecular tunnel to the active site is a relatively fast process. The cleavage of the carbon-halogen bond was previously identified as the rate-limiting step in the wild-type. Here we demonstrate that this step was enhanced in DhaA31 due to a significantly higher number of reactive configurations of the substrate and a decrease of the energy barrier to the S_N2 reaction. C176Y and V245F were identified as the key mutations responsible for most of those improvements. The release of the alcohol product was found to be the rate-limiting step in DhaA31 primarily due to the C176Y mutation. Mutational dissection of DhaA31 and kinetic analysis of the intermediate mutants confirmed the theoretical observations. Overall, our comprehensive computational approach has unveiled mechanistic details of the catalytic cycle which will enable a balanced design of more efficient enzymes. This approach is applicable to deepen the biochemical knowledge of a large number of other systems and may contribute to robust strategies in the development of new biocatalysts.

Ligand diffusion in proteins via enhanced sampling in molecular dynamics

Computational simulations in biophysics describe the dynamics and functions of biological macromolecules at the atomic level. Among motions particularly important for life are the transport processes in heterogeneous media. The process of ligand diffusion inside proteins is an example of a complex rare event that can be modeled using molecular dynamics simulations. The study of physical interactions between a ligand and its biological target is of paramount importance for the design of novel drugs and enzymes. Unfortunately, the process of ligand diffusion is difficult to study experimentally. The need for identifying the ligand egress pathways and understanding how ligands migrate through protein tunnels has spurred the development of several methodological approaches to this problem. The complex topology of protein channels and the transient nature of the ligand passage pose difficulties in the modeling of the ligand entry/escape pathways by canonical molecular dynamics simulations. In this review, we report a methodology involving a reconstruction of the ligand diffusion reaction coordinates and the free-energy profiles along these reaction coordinates using enhanced sampling of conformational space. We illustrate the above methods on several ligand-protein systems, including cytochromes and G-protein-coupled receptors. The methods are general and may be adopted to other transport processes in living matter.

Molecular basis of P450 OleTJE: an investigation of substrate binding mechanism and major pathways

Cytochrome P450 OleT_JE has attracted much attention for its ability to catalyze the decarboxylation of long chain fatty acids to generate alkenes, which are not only biofuel molecule, but also can be used broadly for making lubricants, polymers and detergents. In this study, the molecular basis of the binding mechanism of P450 OleT_JE for arachidic acid, myristic acid, and caprylic acid was investigated by utilizing conventional molecular dynamics simulation and binding free energy calculations. Moreover, random acceleration molecular dynamics (RAMD) simulations were performed to uncover the most probable access/egress channels for different fatty acids. The predicted binding free energy shows an order of arachidic acid < myristic acid < caprylic acid. Key residues interacting with three substrates and residues specifically binding to one of them were identified. The RAMD results suggest the most likely channel for arachidic acid, myristic acid, and caprylic acid are 2e/2b, 2a and 2f/2a, respectively. It is suggested that the reaction is easier to carry out in myristic acid bound system than those in arachidic acid and caprylic acid bound system based on the distance of Hβ atom of substrate relative to P450 OleT_JE Compound I states. This study provided novel insight to understand the substrate preference mechanism of P450 OleT_JE and valuable information for rational enzyme design for short chain fatty acid decarboxylation.

Computational tools for enzyme improvement: why everyone can - and should - use them

This review presents computational methods that experimentalists can readily use to create smart libraries for enzyme engineering and to obtain insights into protein-substrate complexes. Computational tools have the reputation of being hard to use and inaccurate compared to experimental methods in enzyme engineering, yet they are essential to probe datasets of ever-increasing size and complexity. In recent years, bioinformatics groups have made a huge leap forward in providing user-friendly interfaces and accurate algorithms for experimentalists. These methods guide efficient experimental planning and allow the enzyme engineer to rationalize time and resources. Computational tools nevertheless face challenges in the realm of transient modern technology.

Molecular Dynamics and Free Energy Simulations of Phenylacetate and CO2 Release from AMDase and Its G74C/C188S Mutant: A Possible Rationale for the Reduced Activity of the Latter

Arylmalonate decarboxylase (AMDase) catalyzes the decarboxylation of α-aryl-α-methyl malonates to produce optically pure α-arylpropionates of industrial and medicinal importance. Herein, atomistic molecular dynamics simulations have been carried out to delineate the mechanism of the release of product molecules phenylacetate (PAC) and carbon dioxide (CO₂), from the wild-type (WT) and its G74C/C188S mutant enzymes. Both of the product molecules follow a crystallographically characterized solvent-accessible channel to come out of the protein interior. A higher free energy barrier for the release of PAC from G74C/C188S compared to that in the WT is consistent with the experimentally observed compromised efficiency of the mutant. The release of CO₂ precedes that of PAC; free energy barriers for CO₂ and PAC release in the WT enzyme are calculated to be ∼1-2 and ∼23 kcal/mol, respectively. Postdecarboxylation, CO₂ moves toward a hydrophobic pocket formed by Pro 14, Leu 38, Leu 40, Leu 77, and the side chain of Tyr 48 which serves as its temporary "reservoir". CO₂ releases following a channel mainly decorated by apolar residues, unlike in the case of oxalate decarboxylase where polar residues mediate its transport.

A Hidden Active Site in the Potential Drug Target Mycobacterium tuberculosis dUTPase Is Accessible through Small Amplitude Protein Conformational Changes

dUTPases catalyze the hydrolysis of dUTP into dUMP and pyrophosphate to maintain the proper nucleotide pool for DNA metabolism. Recent evidence suggests that dUTPases may also represent a selective drug target in mycobacteria because of the crucial role of these enzymes in maintaining DNA integrity. Nucleotide-hydrolyzing enzymes typically harbor a buried ligand-binding pocket at interdomain or intersubunit clefts, facilitating proper solvent shielding for the catalyzed reaction. The mechanism by which substrate binds this hidden pocket and product is released in dUTPases is unresolved because of conflicting crystallographic and spectroscopic data. We sought to resolve this conflict by using a combination of random acceleration molecular dynamics (RAMD) methodology and structural and biochemical methods to study the dUTPase from Mycobacterium tuberculosis In particular, the RAMD approach used in this study provided invaluable insights into the nucleotide dissociation process that reconciles all previous experimental observations. Specifically, our data suggest that nucleotide binding takes place as a small stretch of amino acids transiently slides away and partially uncovers the active site. The in silico data further revealed a new dUTPase conformation on the pathway to a relatively open active site. To probe this model, we developed the Trp²¹ reporter and collected crystallographic, spectroscopic, and kinetic data that confirmed the interaction of Trp²¹ with the active site shielding C-terminal arm, suggesting that the RAMD method is effective. In summary, our computational simulations and spectroscopic results support the idea that small loop movements in dUTPase allow the shuttlingof the nucleotides between the binding pocket and the solvent.

Product Release Pathways in Human and Plasmodium falciparum Phosphoribosyltransferase

Atomistic molecular dynamics (MD) simulations coupled with the metadynamics technique were carried out to delineate the product (PPi.2Mg and IMP) release mechanisms from the active site of both human (Hs) and Plasmodium falciparum (Pf) hypoxanthine-guanine-(xanthine) phosphoribosyltransferase (HG(X)PRT). An early movement of PPi.2Mg from its binding site has been observed. The swinging motion of the Asp side chain (D134/D145) in the binding pocket facilitates the detachment of IMP, which triggers the opening of flexible loop II, the gateway to the bulk solvent. In PfHGXPRT, PPi.2Mg and IMP are seen to be released via the same path in all of the biased MD simulations. In HsHGPRT too, the product molecules follow similar routes from the active site; however, an alternate but minor escape route for PPi.2Mg has been observed in the human enzyme. Tyr 104 and Phe 186 in HsHGPRT and Tyr 116 and Phe 197 in PfHGXPRT are the key residues that mediate the release of IMP, whereas the motion of PPi.2Mg away from the reaction center is guided by the negatively charged Asp and Glu and a few positively charged residues (Lys and Arg) that line the product release channels. Mutations of a few key residues present in loop II of Trypanosoma cruzi (Tc) HGPRT have been shown to reduce the catalytic efficiency of the enzyme. Herein, in silico mutation of corresponding residues in loop II of HsHGPRT and PfHGXPRT resulted in partial opening of the flexible loop (loop II), thus exposing the active site to bulk water, which offers a rationale for the reduced catalytic activity of these two mutant enzymes. Investigations of the product release from these HsHGPRT and PfHGXPRT mutants delineate the role of these important residues in the enzymatic turnover.

Memetic algorithms for ligand expulsion from protein cavities

Ligand diffusion through a protein interior is a fundamental process governing biological signaling and enzymatic catalysis. A complex topology of channels in proteins leads often to difficulties in modeling ligand escape pathways by classical molecular dynamics simulations. In this paper, two novel memetic methods for searching the exit paths and cavity space exploration are proposed: Memory Enhanced Random Acceleration (MERA) Molecular Dynamics (MD) and Immune Algorithm (IA). In MERA, a pheromone concept is introduced to optimize an expulsion force. In IA, hybrid learning protocols are exploited to predict ligand exit paths. They are tested on three protein channels with increasing complexity: M2 muscarinic G-protein-coupled receptor, enzyme nitrile hydratase, and heme-protein cytochrome P450cam. In these cases, the memetic methods outperform simulated annealing and random acceleration molecular dynamics. The proposed algorithms are general and appropriate in all problems where an accelerated transport of an object through a network of channels is studied.

CAVER: Algorithms for Analyzing Dynamics of Tunnels in Macromolecules

The biological function of a macromolecule often requires that a small molecule or ion is transported through its structure. The transport pathway often leads through void spaces in the structure. The properties of transport pathways change significantly in time; therefore, the analysis of a trajectory from molecular dynamics rather than of a single static structure is needed for understanding the function of pathways. The identification and analysis of transport pathways are challenging because of the high complexity and diversity of macromolecular shapes, the thermal motion of their atoms, and the large amount of conformations needed to properly describe conformational space of protein structure. In this paper, we describe the principles of the CAVER 3.0 algorithms for the identification and analysis of properties of transport pathways both in static and dynamic structures. Moreover, we introduce the improved clustering solution for finding tunnels in macromolecules, which is included in the latest CAVER 3.02 version. Voronoi diagrams are used to identify potential pathways in each snapshot of a molecular dynamics trajectory and clustering is then used to find the correspondence between tunnels from different snapshots. Furthermore, the geometrical properties of pathways and their evolution in time are computed and visualized.

Enhanced sampling simulations to construct free-energy landscape of protein-partner substrate interaction

Molecular dynamics (MD) simulations using all-atom and explicit solvent models provide valuable information on the detailed behavior of protein-partner substrate binding at the atomic level. As the power of computational resources increase, MD simulations are being used more widely and easily. However, it is still difficult to investigate the thermodynamic properties of protein-partner substrate binding and protein folding with conventional MD simulations. Enhanced sampling methods have been developed to sample conformations that reflect equilibrium conditions in a more efficient manner than conventional MD simulations, thereby allowing the construction of accurate free-energy landscapes. In this review, we discuss these enhanced sampling methods using a series of case-by-case examples. In particular, we review enhanced sampling methods conforming to trivial trajectory parallelization, virtual-system coupled multicanonical MD, and adaptive lambda square dynamics. These methods have been recently developed based on the existing method of multicanonical MD simulation. Their applications are reviewed with an emphasis on describing their practical implementation. In our concluding remarks we explore extensions of the enhanced sampling methods that may allow for even more efficient sampling.

Computational study on the unbinding pathways of B-RAF inhibitors and its implication for the difference of residence time: insight from random acceleration and steered molecular dynamics simulations

Random acceleration and steered molecular dynamics simulations reveal the unbinding pathway of B-RAF inhibitors and the difference in the residence time.

AnimoAminoMiner: Exploration of Protein Tunnels and their Properties in Molecular Dynamics

In this paper we propose a novel method for the interactive exploration of protein tunnels. The basic principle of our approach is that we entirely abstract from the 3D/4D space the simulated phenomenon is embedded in. A complex 3D structure and its curvature information is represented only by a straightened tunnel centerline and its width profile. This representation focuses on a key aspect of the studied geometry and frees up graphical estate to key chemical and physical properties represented by surrounding amino acids. The method shows the detailed tunnel profile and its temporal aggregation. The profile is interactively linked with a visual overview of all amino acids which are lining the tunnel over time. In this overview, each amino acid is represented by a set of colored lines depicting the spatial and temporal impact of the amino acid on the corresponding tunnel. This representation clearly shows the importance of amino acids with respect to selected criteria. It helps the biochemists to select the candidate amino acids for mutation which changes the protein function in a desired way. The AnimoAminoMiner was designed in close cooperation with domain experts. Its usefulness is documented by their feedback and a case study, which are included.

Substrate tunnels in enzymes: structure-function relationships and computational methodology

In enzymes, the active site is the location where incoming substrates are chemically converted to products. In some enzymes, this site is deeply buried within the core of the protein, and, in order to access the active site, substrates must pass through the body of the protein via a tunnel. In many systems, these tunnels act as filters and have been found to influence both substrate specificity and catalytic mechanism. Identifying and understanding how these tunnels exert such control has been of growing interest over the past several years because of implications in fields such as protein engineering and drug design. This growing interest has spurred the development of several computational methods to identify and analyze tunnels and how ligands migrate through these tunnels. The goal of this review is to outline how tunnels influence substrate specificity and catalytic efficiency in enzymes with buried active sites and to provide a brief summary of the computational tools used to identify and evaluate these tunnels.

Molecular dynamics, monte carlo simulations, and langevin dynamics: a computational review

Macromolecular structures, such as neuraminidases, hemagglutinins, and monoclonal antibodies, are not rigid entities. Rather, they are characterised by their flexibility, which is the result of the interaction and collective motion of their constituent atoms. This conformational diversity has a significant impact on their physicochemical and biological properties. Among these are their structural stability, the transport of ions through the M2 channel, drug resistance, macromolecular docking, binding energy, and rational epitope design. To assess these properties and to calculate the associated thermodynamical observables, the conformational space must be efficiently sampled and the dynamic of the constituent atoms must be simulated. This paper presents algorithms and techniques that address the abovementioned issues. To this end, a computational review of molecular dynamics, Monte Carlo simulations, Langevin dynamics, and free energy calculation is presented. The exposition is made from first principles to promote a better understanding of the potentialities, limitations, applications, and interrelations of these computational methods.