CAVER: Algorithms for Analyzing Dynamics of Tunnels in Macromolecules

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.

Computational studies on the catalytic potential of the double active site for enzyme engineering

Proteins possessing double active sites have the potential to revolutionise enzyme design strategies. This study extensively explored an enzyme that contains both a natural active site (NAS) and an engineered active site (EAS), focusing on understanding its structural and functional properties. Metadynamics simulations were employed to investigate how substrates interacted with their respective active sites. The results revealed that both the NAS and EAS exhibited similar minimum energy states, indicating comparable binding affinities. However, it became apparent that the EAS had a weaker binding site for the substrate due to its smaller pocket and constrained conformation. Interestingly, the EAS also displayed dynamic behaviour, with the substrate observed to move outside the pocket, suggesting the possibility of substrate translocation. To gain further insights, steered molecular dynamics (SMD) simulations were conducted to study the conformational changes of the substrate and its interactions with catalytic residues. Notably, the substrate adopted distinct conformations, including near-attack conformations, in both the EAS and NAS. Nevertheless, the NAS demonstrated superior binding minima for the substrate compared to the EAS, reinforcing the observation that the engineered active site was less favourable for substrate binding due to its limitations. The QM/MM (Quantum mechanics and molecular mechanics) analyses highlight the energy disparity between NAS and EAS. Specifically, EAS exhibited elevated energy levels due to its engineered active site being located on the surface. This positioning exposes the substrate to solvents and water molecules, adding to the energy challenge. Consequently, the engineered enzyme did not provide a significant advantage in substrate binding over the single active site protein. Further, the investigation of internal channels and tunnels within the protein shed light on the pathways facilitating transport between the two active sites. By unravelling the complex dynamics and functional characteristics of this double-active site protein, this study offers valuable insights into novel strategies of enzyme engineering. These findings establish a solid foundation for future research endeavours aimed at harnessing the potential of double-active site proteins in diverse biotechnological applications.

Meningitis caused by Aeromonas hydrophila in Oreochromis niloticus: Proteomics and druggability of virulence factors

Meningitis caused by Gram-negative bacteria is a serious public health problem, causing morbidity and mortality in both children and adults. Here, we propose a novel experimental model using Nile tilapia (Oreochromis niloticus) to study neuroinflammation. The fish were infected with Aeromonas hydrophila, and the course of infection was monitored in the peripheral blood. Septicemia was obvious in the blood, while in the brain tissue, infection of the meninges was present. The histopathological examination showed suppurative meningitis, and the cellular immune response in the brain tissue during infection was mediated by microglia. These cells were morphologically characterized and phenotyped by MHC class II markers and CD68. The increased production of TNF-α, IL-1β and iNOS supported the infiltration of these cells during the neuroinflammatory process. In the proteomic analysis of A. hydrophila isolated from brain tissue, we found chemotactic and transport proteins, proteolytic enzymes and enzymes associated with the dismutation of nitric oxide (NO), as well as motor proteins and those responsible for cell division. After characterizing the most abundant proteins during the course of infection, we investigated the druggability index of these proteins and identified promising peptide sequences as molecular targets that are similar among bacteria. Thus, these findings deepened the understanding of the pathophysiology of meningitis caused by A. hydrophila. Moreover, through the proteomics analysis, important mechanisms and pathways used by the pathogen to subvert the host response were revealed, providing insights for the development of novel antibiotics and vaccines.

Entrance channels to coproheme in coproporphyrin ferrochelatase probed by exogenous imidazole binding

Iron insertion into porphyrins is an essential step in heme biosynthesis. In the coproporphyrin-dependent pathway, specific to monoderm bacteria, this reaction is catalyzed by the monomeric enzyme coproporphyrin ferrochelatase. In addition to the mechanistic details of the metalation of the porphyrin, the identification of the substrate access channel for ferrous iron to the active site is important to fully understand this enzymatic system. In fact, whether the iron reaches the active site from the distal or the proximal porphyrin side is still under debate. In this study we have thoroughly addressed this question in Listeria monocytogenes coproporphyrin ferrochelatase by X-ray crystallography, steady-state and pre-steady-state imidazole ligand binding studies, together with a detailed spectroscopic characterization using resonance Raman and UV-vis absorption spectroscopies in solution. Analysis of the X-ray structures of coproporphyrin ferrochelatase-coproporphyrin III crystals soaked with ferrous iron shows that iron is present on both sides of the porphyrin. The kinetic and spectroscopic study of imidazole binding to coproporphyrin ferrochelatase‑iron coproporphyrin III clearly indicates the presence of two possible binding sites in this monomeric enzyme that influence each other, which is confirmed by the observed cooperativity at steady-state and a biphasic behavior in the pre-steady-state experiments. The current results are discussed in the context of the entire heme biosynthetic pathway and pave the way for future studies focusing on protein-protein interactions.

Engineering of Rieske dioxygenase variants with improved cis-dihydroxylation activity for benzoates

Rieske dioxygenases have a long history of being utilized as green chemical tools in the organic synthesis of high-value compounds, due to their capacity to perform the cis-dihydroxylation of a wide variety of aromatic substrates. The practical utility of these enzymes has been hampered however by steric and electronic constraints on their substrate scopes, resulting in limited reactivity with certain substrate classes. Herein, we report the engineering of a widely used member of the Rieske dioxygenase class of enzymes, toluene dioxygenase (TDO), to produce improved variants with greatly increased activity for the cis-dihydroxylation of benzoates. Through rational mutagenesis and screening, TDO variants with substantially improved activity over the wild-type enzyme were identified. Homology modeling, docking studies, molecular dynamics simulations, and substrate tunnel analysis were applied in an effort to elucidate how the identified mutations resulted in improved activity for this polar substrate class. These analyses revealed modification of the substrate tunnel as the likely cause of the improved activity observed with the best-performing enzyme variants.

Atomic Insight into the Enzymatic Selectivity of Acetyltransferase for Endogenous Polyamines

The N1-Spermidine/spermine acetyltransferase (SSAT) serves as the rate-limiting enzyme in the polyamine metabolism pathway, specifically catalyzing the acetylation of spermidine, spermine, and other specific polyamines. The source of its enzymatic selectivity remains elusive. Here, we used quantum mechanics and molecular mechanics simulations combined with various technologies to explore the enzymatic mechanism of SSAT for endogenous polyamines from an atomic perspective. The static binding and chemical transformation were considered. The binding affinity was identified to be dependent on protonated state of polyamine. The order of the binding affinity for Spm, Spd, and Put is consistent with the experimental results, which is also verified by the dynamic separation of polyamine and SSAT. Hydrogen bond interactions and salt bridges contribute most, and the common hot residues were identified. In addition, the transfer of acetyl and proton between polyamine and AcCoA was discovered to follow a concert mechanism, and thermodynamic properties are responsible for the catalytic efficiency of SSAT. This work may be helpful for development of polyamine derivatives based on catalysis to regulate polyamine metabolism.

Structural basis of lipid head group entry to the Kennedy pathway by FLVCR1

Phosphatidylcholine and phosphatidylethanolamine, the two most abundant phospholipids in mammalian cells, are synthesized de novo by the Kennedy pathway from choline and ethanolamine, respectively1-6. Despite the essential roles of these lipids, the mechanisms that enable the cellular uptake of choline and ethanolamine remain unknown. Here we show that the protein encoded by FLVCR1, whose mutation leads to the neurodegenerative syndrome posterior column ataxia and retinitis pigmentosa7-9, transports extracellular choline and ethanolamine into cells for phosphorylation by downstream kinases to initiate the Kennedy pathway. Structures of FLVCR1 in the presence of choline and ethanolamine reveal that both metabolites bind to a common binding site comprising aromatic and polar residues. Despite binding to a common site, FLVCR1 interacts in different ways with the larger quaternary amine of choline in and with the primary amine of ethanolamine. Structure-guided mutagenesis identified residues that are crucial for the transport of ethanolamine, but dispensable for choline transport, enabling functional separation of the entry points into the two branches of the Kennedy pathway. Altogether, these studies reveal how FLVCR1 is a high-affinity metabolite transporter that serves as the common origin for phospholipid biosynthesis by two branches of the Kennedy pathway.

Structural insights into the distinct substrate preferences of two bacterial epoxide hydrolases

Epoxide hydrolases (EHs), which catalyze the transformation of epoxides to diols, are present in many eukaryotic and prokaryotic organisms. They have recently drawn considerable attention from organic chemists owing to their application in the semisynthesis of enantiospecific diol compounds. Here, we report the crystal structures of BoEH from Bosea sp. PAMC 26642 and CaEH from Caballeronia sordidicola PAMC 26510 at 1.95 and 2.43 Å resolution, respectively. Structural analysis showed that the overall structures of BoEH and CaEH commonly possess typical α/β hydrolase fold with the same ring-opening residues (Tyr-Tyr) and conserved catalytic triad residues (Asp-Asp-His). However, the two enzymes were found to have significantly different sequence compositions in the cap domain region, which is involved in the formation of the substrate-binding site in both enzymes. Enzyme activity assay results showed that BoEH had the strongest activity toward the linear aliphatic substrates, whereas CaEH had a higher preference for aromatic- and cycloaliphatic substrates. Computational docking simulations and tunnel identification revealed important residues with different substrate-binding preferences. Collectively, structure comparison studies, together with ligand docking simulation results, suggested that the differences in substrate-binding site residues were highly correlated with substrate specificity.

Catalpolaglycone disrupts mitochondrial thermogenesis by specifically binding to a conserved lysine residue of UCP2 on the proton leak tunnel

BACKGROUND

Catalpol (CAT), a naturally occurring iridoid glycoside sourced from the root of Rehmannia glutinosa, affects mitochondrial metabolic functions. However, the mechanism of action of CAT against pyrexia and its plausible targets remain to be fully elucidated.

PURPOSE

This study aimed to identify the specific targets of CAT for blocking mitochondrial thermogenesis and to unveil the unique biological mechanism of action of the orthogonal binding mode between the hemiacetal group and lysine residue on the target protein in vivo.

METHODS

Lipopolysaccharide (LPS)/ carbonyl cyanide 3-chlorophenylhydrazone (CCCP)-induced fever models were established to evaluate the potential antipyretic effects of CAT. An alkenyl-modified CAT probe was designed to identify and capture potential targets. Binding capacity was tested using in-gel imaging and a cellular thermal shift assay. The underlying antipyretic mechanisms were explored using biochemical and molecular biological methods. Catalpolaglycone (CA) was coupled with protein profile identification and molecular docking analysis to evaluate and identify its binding mode to UCP2.

RESULTS

After deglycation of CAT in vivo, the hemiacetal group in CA covalently binds to Lys239 of UCP2 in the mitochondria of the liver via an ɛ-amine nucleophilic addition. This irreversible binding affects proton leakage and improves mitochondrial membrane potential and ADP/ATP transformation efficiency, leading to an antipyretic effect.

CONCLUSION

Our findings highlight the potential role of CA in modulating UCP2 activity or function within the mitochondria and open new avenues for investigating the therapeutic effects of CA on mitochondrial homeostasis.

Fungi Tryptophan Synthases: What Is the Role of the Linker Connecting the α and β Structural Domains in Hemileia vastatrix TRPS? A Molecular Dynamics Investigation

Tryptophan synthase (TRPS) is a complex enzyme responsible for tryptophan biosynthesis. It occurs in bacteria, plants, and fungi as an αββα heterotetramer. Although encoded by independent genes in bacteria and plants, in fungi, TRPS is generated by a single gene that concurrently expresses the α and β entities, which are linked by an elongated peculiar segment. We conducted 1 µs all-atom molecular dynamics simulations on Hemileia vastatrix TRPS to address two questions: (i) the role of the linker segment and (ii) the comparative mode of action. Since there is not an experimental structure, we started our simulations with homology modeling. Based on the results, it seems that TRPS makes use of an already-existing tunnel that can spontaneously move the indole moiety from the α catalytic pocket to the β one. Such behavior was completely disrupted in the simulation without the linker. In light of these results and the αβ dimer's low stability, the full-working TRPS single genes might be the result of a particular evolution. Considering the significant losses that Hemileia vastatrix causes to coffee plantations, our next course of action will be to use the TRPS to look for substances that can block tryptophan production and therefore control the disease.

Active site remodeling in tumor-relevant IDH1 mutants drives distinct kinetic features and potential resistance mechanisms

Mutations in human isocitrate dehydrogenase 1 (IDH1) drive tumor formation in a variety of cancers by replacing its conventional activity with a neomorphic activity that generates an oncometabolite. Little is understood of the mechanistic differences among tumor-driving IDH1 mutants. We previously reported that the R132Q mutant uniquely preserves conventional activity while catalyzing robust oncometabolite production, allowing an opportunity to compare these reaction mechanisms within a single active site. Here, we employed static and dynamic structural methods and found that, compared to R132H, the R132Q active site adopted a conformation primed for catalysis with optimized substrate binding and hydride transfer to drive improved conventional and neomorphic activity over R132H. This active site remodeling revealed a possible mechanism of resistance to selective mutant IDH1 therapeutic inhibitors. This work enhances our understanding of fundamental IDH1 mechanisms while pinpointing regions for improving inhibitor selectivity.

Evolutionary Engineering a Larger Porin Using a Loop-to-Hairpin Mechanism

In protein evolution, diversification is generally driven by genetic duplication. The hallmarks of this mechanism are visible in the repeating topology of various proteins. In outer membrane β-barrels, duplication is visible with β-hairpins as the repeating unit of the barrel. In contrast to the overall use of duplication in diversification, a computational study hypothesized evolutionary mechanisms other than hairpin duplications leading to increases in the number of strands in outer membrane β-barrels. Specifically, the topology of some 16- and 18-stranded β-barrels appear to have evolved through a loop to β-hairpin transition. Here we test this novel evolutionary mechanism by creating a chimeric protein from an 18-stranded β-barrel and an evolutionarily related 16-stranded β-barrel. The chimeric combination of the two was created by replacing loop L3 of the 16-stranded barrel with the sequentially matched transmembrane β-hairpin region of the 18-stranded barrel. We find the resulting chimeric protein is stable and has characteristics of increased strand number. This study provides the first experimental evidence supporting the evolution through a loop to β-hairpin transition.

Structural basis of lipid head group entry to the Kennedy pathway by FLVCR1

Phosphatidylcholine and phosphatidylethanolamine, the two most abundant phospholipids in mammalian cells, are synthesized de novo by the Kennedy pathway from choline and ethanolamine, respectively1-6. Despite the importance of these lipids, the mechanisms that enable the cellular uptake of choline and ethanolamine remain unknown. Here, we show that FLVCR1, whose mutation leads to the neurodegenerative syndrome PCARP7-9, transports extracellular choline and ethanolamine into cells for phosphorylation by downstream kinases to initiate the Kennedy pathway. Structures of FLVCR1 in the presence of choline and ethanolamine reveal that both metabolites bind to a common binding site comprised of aromatic and polar residues. Despite binding to a common site, the larger quaternary amine of choline interacts differently with FLVCR1 than does the primary amine of ethanolamine. Structure-guided mutagenesis identified residues that are critical for the transport of ethanolamine, while being dispensable for choline transport, enabling functional separation of the entry points into the two branches of the Kennedy pathway. Altogether, these studies reveal how FLCVR1 is a high-affinity metabolite transporter that serves as the common origin for phospholipid biosynthesis by two branches of the Kennedy pathway.

Open-chain thiamine analogues as potent inhibitors of thiamine pyrophosphate (TPP)-dependent enzymes

A common approach to studying thiamine pyrophosphate (TPP)-dependent enzymes is by chemical inhibition with thiamine/TPP analogues which feature a neutral aromatic ring in place of the positive thiazolium ring of TPP. These are potent inhibitors but their preparation generally involves multiple synthetic steps to construct the central ring. We report efficient syntheses of novel, open-chain thiamine analogues which potently inhibit TPP-dependent enzymes and are predicted to share the same binding mode as TPP. We also report some open-chain analogues that inhibit pyruvate dehydrogenase E1-subunit (PDH E1) and are predicted to occupy additional pockets in the enzyme other than the TPP-binding pockets. This opens up new possibilities for increasing the affinity and selectivity of the analogues for PDH, which is an established anti-cancer target.

Crystal structures of multidrug efflux transporters from Burkholderia pseudomallei suggest details of transport mechanism

BpeB and BpeF are multidrug efflux transporters from Burkholderia pseudomallei that enable multidrug resistance. Here, we report the crystal structures of BpeB and BpeF at 2.94 Å and 3.0 Å resolution, respectively. BpeB was found as an asymmetric trimer, consistent with the widely-accepted functional rotation mechanism for this type of transporter. One of the monomers has a distinct structure that we interpret as an intermediate along this functional cycle. Additionally, a detergent molecule bound in a previously undescribed binding site provides insights into substrate translocation through the pathway. BpeF shares structural similarities with the crystal structure of OqxB from Klebsiella pneumoniae, where both are symmetric trimers composed of three "binding"-state monomers. The structures of BpeB and BpeF further our understanding of the functional mechanisms of transporters belonging to the HAE1-RND superfamily.

Putting the Puzzle Together To Get the Whole Picture: Molecular Basis of the Affinity of Two Steroid Derivatives to Acetylcholinesterase

Alzheimer's disease (AD) is a progressive neurodegenerative disorder that has no cure because its etiology is still unknown, and its main treatment is the administration of acetylcholinesterase (AChE) inhibitors. The study of the mechanism of action of this family of compounds is critical for the design of new more potent and specific inhibitors. In this work, we study the molecular basis of an uncompetitive inhibitor (compound 1, 2β, 3α-dihydroxy-5α-cholestan-6-one disulfate), which we have proved to be a peripheral anionic site (PAS)-binding AChE inhibitor. The pipeline designed in this work is key to the development of other PAS inhibitors that not only inhibit the esterase action of the enzyme but could also modulate the non-cholinergic functions of AChE linked to the process of amylogenesis. Our studies showed that 1 inhibits the enzyme not simply by blocking the main gate but by an allosteric mechanism. A detailed and careful analysis of the ligand binding position and the protein dynamics, particularly regarding their secondary gates and active site, was necessary to conclude this. The same analysis was executed with an inactive analogue (compound 2, 2β, 3α-dihydroxy-5α-cholestan-6-one). Our first computational results showed no differences in affinity to AChE between both steroids, making further analysis necessary. This work highlights the variables to be considered and develops a refined methodology, for the successful design of new potent dual-action drugs for AD, particularly PAS inhibitors, an attractive strategy to combat AD.

Expanding the Physiological Role of Aryl-Alcohol Flavooxidases as Quinone Reductases

Aryl-alcohol oxidases (AAOs) are members of the glucose-methanol-choline oxidase/dehydrogenase (GMC) superfamily. These extracellular flavoproteins have been described as auxiliary enzymes in the degradation of lignin by several white-rot basidiomycetes. In this context, they oxidize fungal secondary metabolites and lignin-derived compounds using O₂ as an electron acceptor, and supply H₂O₂ to ligninolytic peroxidases. Their substrate specificity, including mechanistic aspects of the oxidation reaction, has been characterized in Pleurotus eryngii AAO, taken as a model enzyme of this GMC superfamily. AAOs show broad reducing-substrate specificity in agreement with their role in lignin degradation, being able to oxidize both nonphenolic and phenolic aryl alcohols (and hydrated aldehydes). In the present work, the AAOs from Pleurotus ostreatus and Bjerkandera adusta were heterologously expressed in Escherichia coli, and their physicochemical properties and oxidizing abilities were compared with those of the well-known recombinant AAO from P. eryngii. In addition, electron acceptors different from O₂, such as p-benzoquinone and the artificial redox dye 2,6-Dichlorophenolindophenol, were also studied. Differences in reducing-substrate specificity were found between the AAO enzymes from B. adusta and the two Pleurotus species. Moreover, the three AAOs oxidized aryl alcohols concomitantly with the reduction of p-benzoquinone, with similar or even higher efficiencies than when using their preferred oxidizing-substrate, O₂. IMPORTANCE In this work, quinone reductase activity is analyzed in three AAO flavooxidases, whose preferred oxidizing-substrate is O₂. The results presented, including reactions in the presence of both oxidizing substrates-benzoquinone and molecular oxygen-suggest that such aryl-alcohol dehydrogenase activity, although less important than its oxidase activity in terms of maximal turnover, may have a physiological role during fungal decay of lignocellulose by the reduction of quinones (and phenoxy radicals) from lignin degradation, preventing repolymerization. Moreover, the resulting hydroquinones would participate in redox-cycling reactions for the production of hydroxyl free radical involved in the oxidative attack of the plant cell-wall. Hydroquinones can also act as mediators for laccases and peroxidases in lignin degradation in the form of semiquinone radicals, as well as activators of lytic polysaccharide monooxygenases in the attack of crystalline cellulose. Moreover, reduction of these, and other phenoxy radicals produced by laccases and peroxidases, promotes lignin degradation by limiting repolymerization reactions. These findings expand the role of AAO in lignin biodegradation.

AsiteDesign: a Semirational Algorithm for an Automated Enzyme Design

With advances in protein structure predictions, the number of available high-quality structures has increased dramatically. In light of these advances, structure-based enzyme engineering is expected to become increasingly important for optimizing biocatalysts for industrial processes. Here, we present AsiteDesign, a Monte Carlo-based protocol for structure-based engineering of active sites. AsiteDesign provides a framework for introducing new catalytic residues in a given binding pocket to either create a new catalytic activity or alter the existing one. AsiteDesign is implemented using pyRosetta and incorporates enhanced sampling techniques to efficiently explore the search space. The protocol was tested by designing an alternative catalytic triad in the active site of Pseudomonas fluorescens esterase (PFE). The designed variant was experimentally verified to be active, demonstrating that AsiteDesign can find alternative catalytic triads. Additionally, the AsiteDesign protocol was employed to enhance the hydrolysis of a bulky chiral substrate (1-phenyl-2-pentyl acetate) by PFE. The experimental verification of the designed variants demonstrated that F158L/F198A and F125A/F158L mutations increased the hydrolysis of 1-phenyl-2-pentyl acetate from 8.9 to 66.7 and 23.4%, respectively, and reversed the enantioselectivity of the enzyme from (R) to (S)-enantiopreference, with 32 and 55% enantiomeric excess (ee), respectively.

Enzyme adaptation to habitat thermal legacy shapes the thermal plasticity of marine microbiomes

Microbial communities respond to temperature with physiological adaptation and compositional turnover. Whether thermal selection of enzymes explains marine microbiome plasticity in response to temperature remains unresolved. By quantifying the thermal behaviour of seven functionally-independent enzyme classes (esterase, extradiol dioxygenase, phosphatase, beta-galactosidase, nuclease, transaminase, and aldo-keto reductase) in native proteomes of marine sediment microbiomes from the Irish Sea to the southern Red Sea, we record a significant effect of the mean annual temperature (MAT) on enzyme response in all cases. Activity and stability profiles of 228 esterases and 5 extradiol dioxygenases from sediment and seawater across 70 locations worldwide validate this thermal pattern. Modelling the esterase phase transition temperature as a measure of structural flexibility confirms the observed relationship with MAT. Furthermore, when considering temperature variability in sites with non-significantly different MATs, the broadest range of enzyme thermal behaviour and the highest growth plasticity of the enriched heterotrophic bacteria occur in samples with the widest annual thermal variability. These results indicate that temperature-driven enzyme selection shapes microbiome thermal plasticity and that thermal variability finely tunes such processes and should be considered alongside MAT in forecasting microbial community thermal response.

Computational Investigation of a Series of Small Molecules as Potential Compounds for Lysyl Hydroxylase-2 (LH2) Inhibition

The catalytic function of lysyl hydroxylase-2 (LH2), a member of the Fe(II)/αKG-dependent oxygenase superfamily, is to catalyze the hydroxylation of lysine to hydroxylysine in collagen, resulting in stable hydroxylysine aldehyde-derived collagen cross-links (HLCCs). Reports show that high amounts of LH2 lead to the accumulation of HLCCs, causing fibrosis and specific types of cancer metastasis. Some members of the Fe(II)/αKG-dependent family have also been reported to have intramolecular O2 tunnels, which aid in transporting one of the required cosubstrates into the active site. While LH2 can be a promising target to combat these diseases, efficacious inhibitors are still lacking. We have used computational simulations to investigate a series of 44 small molecules as lead compounds for LH2 inhibition. Tunneling analyses indicate the existence of several intramolecular tunnels. The lengths of the calculated O2-transporting tunnels in holoenzymes are relatively longer than those in the apoenzyme, suggesting that the ligands may affect the enzyme's structure and possibly block (at least partially) the tunnels. The sequence alignment analysis between LH enzymes from different organisms shows that all of the amino acid residues with the highest occurrence rate in the oxygen tunnels are conserved. Our results suggest that the enolate form of diketone compounds establishes stronger interactions with the Fe(II) in the active site. Branching the enolate compounds with functional groups such as phenyl and pyridinyl enhances the interaction with various residues around the active site. Our results provide information about possible leads for further LH2 inhibition design and development.

Geometry-Based versus Small-Molecule Tracking Method for Tunnel Identification: Benefits and Pitfalls

Different methods for tunnel identification, geometry-based and small-molecule tracking approaches, were compared to provide their benefits and pitfalls. Results obtained for both crystal structures and molecular dynamics (MD) simulations were analyzed to investigate if a more computationally demanding method would be beneficial. Careful examination of the results is essential for the low-diameter tunnel description, and assessment of the tunnel functionality based only on their geometrical parameters is challenging. We showed that the small-molecule tracking approach can provide a detailed description of the system; however, it can also be the most computationally demanding.

Programming xenon diffusion in maltose-binding protein

Protein interiors contain void space that can bind small gas molecules. Determination of gas pathways and kinetics in proteins has been an intriguing and challenging task. Here, we combined computational methods and the hyperpolarized xenon-129 chemical exchange saturation transfer (hyper-CEST) NMR technique to investigate xenon (Xe) exchange kinetics in maltose-binding protein (MBP). A salt bridge ∼9 Å from the Xe-binding site formed upon maltose binding and slowed the Xe exchange rate, leading to a hyper-CEST 129Xe signal from maltose-bound MBP. Xe dissociation occurred faster than dissociation of the salt bridge, as shown by 13C NMR spectroscopy and variable-B1 hyper-CEST experiments. "Xe flooding" molecular dynamics simulations identified a surface hydrophobic site, V23, that has good Xe binding affinity. Mutations at this site confirmed its role as a secondary exchange pathway in modulating Xe diffusion. This shows the possibility for site-specifically controlling xenon protein-solvent exchange. Analysis of the available MBP structures suggests a biological role of MBP's large hydrophobic cavity to accommodate structural changes associated with ligand binding and protein-protein interactions.

Structure prediction-based insights into the patatin family of Mycobacterium tuberculosis

Despite its genome sequencing more than two decades ago, the majority of the genes of Mycobacterium tuberculosis remain functionally uncharacterized. Patatins are one such class of proteins that, despite undergoing an expansion in this pathogenic species compared to their non-pathogenic cousins, remain largely unstudied. Recent advances in protein structure prediction using machine learning tools such as AlphaFold2 have provided high-confidence predicted structures for all M. tuberculosis proteins. Here we present detailed analyses of the patatin family of M. tuberculosis using AlphaFold-predicted structures, providing insights into likely modes of regulation, membrane interaction and substrate binding. Regulatory domains within this family of proteins include cyclic nucleotide binding, lid-like domains and other helical domains. Using structural homologues, we identified the likely membrane localization mechanisms and substrate-binding sites. These analyses reveal diversity in their regulatory capacity, mechanisms of membrane binding and likely length of fatty acid substrates. Together, this analysis suggests unique roles for the eight predicted patatins of M. tuberculosis.

Structural insights into the mechanism of the sodium/iodide symporter

The sodium/iodide symporter (NIS) is the essential plasma membrane protein that mediates active iodide (I-) transport into the thyroid gland, the first step in the biosynthesis of the thyroid hormones-the master regulators of intermediary metabolism. NIS couples the inward translocation of I- against its electrochemical gradient to the inward transport of Na+ down its electrochemical gradient1,2. For nearly 50 years before its molecular identification3, NIS was the molecule at the centre of the single most effective internal radiation cancer therapy: radioiodide (131I-) treatment for thyroid cancer2. Mutations in NIS cause congenital hypothyroidism, which must be treated immediately after birth to prevent stunted growth and cognitive deficiency2. Here we report three structures of rat NIS, determined by single-particle cryo-electron microscopy: one with no substrates bound; one with two Na+ and one I- bound; and one with one Na+ and the oxyanion perrhenate bound. Structural analyses, functional characterization and computational studies show the substrate-binding sites and key residues for transport activity. Our results yield insights into how NIS selects, couples and translocates anions-thereby establishing a framework for understanding NIS function-and how it transports different substrates with different stoichiometries and releases substrates from its substrate-binding cavity into the cytosol.

Evaluating the role of trypsin in silk degumming: An in silico approach

The trypsin being universal enzyme forming family of proteases catalyzes the hydrolysis of proteins into amino acids and regenerates the serine hydroxyl an active site. The trypsin enzyme from D. saccharalis, uses sericin as its preferred substrate. Presence of catalytic triad (serine, aspartic acid and histidine) at the substrate binding site of this enzyme is very important for the catalytic activity. In the current study, the interacting mechanism between the substrate sericin protein and enzyme trypsin protein were explored by integrating various computational approaches including physico-chemical properties, biophysical properties, dynamics, gene ontology, molecular docking, protein - protein interactions, binding free energy calculation and structural motifs were studied. The evolutionary study performed by MEGA X showed that trypsin protein sequence (ALE15212.1) is closely related to cocoonase protein sequence (ADG26770.1) from Antheraea pernyi. 3-D models of trypsin and sericin proteins were predicted using I-TASSER and further validated by PROCHECK, and ProSAweb softwares. The predicted trypsin structure model was assigned E.C. no. 3.4.21.4 which refers hydrolytic mechanism. Gene Ontology predicted by QuickGO showed that trypsin has serine hydrolase activity (GO: 00017171), and part of proteolysis (GO: 0006508) as well as protein metabolic process (GO:0019538) actvity. Molecular docking studies between trypsin and sericin proteins were conducted by the HADDOCK 2.4 having best docked protein complex with Z-score - 1.9. 2D and 3D protein-protein interaction was performed with LIGPLOT⁺ and HAWKDOCK, PDBsum, respectively. The amino acid residues interacting across proteins interface are sericin_chain A representing "Ser133, Tyr214, Thr188, Thr243, Ser225, Ser151, Ser156, His294, Arg293, Gly296″ and trypsin_chain B "Lys120, Tyr246, Asn119, Glu239, Ser62, Tyr194, Ile197, Ser171, Tyr169, Gly170″. Based on our results trypsin shows similarity with cocoonase and presumably trypsin can be used as an alternative source in cocoon degumming.

Structure-Guided Insight into the Specificity and Mechanism of a Parasitic Nucleoside Hydrolase

Trichomonas vaginalis is the causative parasitic protozoan of the disease trichomoniasis, the most prevalent, nonviral sexually transmitted disease in the world. T. vaginalis is a parasite that scavenges nucleosides from the host organism via catalysis by nucleoside hydrolase (NH) enzymes to yield purine and pyrimidine bases. One of the four NH enzymes identified within the genome of T. vaginalis displays unique specificity toward purine nucleosides, adenosine and guanosine, but not inosine, and atypically shares greater sequence similarity to the pyrimidine hydrolases. Bioinformatic analysis of this enzyme, adenosine/guanosine-preferring nucleoside ribohydrolase (AGNH), was incapable of identifying the residues responsible for this uncommon specificity, highlighting the need for structural information. Here, we report the X-ray crystal structures of holo, unliganded AGNH and three additional structures of the enzyme bound to fragment and small-molecule inhibitors. Taken together, these structures facilitated the identification of residue Asp231, which engages in substrate interactions in the absence of those residues that typically support the canonical purine-specific tryptophan-stacking specificity motif. An altered substrate-binding pose is mirrored by repositioning within the protein scaffold of the His80 general acid/base catalyst. The newly defined structure-determined sequence markers allowed the assignment of additional NH orthologs, which are proposed to exhibit the same specificity for adenosine and guanosine alone and further delineate specificity classes for these enzymes.

Transgenic Expression of Haemonchus contortus Cytochrome P450 Hco-cyp-13A11 Decreases Susceptibility to Particular but Not All Macrocyclic Lactones in the Model Organism Caenorhabditis elegans

The number of reported macrocyclic lactones (ML) resistance cases across all livestock hosts is steadily increasing. Different studies in the parasitic nematode Haemonchus contortus assume the participation of cytochrome P450s (Cyps) enzymes in ML resistance. Still, functional data about their individual contribution to resistance or substrate specificity is missing. Via microinjection, transgenic Caenorhabditis elegans expressing HCON_00141052 (transgene-Hco-cyp-13A11) from extrachromosomal arrays were generated. After 24 h of exposure to different concentrations of ivermectin (IVM), ivermectin aglycone (IVMa), selamectin (SEL), doramectin (DRM), eprinomectin (EPR), and moxidectin (MOX), motility assays were performed to determine the impact of the H. contortus Cyp to the susceptibility of the worms against each ML. While transgene-Hco-cyp-13A11 significantly decreased susceptibility to IVM (four-fold), IVMa (2-fold), and SEL (3-fold), a slight effect for DRM and no effect for MOX, and EPR was observed. This substrate specificity of Hco-cyp-13A11 could not be explained by molecular modeling and docking studies. Hco-Cyp-13A11 molecular models were obtained for alleles from isolates with different resistance statuses. Although 14 amino acid polymorphisms were detected, none was resistance specific. In conclusion, Hco-cyp-13A11 decreased IVM, IVMa, and SEL susceptibility to a different extent, but its potential impact on ML resistance is not driven by polymorphisms.

Rhodopsin-bestrophin fusion proteins from unicellular algae form gigantic pentameric ion channels

Many organisms sense light using rhodopsins, photoreceptive proteins containing a retinal chromophore. Here we report the discovery, structure and biophysical characterization of bestrhodopsins, a microbial rhodopsin subfamily from marine unicellular algae, in which one rhodopsin domain of eight transmembrane helices or, more often, two such domains in tandem, are C-terminally fused to a bestrophin channel. Cryo-EM analysis of a rhodopsin-rhodopsin-bestrophin fusion revealed that it forms a pentameric megacomplex (~700 kDa) with five rhodopsin pseudodimers surrounding the channel in the center. Bestrhodopsins are metastable and undergo photoconversion between red- and green-absorbing or green- and UVA-absorbing forms in the different variants. The retinal chromophore, in a unique binding pocket, photoisomerizes from all-trans to 11-cis form. Heterologously expressed bestrhodopsin behaves as a light-modulated anion channel.

Engineering the Entrance of a Flavonoid Glycosyltransferase Promotes the Glycosylation of Etoposide Aglycone

Enzyme entrances, which function as the first molecular filters, influence substrate selectivity and enzymatic activity. Because of low binding affinities, engineering enzyme entrances that recognize non-natural substrates is a major challenge for artificial biocatalyst design. Here, the entrance of flavonoid glycosyltransferase UGT78D2 was engineered to promote the recognition of the aglycone of etoposide, a chemotherapeutic agent. We found that Q258, S446, R444, and R450, the key residues surrounding the substrate entrance, specifically guide the flux of etoposide aglycone, which has a high steric hindrance, into the active site; this activity was inferred to be determined by the entrance size and hydrophobic and electrostatic interactions. Engineering the coordination of Q258 and S446 to increase the entrance size and hydrophobic interaction between UGT78D2 and etoposide aglycone increased the affinity by 10.10-fold and the conversion by 10%. The entrance-engineering strategy applied in this study can improve the design of artificial biocatalysts.

Characterisation of ascocorynin biosynthesis in the purple jellydisc fungus Ascocoryne sarcoides

Background

Non-ribosomal peptide synthetase-like (NRPS-like) enzymes are highly enriched in fungal genomes and can be discriminated into reducing and non-reducing enzymes. Non-reducing NRPS-like enzymes possess a C-terminal thioesterase domain that catalyses the condensation of two identical aromatic α-keto acids under the formation of enzyme-specific substrate-interconnecting core structures such as terphenylquinones, furanones, butyrolactones or dioxolanones. Ascocoryne sarcoides produces large quantities of ascocorynin, which structurally resembles a terphenylquinone produced from the condensation of p-hydroxyphenylpyruvate and phenylpyruvate. Since the parallel use of two different substrates by a non-reducing NRPS-like enzyme appeared as highly unusual, we investigated the biosynthesis of ascocorynin in A. sarcoides.

Results

Here, we searched the genome of A. sarcoides for genes coding for non-reducing NRPS-like enzymes. A single candidate gene was identified that was termed acyN. Heterologous gene expression confirmed that AcyN is involved in ascocorynin production but only produces the non-hydroxylated precursor polyporic acid. Although acyN is embedded in an ascocorynin biosynthesis gene cluster, a gene encoding a monooxygenase required for the hydroxylation of polyporic acid was not present. Expression analyses of all monooxygenase-encoding genes from A. sarcoides identified a single candidate that showed the same expression pattern as acyN. Accordingly, heterologous co-expression of acyN and the monooxygenase gene resulted in the production of ascocorynin. Structural modelling of the monooxygenase suggests that the hydrophobic substrate polyporic acid enters the monooxygenase from a membrane facing entry site and is converted into the more hydrophilic product ascocorynin, which prevents its re-entry for a second round of hydroxylation.

Conclusion

This study characterises the first naturally occurring polyporic acid synthetase from an ascomycete. It confirms the high substrate and product specificity of this non-reducing NRPS-like enzyme and highlights the requirement of a monooxygenase to produce the terphenylquinone ascocorynin.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40694-022-00138-7.

Structure and function of ClpXP, a AAA+ proteolytic machine powered by probabilistic ATP hydrolysis

ClpXP is an archetypical AAA+ protease, consisting of ClpX and ClpP. ClpX is an ATP-dependent protein unfoldase and polypeptide translocase, whereas ClpP is a self-compartmentalized peptidase. ClpXP is currently the only AAA+ protease for which high-resolution structures exist, the molecular basis of recognition for a protein substrate is understood, extensive biochemical and genetic analysis have been performed, and single-molecule optical trapping has allowed direct visualization of the kinetics of substrate unfolding and translocation. In this review, we discuss our current understanding of ClpXP structure and function, evaluate competing sequential and probabilistic mechanisms of ATP hydrolysis, and highlight open questions for future exploration.

Structure and activity of a thermally stable mutant of Acanthamoeba actophorin

Actophorin, which was recently tested for crystallization under microgravity on the International Space Station, was subjected to mutagenesis to identify a construct with improved biophysical properties that were expected to improve the extent of diffraction. First, 20 mutations, including one C-terminal deletion of three residues, were introduced individually into actophorin, resulting in modest increases in thermal stability of between +0.5°C and +2.2°C. All but two of the stabilizing mutants increased both the rates of severing F-actin filaments and of spontaneous polymerization of pyrenyl G-actin in vitro. When the individual mutations were combined into a single actophorin variant, Acto-2, the overall thermal stability was 22°C higher than that of wild-type actophorin. When an inactivating S2P mutation in Acto-2 was restored, Acto-2/P2S was more stable by 20°C but was notably more active than the wild-type protein. The inactivating S2P mutation reaffirms the importance that Ser2 plays in the F-actin-severing reaction. The crystal structure of Acto-2 was solved to 1.7 Å resolution in a monoclinic space group, a first for actophorin. Surprisingly, despite the increase in thermal stability, the extended β-turn region, which is intimately involved in interactions with F-actin, is disordered in one copy of Acto-2 in the asymmetric unit. These observations emphasize the complex interplay among protein thermal stability, function and dynamics.

Molecular understanding of acetylcholinesterase adsorption on functionalized carbon nanotubes for enzymatic biosensors

The immobilization of acetylcholinesterase on different nanomaterials has been widely used in the field of amperometric organophosphorus pesticide (OP) biosensors. However, the molecular adsorption mechanism of acetylcholinesterase on a nanomaterial's surface is still unclear. In this work, multiscale simulations were utilized to study the adsorption behavior of acetylcholinesterase from Torpedo californica (TcAChE) on amino-functionalized carbon nanotube (CNT) (NH₂-CNT), carboxyl-functionalized CNT (COOH-CNT) and pristine CNT surfaces. The simulation results show that the active center and enzyme substrate tunnel of TcAChE are both close to and oriented toward the surface when adsorbed on the positively charged NH₂-CNT, which is beneficial to the direct electron transfer (DET) and accessibility of the substrate molecule. Meanwhile, the NH₂-CNT can also reduce the tunnel cost of the enzyme substrate of TcAChE, thereby further accelerating the transfer rate of the substrate from the surface or solution to the active center. However, for the cases of TcAChE adsorbed on COOH-CNT and pristine CNT, the active center and substrate tunnel are far away from the surface and face toward the solution, which is disadvantageous for the DET and transportation of enzyme substrate. These results indicate that NH₂-CNT is more suitable for the immobilization of TcAChE. This work provides a better molecular understanding of the adsorption mechanism of TcAChE on functionalized CNT, and also provides theoretical guidance for the ordered immobilization of TcAChE and the design, development and improvement of TcAChE-OPs biosensors based on functionalized carbon nanomaterials.

Correlated Motions and Dynamics in Different Domains of Epidermal Growth Factor Receptor With L858R and T790M Mutations

Non-small cell lung cancer with an activating epidermal growth factor receptor (EGFR) mutation responds well to targeted drugs. In most cases, drug resistance appears after about a year. Several studies have been conducted on the kinase domain of EGFR to understand the drug resistance mechanism. Since EGFR is a multi-domain protein, mutation in the kinase domain may affect the other domains as well. In this study, we examine the complete structure of the multi-domain EGFR protein and its mutants. We performed molecular dynamics simulations for wildtype EGFR, EGFR with L858R mutation, and EGFR with L858R and T790M mutations. We applied normal mode analysis and complex network analysis to extract the correlated motions in the domains of EGFR. The normal modes are used to construct the dynamic cross-correlation map (DCCM). Simulation results show different patterns of correlated motions in each domain of EGFR mutants compared to the wildtype. In Domains 1 and 3 of the extracellular region, a small number of weak positively correlated motions are extracted. Domains 2 and 4 show large numbers of both positive and negative motions. However, the negatively correlated motions are stronger in mutant structures compared to the wildtype. In Domain 7, some residues showed a positive correlation around the main diagonal. We also identified different communities, nodes and crucial residues in the domains of the structures, which can be important for the function of EGFR. Moreover, hydrogen bond analysis is performed for the stability analysis. The mutant structures have fewer hydrogen bonds compared to the wildtype. Overall, these findings are useful for understanding the dynamics and communications in EGFR domains.

Structural model of human PORCN illuminates disease-associated variants and drug-binding sites

Wnt signaling is essential for normal development and is a therapeutic target in cancer. The enzyme PORCN, or porcupine, is a membrane-bound O-acyltransferase (MBOAT) that is required for the post-translational modification of all Wnts, adding an essential mono-unsaturated palmitoleic acid to a serine on the tip of Wnt hairpin 2. Inherited mutations in PORCN cause focal dermal hypoplasia, and therapeutic inhibition of PORCN slows the growth of Wnt-dependent cancers. Based on homology to mammalian MBOAT proteins, we developed and validated a structural model of human PORCN. The model accommodates palmitoleoyl-CoA and Wnt hairpin 2 in two tunnels in the conserved catalytic core, shedding light on the catalytic mechanism. The model predicts how previously uncharacterized human variants of uncertain significance can alter PORCN function. Drugs including ETC-159, IWP-L6 and LGK-974 dock in the PORCN catalytic site, providing insights into PORCN pharmacologic inhibition. This structural model enhances our mechanistic understanding of PORCN substrate recognition and catalysis, as well as the inhibition of its enzymatic activity, and can facilitate the development of improved inhibitors and the understanding of disease-relevant PORCN mutants. This article has an associated First Person interview with the joint first authors of the paper.

Proteomic analysis capsule synthesis and redox mechanisms in the intracellular survival of group B Streptococcus in fish microglia

Group B Streptococcus (GBS) causes meningitis in neonates and Nile tilapia (Oreochromis niloticus). The molecular mechanisms regulating the intracellular survival of this pathogen in the host cell are complex and crucial for the progression of infection. Thus, we propose the use of GBS-infected Nile tilapia microglia as an in vitro model system simulating infection caused by homologous bacteria in humans. We used this model to evaluate the phagocytic activity, as well as the functional aspects of the capsular proteins A, B, C, and D and the major redox enzymes, and the synergistic role of mechanisms/proteins involved in blocking phagocytic process. We observed that in the intracellular phase, GBS showed enhanced synthesis of the polysaccharide capsule and used superoxide dismutase, thioredoxin, NADH oxidase, and alkyl hydroperoxide reductase to scavenge reactive oxygen species and reactive nitrogen species produced by the host cell. Furthermore, although these virulence mechanisms were effective during the initial hours of infection, they were not able to subvert microglial responses, which partially neutralized the infection. Altogether, our findings provided important information regarding the intracellular survival mechanisms of GBS and perspectives for the production of new drugs and vaccines, through the druggability analysis of specific proteins. In conclusion, tilapia microglia serve as a potent in vitro experimental model for the study of meningitis.

Computational investigations of selected enzymes from two iron and α-ketoglutarate-dependent families

DNA alkylation is used as the key epigenetic mark in eukaryotes, however, most alkylation in DNA can result in deleterious effects. Therefore, this process needs to be tightly regulated. The enzymes of the AlkB and Ten-Eleven Translocation (TET) families are members of the Fe and alpha-ketoglutarate-dependent superfamily of enzymes that are tasked with dealkylating DNA and RNA in cells. Members of these families span all species and are an integral part of transcriptional regulation. While both families catalyze oxidative dealkylation of various bases, each has specific preference for alkylated base type as well as distinct catalytic mechanisms. This perspective aims to provide an overview of computational work carried out to investigate several members of these enzyme families including AlkB, ALKB Homolog 2, ALKB Homolog 3 and Ten-Eleven Translocate 2. Insights into structural details, mutagenesis studies, reaction path analysis, electronic structure features in the active site, and substrate preferences are presented and discussed.

Catalytic and structural properties of ATP-dependent caprolactamase from Pseudomonas jessenii

Caprolactamase is the first enzyme in the caprolactam degradation pathway of Pseudomonas jessenii. It is composed of two subunits (CapA and CapB) and sequence-related to other ATP-dependent enzymes involved in lactam hydrolysis, like 5-oxoprolinases and hydantoinases. Low sequence similarity also exists with ATP-dependent acetone- and acetophenone carboxylases. The caprolactamase was produced in Escherichia coli, isolated by His-tag affinity chromatography, and subjected to functional and structural studies. Activity toward caprolactam required ATP and was dependent on the presence of bicarbonate in the assay buffer. The hydrolysis product was identified as 6-aminocaproic acid. Quantum mechanical modeling indicated that the hydrolysis of caprolactam was highly disfavored (ΔG0 '= 23 kJ/mol), which explained the ATP dependence. A crystal structure showed that the enzyme exists as an (αβ)2 tetramer and revealed an ATP-binding site in CapA and a Zn-coordinating site in CapB. Mutations in the ATP-binding site of CapA (D11A and D295A) significantly reduced product formation. Mutants with substitutions in the metal binding site of CapB (D41A, H99A, D101A, and H124A) were inactive and less thermostable than the wild-type enzyme. These residues proved to be essential for activity and on basis of the experimental findings we propose possible mechanisms for ATP-dependent lactam hydrolysis.

Functional Effects In Silico Prediction for Androgen Receptor Ligand-Binding Domain Novel I836S Mutation

Over 1000 mutations are described in the androgen receptor (AR) gene. Of those, about 600 were found in androgen insensitivity syndrome (AIS) patients, among which 400 mutations affect the ligand-binding domain (LBD) of the AR protein. Recently, we reported a novel missense mutation c.2507T>G I836S (ClinVarID: 974911) in a patient with complete AIS (CAIS) phenotype. In the present study, we applied a set of computational approaches for the structural analysis of the ligand-binding domains in a wild-type and mutant AR to evaluate the functional impact of the novel I836S mutation. We revealed that the novel I836S substitution leads to a shorter existence time of the ligand's gating tunnel and internal cavity, occurring only in the presence of S836 phosphorylation. Additionally, the analysis of phosphorylation of the 836 mutant residues explained the negative impact on AR homodimerization, since monomer surface changes indirectly impacted the binding site. Our analyses provide evidence that I836S causes disruptions of AR protein functionality and development of CAIS clinical features in patients.

The role of positive charged residue in the proton-transfer mechanism of two-domain laccase from Streptomyces griseoflavus Ac-993

Multi-copper oxidases are capable of coupling the one-electron oxidation of four substrate equivalents to the four-electron reduction of dioxygen to two molecules of water. This process takes place at the trinuclear copper center of the enzymes. Previously, the main catalytic stages for three-domain (3D) laccases have been identified. However, for bacterial small two-domain (2D) laccases several questions remain to be answered. One of them is the nature of the protonation events upon the reductive cleavage of dioxygen to water. In 3D laccases, acidic residues play a key role in the protonation mechanisms. In this study, the role of the Arg240 residue, located within the T2 tunnel of 2D laccase from Streptomyces griseoflavus Ac-993, was investigated. X-ray structural analysis and kinetic characterization of two mutants, R240A and R240H, have provided support for a role of this residue in the protonation events. Communicated by Ramaswamy H. Sarma.

Development of the computational antibiotic screening platform (CLASP) to aid in the discovery of new antibiotics

Bacterial colonization of biotic and abiotic surfaces and antibiotic resistance are grand challenges with paramount societal impacts. However, in the face of increasing bacterial resistance to all known antibiotics, efforts to discover new classes of antibiotics have languished, creating an urgent need to accelerate the antibiotic discovery pipeline. A major deterrent in the discovering of new antibiotics is the limited permeability of molecules across the bacterial envelope. Notably, the Gram-negative bacteria have nutrient specific protein channels (or porins) that restrict the permeability of non-essential molecules, including antibiotics. Here, we have developed the Computational Antibiotic Screening Platform (CLASP) for screening of potential drug molecules through the porins. The CLASP takes advantage of coarse grain (CG) resolution, advanced sampling techniques, and a parallel computing environment to maximize its performance. The CLASP yields comprehensive thermodynamic and kinetic output data of a potential drug molecule within a few hours of wall-clock time. Its output includes the potential of mean force profile, energy barrier, the rate constant, and contact analysis of the molecule with the pore-lining residues, and the orientational analysis of the molecule in the porin channel. In our first CLASP application, we report the transport properties of six carbapenem antibiotics-biapenem, doripenem, ertapenem, imipenem, meropenem, and panipenem-through OccD3, a major channel for carbapenem uptake in Pseudomonas aeruginosa. The CLASP is designed to screen small molecule libraries with a fast turnaround time to yield structure-property relationships to discover antibiotics with high permeability. The CLASP will be freely distributed to enable accelerated antibiotic drug discovery.

Essential role of accessory subunit LYRM6 in the mechanism of mitochondrial complex I

Respiratory complex I catalyzes electron transfer from NADH to ubiquinone (Q) coupled to vectorial proton translocation across the inner mitochondrial membrane. Despite recent progress in structure determination of this very large membrane protein complex, the coupling mechanism is a matter of ongoing debate and the function of accessory subunits surrounding the canonical core subunits is essentially unknown. Concerted rearrangements within a cluster of conserved loops of central subunits NDUFS2 (β1-β2^S2 loop), ND1 (TMH5-6^ND1 loop) and ND3 (TMH1-2^ND3 loop) were suggested to be critical for its proton pumping mechanism. Here, we show that stabilization of the TMH1-2^ND3 loop by accessory subunit LYRM6 (NDUFA6) is pivotal for energy conversion by mitochondrial complex I. We determined the high-resolution structure of inactive mutant F89A^LYRM6 of eukaryotic complex I from the yeast Yarrowia lipolytica and found long-range structural changes affecting the entire loop cluster. In atomistic molecular dynamics simulations of the mutant, we observed conformational transitions in the loop cluster that disrupted a putative pathway for delivery of substrate protons required in Q redox chemistry. Our results elucidate in detail the essential role of accessory subunit LYRM6 for the function of eukaryotic complex I and offer clues on its redox-linked proton pumping mechanism.

Respiratory complex I plays a key role in energy metabolism. Cryo-EM structure of a mutant accessory subunit LYRM6 from the yeast Yarrowia lipolytica and molecular dynamics simulations reveal conformational changes at the interface between LYRM6 and subunit ND3, propagated further into the complex. These findings offer insight into the mechanism of proton pumping by respiratory complex I.

Multi-scale simulation reveals that an amino acid substitution increases photosensitizing reaction inputs in Rhodopsins

Evaluating the availability of molecular oxygen (O₂ ) and energy of excited states in the retinal binding site of rhodopsin is a crucial challenging first step to understand photosensitizing reactions in wild-type (WT) and mutant rhodopsins by absorbing visible light. In the present work, energies of the ground and excited states related to 11-cis-retinal and the O₂ accessibility to the β-ionone ring are evaluated inside WT and human M207R mutant rhodopsins. Putative O₂ pathways within rhodopsins are identified by using molecular dynamics simulations, Voronoi-diagram analysis, and implicit ligand sampling while retinal energetic properties are investigated through density functional theory, and quantum mechanical/molecular mechanical methods. Here, the predictions reveal that an amino acid substitution can lead to enough energy and O₂ accessibility in the core hosting retinal of mutant rhodopsins to favor the photosensitized singlet oxygen generation, which can be useful in understanding retinal degeneration mechanisms and in designing blue-lighting-absorbing proteic photosensitizers.

Cryo-EM Studies of TMEM16F Calcium-Activated Ion Channel Suggest Features Important for Lipid Scrambling

As a Ca²⁺-activated lipid scramblase and ion channel that mediates Ca²⁺ influx, TMEM16F relies on both functions to facilitate extracellular vesicle generation, blood coagulation, and bone formation. How a bona fide ion channel scrambles lipids remains elusive. Our structural analyses revealed the coexistence of an intact channel pore and PIP₂-dependent protein conformation changes leading to membrane distortion. Correlated to the extent of membrane distortion, many tightly bound lipids are slanted. Structure-based mutagenesis studies further reveal that neutralization of some lipid-binding residues or those near membrane distortion specifically alters the onset of lipid scrambling, but not Ca²⁺ influx, thus identifying features outside of channel pore that are important for lipid scrambling. Together, our studies demonstrate that membrane distortion does not require open hydrophilic grooves facing the membrane interior and provide further evidence to suggest separate pathways for lipid scrambling and ion permeation.

TMEM16F is a calcium-activated ion channel and lipid scramblase linked to the bleeding disorder Scott syndrome. Feng et al. examine cryo-EM structures of TMEM16F with or without Ca²⁺ ions and PIP₂ nanodisc supplementation and identify structural features for lipid binding and membrane distortion critical for lipid scrambling activity.