Graphs of dynamic H-bond networks: from model proteins to protein complexes in cell signaling

PyInteraph2 and PyInKnife2 to Analyze Networks in Protein Structural Ensembles

Due to the complex nature of noncovalent interactions and their long-range effects, analyzing protein conformations using network theory can be enlightening. Protein Structure Networks (PSNs) provide a convenient formalism to study protein structures in relation to essential properties such as key residues for structural stability, allosteric communication, and the effects of modifications of the protein. PSNs can be defined according to very different principles, and the available tools have limitations in input formats, supported models, and version control. Other outstanding problems are related to the definition of network cutoffs and the assessment of the stability of the network properties. The protein science community could benefit from a common framework to carry out these analyses and make them easier to reproduce, reuse, and evaluate. We here provide two open-source software packages, PyInteraph2 and PyInKnife2, to implement and analyze PSNs in a reproducible and documented manner. PyInteraph2 interfaces with multiple formats for protein ensembles and incorporates different network models with the possibility of integrating them into a macronetwork and performing various downstream analyses, including hubs, connected components, and several other centrality measures, and visualizes the networks or further analyzes them thanks to compatibility with Cytoscape.PyInKnife2 that supports the network models implemented in PyInteraph2. It employs a jackknife resampling approach to estimate the convergence of network properties and streamline the selection of distance cutoffs. We foresee that the modular structure of the code and the supported version control system will promote the transition to a community-driven effort, boost reproducibility, and establish common protocols in the PSN field. As developers, we will guarantee the introduction of new functionalities and maintenance, assistance, and training of new contributors.

Enhanced Interfacial H-Bond Networks Promote Glycan-Glycan Recognition and Interaction

H-bond networks at heterogeneous interfaces play crucial roles in bioseparation, biocatalysis, biochip array profiling, and functional nanosystem self-assembly, but their precise modulation and enhancement remain challenging. In this study, we have discovered that interfacial hydrophobic hydration significantly enhances H-bond networks at the interface between a glycan-modified adsorbent and a methanol-water-acetonitrile ternary solution. The enhanced H-bond networks greatly promote the adsorbent-solution heterogeneous glycan-glycan recognition and interaction. This novel hydrophobic hydration-enhanced hydrophilic interaction (HEHI) strategy improves the affinity and efficiency of intact glycopeptide enrichment. Compared with the commonly used hydrophilic-interaction enrichment strategy, 23.5 and 48.5% more intact N- and O-glycopeptides are identified, and the enrichment recoveries of half of the glycopeptides are increased >100%. Further, in-depth profiling of both N- and O-glycosylation occurring on SARS-CoV-2 S1 and hACE2 proteins has been achieved with more glycan types and novel O-glycosylation information involved. Interfacial hydrophobic hydration provides a powerful tool for the modulation of hydrophilic interactions in biological systems.

Bioinspired and biomolecular catalysts for energy conversion and storage

Metalloenzymes are remarkable for facilitating challenging redox transformations with high efficiency and selectivity. In the area of alternative energy, scientists aim to capture these properties in bioinspired and engineered biomolecular catalysts for the efficient and fast production of fuels from low-energy feedstocks such as water and carbon dioxide. In this short review, efforts to mimic biological catalysts for proton reduction and carbon dioxide reduction are highlighted. Two important recurring themes are the importance of the microenvironment of the catalyst active site and the key role of proton delivery to the active site in achieving desired reactivity. Perspectives on ongoing and future challenges are also provided.

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.

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.

Deprotonation from an OH on myo-Inositol Promoted by μ2-Bridges with Possible Regioselectivity/Chiral Selectivity

Single-crystal structures of myo-inositol complexes with erbium ([Er₂(C₆H₁₁O₆)₂(H₂O)₅Cl₂]Cl₂(H₂O)₄, denoted ErI hereafter) and strontium (Sr(C₆H₁₂O₆)₂(H₂O)₂Cl₂, denoted SrI hereafter) are described. In ErI, deprotonation occurs on an OH of myo-inositol, although the complex is synthesized in an acidic solution, and the pK_a values of all of the OHs in myo-inositol are larger than 12. The deprotonated OH is involved in a μ₂-bridge. The polarization from two Er³⁺ ions activates the chemically relatively inert OH and promotes deprotonation. In the stable conformation of myo-inositol, there are five equatorial OHs and one axial OH. The deprotonation occurs on the only axial OH, suggesting that the deprotonation possesses characteristics of regioselectivity/chiral selectivity. Two Er³⁺ ions in the μ₂-bridge are stabilized by five-membered rings formed by chelating Er³⁺ with an O-C-C-O moiety. As revealed by the X-ray crystallography study, the absolute values of the O-C-C-O torsion angles decrease from ∼60 to ∼45° upon chelating. Since the O-C-C-O moiety is within a six-membered ring, the variation of the torsion angle may exert distortion of the chair conformation. Quantum chemistry calculation results indicate that an axial OH flanked by two equatorial OHs (double ax-eq motif) is favorable for the formation of a μ₂-bridge, accounting for the selectivity. The double ax-eq motif may be used in a rational design of high-performance catalysts where deprotonation with high regioselectivity/chiral selectivity is carried out.

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.

Activation-Induced Reorganization of Energy Transport Networks in the β2 Adrenergic Receptor

We compute energy exchange networks (EENs) through the β₂ adrenergic receptor (β₂AR), a G-protein coupled receptor (GPCR), in inactive and active states, based on the results of molecular dynamics simulations of this membrane bound protein. We introduce a new definition for the reorganization of EENs upon activation that depends on the relative change in rates of energy transfer across noncovalent contacts throughout the protein. On the basis of the reorganized network that we obtain for β₂AR upon activation, we identify a branched pathway between the agonist binding site and the cytoplasmic region, where a G-protein binds to the receptor when activated. The pathway includes all of the motifs containing molecular switches previously identified as contributing to the allosteric transition of β₂AR upon agonist binding. EENs and their reorganization upon activation are compared with structure-based contact networks computed for the inactive and active states of β₂AR.

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.

A nexus of intrinsic dynamics underlies translocase priming

The cytoplasmic ATPase SecA and the membrane-embedded SecYEG channel assemble to form the Sec translocase. How this interaction primes and catalytically activates the translocase remains unclear. We show that priming exploits a nexus of intrinsic dynamics in SecA. Using atomistic simulations, smFRET, and HDX-MS, we reveal multiple dynamic islands that cross-talk with domain and quaternary motions. These dynamic elements are functionally important and conserved. Central to the nexus is a slender stem through which rotation of the preprotein clamp of SecA is biased by ATPase domain motions between open and closed clamping states. An H-bonded framework covering most of SecA enables multi-tier dynamics and conformational alterations with minimal energy input. As a result, cognate ligands select preexisting conformations and alter local dynamics to regulate catalytic activity and clamp motions. These events prime the translocase for high-affinity reception of non-folded preprotein clients. Dynamics nexuses are likely universal and essential in multi-liganded proteins.

Inhibitor binding influences the protonation states of histidines in SARS-CoV-2 main protease

The main protease (M^pro) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an attractive target for antiviral therapeutics. Recently, many high-resolution apo and inhibitor-bound structures of M^pro, a cysteine protease, have been determined, facilitating structure-based drug design. M^pro plays a central role in the viral life cycle by catalyzing the cleavage of SARS-CoV-2 polyproteins. In addition to the catalytic dyad His41–Cys145, M^pro contains multiple histidines including His163, His164, and His172. The protonation states of these histidines and the catalytic nucleophile Cys145 have been debated in previous studies of SARS-CoV M^pro, but have yet to be investigated for SARS-CoV-2. In this work we have used molecular dynamics simulations to determine the structural stability of SARS-CoV-2 M^pro as a function of the protonation assignments for these residues. We simulated both the apo and inhibitor-bound enzyme and found that the conformational stability of the binding site, bound inhibitors, and the hydrogen bond networks of M^pro are highly sensitive to these assignments. Additionally, the two inhibitors studied, the peptidomimetic N3 and an α-ketoamide, display distinct His41/His164 protonation-state-dependent stabilities. While the apo and the N3-bound systems favored N_δ (HD) and N_ϵ (HE) protonation of His41 and His164, respectively, the α-ketoamide was not stably bound in this state. Our results illustrate the importance of using appropriate histidine protonation states to accurately model the structure and dynamics of SARS-CoV-2 M^pro in both the apo and inhibitor-bound states, a necessary prerequisite for drug-design efforts.

The main protease (M^pro) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an attractive target for antiviral therapeutics.

Tuning Proton Transfer Thermodynamics in SARS-Cov-2 Main Protease: Implications for Catalysis and Inhibitor Design

In this comutational work a hybrid quantum mechanics/molecular mechanics approach, the MD-PMM approach, is used to investigate the proton transfer reaction the activates the catalytic activity of SARS-CoV-2 main protease. The proton transfer thermodynamics is investigated for the apo ensyme (i.e., without any bound substrate or inhibitor) and in the presence of a inhibitor, N3, which was previously shown to covalently bind SARS-CoV-2 main protease.

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.

Inhibitor binding influences the protonation states of histidines in SARS-CoV-2 main protease

The main protease (M ^pro ) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an attractive target for antiviral therapeutics. Recently, many high-resolution apo and inhibitor-bound structures of M ^pro , a cysteine protease, have been determined, facilitating structure-based drug design. M ^pro plays a central role in the viral life cycle by catalyzing the cleavage of SARS-CoV-2 polyproteins. In addition to the catalytic dyad His41-Cys145, M ^pro contains multiple histidines including His163, His164, and His172. The protonation states of these histidines and the catalytic nu-cleophile Cys145 have been debated in previous studies of SARS-CoV M ^pro , but have yet to be investigated for SARS-CoV-2. In this work we have used molecular dynamics simulations to determine the structural stability of SARS-CoV-2 M ^pro as a function of the protonation assignments for these residues. We simulated both the apo and inhibitor-bound enzyme and found that the conformational stability of the binding site, bound inhibitors, and the hydrogen bond networks of M ^pro are highly sensitive to these assignments. Additionally, the two inhibitors studied, the peptidomimetic N3 and an α -ketoamide, display distinct His41/His164 protonation-state-dependent stabilities. While the apo and the N3-bound systems favored N _δ (HD) and N _ϵ (HE) protonation of His41 and His164, respectively, the α -ketoamide was not stably bound in this state. Our results illustrate the importance of using appropriate histidine protonation states to accurately model the structure and dynamics of SARS-CoV-2 M ^pro in both the apo and inhibitor-bound states, a necessary prerequisite for drug-design efforts.