The Increased Amyloidogenicity of Spike RBD and pH-Dependent Binding to ACE2 May Contribute to the Transmissibility and Pathogenic Properties of SARS-CoV-2 Omicron as Suggested by In Silico Study.

SARS-CoV-2 is a rapidly evolving pathogen that has caused a global pandemic characterized by several consecutive waves. Based on epidemiological and NGS data, many different variants of SARS-CoV-2 were described and characterized since the original variant emerged in Wuhan in 2019. Notably, SARS-CoV-2 variants differ in transmissibility and pathogenicity in the human population, although the molecular basis for this difference is still debatable. A significant role is attributed to amino acid changes in the binding surface of the Spike protein to the ACE2 receptor, which may facilitate virus entry into the cell or contribute to immune evasion. We modeled in silico the interaction between Spike RBDs of Wuhan-Hu-1, Delta, and Omicron BA.1 variants and ACE2 at different pHs (pH 5 and pH 7) and showed that the strength of this interaction was higher for the Omicron BA.1 RBD compared to Wuhan-Hu-1 or Delta RBDs and that the effect was more profound at pH 5. This finding is strikingly related to the increased ability of Omicron variants to spread in the population. We also noted that during its spread in the population, SARS-CoV-2 evolved to a more charged, basic composition. We hypothesize that the more basic surface of the Omicron variant may facilitate its spread in the upper respiratory tract but not in the lower respiratory tract, where pH estimates are different. We calculated the amyloidogenic properties of Spike RBDs in different SARS-CoV-2 variants and found eight amyloidogenic regions in the Spike RBDs for each of the variants predicted by the FoldAmyloid program. Although all eight regions were almost identical in the Wuhan to Gamma variants, two of them were significantly longer in both Omicron variants, making the Omicron RBD more amyloidogenic. We discuss how the increased predicted amyloidogenicity of the Omicron variants RBDs may be important for protein stability, influence its interaction with ACE2 and contribute to immune evasion.

Is It Possible to Find an Antimicrobial Peptide That Passes the Membrane Bilayer with Minimal Force Resistance? An Attempt at a Predictive Approach by Molecular Dynamics Simulation.

There is still no answer to the mechanism of penetration of AMP peptides through the membrane bilayer. Several mechanisms for such a process have been proposed. It is necessary to understand whether it is possible, using the molecular dynamics method, to determine the ability of peptides of different compositions and lengths to pass through a membrane bilayer. To explain the passage of a peptide through a membrane bilayer, a method for preparing a membrane phospholipid bilayer was proposed, and 656 steered molecular dynamics calculations were carried out for pulling 7 amyloidogenic peptides with antimicrobial potential, and monopeptides (homo-repeats consisting of 10 residues of the same amino acid: Poly (Ala), Poly (Leu), Poly (Met), Poly (Arg), and Poly (Glu)) with various sequences through the membrane. Among the 15 studied peptides, the peptides exhibiting the least force resistance when passing through the bilayer were found, and the maximum reaction occurred at the boundary of the membrane bilayer entry. We found that the best correlation between the maximum membrane reaction force and the calculated parameters corresponds to the instability index (the correlation coefficient is above 0.9). One of the interesting results of this study is that the 10 residue amyloidogenic peptides and their extended peptides, with nine added residue cell-penetrating peptides and four residue linkers, both with established antimicrobial activity, have the same bilayer resistance force. All calculated data are summarized and posted on the server.

Signal transfer in human protein tyrosine phosphatase PTP1B from allosteric inhibitor P00058.

Protein tyrosine phosphatases constitute a family of cytosolic and receptor-like signal transducing enzymes that catalyze the hydrolysis of phospho-tyrosine residues of phosphorylated proteins. PTP1B, encoded by PTPN1, is a key negative regulator of insulin and leptin receptor signaling, linking it to two widespread diseases: type 2 diabetes mellitus and obesity. Here, we present crystal structures of the PTP1B apo-enzyme and a complex with a newly identified allosteric inhibitor, 2-(2,5-dimethyl-pyrrol-1-yl)-5-hydroxy-benzoic acid, designated as P00058. The inhibitor binding site is located about 18 Å away from the active center. However, the inhibitor causes significant re-arrangements in the active center of enzyme: residues 45-50 of catalytic Tyr-loop are shifted at their Cα-atom positions by 2.6 to 5.8 Å. We have identified an event of allosteric signal transfer from the inhibitor to the catalytic area using molecular dynamic simulation. Analyzing change of complex structure along the fluctuation trajectory we have found the large Cα-atom shifts in external strand, residues 25-40, which occur at the same time with the shifts in adjacent catalytic p-Tyr-loop. Coming of the signal to this loop arises due to dynamic fluctuation of protein structure at about 4.0 nanoseconds after the inhibitor takes up its space. Communicated by Ramaswamy H. Sarma.

Is there an advantageous arrangement of aromatic residues in proteins? Statistical analysis of aromatic interactions in globular proteins.

The aim of this study was to evaluate the favorability of different conformations of aromatic residues in proteins by analysing the occurrence of particular conformations. The clustering of protein structures from the Protein Data Bank (PDB) was performed. Conformations of interacting aromatic residues were analyzed for 511 282 pairs in 35 493 protein structures sharing less than 50% identity. Pairs with a parallel arrangement of aromatic residues made up 6.2% of all possible ones, which was twice as much as expected. Pairs with a perpendicular arrangement of aromatic residues made up 25%. We demonstrate that the most favorable arrangement was at an angle of 60° between the interacting aromatic residues. Among all possible aromatic pairs, the His-His pair was twice as frequent as expected, and the His-Phe pair was less frequent than expected. A server (CARP - Contacts of Aromatic Residues in Proteins) has been created for calculating essential structural features of interacting aromatic residues in proteins: http://bioproteom.protres.ru/arom_q_prog/.

Disordered Residues and Patterns in the Protein Data Bank.

We created a new library of disordered patterns and disordered residues in the Protein Data Bank (PDB). To obtain such datasets, we clustered the PDB and obtained the groups of chains with different identities and marked disordered residues. We elaborated a new procedure for finding disordered patterns and created a new version of the library. This library includes three sets of patterns: unique patterns, patterns consisting of two kinds of amino acids, and homo-repeats. Using this database, the user can: (1) find homologues in the entire Protein Data Bank; (2) perform a statistical analysis of disordered residues in protein structures; (3) search for disordered patterns and homo-repeats; (4) search for disordered regions in different chains of the same protein; (5) download clusters of protein chains with different identity from our database and library of disordered patterns; and (6) observe 3D structure interactively using MView. A new library of disordered patterns will help improve the accuracy of predictions for residues that will be structured or unstructured in a given region.

Comparative mechanical unfolding studies of spectrin domains R15, R16 and R17.

Spectrins belong to repetitive three-helix bundle proteins that have vital functions in multicellular organisms and are of potential value in nanotechnology. To reveal the unique physical features of repeat proteins we have studied the structural and mechanical properties of three repeats of chicken brain α-spectrin (R15, R16 and R17) at the atomic level under stretching at constant velocities (0.01, 0.05 and 0.1 ŷps-1) and constant forces (700 and 900 pN) using molecular dynamics (MD) simulations at T = 300 K. 114 independent MD simulations were performed and their analysis has been done. Despite structural similarity of these domains we have found that R15 is less mechanically stable than R16, which is less stable than R17. This result is in agreement with the thermal unfolding rates. Moreover, we have observed the relationship between mechanical stability, flexibility of the domains and the number of aromatic residues involved in aromatic clusters.

Mechanical stability analysis of the protein L immunoglobulin-binding domain by full alanine screening using molecular dynamics simulations.

This article is the first to study the mechanical properties of the immunoglobulin-binding domain of protein L (referred to as protein L) and its mutants at the atomic level. In the structure of protein L, each amino acid residue (except for alanines and glycines) was replaced sequentially by alanine. Thus, 49 mutants of protein L were obtained. The proteins were stretched at their termini at constant velocity using molecular dynamics simulations in water, i.e. by forced unfolding. 19 out of 49 mutations resulted in a large decrease of mechanical protein stability. These amino acids were affecting either the secondary structure (11 mutations) or loop structures (8 mutations) of protein L. Analysis of mechanical unfolding of the generated protein that has the same topology as protein L but consists of only alanines and glycines allows us to suggest that the mechanical stability of proteins, and specifically protein L, is determined by interactions between certain amino acid residues, although the unfolding pathway depends on the protein topology. This insight can now be used to modulate the mechanical properties of proteins and their unfolding pathways in the desired direction for using them in various biochips, biosensors and biomaterials for medicine, industry, and household purposes.

Right- and left-handed three-helix proteins. II. Similarity and differences in mechanical unfolding of proteins.

Here, we study mechanical properties of eight 3-helix proteins (four right-handed and four left-handed ones), which are similar in size under stretching at a constant speed and at a constant force on the atomic level using molecular dynamics simulations. The analysis of 256 trajectories from molecular dynamics simulations with explicit water showed that the right-handed three-helix domains are more mechanically resistant than the left-handed domains. Such results are observed at different extension velocities studied (192 trajectories obtained at the following conditions: v = 0.1, 0.05, and 0.01 Å ps(-1) , T = 300 K) and under constant stretching force (64 trajectories, F = 800 pN, T = 300 K). We can explain this by the fact, at least in part, that the right-handed domains have a larger number of contacts per residue and the radius of cross section than the left-handed domains.