Bridge: A Graph-Based Algorithm to Analyze Dynamic H-Bond Networks in Membrane Proteins

Natural Products in Precision Neurological Disease (Cryptococcal Meningitis): Structure-Based Phytochemical Screening of Glycyrrhiza glabra Plant Against Cryptococcus neoformans Farnesyltransferase (FTase)

Cryptococcus neoformans causes cryptococcal meningitis, which is lethal to immune-compromised people, especially AIDS patients. This study employed diverse in silico techniques to find the best phytochemical to block farnesyltransferase (FTase). Based on molecular docking, the top two compounds selected from a screening of 5807 phytochemical compounds from 29 medicinal plants were CID_8299 (hydroxyacetone) and CID_71346280 (1,7-bis (4-hydroxyphenyl)-1,4,6-heptatrien-3-one), with docking scores of -5.786 and -0.078 kcal/mol, respectively, indicating stronger binding affinities than the control CID_3365 (fluconazole), which scored -4.2 kcal/mol. The control and lead compounds bind at the common active site of protein by interacting with common amino acid residues (HIS97, GLN408, PHE93, and TRP94). Post-docking MM-GBSA verified docking score where CID_8299 and CID_71346280 had negative binding free energies of -19.81 and -0.27 kcal/mol, respectively. These two lead compounds were reassessed through molecular dynamics simulation (100 ns), and several post-dynamics analyses were conducted. CID_71346280, 8299, and 3365 (control) showed average RSMD values of 3.17, 1.904, and 2.08; average root mean square fluctuation values of 1.167, 0.886, and 1.028 Å; average radius of gyration values of 5.13, 1.58, and 3.54 Å; average solvent accessible surface area values of 121.16, 3.51, and 183.81 Å2; average H-bond values of 466.05, 470.84, and 456.84 Å, respectively. The results revealed that CID_8299 had the highest stability and consistent interaction with the target protein throughout the simulation period. According to the toxicity analysis, CID_8299, which is found in the Glycyrrhiza glabra plant, can also cross the BBB, which makes it unbeatable in treating neuro-disease caused by C. neoformans and may potentially block FTase protein's activity inhibiting post-translational lipidation of essential signal transduction protein.

The PH domain in the ArfGAP ASAP1 drives catalytic activation through an unprecedented allosteric mechanism

ASAP1 is a multidomain Arf GTPase-activating protein (ArfGAP) that catalyzes GTP hydrolysis on the small GTPase Arf1 and is implicated in cancer progression. The PH domain of ASAP1 enhances its activity greater than 7 orders of magnitude but the underlying mechanisms remain poorly understood. Here, we combined Nuclear Magnetic Resonance (NMR), Molecular Dynamic (MD) simulations and mathematical modeling of functional data to build a comprehensive structural-mechanistic model of the complex of Arf1 and the ASAP1 PH domain on a membrane surface. Our results support a new conceptual model in which the PH domain contributes to efficient catalysis not only by membrane recruitment but by acting as a critical component of the catalytic interface, binding Arf·GTP and allosterically driving it towards the catalytic transition state. We discuss the biological implications of these results and how they may apply more broadly to poorly understood membrane-dependent regulatory mechanisms controlling catalysis of the ArfGAP superfamily as well as other peripheral membrane enzymes.

A conserved H-bond network in human aquaporin-1 is necessary for native folding and oligomerization

Aquaporins (AQPs) are α-helical transmembrane proteins that conduct water through membranes with high selectivity and permeability. For human AQP1, in addition to the functional Asn-Pro-Ala motifs and the aromatic/Arg selectivity filter within the pore, there are several highly conserved residues that form an expansive hydrogen-bonding network. Previous solid-state nuclear magnetic resonance studies and structural conservation analysis have detailed which residues may be involved in this network. We explored this network by mutating the side chains or backbones involved in hydrogen-bonding, generating the following mutants: N127A, V133P, E142A, T187A, R195A, and S196A. The fold and stability of these mutants were assessed with attenuated total reflection Fourier transform infrared spectroscopy coupled with hydrogen/deuterium exchange upon increasing temperature. We found that replacement of any of the chosen residues to alanine leads to either partial instability or outright misfolding at room temperature, with the latter being most pronounced for the N127A, V133P, T187A, and R195A mutants. Deconvolution analysis of the amide I band revealed considerable secondary structure deviations, with some mutants exhibiting new random coil and β sheet structures. We also found that some of these mutations potentially disrupt the oligomerization of human AQP1. BN-PAGE and DLS data provide evidence toward the loss of tetramers within most of the mutants, meanwhile only the S196A mutant retains tetrameric organization. The molecular dynamics simulation of the wild-type, and the N127A, E142A, and T187A mutants show that these mutations result in major rearrangements of intra- and intermonomer hydrogen-bond networks. Overall, we show that specific point mutations that perturb hydrogen-bonding clusters result in severe misfolding in hAQP1 and disruption of its oligomerization. These data provide valuable insight into the structural stability of human aquaporin-1 and have implications toward other members of the AQP family, as these networks are largely conserved among a variety of human and nonmammalian AQP homologs.

Exploring TAS2R46 biomechanics through molecular dynamics and network analysis

Understanding the intricate interplay between structural features and signal-processing events is crucial for unravelling the mechanisms of biomolecular systems. G protein-coupled receptors (GPCRs), a pervasive protein family in humans, serve a wide spectrum of vital functions. TAS2Rs, a subfamily of GPCRs, play a primary role in recognizing bitter molecules and triggering events leading to the perception of bitterness, a crucial defence mechanism against spoiled or poisonous food. Beyond taste, TAS2Rs function is associated with many diseases as they are expressed in several extra-oral tissues. Given that the precise functioning mechanisms of TAS2R remain poorly understood, this study employed molecular dynamics simulations combined with network-based analysis to investigate local conformational changes and global structural correlations in different states of the receptor. The focus was on the human TAS2R46 bitter taste receptor, recently resolved experimentally, both in the presence and absence of strychnine, a known bitter agonist. The results showed that the ligand-bound state of the receptor exhibited more correlated dynamics compared to the apo state, and the presence of the agonist mediated the allosteric network between two helices (TM3 and TM6) which mainly convey the signal transferring from the extracellular to the intracellular region. By elucidating the hallmarks of the conformational changes and allosteric network of TAS2R46 under varying conditions, this study has enabled the identification of the unique structural and dynamics features of this receptor, thereby establishing a foundation for a more profound characterisation of this intriguing class of receptors.

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.

Molecular Mechanisms Underlying the Spectral Shift in Zebrafish Cone Opsins

Visual pigments are essential for converting light into electrical signals during vision. Composed of an opsin protein and a retinal-based chromophore, pigments in vertebrate rods (Rh1) and cones (Rh2) have different spectral sensitivities, with distinct peak absorption wavelengths determined by the shape and composition of the chromophore binding pocket. Despite advances in understanding Rh1 pigments such as bovine rhodopsin, the molecular basis of spectral shifts in Rh2 cone opsins has been less studied, particularly the E122Q mutation, which accounts for about half of the observed spectral shift in these pigments. In this study, we employed molecular modeling and quantum mechanical techniques to investigate the molecular mechanisms behind the spectral difference in blue-shifted Rh2-1 (absorption peak = 467 nm, 122Q) and green-shifted Rh2-4 (absorption peak = 505 nm, 122E) zebrafish cone opsins. We modeled the pigments 3D structures based on their sequences and conducted all-atom molecular dynamics simulations totaling 2 microseconds. Distance analysis of the trajectories identified three key sites: E113, E181, and E122. The E122Q mutation, previously known, validates our findings, while E181 and E113 are newly identified contributors. Structural analysis revealed key features with differing values that explain the divergent spectral sensitivities of Rh2-1 and Rh2-4: 1) chromophore atom fluctuations and C5-C6 torsion angle, 2) binding pocket volume, 3) hydration patterns, and 4) E113-chromophore interaction stability. Quantum mechanics further confirms the critical role of residue E181 in Rh2-1 and E122 in Rh2-4 for their spectral behavior. Our study provides new insights into the molecular determinants of spectral shifts in cone opsins, and we anticipate that it will serve as a starting point for a broader understanding of the functional diversity of visual pigments.

The retinal chromophore environment in an inward light-driven proton pump studied by solid-state NMR and hydrogen-bond network analysis

Inward proton pumping is a relatively new function for microbial rhodopsins, retinal-binding light-driven membrane proteins. So far, it has been demonstrated for two unrelated subgroups of microbial rhodopsins, xenorhodopsins and schizorhodopsins. A number of recent studies suggest unique retinal-protein interactions as being responsible for the reversed direction of proton transport in the latter group. Here, we use solid-state NMR to analyze the retinal chromophore environment and configuration in an inward proton-pumping Antarctic schizorhodopsin. Using fully 13C-labeled retinal, we have assigned chemical shifts for every carbon atom and, assisted by structure modelling and molecular dynamics simulations, made a comparison with well-studied outward proton pumps, identifying locations of the unique protein-chromophore interactions for this functional subclass of microbial rhodopsins. Both the NMR results and molecular dynamics simulations point to the distinctive polar environment in the proximal part of the retinal, which may result in a hydration pattern dramatically different from that of the outward proton pumps, causing the reversed proton transport.

In silico prediction of the interaction of legacy and novel per- and poly-fluoroalkyl substances (PFAS) with selected human transporters and of their possible accumulation in the human body

Per- and poly-fluorinated compounds constitute a wide group of fluorocarbon chemicals with widespread industrial applications, ranging from non-stick coating in cookware to water surfactants, from fire-fighting foams to water-repellent coatings on textiles. Presently, over 12,000 PFAS are known worldwide. In recent years, extensive research has focused on investigating the biological effects of these molecules on various organisms, including humans. Here, we conducted in silico simulations to examine the potential binding of a representative selection of PFAS to various human proteins known to be involved in chemical transportation and accumulation processes. Specifically, we targeted human serum albumin (HSA), transthyretin (TTR), thyroxine binding protein (TBG), fatty acid binding proteins (FABPs), organic anion transporters (OATs), aiming to assess the potential for bioaccumulation. Molecular docking simulations were employed for this purpose, supplemented by molecular dynamics (MD) simulations to account for protein flexibility, when necessary. Our findings indicate that so-called "legacy PFAS" such as PFOA or PFOS exhibit a higher propensity for interaction with the analysed human protein targets compared to newly formulated PFAS, characterised by higher branching and hydrophilicity, and possibly a higher accumulation in the human body.

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.

Structural Determinants of Peptide Nanopore Formation

We have evolved the nanopore-forming macrolittin peptides from the bee venom peptide melittin using successive generations of synthetic molecular evolution. Despite their sequence similarity to the broadly membrane permeabilizing cytolytic melittin, the macrolittins have potent membrane selectivity. They form nanopores in synthetic bilayers made from 1-palmitoyl, 2-oleoyl-phosphatidylcholine (POPC) at extremely low peptide concentrations and yet have essentially no cytolytic activity against any cell membrane, even at high concentration. Here, we explore the structural determinants of macrolittin nanopore stability in POPC bilayers using atomistic molecular dynamics simulations and experiments on macrolittins and single-site variants. Simulations of macrolittin nanopores in POPC bilayers show that they are stabilized by an extensive, cooperative hydrogen bond network comprised of the many charged and polar side chains interacting with each other via bridges of water molecules and lipid headgroups. Lipid molecules with unusual conformations participate in the H-bond network and are an integral part of the nanopore structure. To explore the role of this H-bond network on membrane selectivity, we swapped three critical polar residues with the nonpolar residues found in melittin. All variants have potency, membrane selectivity, and cytotoxicity that were intermediate between a cytotoxic melittin variant called MelP5 and the macrolittins. Simulations showed that the variants had less organized H-bond networks of waters and lipids with unusual structures. The membrane-spanning, cooperative H-bond network is a critical determinant of macrolittin nanopore stability and membrane selectivity. The results described here will help guide the future design and optimization of peptide nanopore-based applications.

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.

A Proteorhodopsin-Related Photosensor Expands the Repertoire of Structural Motifs Employed by Sensory Rhodopsins

Microbial rhodopsins are light-activated retinal-binding membrane proteins that perform a variety of ion transport and photosensory functions. They display several cases of convergent evolution where the same function is present in unrelated or very distant protein groups. Here we report another possible case of such convergent evolution, describing the biophysical properties of a new group of sensory rhodopsins. The first representative of this group was identified in 2004 but none of the members had been expressed and characterized. The well-studied haloarchaeal sensory rhodopsins interacting with methyl-accepting Htr transducers are close relatives of the halobacterial proton pump bacteriorhodopsin. In contrast, the sensory rhodopsins we describe here are relatives of proteobacterial proton pumps, proteorhodopsins, but appear to interact with Htr-like transducers likewise, even though they do not conserve the residues important for the interaction of haloarchaeal sensory rhodopsins with their transducers. The new sensory rhodopsins display many unusual amino acid residues, including those around the retinal chromophore; most strikingly, a tyrosine in place of a carboxyl counterion of the retinal Schiff base on helix C. To characterize their unique sequence motifs, we augment the spectroscopy and biochemistry data by structural modeling of the wild-type and three mutants. Taken together, the experimental data, bioinformatics sequence analyses, and structural modeling suggest that the tyrosine/aspartate complex counterion contributes to a complex water-mediated hydrogen-bonding network that couples the protonated retinal Schiff base to an extracellular carboxylic dyad.

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.

Graph-Based Analyses of Dynamic Water-Mediated Hydrogen-Bond Networks in Phosphatidylserine: Cholesterol Membranes

Phosphatidylserine lipids are anionic molecules present in eukaryotic plasma membranes, where they have essential physiological roles. The altered distribution of phosphatidylserine in cells such as apoptotic cancer cells, which, unlike healthy cells, expose phosphatidylserine, is of direct interest for the development of biomarkers. We present here applications of a recently implemented Depth-First-Search graph algorithm to dissect the dynamics of transient water-mediated lipid clusters at the interface of a model bilayer composed of 1-palmytoyl-2-oleoyl-sn-glycero-2-phosphatidylserine (POPS) and cholesterol. Relative to a reference POPS bilayer without cholesterol, in the POPS:cholesterol bilayer there is a somewhat less frequent sampling of relatively complex and extended water-mediated hydrogen-bond networks of POPS headgroups. The analysis protocol used here is more generally applicable to other lipid:cholesterol bilayers.

Structures of channelrhodopsin paralogs in peptidiscs explain their contrasting K+ and Na+ selectivities

Kalium channelrhodopsin 1 from Hyphochytrium catenoides (HcKCR1) is a light-gated channel used for optogenetic silencing of mammalian neurons. It selects K+ over Na+ in the absence of the canonical tetrameric K+ selectivity filter found universally in voltage- and ligand-gated channels. The genome of H. catenoides also encodes a highly homologous cation channelrhodopsin (HcCCR), a Na+ channel with >100-fold larger Na+ to K+ permeability ratio. Here, we use cryo-electron microscopy to determine atomic structures of these two channels embedded in peptidiscs to elucidate structural foundations of their dramatically different cation selectivity. Together with structure-guided mutagenesis, we show that K+ versus Na+ selectivity is determined at two distinct sites on the putative ion conduction pathway: in a patch of critical residues in the intracellular segment (Leu69/Phe69, Ile73/Ser73 and Asp116) and within a cluster of aromatic residues in the extracellular segment (primarily, Trp102 and Tyr222). The two filters are on the opposite sides of the photoactive site involved in channel gating.

Network analysis of chromophore binding site in LOV domain

Photoreceptor proteins are versatile toolbox for developing biosensors for optogenetic applications. These molecular tools get activated upon illumination of blue light, which in turn offers a non-invasive method for gaining high spatiotemporal resolution and precise control of cellular signal transduction. The Light-Oxygen-Voltage (LOV) domain family of proteins is a well-recognized system for constructing optogenetic devices. Translation of these proteins into efficient cellular sensors is possible by tuning their photochemistry lifetime. However, the bottleneck is the need for more understanding of the relationship between the protein environment and photocycle kinetics. Significantly, the effect of the local environment also modulates the electronic structure of chromophore, which perturbs the electrostatic and hydrophobic interaction within the binding site. This work highlights the critical factors hidden in the protein networks, linking with their experimental photocycle kinetics. It presents an opportunity to quantitatively examine the alternation in chromophore's equilibrium geometry and identify details which have substantial implications in designing synthetic LOV constructs with desirable photocycle efficiency.

Starch-Derived Superabsorbent Polymer in Remediation of Solid Waste Sludge Based on Water–Polymer Interaction

The purpose of this study is to assess water–polymer interaction in synthesized starch-derived superabsorbent polymer (S-SAP) for the treatment of solid waste sludge. While S-SAP for solid waste sludge treatment is still rare, it offers a lower cost for the safe disposal of sludge into the environment and recycling of treated solid as crop fertilizer. For that to be possible, the water–polymer interaction on S-SAP must first be fully comprehended. In this study, the S-SAP was prepared through graft polymerization of poly (methacrylic acid-co-sodium methacrylate) on the starch backbone. By analyzing the amylose unit, it was possible to avoid the complexity of polymer networks when considering S-SAP using molecular dynamics (MD) simulations and density functional theory (DFT). Through the simulations, formation of hydrogen bonding between starch and water on the H06 of amylose was assessed for its flexibility and less steric hindrance. Meanwhile, water penetration into S-SAP was recorded by the specific radial distribution function (RDF) of atom–molecule interaction in the amylose. The experimental evaluation of S-SAP correlated with high water capacity by measuring up to 500% of distilled water within 80 min and more than 195% of the water from solid waste sludge for 7 days. In addition, the S-SAP swelling showed a notable performance of a 77 g/g swelling ratio within 160 min, while a water retention test showed that S-SAP was capable of retaining more than 50% of the absorbed water within 5 h of heating at 60 °C. The water retention of S-SAP adheres to pseudo-second-order kinetics for chemisorption reactions. Therefore, the prepared S-SAP might have potential applications as a natural superabsorbent, especially for the development of sludge water removal technology.

Mechanistic Insights of Polyphenolic Compounds from Rosemary Bound to Their Protein Targets Obtained by Molecular Dynamics Simulations and Free-Energy Calculations

Rosemary represents an important medicinal plant that has been attributed with various health-promoting properties, especially antioxidative, anti-inflammatory, and anticarcinogenic activities. Carnosic acid, carnosol, and rosmanol, as well as the phenolic acid ester rosmarinic acid, are the main compounds responsible for these actions. In our earlier research, we carried out an inverse molecular docking at the proteome scale to determine possible protein targets of the mentioned compounds. Here, we subjected the previously identified ligand-protein complexes with HIV-1 protease, K-RAS, and factor X to molecular dynamics simulations coupled with free-energy calculations. We observed that carnosic acid and rosmanol act as viable binders of the HIV-1 protease. In addition, carnosol represents a potential binder of the oncogene protein K-RAS. On the other hand, rosmarinic acid was characterized as a weak binder of factor X. We also emphasized the importance of water-mediated hydrogen-bond networks in stabilizing the binding conformation of the studied polyphenols, as well as in mechanistically explaining their promiscuous nature.

Interplay between local protein interactions and water bridging of a proton antenna carboxylate cluster

Proteins that bind protons at cell membrane interfaces often expose to the bulk clusters of carboxylate and histidine sidechains that capture protons transiently and, in proton transporters, deliver protons to an internal site. The protonation-coupled dynamics of bulk-exposed carboxylate clusters, also known as proton antennas, is poorly described. An essential open question is how water-mediated bridges between sidechains of the cluster respond to protonation change and facilitate transient proton storage. To address this question, here I studied the protonation-coupled dynamics at the proton-binding antenna of PsbO, a small extrinsinc subunit of the photosystem II complex, with atomistic molecular dynamics simulations and systematic graph-based analyses of dynamic protein and protein-water hydrogen-bond networks. The protonation of specific carboxylate groups is found to impact the dynamics of their local protein-water hydrogen-bond clusters. Regardless of the protonation state considered for PsbO, carboxylate pairs that can sample direct hydrogen bonding, or bridge via short hydrogen-bonded water chains, anchor to nearby basic or polar protein sidechains. As a result, carboxylic sidechains of the hypothesized antenna cluster are part of dynamic hydrogen bond networks that may rearrange rapidly when the protonation changes.

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.

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.

Algorithm to catalogue topologies of dynamic lipid hydrogen-bond networks

Lipid membrane interfaces host reactions essential for the functioning of cells. The hydrogen-bonding environment at the membrane interface is particularly important for binding of proteins, drug molecules, and ions. We present here the implementation and applications of a depth-first search algorithm that analyzes dynamic lipid interaction networks. Lipid hydrogen-bond networks sampled transiently during simulations of lipid bilayers are clustered according to main types of topologies that characterize three-dimensional arrangements of lipids connected to each other via short water bridges. We characterize the dynamics of hydrogen-bonded lipid clusters in simulations of model POPE and POPE:POPG membranes that are often used for bacterial membrane proteins, in a model of the Escherichia coli membrane with six different lipid types, and in POPS membranes. We find that all lipids sample dynamic hydrogen-bonded networks with linear, star, or circular arrangements of the lipid headgroups, and larger networks with combinations of these three types of topologies. Overall, linear lipid-water bridges tend to be short. Water-mediated lipid clusters in all membranes with PE lipids tend to be somewhat small, with about four lipids in all membranes studied here. POPS membranes allow circular arrangements of three POPS lipids to be sampled frequently, and complex arrangements of linear, star, and circular paths may also be sampled. These findings suggest a molecular picture of the membrane interface whereby lipid molecules transiently connect in clusters with somewhat small spatial extension.

Characterizing Protein Protonation Microstates Using Monte Carlo Sampling

Proteins are polyelectrolytes with acidic and basic amino acids Asp, Glu, Arg, Lys, and His, making up ≈25% of the residues. The protonation state of residues, cofactors, and ligands defines a "protonation microstate". In an ensemble of proteins some residues will be ionized and others neutral, leading to a mixture of protonation microstates rather than in a single one as is often assumed. The microstate distribution changes with pH. The protein environment also modifies residue proton affinity so microstate distributions change in different reaction intermediates or as ligands are bound. Particular protonation microstates may be required for function, while others exist simply because there are many states with similar energy. Here, the protonation microstates generated in Monte Carlo sampling in MCCE are characterized in HEW lysozyme as a function of pH and bacterial photosynthetic reaction centers (RCs) in different reaction intermediates. The lowest energy and highest probability microstates are compared. The ΔG, ΔH, and ΔS between the four protonation states of Glu35 and Asp52 in lysozyme are shown to be calculated with reasonable precision. At pH 7 the lysozyme charge ranges from 6 to 10, with 24 accepted protonation microstates, while RCs have ≈50,000. A weighted Pearson correlation analysis shows coupling between residue protonation states in RCs and how they change when the quinone in the Q_B site is reduced. Protonation microstates can be used to define input MD parameters and provide insight into the motion of protons coupled to reactions.

Mechanisms of long-distance allosteric couplings in proton-binding membrane transporters

Membrane transporters that use proton binding and proton transfer for function couple local protonation change with changes in protein conformation and water dynamics. Changes of protein conformation might be required to allow transient formation of hydrogen-bond networks that bridge proton donor and acceptor pairs separated by long distances. Inter-helical hydrogen-bond networks adjust rapidly to protonation change, and ensure rapid response of the protein structure and dynamics. Membrane transporters with known three-dimensional structures and proton-binding groups inform on general principles of protonation-coupled protein conformational dynamics. Inter-helical hydrogen bond motifs between proton-binding carboxylate groups and a polar sidechain are observed in unrelated membrane transporters, suggesting common principles of coupling protonation change with protein conformational dynamics.

Oxa-376 and Oxa-530 variants of β-lactamase: computational study uncovers potential therapeutic targets of Acinetobacter baumannii

Antimicrobial resistance is a major global health crisis, resulting in thousands of deaths each year. Antibiotics' effectiveness against microorganisms deteriorates over time as multidrug resistance (MDR) develops, which is exacerbated by irregular antibiotic use, poor disease management, and the evasive nature of bacteria. The World Health Organization has recognized multidrug resistance as a critical public health concern, and Acinetobacter baumannii has been at the center of attention due to its ability to develop multidrug resistance (MDR). It generally produces carbapenem-hydrolyzing oxacillinase, which has been identified as the primary source of beta-lactam resistance in MDR bacteria. Recently, point mutations in A. baumannii have been identified as a key factor of multidrug resistance, making them a prime concern for researchers. The goal of the current work was to establish a unique way of finding multidrug-resistant variants and identify the most damaging mutations in the existing databases. We characterized the deleterious variants of oxacillinases using several computational tools. Following a thorough analysis, Oxa-376 and Oxa-530 were found to be more damaging when compared with the wild-type Oxa-51. The mutants' 3D structures were then prepared and refined with RaptorX, GalaxyRefine, and SAVES servers. Our research incorporates seven antimicrobial agents to illustrate the resistance capability of the variants of oxacillinase by evaluating binding affinity in Autodock-vina and Schrodinger software. RMSD, RMSF, Radius of gyration analysis, the solvent-accessible surface area (SASA), hydrogen bonding analysis and MM-GBSA from Molecular Dynamics Simulation revealed the dynamic nature and stability of wild-type and Oxa-376 and Oxa-530 variants. Our findings will benefit researchers looking for the deleterious mutations of Acinetobacter baumannii and new therapeutics to combat those variants. However, further studies are necessary to evaluate the mechanism of hydrolyzing activity and antibiotic resistance of these variants.

Determining novel therapeutic targets of Acinetobacter baumannii. Deleterious variants, causing antibiotic resistance, were identified by molecular docking and molecular dynamics simulation suggesting new therapeutic targets Oxa-376 and Oxa-530.

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.

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.

Frontiers in Multiscale Modeling of Photoreceptor Proteins

This perspective article highlights the challenges in the theoretical description of photoreceptor proteins using multiscale modeling, as discussed at the CECAM workshop in Tel Aviv, Israel. The participants have identified grand challenges and discussed the development of new tools to address them. Recent progress in understanding representative proteins such as green fluorescent protein, photoactive yellow protein, phytochrome, and rhodopsin is presented, along with methodological developments.

Phosphatidylglyerol Lipid Binding at the Active Site of an Intramembrane Protease

Transmembrane substrate cleavage by the small Escherichia coli rhomboid protease GlpG informs on mechanisms by which lipid interactions shape reaction coordinates of membrane-embedded enzymes. Here, I review and discuss new work on the molecular picture of protein-lipid interactions that might govern the formation of the substrate-enzyme complex in fluid lipid membranes. Negatively charged PG-type lipids are of particular interest, because they are a major component of bacterial membranes. Atomistic computer simulations indicate POPG and DOPG lipids bridge remote parts of GlpG and might pre-occupy the substrate-docking site. Inhibition of catalytic activity by PG lipids could arise from ligand-like lipid binding at the active site, which could delay or prevent substrate docking. Dynamic protein-lipid H-bond networks, water access to the active site, and fluctuations in the orientation of GlpG suggest that GlpG has lipid-coupled dynamics that could shape the energy landscape of transmembrane substrate docking.

A graph-based approach identifies dynamic H-bond communication networks in spike protein S of SARS-CoV-2

Corona virus spike protein S is a large homo-trimeric protein anchored in the membrane of the virion particle. Protein S binds to angiotensin-converting-enzyme 2, ACE2, of the host cell, followed by proteolysis of the spike protein, drastic protein conformational change with exposure of the fusion peptide of the virus, and entry of the virion into the host cell. The structural elements that govern conformational plasticity of the spike protein are largely unknown. Here, we present a methodology that relies upon graph and centrality analyses, augmented by bioinformatics, to identify and characterize large H-bond clusters in protein structures. We apply this methodology to protein S ectodomain and find that, in the closed conformation, the three protomers of protein S bring the same contribution to an extensive central network of H-bonds, and contribute symmetrically to a relatively large H-bond cluster at the receptor binding domain, and to a cluster near a protease cleavage site. Markedly different H-bonding at these three clusters in open and pre-fusion conformations suggest dynamic H-bond clusters could facilitate structural plasticity and selection of a protein S protomer for binding to the host receptor, and proteolytic cleavage. From analyses of spike protein sequences we identify patches of histidine and carboxylate groups that could be involved in transient proton binding.

Molecular movie of nucleotide binding to a motor protein

BACKGROUND

The SecA DEAD (Asp-Glu-Ala-Asp) motor protein uses binding and hydrolysis of adenosine triphosphate (ATP) to push secretory proteins across the plasma membrane of bacteria. The reaction coordinate of nucleotide exchange is unclear at the atomic level of detail.

METHODS

We performed multiple atomistic computations of the DEAD motor domain of SecA with different occupancies of the nucleotide and magnesium ion sites, for a total of ~1.7 μs simulation time. To characterize dynamics at the active site we analyzed hydrogen-bond networks.

RESULTS

ATP and ADP can bind spontaneously at the interface between the nucleotide binding domains, albeit at an intermediate binding site distinct from the native site. Binding of the nucleotide is facilitated by the presence of a magnesium ion close to the glutamic group of the conserved DEAD motif. In the absence of the magnesium ion, protein interactions of the ADP molecule are perturbed.

CONCLUSIONS

A protein hydrogen-bond network whose dynamics couples to the occupancy of the magnesium ion site helps guide the nucleotide along the nucleotide exchange path. In SecA, release of magnesium might be required to destabilize the ADP binding site prior to release of the nucleotide.

GENERAL SIGNIFICANCE

We identified dynamic hydrogen-bond networks that help control nucleotide exchange in SecA, and stabilize ADP at an intermediate site that could explain slow release. The reaction coordinate of the protein motor involves complex rearrangements of a hydrogen-bond network at the active site, with perturbation of the magnesium ion site likely occurring prior to the release of ADP.