On the thermodynamics of DNA in the Peyrard-Bishop-Dauxois model

Based on the direct method of molecular dynamics, the temperature dependence of energy, partition function, entropy and free energy, describing the first-order phase transition, is calculated for the classical Peyrard-Bishop-Dauxois model of the DNA molecule. It is shown that without taking into account the quantum freezing of the degrees of freedom at low temperatures, the relationship between these quantities is not determined. An estimation was given for the number of bubbles in homogeneous DNA chains: PolyA/PolyT and PolyG/PolyC.Comment: 4 pages, 2 figure

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

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

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 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 par allel 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 protein