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

We apply graph-based approaches to identify H-bond clusters in protein complexes. Three conformations of spike protein S have distinct H-bond clusters at key sites. Hydrogen-bond clusters could govern structural plasticity of spike protein S. Protein S binds to ACE2 receptor via H-bond clusters extending deep across interface.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.PSI COVID19 Emergency Science FundSpanish Ministry of Science, Innovation and Universities RTI2018-098983-B-I00Excellence Initiative of the German Federal and State Governments via the Freie Universitat BerlinGerman Research Foundation (DFG) SFB 107

Biophys. J.

The dynamics of proton transfer between the surface of purple membrane and the aqueous bulk have recently been investigated by the Laser Induced Proton Pulse Method. Following a Delta- function release of protons to the bulk, the system was seen to regain its state of equilibrium within a few hundreds of microseconds. These measurements set the time frame for the relaxation of any state of acid-base disequilibrium between the bacteriorhodopsin's surface and the bulk. It was also deduced that the released protons react with the various proton binding within less than 10 mus. In the present study, we monitored the photocycle and the proton-cycle of photo-excited bacteriorhodopsin, in the absence of added buffer, and calculated the proton balance between the Schiff base and the bulk phase in a time-resolved mode. It was noticed that the late phase of the M decay (beyond 1 ms) is characterized by a slow (subsecond) relaxation of disequilibrium, where the Schiff base is already reprotonated but the pyranine still retains protons. Thus, it appears that the protonation of D96 is a slow rate-limiting process that generates a "proton hole" in the cytoplasmic section of the protein. The velocity of the hole propagation is modulated by the ionic strength of the solution and by selective replacements of charged residues on the interhelical loops of the protein, at domains that seems to be remote from the intraprotein proton conduction trajectory

Determination of a Unique Solution to Parallel Proton Transfer Reactions Using the Genetic Algorithm

Kinetic analysis of the dynamics as measured in multiequilibria systems is readily attained by curve-fitting methodologies, a treatment that can accurately retrace the shape of the measured signal. Still, these reconstructions are not related to the detailed mechanism of the process. In this study we subjected multiple proton transfer reactions to rigorous kinetic analysis, which consists of solving a set of coupled-nonlinear differential rate equations. The manual analysis of such systems can be biased by the operator; thus the analysis calls for impartial corroboration. What is more, there is no assurance that such a complex system has a unique solution. In this study, we used the Genetic Algorithm to investigate whether the solution of the system will converge into a single global minimum in the multidimensional parameter space. The experimental system consisted of proton transfer between four proton-binding sites with seven independent adjustable parameters. The results of the search indicate that the solution is unique and all adjustable parameters converge into a single minimum in the multidimensional parameter space, thus corroborating the accuracy of the manual analysis

Subsecond Proton-Hole Propagation in Bacteriorhodopsin

The dynamics of proton transfer between the surface of purple membrane and the aqueous bulk have recently been investigated by the Laser Induced Proton Pulse Method. Following a Δ-function release of protons to the bulk, the system was seen to regain its state of equilibrium within a few hundreds of microseconds. These measurements set the time frame for the relaxation of any state of acid-base disequilibrium between the bacteriorhodopsin's surface and the bulk. It was also deduced that the released protons react with the various proton binding within less than 10 μs. In the present study, we monitored the photocycle and the proton-cycle of photo-excited bacteriorhodopsin, in the absence of added buffer, and calculated the proton balance between the Schiff base and the bulk phase in a time-resolved mode. It was noticed that the late phase of the M decay (beyond 1 ms) is characterized by a slow (subsecond) relaxation of disequilibrium, where the Schiff base is already reprotonated but the pyranine still retains protons. Thus, it appears that the protonation of D96 is a slow rate-limiting process that generates a “proton hole” in the cytoplasmic section of the protein. The velocity of the hole propagation is modulated by the ionic strength of the solution and by selective replacements of charged residues on the interhelical loops of the protein, at domains that seems to be remote from the intraprotein proton conduction trajectory