Active Site Water Research Articles

Molecular, Genetic, and Biochemical Characterization of OXA-484 Carbapenemase, a Difficult-to-Detect R214G Variant of OXA-181.

OXA-244, an R214G variant of OXA-48, is silently spreading worldwide likely because of difficulties in detection using classical screening media. Here, we characterized two clinical isolates of Escherichia coli and Citrobacter youngae that displayed reduced susceptibility to carbapenems but were lacking significant carbapenemase activity as revealed by negative Carba NP test results. However, positive test results were seen for OXA-48-like enzymes by lateral flow immunoassays. WGS revealed the presence of a blaOXA-181-like gene that codes for OXA-484, an R214G variant of OXA-181. BlaOXA-484 gene was located on a 58.4-kb IncP1-like plasmid (pN-OXA-484), that upon transfer into E. coli HB4 with impaired permeability, conferred carbapenem and temocillin resistance (MICs > 32 mg/L). E. coli TOP10 (pTOPO-OXA-484) revealed reduced MICs in most substrates as compared to E. coli TOP10 (pTOPO-OXA-181), especially for imipenem (0.25 mg/L versus 0.75 mg/L) and temocillin (16 mg/L versus 1028 mg/L). Catalytic efficiencies of OXA-484 were reduced as compared to OXA-181 for most ß-lactams including imipenem and temocillin with 27.5- and 21.7-fold reduction, respectively. Molecular modeling confirmed that the salt bridges between R214, D159, and the R1 substituent's carboxylate group of temocillin were not possible with G214 in OXA-484, explaining the reduced affinity for temocillin. In addition, changes in active site's water network may explain the decrease in hydrolysis rate of carbapenems. OXA-484 has weak imipenem and temocillin hydrolytic activities, which may lead to silent spread due to underdetection using selective screening media or biochemical imipenem hydrolysis confirmatory tests.

Simulation-Guided Engineering Enables a Functional Switch in Selinadiene Synthase toward Hydroxylation.

Engineering sesquiterpene synthases to form predefined alternative products is a major challenge due to their diversity in cyclization mechanisms and our limited understanding of how amino acid changes affect the steering of these mechanisms. Here, we use a combination of atomistic simulation and site-directed mutagenesis to engineer a selina-4(15),7(11)-diene synthase (SdS) such that its final reactive carbocation is quenched by trapped active site water, resulting in the formation of a complex hydroxylated sesquiterpene (selin-7(11)-en-4-ol). Initially, the SdS G305E variant produced 20% selin-7(11)-en-4-ol. As suggested by modeling of the enzyme-carbocation complex, selin-7(11)-en-4-ol production could be further improved by varying the pH, resulting in selin-7(11)-en-4-ol becoming the major product (48%) at pH 6.0. We incorporated the SdS G305E variant along with genes from the mevalonate pathway into bacterial BL21(DE3) cells and demonstrated the production of selin-7(11)-en-4-ol at a scale of 10 mg/L in batch fermentation. These results highlight opportunities for the simulation-guided engineering of terpene synthases to produce predefined complex hydroxylated sesquiterpenes.

The β-subunit of tryptophan synthase is a latent tyrosine synthase.

Aromatic amino acids and their derivatives are diverse primary and secondary metabolites with critical roles in protein synthesis, cell structure and integrity, defense and signaling. All de novo aromatic amino acid production relies on a set of ancient and highly conserved chemistries. Here we introduce a new enzymatic transformation for L-tyrosine synthesis by demonstrating that the β-subunit of tryptophan synthase-which natively couples indole and L-serine to form L-tryptophan-can act as a latent 'tyrosine synthase'. A single substitution of a near-universally conserved catalytic residue unlocks activity toward simple phenol analogs and yields exclusive para carbon-carbon bond formation to furnish L-tyrosines. Structural and mechanistic studies show how a new active-site water molecule orients phenols for a nonnative mechanism of alkylation, with additional directed evolution resulting in a net >30,000-fold rate enhancement. This new biocatalyst can be used to efficiently prepare valuable L-tyrosine analogs at gram scales and provides the missing chemistry for a conceptually different pathway to L-tyrosine.

Structural and biochemical characterization of a nucleotide hydrolase from Streptococcus pneumonia

In this report, we structurally and biochemically characterized the unknown gene product SP1746 from Streptococcus pneumoniae serotype 4. Various crystal structures of SP1746 in the apo form and in complex with different nucleotides were determined. SP1746 is a globular protein, which belongs to the histidine-aspartate (HD) domain superfamily with two Fe3+ ions in the active site that are coordinated by key active site residues and water molecules. All nucleotides bind in a similar orientation in the active site with their phosphate groups anchored to the diiron cluster. Biochemically, SP1746 hydrolyzes different nucleotide substrates. SP1746 most effectively hydrolyzes diadenosine tetraphosphate (Ap4A) to two ADPs. Based on the aforementioned data, we annotated SP1746 as an Ap4A hydrolase, belonging to the YqeK family. Our invitro data indicate a potential role for SP1746 in regulating Ap4A homeostasis, which requires validation with invivo experiments in bacteria in the future.

Abstract 1168 Water as structure in a cytochrome P450 enzyme: CYP121 from Mycobacterium tuberculosis relies on an ordered active site water network for its structure and function

Abstract 1168 Water as structure in a cytochrome P450 enzyme: CYP121 from Mycobacterium tuberculosis relies on an ordered active site water network for its structure and function

Asparagine-85 Stabilizes a Structural Active Site Water Network in CYP121A1 of Mycobacterium tuberculosis.

The cytochrome P450 enzyme CYP121A1 endogenously catalyzes the formation of a carbon-carbon bond between the two phenol groups of dicyclotyrosine (cYY) in Mycobacterium tuberculosis (Mtb). One of 20 CYP enzymes in Mtb, CYP121A1 continues to garner significant interest as a potential drug target. The accompanying reports the use of 19F NMR spectroscopy, reconstituted activity assays, and molecular dynamics simulations to investigate the significance of hydrogen bonding interactions that were theorized to stabilize a static active site water network. The active site residue Asn-85, whose hydrogen bonds with the diketopiperazine ring of cYY contributes to a contiguous active site water network in the absence of cYY, was mutated to a serine (N85S) and to a glutamine (N85Q). These conservative changes in the hydrogen bond donor side chain result in inactivation of the enzyme. Moreover, the N85S mutation induces reverse type-I binding as measured by absorbance difference spectra. NMR spectra monitoring the ligand-adaptive FG-loop and the active site Trp-182 side chain confirm that disruption of the active site water network also significantly alters the structure of the active site. These data were consistent with dynamics simulations of N85S and N85Q that demonstrate that a compromised water network is responsible for remodeling of the active site B-helix and a repositioning of cYY toward the heme. These findings implicate a slowly exchanging water network as a critical factor in CYP121A1 function and a likely contributor to the unusual rigidity of the structure.

PH-Dependent Structure and Dynamics of the Catalytic Domains of Human Carbonic Anhydrase II and IX.

Extensive computer simulation studies have been carried out to probe the pH-dependent structure and dynamics of the two most efficient isoenzymes II and IX of human carbonic anhydrase (HCA) that control the pH in the human body. The equilibrium structure and hydration of their catalytic domains are found to be largely unaffected by the variation of pH in the range studied, in close agreement with the known experimental results. In contrast, a significant effect of the change in pH is observed for the first time on the local electrostatic potential of the active site walls and the dynamics of active site water molecules. We also report for the first time the free energy and kinetics of coupled fluctuations of orientation and protonation states of the well-known His-mediated proton shuttle (His-64) in both isozymes at pH 7 and 8. The transitions between different tautomers of in or out conformations of His-64 side chain range between 109 and 106 s-1 depending on pH. Possible implications of these results on conformation-dependent pKa of His-64 side chain and its role in driving the catalysis toward hydration of CO2 or dehydration of HCO3- with varying pH are discussed.

Computational Insights on the Hydride and Proton Transfer Mechanisms of D-Arginine Dehydrogenase.

D-Arginine dehydrogenase from Pseudomonas aeruginosa (PaDADH) is an amine oxidase which catalyzes the conversion of D-arginine into iminoarginine. It contains a non-covalent FAD cofactor that is involved in the oxidation mechanism. Based on substrate, solvent, and multiple kinetic isotope effects studies, a stepwise hydride transfer mechanism is proposed. It was shown that D-arginine binds to the active site of enzyme as α-amino group protonated, and it is deprotonated before a hydride ion is transferred from its α-C to FAD. Based on a mutagenesis study, it was concluded that a water molecule is the most likely catalytic base responsible from the deprotonation of α-amino group. In this study, we formulated computational models based on ONIOM method to elucidate the oxidation mechanism of D-arginine into iminoarginine using the crystal structure of enzyme complexed with iminoarginine. The calculations showed that Arg222, Arg305, Tyr249, Glu87, His 48, and two active site water molecules play key roles in binding and catalysis. Model systems showed that the deprotonation step occurs prior to hydride transfer step, and active site water molecule(s) may have participated in the deprotonation process.

Porosity as a Design Element for Developing Catalytic Molecular Materials for Electrochemical and Photochemical Carbon Dioxide Reduction.

The catalytic reduction of carbon dioxide (CO2 ) using sustainable energy inputs is a promising strategy for upcycling of atmospheric carbon into value-added chemical products. This goal has inspired the development of catalysts for selective and efficient CO2 conversion using electrochemical and photochemical methods. Among the diverse array of catalyst systems designed for this purpose, two- and three-dimensional platforms that feature porosity offer the potential to combine carbon capture and conversion. Included are covalent organic frameworks (COFs), metal-organic frameworks (MOFs), porous molecular cages, and other hybrid molecular materials developed to increase active site exposure, stability, and water compatibility while maintaining precise molecular tunability. This mini-review showcases catalysts for the CO2 reduction reaction (CO2 RR) that incorporate well-defined molecular elements integrated into porous materials structures. Selected examples provide insights into how different approaches to this overall design strategy can augment their electrocatalytic and/or photocatalytic CO2 reduction activity. This article is protected by copyright. All rights reserved.

Thiolase: A Versatile Biocatalyst Employing CoA-Thioester Chemistry for Making and Breaking C-C Bonds.

Thiolases are CoA-dependent enzymes that catalyze the thiolytic cleavage of 3-ketoacyl-CoA, as well as its reverse reaction, which is the thioester-dependent Claisen condensation reaction. Thiolases are dimers or tetramers (dimers of dimers). All thiolases have two reactive cysteines: (a) a nucleophilic cysteine, which forms a covalent intermediate, and (b) an acid/base cysteine. The best characterized thiolase is the Zoogloea ramigera thiolase, which is a bacterial biosynthetic thiolase belonging to the CT-thiolase subfamily. The thiolase active site is also characterized by two oxyanion holes, two active site waters, and four catalytic loops with characteristic amino acid sequence fingerprints. Three thiolase subfamilies can be identified, each characterized by a unique sequence fingerprint for one of their catalytic loops, which causes unique active site properties. Recent insights concerning the thiolase reaction mechanism, as obtained from recent structural studies, as well as from classical and recent enzymological studies, are addressed, and open questions are discussed. Expected final online publication date for the Annual Review of Biochemistry, Volume 92 is June 2023. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Structural insights into the substrate recognition of serine palmitoyltransferase from Sphingobacterium multivorum

Serine palmitoyltransferase (SPT) is a key enzyme of sphingolipid biosynthesis, which catalyzes the pyridoxal-5'-phosphate-dependent decarboxylative condensation reaction of L-serine (L-Ser) and palmitoyl-CoA (PalCoA) to form 3-ketodihydrosphingosine called long chain base (LCB). SPT is also able to metabolize L-alanine (L-Ala) and glycine (Gly), albeit with much lower efficiency. Human SPT is a membrane-bound large protein complex containing SPTLC1/SPTLC2 heterodimer as the core subunits, and it is known that mutations of the SPTLC1/SPTLC2 genes increase the formation of deoxy-type of LCBs derived from L-Ala and Gly to cause some neurodegenerative diseases. In order to study the substrate recognition of SPT, we examined the reactivity of Sphingobacterium multivorum SPT on various amino acids in the presence of PalCoA. The S. multivorum SPT could convert not only L-Ala and Gly but also L-homoserine, in addition to L-Ser, into the corresponding LCBs. Furthermore, we obtained high-quality crystals of the ligand-free form and the binary complexes with a series of amino acids, including a nonproductive amino acid, L-threonine, and determined the structures at 1.40-1.55 Å resolutions. The S. multivorum SPT accommodated various amino acid substrates through subtle rearrangements of the active-site amino acid residues and water molecules. It was also suggested that non-active-site residues mutated in the human SPT genes might indirectly influence the substrate specificity by affecting the hydrogen-bonding networks involving the bound substrate, water molecules, and amino acid residues in the active site of this enzyme. Collectively, our results highlight SPT structural features affecting substrate specificity for this stage of sphingolipid biosynthesis.

Mutagenic analysis of actin reveals the mechanism of His161 flipping that triggers ATP hydrolysis.

The dynamic assembly of actin is controlled by the hydrolysis of ATP, bound to the center of the molecule. Upon polymerization, actin undergoes a conformational change from the monomeric G-form to the fibrous F-form, which is associated with the flipping of the side chain of His161 toward ATP. His161 flipping from the gauche-minus to gauche-plus conformation leads to a rearrangement of the active site water molecules, including ATP attacking water (W1), into an orientation capable of hydrolysis. We previously showed that by using a human cardiac muscle α-actin expression system, mutations in the Pro-rich loop residues (A108G and P109A) and in a residue that was hydrogen-bonded to W1 (Q137A) affect the rate of polymerization and ATP hydrolysis. Here, we report the crystal structures of the three mutant actins bound to AMPPNP or ADP-Pi determined at a resolution of 1.35-1.55Å, which are stabilized in the F-form conformation with the aid of the fragmin F1 domain. In A108G, His161 remained non-flipped despite the global actin conformation adopting the F-form, demonstrating that the side chain of His161 is flipped to avoid a steric clash with the methyl group of A108. Because of the non-flipped His161, W1 was located away from ATP, similar to G-actin, which was accompanied by incomplete hydrolysis. In P109A, the absence of the bulky proline ring allowed His161 to be positioned near the Pro-rich loop, with a minor influence on ATPase activity. In Q137A, two water molecules replaced the side-chain oxygen and nitrogen of Gln137 almost exactly at their positions; consequently, the active site structure, including the W1 position, is essentially conserved. This seemingly contradictory observation to the reported low ATPase activity of the Q137A filament could be attributed to a high fluctuation of the active site water. Together, our results suggest that the elaborate structural design of the active site residues ensures the precise control of the ATPase activity of actin.

A Multiple Proton Transfer Mechanism for the Charging Step of the Aminoacylation Reaction at the Active Site of Aspartyl tRNA Synthetase.

Aspartyl-tRNA synthetase catalyzes the attachment of aspartic acid to its cognate tRNA by the aminoacylation reaction during the initiation of the protein biosynthesis process. In the second step of the aminoacylation reaction, known as the charging step, the aspartate moiety is transferred from aspartyl-adenylate to the 3'-OH of A76 of tRNA through a proton transfer process. We have investigated different pathways for the charging step through three separate QM/MM simulations combined with the enhanced sampling method of well-sliced metadynamics and found out the most feasible pathway for the reaction at the active site of the enzyme. In the charging reaction, both the phosphate group and the ammonium group after deprotonation can potentially act as a base for proton transfer in the substrate-assisted mechanism. We have considered three possible mechanisms involving different pathways of proton transfer, and only one of them is determined to be enzymatically feasible. The free energy landscape along reaction coordinates where the phosphate group acts as the general base showed that, in the absence of water, the barrier height is 52.6 kcal/mol. The free energy barrier is reduced to 39.7 kcal/mol when the active site water molecules are also treated quantum mechanically, thus allowing a water mediated proton transfer. The charging reaction involving the ammonium group of the aspartyl adenylate is found to follow a path where first a proton from the ammonium group moves to a water in the vicinity forming a hydronium ion (H3O+) and NH2 group. The hydronium ion subsequently passes the proton to the Asp233 residue, thus minimizing the chance of back proton transfer from hydronium to the NH2 group. The neutral NH2 group subsequently takes the proton from the O3' of A76 with a free energy barrier of 10.7 kcal/mol. In the next step, the deprotonated O3' makes a nucleophilic attack to the carbonyl carbon forming a tetrahedral transition state with a free energy barrier of 24.8 kcal/mol. Thus, the present work shows that the charging step proceeds through a multiple proton transfer mechanism where the amino group formed after deprotonation acts as the base to capture a proton from O3' of A76 rather than the phosphate group. The current study also shows the important role played by Asp233 in the proton transfer process.

Insight into the Role of Active Site Protonation States and Water Molecules in the Catalytic Inhibition of DPP4 by Vildagliptin.

Vildagliptin (VIL) is an antidiabetic drug that inhibits dipeptidyl peptidase-4 (DPP4) through a covalent mechanism. The molecular bases for this inhibitory process have been addressed experimentally and computationally. Nevertheless, relevant issues remain unknown such as the roles of active site protonation states and conserved water molecules nearby the catalytic center. In this work, molecular dynamics simulations were applied to examine the structures of 12 noncovalent VIL-DPP4 complexes encompassing all possible protonation states of three noncatalytic residues (His126, Asp663, Asp709) that were inconclusively predicted by different computational tools. A catalytically competent complex structure was only achieved in the system with His126 in its ε-form and nonconventional neutral states for Asp663/Asp709. This complex suggested the involvement of one water molecule in the catalytic process of His740/Ser630 activation, which was confirmed by QM/MM simulations. Our findings support the suitability of a novel water-mediated mechanism in which His740/Ser630 activation occurs concertedly with the nucleophilic attack on VIL and the imidate protonation by Tyr547. Then, the restoration of His740/ Tyr547 protonation states occurs via a two-water hydrogen bonding network in a low-barrier process, thus describing the final step of the catalytic cycle for the first time. Additionally, two hydrolytic mechanisms were proposed based on the hydrogen bonding networks formed by water molecules and the catalytic residues along the inhibitory mechanism. These findings are valuable to unveil the molecular features of the covalent inhibition of DPP4 by VIL and support the future development of novel derivatives with improved structural or mechanistic profiles.

Integrated interfacial design of covalent organic framework photocatalysts to promote hydrogen evolution from water

Attempts to develop photocatalysts for hydrogen production from water usually result in low efficiency. Here we report the finding of photocatalysts by integrated interfacial design of stable covalent organic frameworks. We predesigned and constructed different molecular interfaces by fabricating ordered or amorphous π skeletons, installing ligating or non-ligating walls and engineering hydrophobic or hydrophilic pores. This systematic interfacial control over electron transfer, active site immobilisation and water transport enables to identify their distinct roles in the photocatalytic process. The frameworks, combined ordered π skeletons, ligating walls and hydrophilic channels, work under 300–1000 nm with non-noble metal co-catalyst and achieve a hydrogen evolution rate over 11 mmol g–1 h–1, a quantum yield of 3.6% at 600 nm and a three-order-of-magnitude-increased turnover frequency of 18.8 h–1 compared to those obtained with hydrophobic networks. This integrated interfacial design approach is a step towards designing solar-to-chemical energy conversion systems.

Probing the general base for DNA polymerization in telomerase: a molecular dynamics investigation.

Telomerase is an RNA-dependent DNA polymerase that plays a role in the maintenance of the 3' end of the eukaryotic chromosome, known as a telomere, by catalyzing the DNA polymerization reaction in cancer and embryonic stem cells. The detailed molecular details of the DNA polymerization by telomerase, especially the general base for deprotonating the terminal 3'-hydroxyl, which triggers the chemical reaction, remain elusive. We conducted a computational investigation using hybrid quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulations to probe the detailed mechanism of the reaction. Our simulations started with the telomerase:RNA:DNA:dNTP ternary complex, and by using enhanced sampling QM/MM MD simulations, we probed the general base involved directly in the polymerization. We report the participation of an aspartate (Asp344) coordinated to Mg and an active site water molecule, jointly acting as a base during nucleic acid addition. The Asp344 residue remains transiently protonated during the course of the reaction, and later it deprotonates by transferring its proton to the water at the end of the reaction.

Structural and Energetic Origin of Different Product Specificities and Activities for SETD3 and Its Mutants on the Methylation of the β-Actin H73K Peptide: Insights from a QM/MM Study.

The methylation of the lysine residue can affect some fundamental biological processes, and specific biological effects of the methylations are often related to product specificity of methyltransferases. The question remains concerning how active-site structural features and dynamics control the activity as well as the number (1, 2, or 3) of methyl groups on methyl lysine products. SET domain containing protein 3 (SETD3) has been identified recently as the β-actin histidine73-N3 methyltransferase, and also, it has a weak methylation activity on the H73K β-actin peptide for which the target H73 residue is mutated into K73. Interestingly, the K73 methylation activity of SETD3 increases significantly as a result of the N255 → A or N255 → F/W273 → A mutation, and the N255A product specificity also differs from that of wild-type. Here, we performed QM/MM molecular dynamics and potential of mean force (PMF) simulations for SETD3 and its mutants (N255A and N255F/W273A) to study how SETD3 and its mutants could have different product specificities and activities for the K73 methylation. The PMF simulations show that the barrier for the first methylation of K73 is higher compared to the barrier of the H73 methylation in SETD3. Moreover, the second methylation of K73 has been found to have a barrier from the free energy simulation that is higher by 2.2 kcal/mol compared to the barrier of the first methyl transfer to K73, agreeing with the suggestion that SETD3 is a monomethylase. For the first, second, and third methylations of K73 in the N255A mutant, the barriers obtained from the PMF simulations for transferring the second and third methyl groups are found to be lower relative to the barrier for the first methyl transfer. Thus, N255A can be considered as a trimethyl lysine methyltransferase. In addition, for the first K73 methylation, the activities from the PMF simulations follow the order of N255F/W273A > N255A > WT, in agreement with experiments. The examination of the structural and dynamic results at the active sites provides better understanding of different product specificities and activities for the K73 methylations in SETD3 and its mutants. It is demonstrated that the existence of well-balanced interactions at the active site leading to the near attack conformation is of crucial importance for the efficient methyl transfers. Moreover, the presence of potential interactions (e.g., the C-H···O and cation-π interactions) that are strengthening at the transition state can also be important. Furthermore, the activity as well as product specificity of the K73 methylation also seems to be controlled by certain active-site water molecules which may be released to provide extra space for the addition of more methyl groups on K73.

Molecular dynamics simulation study on Bacillus subtilis EngA: the presence of Mg2+ at the active-sites promotes the functionally important conformation

EngA, a GTPase contains two GTP binding domains [GD1, GD2], and the C-terminal KH domain shown to be involved in the later stages of ribosome maturation. Association of EngA to the ribosomal subunit in the intermediate stage of maturation is essential for complete ribosome maturation. However, this association was shown to be dependent on the nucleotide bound combinations. This nucleotide dependent association tendency is attributed to the conformational changes that occur among different nucleotide bound combinations. Therefore, to explore the conformational changes, all-atom molecular dynamics simulations for Bacillus subtilis EngA in different nucleotide bound combinations along with the presence or absence of Mg2+ in the active-sites were carried out. The presence of Mg2+ along with the bound nucleotide at the GD2 active-site dictates the GD2-Sw-II mobility, but the GD1-Sw-II mobility has not shown any nucleotide or Mg2+ dependent movement. However, the GD1-Sw-II secondary conformations are shown to be influenced by the GD2 nucleotide bound state. This allosteric connection between the GD2 active-site and the GD1-Sw-II is also observed through the dynamic network analysis. Further, the exploration of the GD1-KH interface interactions exhibited a more attractive tendency when GD1 is bound to GTP-Mg2+. In addition, the presence of Mg2+ stabilizes active-site water and also increases the distances between the α- and γ- phosphates of the bound GTP. Curiously, three water molecules in the GD1 active-site and only one water molecule in the GD2 active-site are stabilized. This indicates that the probability of GTP hydrolysis is more in GD1 compared to GD2. Communicated by Ramaswamy H. Sarma

Monofunctional Heme-Catalases.

The review focuses on four issues that are critical for the understanding of monofunctional catalases. How hydrogen peroxide (H2O2) reaches the active site and outcompetes water molecules to be able to function at a very high rate is one of the issues examined. Part of the answer is a gate valve system that is instrumental to drive out solvent molecules from the final section of the main channel. A second issue relates to how the enzyme deals with an unproductive reactive compound I (Cpd I) intermediate. Peroxidatic two and one electron donors and the transfer of electrons to the active site from NADPH and other compounds are reviewed. The new ascribed catalase reactions are revised, indicating possible measurement pitfalls. A third issue concerns the heme b to heme d oxidation, why this reaction occurs only in some large-size subunit catalases (LSCs), and the possible role of singlet oxygen in this and other modifications. The formation of a covalent bond between the proximal tyrosine with the vicinal residue is analyzed. The last issue refers to the origin and function of the additional C-terminal domain (TD) of LSCs. The TD has a molecular chaperone activity that is traced to a gene fusion between a Hsp31-type chaperone and a small-size subunit catalase (SSC).

Absolute Binding Free Energy Calculations for Buried Water Molecules.

Water often plays a key role in mediating protein-ligand interactions. Understanding contributions from active-site water molecules to binding thermodynamics of a ligand is important in predicting binding free energies for ligand optimization. In this work, we tested a non-equilibrium switching method for absolute binding free energy calculations on water molecules in binding sites of 13 systems. We discuss the lessons we learned about identified issues that affected our calculations and ways to address them. This work fits with our larger focus on how to do accurate ligand binding free energy calculations when water rearrangements are very slow, such as rearrangements due to ligand modification (as in relative free energy calculations) or ligand binding (as in absolute free energy calculations). The method studied in this work can potentially be used to account for limited water sampling via providing endpoint corrections to free energy calculations using our calculated binding free energy of water.