Protein-water hydrogen-bond networks of G protein-coupled receptors: Graph-based analyses of static structures and molecular dynamics

Rational Design of Formate Dehydrogenase for Enhanced Thermal Stability and Catalytic Activity in Bioelectrocatalysis

Formate dehydrogenase can be utilized as a biocatalyst in the bioelectrocatalysis of converting CO2 into formic acid. However, its industrial application has been hindered by limited thermal stability. This study successfully obtained a mutant (D533S/E684I) with enhanced thermal stability and catalytic activity through the rational design of flexible regions. The mutant exhibited a half-life (t1/2) 1.5 times longer than the wild type (WT) at 35 °C, along with a specific enzyme activity 7.46 times higher than that of the WT. Additionally, the catalytic efficiency (kcat/Km value) of the mutant toward the substrate was 2.72 s-1·mM-1, representing a 19.4-fold increase compared to the WT (0.14 s-1·mM-1). Formic acid production reached 53.4 mM through bioelectrocatalysis after 10 h, utilizing the mutant as the biocatalyst. Molecular dynamics simulations and structural analysis were employed to investigate the molecular mechanisms behind the enhanced thermal stability and activity. The displacement of a highly flexible region in the mutant may counteract the stability-activity trade-off. This study proposed a method for improving both thermal stability and activity in enzyme evolution.

Effect of Direct and Water-Mediated Interactions on the Identification of Hotspots in Biomolecular Complexes with Multiple Subsystems

Identification of important residues in biochemical complexes is often a crucial step for many problems in molecular biology and biochemistry. A method is proposed to identify hotspots in biomolecular complexes based on a new metric, derived from networks representing molecular subunits (residues, bridging solvent molecules, ligands etc.) connected by interactions. A singular value decomposition of the weighted adjacency matrix is used to construct a scalar rank for each subunit that reflects its importance in the residue interaction network. This metric is called the singular value centrality. In addition, a new formalism is proposed to account for water-mediated interactions in the ranking of residues. Interactions for a residue network can be provided by various computational methods. In this work interactions are obtained from full quantum-mechanical calculations of protein-protein complexes using the fragment molecular orbital method. The ranking results are shown to be in good agreement with earlier computational and experimental studies. The developed method can be used to gain a deeper insight into the role of subunits in complex biomolecular systems.

Structural basis of frizzled 7 activation and allosteric regulation

Frizzleds (ten paralogs: FZD1-10) belong to the class F of G protein-coupled receptors (GPCRs), which remains poorly understood despite its crucial role in multiple key biological functions including embryonic development, stem cell regulation, and homeostasis in the adult. FZD7, one of the most studied members of the family, is more specifically involved in the migration of mesendoderm cells during the development and renewal of intestinal stem cells in adults. Moreover, FZD7 has been highlighted for its involvement in tumor development predominantly in the gastrointestinal tract. This study reports the structure of inactive FZD7, without any stabilizing mutations, determined by cryo-electron microscopy (cryo-EM) at 1.9 Å resolution. We characterize a fluctuating water pocket in the core of the receptor important for FZD7 dynamics. Molecular dynamics simulations are used to investigate the temporal distribution of those water molecules and their importance for potential conformational changes in FZD7. Moreover, we identify lipids interacting with the receptor core and a conserved cholesterol-binding site, which displays a key role in FZD7 association with a transducer protein, Disheveled (DVL), and initiation of downstream signaling and signalosome formation.

Searching for proton transfer channels in respiratory complex I

We have explored a strategy to identify potential proton transfer channels using computational analysis of a protein structure based on Voronoi partitioning and applied it for the analysis of proton transfer pathways in redox-driven proton-pumping respiratory complex I. The analysis results in a network of connected voids/channels, which represent the dual structure of the protein; we then hydrated the identified channels using our water placement program Dowser++. Many theoretical water molecules found in the channels perfectly match the observed experimental water molecules in the structure; some other predicted water molecules have not been resolved in the experiments. The channels are of varying cross sections. Some channels are big enough to accommodate water molecules that are suitable to conduct protons; others are too narrow to hold water but require only minor conformational changes to accommodate proton transfer. We provide a preliminary analysis of the proton conductivity of the network channels, classifying the proton transfer channels as open, closed, and partially open, and discuss possible conformational changes that can modulate, i.e., open and close, the channels.

Hydrogen-Bonding and Hydrophobic Interaction Networks as Structural Determinants of Microbial Rhodopsin Function

Microbial pump rhodopsins are highly versatile light-driven membrane proteins that couple protein conformational dynamics with ion translocation across the cell membranes. Understanding how microbial pump rhodopsins use specific amino acid residues at key functional sites to control ion selectivity and ion pumping direction is of general interest for membrane transporters, and could guide site-directed mutagenesis for optogenetics applications. To enable direct comparisons between proteins with different sequences we implement, for the first time, a unique numbering scheme for the microbial pump rhodopsin residues, NS-mrho. We use NS-mrho to show that distinct microbial pump rhodopsins typically have hydrogen-bond networks that are less conserved than anticipated from the amino acid residue conservation, whereas their hydrophobic interaction networks are largely conserved. To illustrate the role of the hydrogen-bond networks as structural elements that determine the functionality of microbial pump rhodopsins, we performed experiments, atomic-level simulations, and hydrogen bond network analyses on GR, the outward proton pump from Gloeobacter violaceus, and KR2, the outward sodium pump from Krokinobacter eikastus. The experiments indicate that multiple mutations that recover KR2 amino acid residues in GR not only fail to convert it into a sodium pump, but completely inactivate GR by abolishing photoisomerization of the retinal chromophore. This observation could be attributed to the drastically altered hydrogen-bond interaction network identified with simulations and network analyses. Taken together, our findings suggest that functional specificity could be encoded in the collective hydrogen-bond network of microbial pump rhodopsins.

Graph-based algorithms to dissect long-distance water-mediated H-bond networks for conformational couplings in GPCRs

Changes in structure and dynamics elicited by agonist ligand binding at the extracellular side of G protein coupled receptors (GPCRs) must be relayed to the cytoplasmic G protein binding side of the receptors. To decipher the role of water-mediated hydrogen-bond networks in this relay mechanism, we have developed graph-based algorithms and analysis methodologies applicable to datasets of static structures of distinct GPCRs. For a reference dataset of static structures of bovine rhodopsin solved at the same resolution, we show that graph analyses capture the internal protein-water hydrogen-bond network. The extended analyses of static structures of rhodopsins and opioid receptors suggest a relay mechanism whereby inactive receptors have in place much of the internal core hydrogen-bond network required for long-distance relay of structural change, with extensive local H-bond clusters observed in structures solved at high resolution and with internal water molecules.

Improved Protein Dynamics and Hydration in the Martini3 Coarse-Grain Model

The Martini coarse-grain force-field has emerged as an important framework to probe cellular processes at experimentally relevant time- and length-scales. However, the recently developed version, the Martini3 force-field with the implemented Go̅ model (Martini3Go̅), as well as previous variants of the Martini model have not been benchmarked and rigorously tested for globular proteins. In this study, we consider three globular proteins, ubiquitin, lysozyme, and cofilin, and compare protein dynamics and hydration with observables from experiments and all-atom simulations. We show that the Martini3Go̅ model is able to accurately model the structural and dynamic features of small globular proteins. Overall, the structural integrity of the proteins is maintained, as validated by contact maps, radii of gyration (Rg), and SAXS profiles. The chemical shifts predicted from the ensemble sampled in the simulations are consistent with the experimental data. Further, a good match is observed in the protein-water interaction energetics, and the hydration levels of the residues are similar to atomistic simulations. However, the protein-water interaction dynamics is not accurately represented and appears to depend on the protein structural complexity, residue specificity, and water dynamics. Our work is a step toward testing and assessing the Martini3Go̅ model and provides insights into future efforts to refine Martini models with improved solvation effects and better correspondence to the underlying all-atom systems.

Intra-protein interacting collective modes in the terahertz frequency region

Understanding how proteins work requires a thorough understanding of their internal dynamics. Proteins support a wide range of motions, from the femtoseconds to seconds time scale, relevant to crucial biological functions. In this context, the term "protein collective dynamics" refers to the complex patterns of coordinated motions of numerous atoms throughout the protein in the sub-picosecond time scale (terahertz frequency region). It is hypothesized that these dynamics have a substantial impact on the regulation of functional dynamical mechanisms, including ligand binding and allosteric signalling, charge transport direction, and the regulation of thermodynamic and thermal transport properties. Using the theoretical framework of hydrodynamics, the collective dynamics of proteins had previously been described in a manner akin to that of simple liquids, i.e. in terms of a single acoustic-like excitation, related to intra-protein vibrational motions. Here, we employ an interacting-mode model to analyse the results from molecular dynamics simulations and we unveil that the vibrational landscape of proteins is populated by multiple acoustic-like and low-frequency optic-like modes, with mixed symmetry and interfering with each other. We propose an interpretation at the molecular level of the observed scenario that we relate to the side-chains and the hydrogen-bonded networks dynamics. The present insights provide a perspective for understanding the molecular mechanisms underlying the energy redistribution processes in the interior of proteins.

Quantum Effects of Hydrogen Nuclei on the Nuclear Magnetic Shielding Tensor of Ice Ih

Ice is the most fundamental hydrogen-bonded system in which the hydrogen nuclear quantum effect significantly impacts the structure and relevant thermochemical and spectroscopic properties. While ice was experimentally investigated using proton nuclear magnetic resonance spectroscopy more than 40 years ago, the corresponding theoretical investigations have been rarely reported due to the difficulty in evaluating how the proton nuclear quantum effect influences the spectral characteristics of such a condensed material. In this study, we applied a combination of the ONIOM and multicomponent molecular orbital (MC_MO) methods for calculating the anisotropic and isotropic components of the nuclear magnetic shielding tensor of the hexagonal ice crystal to quantify the effects of nuclear quantum fluctuations on the spectroscopic properties of ice. The nuclear magnetic shielding values computed by incorporating the hydrogen nuclear quantum effect reasonably agree with the experimental values. The nuclear quantum effects were found to increase the anisotropic component of the magnetic shielding tensor while decreasing the isotropic component. Such a difference can be explained by their distinct dependence on the electrostatic field and hydrogen-bonding structural parameters.

Dynamics of Hydrogen Bonds between Water and Intrinsically Disordered and Structured Regions of Proteins

Recent studies indicate more restricted dynamics of water around intrinsically disordered proteins (IDPs) than structured proteins. We examine here the dynamics of hydrogen bonds between water molecules and two proteins, small ubiquitin-related modifier-1 (SUMO-1) and ubiquitin-conjugating enzyme E2I (UBC9), which we compare around intrinsically disordered regions (IDRs) and structured regions of these proteins. It has been recognized since some time that excluded volume effects, which influence access of water molecules to hydrogen-bonding sites, and the strength of hydrogen bonds between water and protein affect hydrogen bond lifetimes. While we find those two properties to mediate lifetimes of hydrogen bonds between water and protein residues in this study, we also find that the lifetimes are affected by the concentration of charged groups on other nearby residues. These factors are more important in determining the hydrogen bond lifetimes than whether a residue hydrogen bonding with water belongs to an IDR or to a structured region.

Protons taken hostage: Dynamic H-bond networks of the pH-sensing GPR68

Proton-sensing G Protein Coupled Receptors (GPCRs) sense changes in the extracellular pH to effect cell signaling for cellular homeostasis. They tend to be overexpressed in solid tumors associated with acidic extracellular pH, and are of direct interest as drug targets. How proton-sensing GPCRs sense extracellular acidification and activate upon protonation change is important to understand, because it may guide the design of therapeutics. Lack of publicly available experimental structures make it challenging to discriminate between conflicting mechanisms proposed for proton-binding, as main roles have been assigned to either an extracellular histidine cluster or to an internal carboxylic triad. Here we present a protocol to derive and evaluate structural models of the proton-sensing GPR68. This approach integrates state-of-the-art homology modeling with microsecond-timescale atomistic simulations, and with a detailed assessment of the compatibility of the structural models with known structural features of class A GPCRs. To decipher structural elements of potential interest for protonation-coupled conformational changes of GPR68, we used the best-compatible model as a starting point for independent atomistic simulations of GPR68 with different protonation states, and graph computations to characterize the response of GPR68 to changes in protonation. We found that GPR68 hosts an extended hydrogen-bond network that inter-connects the extracellular histidine cluster to the internal carboxylic triad, and which can even reach groups at the cytoplasmic G-protein binding site. Taken together, results suggest that GPR68 relies on dynamic, hydrogen-bond networks to inter-connect extracellular and internal proton-binding sites, and to elicit conformational changes at the cytoplasmic G-protein binding site.

Fentanyl and the Fluorinated Fentanyl Derivative NFEPP Elicit Distinct Hydrogen-Bond Dynamics of the Opioid Receptor

The development of safe therapeutics to manage pain is of central interest for biomedical applications. The fluorinated fentanyl derivative N-(3-fluoro-1-phenethylpiperidin-4-yl)-N-phenylpropionamide (NFEPP) is potentially a safer alternative to fentanyl because unlike fentanyl─which binds to the μ-opioid receptor (MOR) at both physiological and acidic pH─NFEPP might bind to the MOR only at acidic pH typical of inflamed tissue. Knowledge of the protonation-coupled dynamics of the receptor-drug interactions is thus required to understand the molecular mechanism by which receptor activation initiates cell signaling to silence pain. To this end, here we have carried out extensive atomistic simulations of the MOR in different protonation states, in the absence of opioid drugs, and in the presence of fentanyl vs NFEPP. We used graph-based analyses to characterize internal hydrogen-bond networks that could contribute to the activation of the MOR. We find that fentanyl and NFEPP prefer distinct binding poses and that, in their binding poses, fentanyl and NFEPP partake in distinct internal hydrogen-bond networks, leading to the cytoplasmic G-protein-binding region. Moreover, the protonation state of functionally important aspartic and histidine side chains impacts hydrogen-bond networks that extend throughout the receptor, such that the ligand-bound MOR presents at its cytoplasmic G-protein-binding side, a hydrogen-bonding environment where dynamics depend on whether fentanyl or NFEPP is bound, and on the protonation state of specific MOR groups. The exquisite sensitivity of the internal protein-water hydrogen-bond network to the protonation state and to details of the drug binding could enable the MOR to elicit distinct pH- and opioid-dependent responses at its cytoplasmic G-protein-binding site.

Computational investigation of functional water molecules in GPCRs bound to G protein or arrestin

G protein-coupled receptors (GPCRs) are membrane proteins constituting the largest family of drug targets. The activated GPCR binds either the heterotrimeric G proteins or arrestin through its activation cycle. Water molecules have been reported to play a role in GPCR activation. Nevertheless, reported studies are focused on the hydrophobic helical bundle region. How water molecules function in GPCR bound either G protein or arrestin is rarely studied. To address this issue, we carried out computational studies on water molecules in both GPCR/G protein complexes and GPCR/arrestin complexes. Using inhomogeneous fluid theory (IFT), we locate all possible hydration sites in GPCRs binding either to G protein or arrestin. We observe that the number of water molecules on the interaction surface between GPCRs and signal proteins are correlated with the insertion depths of the α5-helix from G-protein or "finger loop" from arrestin in GPCRs. In three out of the four simulation pairs, the interfaces of Rhodopsin, M₂R and NTSR1 in the G protein-associated systems show more water-mediated hydrogen-bond networks when compared to these in arrestin-associated systems. This reflects that more functionally relevant water molecules may probably be attracted in G protein-associated structures than that in arrestin-associated structures. Moreover, we find the water-mediated interaction networks throughout the NPxxY region and the orthosteric pocket, which may be a key for GPCR activation. Reported studies show that non-biased agonist, which can trigger both GPCR-G protein and GPCR-arrestin activation signal, can result in pharmacologically toxicities. Our comprehensive studies of the hydration sites in GPCR/G protein complexes and GPCR/arrestin complexes may provide important insights in the design of G-protein biased agonists.

Structural elucidation of ETHR-A and ETHR-B from Plutella xylostella and insight into non-conservative mutations in transmembrane helix-6

The development of Diamondback moth (DBM) depends on the ecdysis triggering hormone receptor (ETHR); a neuronal membrane G-protein coupled receptor (GPCR) connected to the metamorphosis cascade. Lepidopteran insect DBM is an infamous pest of cruciferous plants. This study examined the full-length coding sequences (CDS) of PxETHR-A and PxETHR-B from the DBM genome. The three-dimensional (3 D) models of both receptors were generated in an inactive state. The behaviour and stability of receptors were examined using molecular dynamics simulations in a lipid membrane system for 300 ns and established a GPCR family-based view. Secondary interactions within receptors were studied to know more about factors contributing to their stability. Multiple sequence alignment revealed conserved features of insect ETHRs those compared to the GPCR family proteins. These features were helpful during the evaluation of the molecular models of both receptors. Side-chain orientation of conserved residues, non-conserved and conserved hydrogen-bond networks (HBN) and hydrophobic clusters were examined in the structures of both receptors. The non-conserved residues L6.35, T6.39, C/S6.43, and L6.48, are present in a conserved position on the transmembrane helix-6 (TM6) of ETHRs. In TM6, PxETHR-A and PxETHR-B differ at positions C/S6.43 and Y/F6.51, both being part of the HBN.Communicated by Ramaswamy H. Sarma.

Graphs of protein-water hydrogen bond networks to dissect structural movies of ion-transfer microbial rhodopsins

Microbial rhodopsins are membrane proteins that use the energy absorbed by the covalently bound retinal chromophore to initiate reaction cycles resulting in ion transport or signal transduction. Thousands of distinct microbial rhodopsins are known and, for many rhodopsins, three-dimensional structures have been solved with structural biology, including as entire sets of structures solved with serial femtosecond crystallography. This sets the stage for comprehensive studies of large datasets of static protein structures to dissect structural elements that provide functional specificity to the various microbial rhodopsins. A challenge, however, is how to analyze efficiently intra-molecular interactions based on large datasets of static protein structures. Our perspective discusses the usefulness of graph-based approaches to dissect structural movies of microbial rhodopsins solved with time-resolved crystallography.

Energy Transport in Class B GPCRs: Role of Protein-Water Dynamics and Activation

We compute energy exchange networks (EENs) through glucagon-like peptide-1 receptor (GLP-1R), a class B G-protein-coupled receptor (GPCR), in inactive and two active states, one activated by a peptide ligand and the other by a small molecule agonist, from results of molecular dynamics simulations. The reorganized network upon activation contains contributions from structural as well as from dynamic changes and corresponding entropic contributions to the free energy of activation, which are estimated in terms of the change in rates of energy transfer across non-covalent contacts. The role of water in the EENs and in activation of GLP-1R is also investigated. The dynamics of water in contact with the central polar network of the transmembrane region is found to be significantly slower for both activated states compared to the inactive state. This result is consistent with the contribution of water molecules to activation of GLP-1R previously suggested and resembles water dynamics in parts of the transmembrane region found in earlier studies of rhodopsin-like GPCRs.

Structural basis for receptor selectivity and inverse agonism in S1P5 receptors

The bioactive lysophospholipid sphingosine-1-phosphate (S1P) acts via five different subtypes of S1P receptors (S1PRs) - S1P1-5. S1P5 is predominantly expressed in nervous and immune systems, regulating the egress of natural killer cells from lymph nodes and playing a role in immune and neurodegenerative disorders, as well as carcinogenesis. Several S1PR therapeutic drugs have been developed to treat these diseases; however, they lack receptor subtype selectivity, which leads to side effects. In this article, we describe a 2.2 Å resolution room temperature crystal structure of the human S1P5 receptor in complex with a selective inverse agonist determined by serial femtosecond crystallography (SFX) at the Pohang Accelerator Laboratory X-Ray Free Electron Laser (PAL-XFEL) and analyze its structure-activity relationship data. The structure demonstrates a unique ligand-binding mode, involving an allosteric sub-pocket, which clarifies the receptor subtype selectivity and provides a template for structure-based drug design. Together with previously published S1PR structures in complex with antagonists and agonists, our structure with S1P5-inverse agonist sheds light on the activation mechanism and reveals structural determinants of the inverse agonism in the S1PR family.

Development of fluorine-substituted NH2-biphenyl-diarylpyrimidines as highly potent non-nucleoside reverse transcriptase inhibitors: Boosting the safety and metabolic stability

Our recent studies for nonnucleoside reverse transcriptase inhibitors identified a highly potent compound JK-4b against WT HIV-1 (EC₅₀ = 1.0 nmol/L), but the poor metabolic stability in human liver microsomes (t _1/2 = 14.6 min) and insufficient selectivity (SI = 2059) with high cytotoxicity (CC₅₀ = 2.08 μmol/L) remained major issues associated with JK-4b. The present efforts were devoted to the introduction of fluorine into the biphenyl ring of JK-4b, leading to the discovery of a novel series of fluorine-substituted NH₂-biphenyl-diarylpyrimidines with noticeable inhibitory activity toward WT HIV-1 strain (EC₅₀ = 1.8-349 nmol/L). The best compound 5t in this collection (EC₅₀ = 1.8 nmol/L, CC₅₀ = 117 μmol/L) was 32-fold in selectivity (SI = 66,443) compared to JK-4b and showed remarkable potency toward clinically multiple mutant strains, such as L100I, K103N, E138K, and Y181C. The metabolic stability of 5t was also significantly improved (t _1/2 = 74.52 min), approximately 5-fold higher than JK-4b in human liver microsomes (t _1/2 = 14.6 min). Also, 5t possessed good stability in both human and monkey plasma. No significant in vitro inhibition effect toward CYP enzyme and hERG was observed. The single-dose acute toxicity test did not induce mice death or obvious pathological damage. These findings pave the way for further development of 5t as a drug candidate.

Roles of interfacial water states on advanced biomedical material design

When biomedical materials come into contact with body fluids, the first reaction that occurs on the material surface is hydration; proteins are then adsorbed and denatured on the hydrated material surface. The amount and degree of denaturation of adsorbed proteins affect subsequent cell behavior, including cell adhesion, migration, proliferation, and differentiation. Biomolecules are important for understanding the interactions and biological reactions of biomedical materials to elucidate the role of hydration in biomedical materials and their interaction partners. Analysis of the water states of hydrated materials is complicated and remains controversial; however, knowledge about interfacial water is useful for the design and development of advanced biomaterials. Herein, we summarize recent findings on the hydration of synthetic polymers, supramolecular materials, inorganic materials, proteins, and lipid membranes. Furthermore, we present recent advances in our understanding of the classification of interfacial water and advanced polymer biomaterials, based on the intermediate water concept.

Graphs of Hydrogen-Bond Networks to Dissect Protein Conformational Dynamics

Dynamic hydrogen bonds and hydrogen-bond networks are ubiquitous in proteins and protein complexes. Functional roles that have been assigned to hydrogen-bond networks include structural plasticity for protein function, allosteric conformational coupling, long-distance proton transfers, and transient storage of protons. Advances in structural biology provide invaluable insights into architectures of large proteins and protein complexes of direct interest to human physiology and disease, including G Protein Coupled Receptors (GPCRs) and the SARS-Covid-19 spike protein S, and give rise to the challenge of how to identify those interactions that are more likely to govern protein dynamics. This Perspective discusses applications of graph-based algorithms to dissect dynamical hydrogen-bond networks of protein complexes, with illustrations for GPCRs and spike protein S. H-bond graphs provide an overview of sites in GPCR structures where hydrogen-bond dynamics would be required to assemble longer-distance networks between functionally important motifs. In the case of spike protein S, graphs identify regions of the protein where hydrogen bonds rearrange during the reaction cycle and where local hydrogen-bond networks likely change in a virus variant of concern.

Conserved hydrogen-bond motifs of membrane transporters and receptors

Membrane transporters and receptors often rely on conserved hydrogen bonds to assemble transient paths for ion transfer or long-distance conformational couplings. For transporters and receptors that use proton binding and proton transfer for function, inter-helical hydrogen bonds of titratable protein sidechains that could change protonation are of central interest to formulate hypotheses about reaction mechanisms. Knowledge of hydrogen bonds common at sites of potential interest for proton binding could thus inform and guide studies on functional mechanisms of protonation-coupled membrane proteins. Here we apply graph-theory approaches to identify hydrogen-bond motifs of carboxylate and histidine sidechains in a large data set of static membrane protein structures. We find that carboxylate-hydroxyl hydrogen bonds are present in numerous structures of the dataset, and can be part of more extended H-bond clusters that could be relevant to conformational coupling. Carboxylate-carboxyamide and imidazole-imidazole hydrogen bonds are represented in comparably fewer protein structures of the dataset. Atomistic simulations on two membrane transporters in lipid membranes suggest that many of the hydrogen bond motifs present in static protein structures tend to be robust, and can be part of larger hydrogen-bond clusters that recruit additional hydrogen bonds.

True-atomic-resolution insights into the structure and functional role of linear chains and low-barrier hydrogen bonds in proteins

Hydrogen bonds are fundamental to the structure and function of biological macromolecules and have been explored in detail. The chains of hydrogen bonds (CHBs) and low-barrier hydrogen bonds (LBHBs) were proposed to play essential roles in enzyme catalysis and proton transport. However, high-resolution structural data from CHBs and LBHBs is limited. The challenge is that their 'visualization' requires ultrahigh-resolution structures of the ground and functionally important intermediate states to identify proton translocation events and perform their structural assignment. Our true-atomic-resolution structures of the light-driven proton pump bacteriorhodopsin, a model in studies of proton transport, show that CHBs and LBHBs not only serve as proton pathways, but also are indispensable for long-range communications, signaling and proton storage in proteins. The complete picture of CHBs and LBHBs discloses their multifunctional roles in providing protein functions and presents a consistent picture of proton transport and storage resolving long-standing debates and controversies.

Comparison of proton transfer paths to the QA and QB sites of the Rb. sphaeroides photosynthetic reaction centers

The photosynthetic bacterial reaction centers from purple non-sulfur bacteria use light energy to drive the transfer of electrons from cytochrome c to ubiquinone. Ubiquinone bound in the Q_A site cycles between quinone, Q_A, and anionic semiquinone, Q_A^·-, being reduced once and never binding protons. In the Q_B site, ubiquinone is reduced twice by Q_A^·-, binds two protons and is released into the membrane as the quinol, QH₂. The network of hydrogen bonds formed in a molecular dynamics trajectory was drawn to investigate proton transfer pathways from the cytoplasm to each quinone binding site. Q_A is isolated with no path for protons to enter from the surface. In contrast, there is a complex and tangled network requiring residues and waters that can bring protons to Q_B. There are three entries from clusters of surface residues centered around HisH126, GluH224, and HisH68. The network is in good agreement with earlier studies, Mutation of key nodes in the network, such as SerL223, were previously shown to slow proton delivery. Mutational studies had also shown that double mutations of residues such as AspM17 and AspL210 along multiple paths in the network presented here slow the reaction, while single mutations do not. Likewise, mutation of both HisH126 and HisH128, which are at the entry to two paths reduce the rate of proton uptake.

Investigating the binding mechanism of topiramate with bovine serum albumin using spectroscopic and computational methods

Various spectroscopic techniques involving fluorescence spectroscopy, circular dichroism (CD), and computational approaches were used to elucidate the molecular aspects of interaction between the antiepileptic drug topiramate and the multifunctional transport protein bovine serum albumin (BSA) under physiological conditions. Topiramate quenched BSA fluorescence in a static quenching mode, according to the Stern-Volmer quenching constant (K_sv ) data derived from fluorescence spectroscopy for the topiramate-BSA complex. The binding constant was also used to calculate the binding affinity for the topiramate-BSA interaction. Fluorescence and circular dichroism experiments demonstrate that the protein's tertiary structure is affected by the microenvironmental alterations generated by topiramate binding to BSA. To establish the exact binding site, interacting residues, and interaction forces involved in the binding of topiramate to BSA, molecular modeling and simulation approaches were used. According to the MMPBSA calculations, the average binding energy between topiramate and BSA is -421.05 kJ/mol. Topiramate was discovered to have substantial interactions with BSA, changing the structural dynamic and Gibbs free energy landscape patterns.

Opioid Receptors and Protonation-Coupled Binding of Opioid Drugs

Opioid receptors are G-protein-coupled receptors (GPCRs) part of cell signaling paths of direct interest to treat pain. Pain may associate with inflamed tissue characterized by acidic pH. The potentially low pH at tissue targeted by opioid drugs in pain management could impact drug binding to the opioid receptor, because opioid drugs typically have a protonated amino group that contributes to receptor binding, and the functioning of GPCRs may involve protonation change. In this review, we discuss the relationship between structure, function, and dynamics of opioid receptors from the perspective of the usefulness of computational studies to evaluate protonation-coupled opioid-receptor interactions.

C-Graphs Tool with Graphical User Interface to Dissect Conserved Hydrogen-Bond Networks: Applications to Visual Rhodopsins

Dynamic hydrogen-bond networks provide proteins with structural plasticity required to translate signals such as ligand binding into a cellular response or to transport ions and larger solutes across membranes and, thus, are of central interest to understand protein reaction mechanisms. Here, we present C-Graphs, an efficient tool with graphical user interface that analyzes data sets of static protein structures or of independent numerical simulations to identify conserved, vs unique, hydrogen bonds and hydrogen-bond networks. For static structures, which may belong to the same protein or to proteins with different sequences, C-Graphs uses a clustering algorithm to identify sites of the hydrogen-bond network where waters are conserved among the structures. Using C-Graphs, we identify an internal protein-water hydrogen-bond network common to static structures of visual rhodopsins and adenosine A2A G protein-coupled receptors (GPCRs). Molecular dynamics simulations of a visual rhodopsin indicate that the conserved hydrogen-bond network from static structure can recruit dynamic hydrogen bonds and extend throughout most of the receptor. We release with this work the code for C-Graphs and its graphical user interface.

Molecular Dynamics Simulations of the Aptamer Domain of Guanidinium Ion Binding Riboswitch ykkC-III: Structural Insights into the Discrimination of Cognate and Alternate Ligands

Guanidinium ion is a toxic cellular metabolite. The ykkC-III riboswitch, an mRNA stretch, regulates the gene expression by undergoing a conformational change in response to the binding of a free guanidinium ion and thereby plays a potentially important role in alleviating guanidinium toxicity in cells. An experimental crystal structure of the guanidinium-bound aptamer domain of the riboswitch from Thermobifida Fusca revealed the overall RNA architecture and mapped the specific noncovalent interactions that stabilize the ligand within the binding pocket aptamer. However, details of how the aptamer domain discriminates the cognate ligand from its closest structurally analogous physiological metabolites (arginine and urea), and how the binding of cognate ligand arrays information from the aptamer domain to the expression platform for regulating the gene expression, are not well understood. To fill this void, we perform a cumulative of 2 μs all-atom explicit-solvent molecular dynamics (MD) simulations on the full aptamer domain, augmented with quantum-chemical calculations on the ligand-binding pocket, to compare the structural and dynamical details of the guanidinium-bound state with the arginine or urea bound states, as well as the unbound (open) state. Analysis of the ligand-binding pocket reveals that due to unfavorable interactions with the binding-pocket residues, urea cannot bind the aptamer domain and thereby cannot alter the gene expression. Although interaction of the guanidyl moiety of arginine within the binding pocket is either comparable or stronger than the guanidinium ion, additional non-native hydrogen-bonding networks, as well as differences in the dynamical details of the arginine-bound state, explain why arginine cannot transmit the information from the aptamer domain to the expression platform. Based on our simulations, we propose a mechanism of how the aptamer domain communicates with the expression platform. Overall, our work provides interesting insights into the ligand recognition by a specific class of riboswitches and may hopefully inspire future studies to further understand the gene regulation by riboswitches.

Hydration Dynamics in Biological Membranes: Emerging Applications of Terahertz Spectroscopy

Water drives the spontaneous self-assembly of lipids and proteins into quasi two-dimensional biological membranes that act as catalytic scaffolds for numerous processes central to life. However, the functional relevance of hydration in membrane biology is only beginning to be addressed, predominantly because of challenges associated with direct measurements of hydration microstructure and dynamics in a biological milieu. Our recent work on the novel interplay of membrane electrostatics and crowding in shaping membrane hydration dynamics utilizing terahertz (THz) spectroscopy represents an important step in this context. In this Perspective, we provide a glimpse into the ever-broadening functional landscape of hydration dynamics in biological membranes in the backdrop of the unique physical chemistry of water molecules. We further highlight the immense (and largely untapped) potential of the THz toolbox in addressing contemporary problems in membrane biology, while emphasizing the adaptability of the analytical framework reported recently by us to such studies.

Concerted conformational dynamics and water movements in the ghrelin G protein-coupled receptor

There is increasing support for water molecules playing a role in signal propagation through G protein-coupled receptors (GPCRs). However, exploration of the hydration features of GPCRs is still in its infancy. Here, we combined site-specific labeling with unnatural amino acids to molecular dynamics to delineate how local hydration of the ghrelin receptor growth hormone secretagogue receptor (GHSR) is rearranged upon activation. We found that GHSR is characterized by a specific hydration pattern that is selectively remodeled by pharmacologically distinct ligands and by the lipid environment. This process is directly related to the concerted movements of the transmembrane domains of the receptor. These results demonstrate that the conformational dynamics of GHSR are tightly coupled to the movements of internal water molecules, further enhancing our understanding of the molecular bases of GPCR-mediated signaling.

CHARMM Force-Field Parameters for Morphine, Heroin, and Oliceridine, and Conformational Dynamics of Opioid Drugs

Opioid drug binding to specialized G protein-coupled receptors (GPCRs) can lead to analgesia upon activation via downstream G_i protein signaling and to severe side effects via activation of the β-arrestin signaling pathway. Knowledge of how different opioid drugs interact with receptors is essential, as it can inform and guide the design of safer therapeutics. We performed quantum and classical mechanical computations to explore the potential energy landscape of four opioid drugs: morphine and its derivatives heroin and fentanyl and for the unrelated oliceridine. From potential energy profiles for bond twists and from interactions between opioids and water, we derived a set of force-field parameters that allow a good description of structural properties and intermolecular interactions of the opioids. Potential of mean force profiles computed from molecular dynamics simulations indicate that fentanyl and oliceridine have complex energy landscapes with relatively small energy penalties, suggesting that interactions with the receptor could select different binding poses of the drugs.

Proton-Binding Motifs of Membrane-Bound Proteins: From Bacteriorhodopsin to Spike Protein S

Membrane-bound proteins that change protonation during function use specific protein groups to bind and transfer protons. Knowledge of the identity of the proton-binding groups is of paramount importance to decipher the reaction mechanism of the protein, and protonation states of prominent are studied extensively using experimental and computational approaches. Analyses of model transporters and receptors from different organisms, and with widely different biological functions, indicate common structure-sequence motifs at internal proton-binding sites. Proton-binding dynamic hydrogen-bond networks that are exposed to the bulk might provide alternative proton-binding sites and proton-binding pathways. In this perspective article I discuss protonation coupling and proton binding at internal and external carboxylate sites of proteins that use proton transfer for function. An inter-helical carboxylate-hydroxyl hydrogen-bond motif is present at functionally important sites of membrane proteins from archaea to the brain. External carboxylate-containing H-bond clusters are observed at putative proton-binding sites of protonation-coupled model proteins, raising the question of similar functionality in spike protein S.

Interactive Interface for Graph-Based Analyses of Dynamic H-Bond Networks: Application to Spike Protein S

Dynamic hydrogen-bond networks are key determinants of protein conformational dynamics. In the case of macromolecular protein complexes, which can have a large number of hydrogen bonds giving rise to extensive hydrogen-bond networks, efficient algorithms are required to analyze interactions that could be important for the dynamics and biological function of the complex. We present here a highly efficient, standalone interface designed for analyses of dynamical hydrogen-bond networks of biomolecules and macromolecular complexes. To facilitate a comprehensive description of protein dynamics, the interface includes analyses of hydrophobic interactions. We illustrate the usefulness and workflow of the interface by dissecting the dynamics of the ectodomain of SARS-CoV-2 protein S in its closed conformation. We find that protein S contains numerous local clusters of dynamic hydrogen bonds and identify hydrogen bonds that are sampled persistently. The receptor binding domain of the spike protein hosts only a handful of persistent hydrogen-bond clusters, suggesting structural plasticity. Our data analysis interface is released here for open use.

Molecular elucidation of pancreatic elastase inhibition by baicalein

The serine protease, elastase exists in various forms and plays diverse roles in the human body. Pharmacological inhibition of elastase has been investigated for its therapeutic role in managing conditions such as diabetes, pneumonia and arthritis. Sivelestat, a synthetic molecule, is the only elastase inhibitor to have been approved by any major drug regulatory authority (PMDA, in this case) - but still has failed to attain widespread clinical usage owing to its high price, cumbersome administration and obscure long-term safety profile. In order to find a relatively better-suited alternative, screening was conducted using plant flavonoids, which yielded baicalein, a molecule that showed robust inhibition against Pancreatic Elastase inhibition (IC₅₀: 3.53 μM). Other than having a considerably lower IC₅₀than sivelestat, baicalein is also cheaper, safer and easier to administer. While MicroScale Thermophoresis validated baicalein-elastase interaction, enzyme-kinetic studies, molecular docking and molecular dynamic simulation revealed the mode of inhibition to be non-competitive. Baicalein exhibited binding to a distinct allosteric site on the enzyme. The current study demonstrates the elastase inhibition properties of baicalein in an in-vitro and in-silico environment.Communicated by Ramaswamy H. Sarma.

Biochemical characterization and molecular docking analysis of novel esterases from Sphingobium chungbukense DJ77

We identified three novel microbial esterase (Est1, Est2, and Est3) from Sphingobium chungbukense DJ77. Multiple sequence alignment showed the Est1 and Est3 have distinct motifs, such as tetrapeptide motif HGGG, a pentapeptide sequence motif GXSXG, and catalytic triad residues Ser-Asp-His, indicating that the identified enzymes belong to family IV esterases. Interestingly, Est1 exhibited strong activity toward classical esterase substrates, p-nitrophenyl ester of short-chain fatty acids and long-chain. However, Est3 did not exhibit any activity despite having high sequence similarity and sharing the identical catalytic active residues with Est1. Est3 only showed hydrolytic degradation activity to polycaprolactone (PCL). MOE-docking prediction also provided the parameters consisting of binding energy, molecular docking score, and molecular distance between substrate and catalytic nucleophilic residue, serine. The engineered mutEst3 has hydrolytic activity for a variety of esters ranging from p-nitrophenyl esters to PCL. In the present study, we demonstrated that MOE-docking simulation provides a valuable insight for facilitating biocatalytic performance.