Evolutionary trajectory of receptor binding specificity and promiscuity of the spike protein of SARS-CoV-2.

SARS-CoV-2 infects cells by attachment to its receptor-the angiotensin converting enzyme 2 (ACE2). Regardless of the wealth of structural data, little is known about the physicochemical mechanism of interactions of the viral spike (S) protein with ACE2 and how this mechanism has evolved during the pandemic. Here, we applied experimental and computational approaches to characterize the molecular interaction of S proteins from SARS-CoV-2 variants of concern (VOC). Data on kinetics, activation-, and equilibrium thermodynamics of binding of the receptor binding domain (RBD) from VOC with ACE2 as well as data from computational protein electrostatics revealed a profound remodeling of the physicochemical characteristics of the interaction during the evolution. Thus, as compared to RBDs from Wuhan strain and other VOC, Omicron RBD presented as a unique protein in terms of conformational dynamics and types of non-covalent forces driving the complex formation with ACE2. Viral evolution resulted in a restriction of the RBD structural dynamics, and a shift to a major role of polar forces for ACE2 binding. Further, we investigated how the reshaping of the physicochemical characteristics of interaction affects the binding specificity of S proteins. Data from various binding assays revealed that SARS-CoV-2 Wuhan and Omicron RBDs manifest capacity for promiscuous recognition of unrelated human proteins, but they harbor distinct reactivity patterns. These findings might contribute for mechanistic understanding of the viral tropism and capacity to evade immune responses during evolution.

The key role of E418 carboxyl group in the formation of Nt.BspD6I nickase active site: Structural and functional properties of Nt.BspD6I E418A mutant.

The mutated nickase Nt.BspD6I E418A has been obtained by site-directed mutagenesis. The purified protein has been crystallized, and its spatial structure has been determined at 2.45 Å resolution. An analysis of the crystal structures of the wild-type and mutated nickase have shown that the elimination of a carboxyl group due to the E418A mutation initiates marked conformational changes in both the N-terminal recognition domain and the C-terminal catalytic domain of nickase and insignificantly affects its linker domain. This is supported by changes in the functional properties of mutated nickase: an increase in the oligomerization capacity in the presence of a substrate, a reduction in the capacity to bind a substrate, and complete loss of catalytic activity.

Structural analysis and molecular modelling of the Cu/Zn-SOD from fungal strain Humicola lutea 103.

The native form of Cu/Zn-superoxide dismutase, isolated from fungal strain Humicola lutea 103 is a homodimer that coordinates one Cu(2+) and one Zn(2+) per monomer. Cu(2+) and Zn(2+) ions play crucial roles in enzyme activity and structural stability, respectively. It was established that HLSOD shows high pH and temperature stability. Thermostability of the glycosylated enzyme Cu/Zn-SOD, isolated from fungal strain H. lutea 103, was determined by CD spectroscopy. Determination of reversibility toward thermal denaturation for HLSOD allowed several thermodynamic parameters to be calculated. In this communication we report the conditions under which reversible denaturation of HLSOD exists. The narrow range over which the system is reversible has been determined using the strongest test of two important thermodynamic independent variables (T and pH). Combining both these variables, the "phase diagram" was determined, as a result of which the real thermodynamic parameters (ΔC(p), ΔH(exp)°, and ΔG(exp)°) was established. Because very narrow pH-interval of transitions we assume they are as result of overlapping of two simple transitions. It was found that ΔH(o) is independent from pH with a value of 1.3 kcal/mol and 2.8 kcal/mol for the first and the second transition, respectively. ΔG(o) was pH-dependent in all studied pH-interval. This means that the transitions are entropically driven, these. Based on this, these processes can be described as hydrophobic rearrangement of the quaternary structure. It was also found that glycosylation does not influence the stability of the enzyme because the carbohydrate chain is exposed on the surface of the molecule.

Heme interacts with c1q and inhibits the classical complement pathway.

C1q is the recognition subunit of the first component of the classical complement pathway. It participates in clearance of immune complexes and apoptotic cells as well as in defense against pathogens. Inappropriate activation of the complement contributes to cellular and tissue damage in different pathologies, urging the need for the development of therapeutic agents that are able to inhibit the complement system. In this study, we report heme as an inhibitor of C1q. Exposure of C1q to heme significantly reduced the activation of the classical complement pathway, mediated by C-reactive protein (CRP) and IgG. Interaction analyses revealed that heme reduces the binding of C1q to CRP and IgG. Furthermore, we demonstrated that the inhibition of C1q interactions results from a direct binding of heme to C1q. Formation of complex of heme with C1q caused changes in the mechanism of recognition of IgG and CRP. Taken together, our data suggest that heme is a natural negative regulator of the classical complement pathway at the level of C1q. Heme may play a role at sites of excessive tissue damage and hemolysis where large amounts of free heme are released.

Disoxaril mutants of Coxsackievirus B1: phenotypic characteristics and analysis of the target VP1 gene.

Disoxaril inhibits enterovirus replication by binding to the hydrophobic pocket within the VP1 coat protein, thus stabilizing the virion and blocking its uncoating. Disoxaril-resistant (RES) mutants of the Coxsackievirus B1 (CVB1/RES) were derived from the wild disoxaril-sensitive (SOF) strain (CVB1/SOF) using a selection approach. A disoxaril-dependent (DEP) mutant (CVB1/DEP) was obtained following nine consecutive passages of the disoxaril-resistant mutant in the presence of disoxaril. Phenotypic characteristics of the disoxaril mutants were investigated. A timing-of-addition study of the CVB1/DEP replication demonstrated that in the absence of disoxaril the virus particle assembly stopped. VP1 RNA sequences of disoxaril mutants were compared with the existing Gen Bank CVB1 reference structure. The amino acid sequence of a large VP1 196-258 peptide (disoxaril-binding region) of CVB1/RES was significantly different from that of the CVB1/SOF. Crucially important changes in CVB1/RES were two point mutations, M213H and F237L, both in the ligand-binding pocket. The sequence analysis of the CVB1/DEP showed some reversion to CVB1/SOF. The amino acid sequences of the three VP1 proteins are presented.

Glycan structures and antiviral effect of the structural subunit RvH2 of Rapana hemocyanin.

Molluscan hemocyanins are very large biological macromolecules and they act as oxygen-transporting glycoproteins. Most of them are glycoproteins with molecular mass around 9000 kDa. The oligosaccharide structures of the structural subunit RvH2 of Rapana venosa hemocyanin (RvH) were studied by sequence analysis of glycans using MALDI-TOF-MS and tandem mass spectrometry on a Q-Trap mass spectrometer after enzymatical liberation of the N-glycans from the polypeptides. Our study revealed a highly heterogeneous mixture of glycans of the compositions Hex(0-9) HexNAc(2-4) Hex(0-3) Pent(0-3) Fuc(0-3). A novel type of N-glycan, with an internal fucose residue connecting one GalNAc(ß1-2) and one hexuronic acid, was detected, as also occurs in subunit RvH1. A glycan with the same structure but with two deoxyhexose residues was observed as a doubly charged ion. Antiviral effects of the native molecules of RvH and also of Helix lucorum hemocyanin (HlH), of their structural subunits, and of the glycosylated functional unit RvH2-e and the non-glycosylated unit RvH2-c on HSV virus type 1 were investigated. Only glycosylated FU RvH2-e exhibits this antiviral activity. The carbohydrate chains of the FU are likely to interact with specific regions of glycoproteins of HSV, through van der Waals interactions in general or with certain amino acid residues in particular. Several clusters of these residues can be identified on the surface of RvH2-e.

Structural analysis and molecular modeling of the RvH2-e functional unit of Rapana venosa hemocyanin.

Rapana venosa hemocyanin (RvH), a circulating glycoprotein of the marine snail, has a complex structure. To provide details on the stability of the protein, one functional unit, RvH2-e, was compared with the native molecule and the structural subunits, RvH1 and RvH2, via pH-T diagrams, typical phase portraits for stability and denaturation reversibility. By analyzing the T transition curves of RvH2-e at different pH values, several parameters of the thermodynamic functions were obtained. Increasing the temperature from 25°C to 55°C, the reversibility of the molecule of protein also increases, opening a reversibility window within the range of pH 4.0-8.0. On analyzing the pH transition curves, the start of the acid denaturation (below pH 6) and alkaline denaturation (above pH 9) was determined to be between 20°C and 35°C. For this range, the thermodynamic functions ΔH° and ΔG° for a standard temperature of 25°C were calculated.

Kinetics and thermodynamics of interaction of coagulation factor VIII with a pathogenic human antibody.

Replacement therapy in hemophilia A with exogenous coagulation factor VIII (FVIII) often results in the development of FVIII-neutralizing antibodies, referred to as inhibitors. Despite of large number of studies on the functional properties of FVIII inhibitors, detailed physicochemical characterization of their interactions is not available. Here we studied the biophysical mechanism of the interaction between a human pathogenic antibody--BO2C11 and its target antigen--FVIII. Kinetic and thermodynamic analyses implied that this interaction is not accompanied by significant conformational changes in the proteins. The data also suggested that association of BO2C11 to FVIII is driven mainly by a hydrophobic effect. The protein electrostatics however played a decisive role in this association. Thus, a gradual increase in ionic strength resulted in a considerable increase in the association rate of binding of BO2C11 to FVIII. Such an ionic strength-dependency is uncommon for other antibody-antigen interactions. Our data suggest that electrostatic effects observed for BO2C11-FVIII association may arise from high-energy penalty of desolvation of the charged residues at the binding interfaces. We hypothesize that untypical ionic strength dependence of association of BO2C11 to FVIII reflects the nature of the recognized epitope, namely a molecular surface involved in the binding of FVIII to phospholipids. The presented data provide mechanistic information about FVIII neutralization by an inhibitory antibody and also contribute to the understanding of the general mechanisms of antibody-antigen interactions.

PHEMTO: protein pH-dependent electric moment tools.

PHEMTO (protein pH-dependent electric moment tools) is released in response to the high demand in protein science community for evaluation of electrostatic characteristics in relations to molecular recognition. PHEMTO will serve protein scientists with new advanced features for analysis of protein molecular interactions: Electric/dipole moments, their pH-dependence and in silico charge mutagenesis effects on these properties as well as alternative algorithms for electric/dipole moment computation--Singular value decomposition of electrostatic potential (EP) to account for reaction field. The implementation is based on long-term experience--PHEI mean field electrostatics and PHEPS server for evaluation of global and local pH-dependent properties. However, PHEMTO is not just an update of our PHEPS server. Besides standard electrostatics, we offer new, advanced and useful features for analysis of protein molecular interactions. In addition our algorithms are very fast. Special emphasis is given to the interface--intuitive and user-friendly. The input is comprised of the atomic coordinate file in Protein Data Bank format. The advanced user is provided with a special input section for addition of non-polypeptide charges. The output covers actually full electrostatic characteristics but special emphasis is given to electric/dipole moments and their interactive visualization. PHEMTO server can be accessed at http://phemto.orgchm.bas.bg/.

Reversibility and "pH-T phase diagrams" of Rapana venosa hemocyanin and its structural subunits.

We have studied the stability and reassociation behaviour of native molecules of Rapana venosa hemocyanin and its two subunits, termed RvH1 and RvH2. In the presence of different concentrations of Ca(2+) and Mg(2+) ions and pH values, the subunits differ not only in their reassociation behaviour, but also in their formation of helical tubules and multidecamers. RvH1 revealed a greater stability at higher pH values compared to RvH2. Overall, the stability of reassociated RvH and its structural subunits was found to be pH-dependent. The increasing stability of native Hc and its subunits, shown by pH-induced CD transitions (acid and alkaline denaturation), can be explained with the formation of quaternary structure. The absence of a Cotton effect at temperatures 20-40 degrees C in the pH-transition curves of RvH2 indicates that this subunit is stabilized by additional "factors", e.g.: non-ionic/hydrophobic stabilization and interactions of carbohydrate moieties. A similar behaviour was observed for the T-transition curves in a wide pH interval for RvH and its structural subunits. At higher temperatures, many of the secondary structural elements are preserved especially at neutral pH, even at extreme high temperatures above 90 degrees C the protein structures resemble a "globule state".

Complement C1q-target proteins recognition is inhibited by electric moment effectors.

Classical complement pathway is an important innate immune mechanism, which is usually triggered by binding of C1q to immunoglobulins, pentraxins and other target molecules. Although the activation of the classical pathway is crucial in the host defence, its undesirable and uncontrolled activation can lead to tissue damage. Thus, understanding the molecular basis of complement activation and its inhibition are of great biomedical importance. Recently, we proposed a mechanism for target recognition and classical pathway activation by C1q, which is likely governed by calcium-controlled reorientation of macromolecular electric moment vectors. Here we sought to define the mechanism of C1q inhibition by low molecular weight disulphate compounds that bind to the globular (gC1q) domain, using experimental, computational docking and theoretical modelling approaches. Our experimental results demonstrate that betulin disulphate (B2S) and 9,9-bis(4'-hydroxyphenyl)fluorene disulphate (F2S) inhibit the interaction of C1q and its recombinant globular modules with target molecules IgG1, C-reactive protein (CRP) and long pentraxin 3 (PTX3). In most C1q-inhibitor docked complexes, there is a reduction of electric moment scalar values and similarly altered direction of electric/dipole moment vectors. This could explain the inhibitory effect by impaired electrostatic steering, lacking optimal target recognition and formation of functional complex. In the presence of the inhibitor, the tilt of gC1q domains is likely to be blocked by the altered direction of the electric moment vector. Thus, the transition from the inactive (closed) towards the active (open) conformation of C1q (i.e. the complement activation signal transmission) will be impaired and the cascade initiation disrupted. These results could serve as a starting point for the exploration of a new form of 'electric moment inhibitors/effectors'.

PHEPS: web-based pH-dependent Protein Electrostatics Server.

PHEPS (pH-dependent Protein Electrostatics Server) is a web service for fast prediction and experiment planning support, as well as for correlation and analysis of experimentally obtained results, reflecting charge-dependent phenomena in globular proteins. Its implementation is based on long-term experience (PHEI package) and the need to explain measured physicochemical characteristics at the level of protein atomic structure. The approach is semi-empirical and based on a mean field scheme for description and evaluation of global and local pH-dependent electrostatic properties: protein proton binding; ionic sites proton population; free energy electrostatic term; ionic groups proton affinities (pK(a,i)) and their Coulomb interaction with whole charge multipole; electrostatic potential of whole molecule at fixed pH and pH-dependent local electrostatic potentials at user-defined set of points. The speed of calculation is based on fast determination of distance-dependent pair charge-charge interactions as empirical three exponential function that covers charge-charge, charge-dipole and dipole-dipole contributions. After atomic coordinates input, all standard parameters are used as defaults to facilitate non-experienced users. Special attention was given to interactive addition of non-polypeptide charges, extra ionizable groups with intrinsic pK(a)s or fixed ions. The output information is given as plain-text, readable by 'RasMol', 'Origin' and the like. The PHEPS server is accessible at http://pheps.orgchm.bas.bg/home.html.

Role of Ca2+ in the electrostatic stability and the functional activity of the globular domain of human C1q.

C1q is the recognition subunit of the classical pathway of the complement system and a major connecting link between classical pathway-driven innate immunity and IgG- or IgM-mediated acquired immunity. The basic structural subunit of C1q is composed of an N-terminal triple-helical collagen-like region and a C-terminal heterotrimeric globular head domain (gC1q) that is made up of individual A, B, and C chains. Recent crystallographic studies have revealed that the gC1q domain, which is the main target-binding region of C1q, has a compact and spherical heterotrimeric assembly, held together by both electrostatic and nonpolar interactions, with quasi-3-fold symmetry. A characteristic feature of the gC1q domain is the presence of a exposed Ca(2+) located near the apex. We have investigated, using theoretical and experimental approaches, the role of Ca(2+) in the electrostatic stability and target-binding properties of the native C1q as well as recombinant monomeric forms of the C-terminal regions of the A, B, and C chains. Here, we report that Ca(2+) primarily influences the target recognition properties of C1q toward IgG, IgM, C-reactive protein, and pentraxin 3. At pH 7.4, the loss of Ca(2+) leads to changes in the direction of electric moment from coaxial (where the putative C-reactive protein-binding site is located) to perpendicular to the molecular axis (toward the most likely IgG-binding site), which appears important for target recognition by C1q and subsequent complement activation.

Structure and stability of arthropodan hemocyanin Limulus polyphemus.

In the hemolymph of many arthropodan species, respiratory copper proteins of high molecular weight, termed hemocyanins (Hcs) are dissolved. In this communication, we report on the protein stability of different hemocyanin species (Crustacea and Chelicerata) using fluorescence spectroscopy. Five to seven major electrophoretically separable protein chains (structural subunits) were purified by fast protein liquid chromatography (FPLC) ion exchange chromatography from different hemocyanins with very high sequence homology of the active site regions binding copper ions (CuA and CuB), and especially the relative sequence positions of histidine (His) and tryptophan (Trp) residues of these protein segments are in all cases identical. The conformational stabilities of the native dodecameric aggregates and their isolated structural subunits towards various denaturants (pH and guanidine hydrochloride (Gdn.HCl)) indicate that the quaternary structure is stabilized by hydrophilic and polar forces, whereby both, the oxy- and apo-forms of the protein are considered. These two classes of Crustacea and Chelicerata Hcs have the similar Trp-fluorescence quantum yields, but different values of lambda(max) emission (about 325 and 337 nm, respectively). Differences in the quantum yields are observed of the oxy- and apo-forms, which must be attributed to the fluorescence quenching effect of the two copper ions (CuA and CuB) in the active site. The position of emission maximum indicates tryptophan side chains are situated in a non-polar environment. Denaturation studies of Hcs by Gdn.HCl indicate that the denaturation process consists of two steps: dissociation of the native molecule into its structural subunits and denaturation of the subunits at concentrations >1.5M Gdn.HCl. Two steps of denaturation are also observed after keeping the protein in buffer solutions at different pH values with different pH-stability for holo-oxy and apo-Hc forms.

The role of the salt concentration, proton, and phosphate binding on the thermal stability of wild and cloned DNA-binding protein Sso7d from Sulfolobus solfataricus.

The acidic pH (1.5-7.0) and ionic strength (0.005-0.2M) dependence of thermodynamic functions of protein Sso7d from Sulfolobus solfataricus, cloned (c-Sso7d) and N-heptapeptide deleted [c-des(1-7)Sso7d] in glycine, and phosphate buffers was studied by means of adiabatic scanning calorimetry. The difference of proton binding was estimated from deltaHcal(pH), Td(pH), and (deltaTd/deltapH). It was found that a single group non-co-operative ionization with apparent pKa = 3.25 for both cloned and deleted proteins govern the thermal unfolding of two different (protonated and unprotonated) forms. deltaH degrees is found to be pH-independent and the changes in stability (deltaG degrees ) originate from changes in entropy terms. The apparent pKa measured at high salt concentrations decreases with 0.5 pH units from glycine to phosphate and the free energy of transfer at high ionic strength is 0.7 kcal/mol. The ionic strength dependence for the pH-dependent D-states is very different at pH 6.0 and 1.5. This is consistent with the property of denatured state to be more compacted or "closed" (Dc) at neutral or weak acidic pH and more random or "open" (Do) at acidic pH. From the Bjerrum's relation was found the number of screened charges important for the unfolding process. The main conclusions are: (1) the thermal stability of Sso7d has prominently entropic nature; (2) a single non-co-operative ionization controls the conformations in the D-state; and (3) pH-dependent conformational equilibrium could be functionally important in Sso7d-DNA recognition.

Mutational analyses of the recombinant globular regions of human C1q A, B, and C chains suggest an essential role for arginine and histidine residues in the C1q-IgG interaction.

The first step in the activation of the classical complement pathway by immune complexes involves the binding of the globular domain (gC1q) of C1q to the Fc regions of aggregated IgG or IgM. Each gC1q domain is a heterotrimer of the C-terminal halves of one A (ghA), one B (ghB), and one C (ghC) chain. Our recent studies have suggested a modular organization of gC1q, consistent with the view that ghA, ghB, and ghC are functionally autonomous modules and have distinct and differential ligand-binding properties. Although C1q binding sites on IgG have been previously identified, the complementary interacting sites on the gC1q domain have not been precisely defined. The availability of the recombinant constructs expressing ghA, ghB, and ghC has allowed us, for the first time, to engineer single-residue substitution mutations and identify residues on the gC1q domain, which are involved in the interaction between C1q and IgG. Because C1q is a charge pattern recognition molecule, we have sequentially targeted arginine and histidine residues in each chain. Consistent with previous chemical modification studies and the recent crystal structure of gC1q, our results support a central role for arginine and histidine residues, especially Arg(114) and Arg(129) of the ghB module, in the C1q-IgG interaction.

pH-dependent quenching of the fluorescence of tryptophan residues in class A beta-lactamase from E. coli (TEM-1).

We performed an investigation of the pH-dependent quenching of the fluorescence of tryptophan residues of TEM-1 beta-lactamase from E. coli by uncharged and charged quenchers. pH-dependent Stern-Volmer constants (Ksv/pH) of tryptophan residues allowed us to determine subtle but discrete structurally and functionally important processes.