Colabind: A Cloud-Based Approach for Prediction of Binding Sites Using Coarse-Grained Simulations with Molecular Probes.

Binding site prediction is a crucial step in understanding protein-ligand and protein-protein interactions (PPIs) with broad implications in drug discovery and bioinformatics. This study introduces Colabind, a robust, versatile, and user-friendly cloud-based approach that employs coarse-grained molecular dynamics simulations in the presence of molecular probes, mimicking fragments of drug-like compounds. Our method has demonstrated high effectiveness when validated across a diverse range of biological targets spanning various protein classes, successfully identifying orthosteric binding sites, as well as known druggable allosteric or PPI sites, in both experimentally determined and AI-predicted protein structures, consistently placing them among the top-ranked sites. Furthermore, we suggest that careful inspection of the identified regions with a high affinity for specific probes can provide valuable insights for the development of pharmacophore hypotheses. The approach is available at https://github.com/porekhov/CG_probeMD.

Immune Escape Associated with RBD Omicron Mutations and SARS-CoV-2 Evolution Dynamics.

The evolution and the emergence of new mutations of viruses affect their transmissibility and/or pathogenicity features, depending on different evolutionary scenarios of virus adaptation to the host. A typical trade-off scenario of SARS-CoV-2 evolution has been proposed, which leads to the appearance of an Omicron strain with lowered lethality, yet enhanced transmissibility. This direction of evolution might be partly explained by virus adaptation to therapeutic agents and enhanced escape from vaccine-induced and natural immunity formed by other SARS-CoV-2 strains. Omicron's high mutation rate in the Spike protein, as well as its previously described high genome mutation rate (Kandeel et al., 2021), revealed a gap between it and other SARS-CoV-2 strains, indicating the absence of a transitional evolutionary form to the Omicron strain. Therefore, Omicron has emerged as a new serotype divergent from the evolutionary lineage of other SARS-CoV-2 strains. Omicron is a rapidly evolving variant of high concern, whose new subvariants continue to manifest. Its further understanding and the further monitoring of key mutations that provide virus immune escape and/or high affinity towards the receptor could be useful for vaccine and therapeutic development in order to control the evolutionary direction of the COVID-19 pandemic.

Mechanisms of Formation, Structure, and Dynamics of Lipoprotein Discs Stabilized by Amphiphilic Copolymers: A Comprehensive Review.

Amphiphilic copolymers consisting of alternating hydrophilic and hydrophobic units account for a major recent methodical breakthrough in the investigations of membrane proteins. Styrene-maleic acid (SMA), diisobutylene-maleic acid (DIBMA), and related copolymers have been shown to extract membrane proteins directly from lipid membranes without the need for classical detergents. Within the particular experimental setup, they form disc-shaped nanoparticles with a narrow size distribution, which serve as a suitable platform for diverse kinds of spectroscopy and other biophysical techniques that require relatively small, homogeneous, water-soluble particles of separate membrane proteins in their native lipid environment. In recent years, copolymer-encased nanolipoparticles have been proven as suitable protein carriers for various structural biology applications, including cryo-electron microscopy (cryo-EM), small-angle scattering, and conventional and single-molecule X-ray diffraction experiments. Here, we review the current understanding of how such nanolipoparticles are formed and organized at the molecular level with an emphasis on their chemical diversity and factors affecting their size and solubilization efficiency.

Computational Analysis of Mutations in the Receptor-Binding Domain of SARS-CoV-2 Spike and Their Effects on Antibody Binding.

Currently, SARS-CoV-2 causing coronavirus disease 2019 (COVID-19) is responsible for one of the most deleterious pandemics of our time. The interaction between the ACE2 receptors at the surface of human cells and the viral Spike (S) protein triggers the infection, making the receptor-binding domain (RBD) of the SARS-CoV-2 S-protein a focal target for the neutralizing antibodies (Abs). Despite the recent progress in the development and deployment of vaccines, the emergence of novel variants of SARS-CoV-2 insensitive to Abs produced in response to the vaccine administration and/or monoclonal ones represent a potential danger. Here, we analyzed the diversity of neutralizing Ab epitopes and assessed the possible effects of single and multiple mutations in the RBD of SARS-CoV-2 S-protein on its binding affinity to various antibodies and the human ACE2 receptor using bioinformatics approaches. The RBD-Ab complexes with experimentally resolved structures were grouped into four clusters with distinct features at sequence and structure level. The performed computational analysis indicates that while single amino acid replacements in RBD may only cause partial impairment of the Abs binding, moreover, limited to specific epitopes, the variants of SARS-CoV-2 with multiple mutations, including some which were already detected in the population, may potentially result in a much broader antigenic escape. Further analysis of the existing RBD variants pointed to the trade-off between ACE2 binding and antigenic escape as a key limiting factor for the emergence of novel SAR-CoV-2 strains, as the naturally occurring mutations in RBD tend to reduce its binding affinity to Abs but not to ACE2. The results provide guidelines for further experimental studies aiming to identify high-risk RBD mutations that allow for an antigenic escape.

Synergistic Effect of Chemical Penetration Enhancers on Lidocaine Permeability Revealed by Coarse-Grained Molecular Dynamics Simulations.

The search for new formulations for transdermal drug delivery (TDD) is an important field in medicine and cosmetology. Molecules with specific physicochemical properties which can increase the permeability of active ingredients across the stratum corneum (SC) are called chemical penetration enhancers (CPEs), and it was shown that some CPEs can act synergistically. In this study, we performed coarse-grained (CG) molecular dynamics (MD) simulations of the lidocaine delivery facilitated by two CPEs-linoleic acid (LA) and ethanol-through the SC model membrane containing cholesterol, N-Stearoylsphingosine (DCPE), and behenic acid. In our simulations, we probed the effects of individual CPEs as well as their combination on various properties of the SC membrane and the lidocaine penetration across it. We demonstrated that the addition of both CPEs decreases the membrane thickness and the order parameters of the DPCE hydrocarbon chains. Moreover, LA also enhances diffusion of the SC membrane components, especially cholesterol. The estimated potential of mean force (PMF) profiles for the lidocaine translocation across SC in the presence/absence of two individual CPEs and their combination demonstrated that while ethanol lowers the free energy barrier for lidocaine to enter SC, LA decreases the depth of the free energy minima for lidocaine inside SC. These two effects supposedly result in synergistic penetration enhancement of drugs. Altogether, the present simulations provide a detailed molecular picture of CPEs' action and their synergistic effect on the penetration of small molecular weight therapeutics that can be beneficial for the design of novel drug and cosmetics formulations.

Strategic aromatic residues in the catalytic cleft of the xyloglucanase MtXgh74 modifying thermostability, mode of enzyme action, and viscosity reduction ability.

The thermostable endo-processive xyloglucanase MtXgh74 from Myceliophthora thermophila was used to study the influence of aromatic amino acids in the catalytic cleft on the mode of action and the ability of enzyme to reduce xyloglucan viscosity. The enzyme derivative Mut I with mutations W64A/W67A in the "negative" subsites of the catalytic cleft resulted in a 5.5-fold increase of the Km value. Mut I produced oligosaccharides of various lengths in addition to xyloglucan building blocks. The W320A/W321A substitutions in the "positive" subsites of the mutated enzyme Mut II catalytic cleft increased the Km value 54-fold and resulted in an endo-dissociative mode of action. The ability of Mut II to reduce the viscosity of xyloglucan at 50 °C was much better than that of other MtXgh74 variants. Besides, Mut II efficiently reduced viscosity of a natural substrate, the pulp of xyloglucan-containing tamarind seed flour. The Km, Vmax, and kcat values and viscosity reduction ability of the enzyme derivative Mut III (W320A/W321A/G446Y) returned to levels close to that of MtXgh74. The pattern of xyloglucan hydrolysis by Mut III was typical for endo-processive xyloglucanases. The thermostability of Mut I and Mut II at 60 °C decreased significantly compared to the wild type, whereas the thermostability of Mut III at 60 °C restored almost to the MtXgh74-wt value. All mutants lost the ability to cleave the backbone of xyloglucan building blocks which was a characteristic of MtXgh74. Instead they acquired a low branch removing activity. Molecular dynamics simulations revealed the role of mutated amino acids in the complex action mechanism of GH74 enzymes. KEY POINTS: • Endo-processive mode of action of the xyloglucanase MtXgh74 was altered by rational design. • The endo-dissociative mutant Mut II (W320A/W321A) efficiently reduced XyG viscosity. • The substitutions W320A/W321A/G446Y in Mut III recovered the endo-processive mode. • Mut II can be used to reduce the viscosity of biomass slurries containing tamarind seed flour.

CPE-DB: An Open Database of Chemical Penetration Enhancers.

The cutaneous delivery route currently accounts for almost 10% of all administered drugs and it is becoming more common. Chemical penetration enhancers (CPEs) increase the transport of drugs across skin layers by different mechanisms that depend on the chemical nature of the penetration enhancers. In our work, we created a chemical penetration enhancer database (CPE-DB) that is, to the best of our knowledge, the first CPE database. We collected information about known enhancers and their derivatives in a single database, and classified and characterized their molecular diversity in terms of scaffold content, key chemical moieties, molecular descriptors, etc. CPE-DB can be used for virtual screening and similarity search to identify new potent and safe enhancers, building quantitative structure-activity relationship (QSAR) and quantitative structure-property relationship (QSPR) models, and other machine-learning (ML) applications for the prediction of biological activity.

Using a viral 2A peptide-based strategy to reconstruct the bovine P450scc steroidogenic system in S. cerevisiae: Bovine P450scc system expression using 2A peptides.

Cytochrome P450scc system performs the first rate-limiting stage of steroidogenesis in mammals. The bovine P450scc system was reconstructed in Saccharomyces cerevisiae, using a foot-and-mouth disease virus 2A peptide (F2A)-based construct, to co-express cytochrome P450scc, adrenodoxin (Adx), and adrenodoxin reductase (AdR). During the translation of the self-processing fusion protein P450scc-F2A-Adx-F2A-AdR, the first and the second linkers are cleaved with different efficiencies (96 % and 11 %, respectively), resulting in the unbalanced expression of individual proteins. The low cleavage efficiency and the relative Adx and AdR protein levels were increased through replacing the second F2A peptide with different sequences and changing the order of Adx and AdR. The P450scc, AdR, and Adx sequences located upstream of the F2A affected F2A processing, to various degrees. Moreover, using molecular dynamics (MD) simulations, we showed that the 2A peptide fused to the C-terminus of Adx formed the steric hindrance during enzymatic complex formation, resulting in the reduction of catalytic activity. Thus, the functional activity of the reconstructed P450scc system was determined not only by the efficiency of 2A peptides but also by the overall sequence of the expressed 2A-polyprotein. Our results can be applied to the development of 2A-based co-translation strategies, to produce other multicomponent protein systems.

Radioprotectors.org: an open database of known and predicted radioprotectors.

The search for radioprotectors is an ambitious goal with many practical applications. Particularly, the improvement of human radioresistance for space is an important task, which comes into view with the recent successes in the space industry. Currently, all radioprotective drugs can be divided into two large groups differing in their effectiveness depending on the type of exposure. The first of these is radioprotectors, highly effective for pulsed, and some types of relatively short exposure to irradiation. The second group consists of long-acting radioprotectors. These drugs are effective for prolonged and fractionated irradiation. They also protect against impulse exposure to ionizing radiation, but to a lesser extent than short-acting radioprotectors. Creating a database on radioprotectors is a necessity dictated by the modern development of science and technology. We have created an open database, Radioprotectors.org, containing an up-to-date list of substances with proven radioprotective properties. All radioprotectors are annotated with relevant chemical and biological information, including transcriptomic data, and can be filtered according to their properties. Additionally, the performed transcriptomics analysis has revealed specific transcriptomic profiles of radioprotectors, which should facilitate the search for potent radioprotectors.

Lipid dynamics in nanoparticles formed by maleic acid-containing copolymers: EPR spectroscopy and molecular dynamics simulations.

Amphiphilic maleic acid-containing copolymers account for a recent methodical breakthrough in the study of membrane proteins. Their application enables a detergent-free extraction of membrane proteins from lipid bilayers, yielding stable water-soluble, discoidal lipid bilayer particles with incorporated proteins, which are wrapped with copolymers. Although many studies confirm the potential of this approach for membrane protein research, the interactions between the maleic acid-containing copolymers and extracted lipids, as well as possible effects of the copolymers on lipid-embedded proteins deserve further scrutinization. Here, we combine electron paramagnetic resonance spectroscopy and coarse-grain molecular dynamics simulations to compare the distribution and dynamics of lipids in lipid particles of phospholipid bilayers encased either by an aliphatic diisobutylene/maleic acid copolymer (DIBMALPs) or by an aromatic styrene/maleic acid copolymer (SMALPs). Nitroxides located at the 5th, 12th or 16th carbon atom positions in phosphatidylcholine-based spin labels experience restrictions of their reorientational motion depending on the type of encasing copolymer. The dynamics of the lipids was less constrained in DIBMALPs than in SMALPs with the affinity of spin labeled lipids to the polymeric rim being more pronounced in SMALPs.

Replicative and radiation-induced aging: a comparison of gene expression profiles.

All living organisms are subject to the aging process and experience the effect of ionizing radiation throughout their life. There have been a number of studies that linked ionizing radiation process to accelerated aging, but comprehensive signalome analysis of both processes was rarely conducted. Here we present a comparative signaling pathway based analysis of the transcriptomes of fibroblasts irradiated with different doses of ionizing radiation, replicatively aged fibroblasts and fibroblasts collected from young, middle age and old patients. We demonstrate a significant concordance between irradiation-induced and replicative senescence signalome signatures of fibroblasts. Additionally, significant differences in transcriptional response were also observed between fibroblasts irradiated with high and low dose. Our data shows that the transcriptome of replicatively aged fibroblasts is more similar to the transcriptome of the cells irradiated with 2 Gy, than with 5 сGy.This work revealed a number of signaling pathways that are shared between senescence and irradiation processes and can potentially be targeted by the new generation of gero- and radioprotectors.

Proton leakage across lipid bilayers: Oxygen atoms of phospholipid ester linkers align water molecules into transmembrane water wires.

Up to half of the cellular energy gets lost owing to membrane proton leakage. The permeability of lipid bilayers to protons is by several orders of magnitude higher than to other cations, which implies efficient proton-specific passages. The nature of these passages remains obscure. By combining experimental measurements of proton flow across phosphatidylcholine vesicles, steered molecular dynamics (MD) simulations of phosphatidylcholine bilayers and kinetic modelling, we have analyzed whether protons could pass between opposite phospholipid molecules when they sporadically converge. The MD simulations showed that each time, when the phosphorus atoms of the two phosphatidylcholine molecules got closer than 1.6 nm, the eight oxygen atoms of their ester linkages could form a transmembrane 'oxygen passage' along which several water molecules aligned into a water wire. Proton permeability along such water wires would be limited by rearrangement of oxygen atoms, which could explain the experimentally shown independence of the proton permeability of pH, H2O/D2O substitution, and membrane dipole potential. We suggest that protons can cross lipid bilayers by moving along short, self-sustaining water wires supported by oxygen atoms of lipid ester linkages.

Styrene/Maleic Acid Copolymers Form SMALPs by Pulling Lipid Patches out of the Lipid Bilayer.

Amphiphilic copolymers composed of styrene and maleic acid (SMA) monomers caused a major methodical breakthrough in the study of membrane proteins. They were found to directly release phospholipids and membrane proteins both from artificial and natural lipid bilayers, yielding stable water-soluble discoidal SMA/lipid particles (SMALPs) of uniform size. Although many empirical studies indicate the great potency of SMALPs for membrane protein research, the mechanisms of their formation remain obscure. It is unknown which factors account for the very assembly of SMALPs and govern their uniform size. We have developed a coarse-grained (CG) molecular model of SMA copolymers based on the MARTINI CG force field and used it to probe the behavior of SMA copolymers with varying composition/charge/concentration in solution as well as their interaction with lipid membranes. First, we found that SMA copolymers tend to aggregate in solution into clusters, which could account for the uniform size of SMALPs. Next, molecular dynamics (MD) simulations showed that periodic SMA copolymers with styrene/maleic acid ratios of 2:1 ([SSM] n) and 3:1 ([SSSM] n) differently interacted with lipid bilayers. While clusters of 2:1 SMA copolymers induced membrane poration, the clusters of 3:1 SMA copolymers extracted lipid patches from the membrane yielding SMALP-like structures. Extraction of lipid patches was also observed when we simulated the behavior of 3:1 copolymers with varying lengths and statistical distribution of styrene and MA units. Analysis of MD simulation trajectories and comparison with experimental data indicate that the formation of SMALPs requires copolymer molecules with a sufficient number of units made of more than two sequential styrene monomers.

Molecular Mechanism of Uptake of Cationic Photoantimicrobial Phthalocyanine across Bacterial Membranes Revealed by Molecular Dynamics Simulations.

Phthalocyanines are aromatic macrocyclic compounds, which are structurally related to porphyrins. In clinical practice, phthalocyanines are used in fluorescence imaging and photodynamic therapy of cancer and noncancer lesions. Certain forms of the substituted polycationic metallophthalocyanines have been previously shown to be active in photodynamic inactivation of both Gram-negative and Gram-positive bacteria; one of them is zinc octakis(cholinyl)phthalocyanine (ZnPcChol8+). However, the molecular details of how these compounds translocate across bacterial membranes still remain unclear. In the present work, we have developed a coarse-grained (CG) molecular model of ZnPcChol8+ within the framework of the popular MARTINI CG force field. The obtained model was used to probe the solvation behavior of phthalocyanine molecules, which agreed with experimental results. Subsequently, it was used to investigate the molecular details of interactions between phthalocyanines and membranes of various compositions. The results demonstrate that ZnPcChol8+ has high affinity to both the inner and the outer model membranes of Gram-negative bacteria, although this species does not show noticeable affinity to the 1-palmitoyl-2-oleoyl- sn-glycero-3-phosphatidylcholine membrane. Furthermore, we found out that the process of ZnPcChol8+ penetration toward the center of the outer bacterial membrane is energetically favorable and leads to its overall disturbance and formation of the aqueous pore. Such intramembrane localization of ZnPcChol8+ suggests their twofold cytotoxic effect on bacterial cells: (1) via induction of lipid peroxidation by enhanced production of reactive oxygen species (i.e., photodynamic toxicity); (2) via rendering the bacterial membrane more permeable for additional Pc molecules as well as other compounds. We also found that the kinetics of penetration depends on the presence of phospholipid defects in the lipopolysaccharide leaflet of the outer membrane and the type of counterions, which stabilize it. Thus, the results of our simulations provide a detailed molecular view of ZnPcChol8+ "self-promoted uptake", the pathway previously proposed for some small molecules crossing the outer bacterial membrane.

Elastic Membrane Deformations Govern Interleaflet Coupling of Lipid-Ordered Domains.

The mechanism responsible for domain registration in two membrane leaflets has thus far remained enigmatic. Using continuum elasticity theory, we show that minimum line tension is achieved along the rim between thicker (ordered) and thinner (disordered) domains by shifting the rims in opposing leaflets by a few nanometers relative to each other. Increasing surface tension yields an increase in line tension, resulting in larger domains. Because domain registration is driven by lipid deformation energy, it does not require special lipid components or interactions at the membrane midplane.

Comparative computational study of interaction of C60-fullerene and tris-malonyl-C60-fullerene isomers with lipid bilayer: relation to their antioxidant effect.

Oxidative stress induced by excessive production of reactive oxygen species (ROS) has been implicated in the etiology of many human diseases. It has been reported that fullerenes and some of their derivatives-carboxyfullerenes-exhibits a strong free radical scavenging capacity. The permeation of C60-fullerene and its amphiphilic derivatives-C3-tris-malonic-C60-fullerene (C3) and D3-tris-malonyl-C60-fullerene (D3)-through a lipid bilayer mimicking the eukaryotic cell membrane was studied using molecular dynamics (MD) simulations. The free energy profiles along the normal to the bilayer composed of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) for C60, C3 and D3 were calculated. We found that C60 molecules alone or in clusters spontaneously translocate to the hydrophobic core of the membrane and stay inside the bilayer during the whole period of simulation time. The incorporation of cluster of fullerenes inside the bilayer changes properties of the bilayer and leads to its deformation. In simulations of the tris-malonic fullerenes we discovered that both isomers, C3 and D3, adsorb at the surface of the bilayer but only C3 tends to be buried in the area of the lipid headgroups forming hydrophobic contacts with the lipid tails. We hypothesize that such position has implications for ROS scavenging mechanism in the specific cell compartments.