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.

Nearly complete structure of bacteriophage DT57C reveals architecture of head-to-tail interface and lateral tail fibers.

The T5 family of viruses are tailed bacteriophages characterized by a long non-contractile tail. The bacteriophage DT57C is closely related to the paradigmal T5 phage, though it recognizes a different receptor (BtuB) and features highly divergent lateral tail fibers (LTF). Considerable portions of T5-like phages remain structurally uncharacterized. Here, we present the structure of DT57C determined by cryo-EM, and an atomic model of the virus, which was further explored using all-atom molecular dynamics simulations. The structure revealed a unique way of LTF attachment assisted by a dodecameric collar protein LtfC, and an unusual composition of the phage neck constructed of three protein rings. The tape measure protein (TMP) is organized within the tail tube in a three-stranded parallel α-helical coiled coil which makes direct contact with the genomic DNA. The presence of the C-terminal fragment of the TMP that remains within the tail tip suggests that the tail tip complex returns to its original state after DNA ejection. Our results provide a complete atomic structure of a T5-like phage, provide insights into the process of DNA ejection as well as a structural basis for the design of engineered phages and future mechanistic studies.

Sub-millisecond conformational dynamics of the A2A adenosine receptor revealed by single-molecule FRET.

The complex pharmacology of G-protein-coupled receptors (GPCRs) is defined by their multi-state conformational dynamics. Single-molecule Förster Resonance Energy Transfer (smFRET) is well suited to quantify dynamics for individual protein molecules; however, its application to GPCRs is challenging. Therefore, smFRET has been limited to studies of inter-receptor interactions in cellular membranes and receptors in detergent environments. Here, we performed smFRET experiments on functionally active human A2A adenosine receptor (A2AAR) molecules embedded in freely diffusing lipid nanodiscs to study their intramolecular conformational dynamics. We propose a dynamic model of A2AAR activation that involves a slow (>2 ms) exchange between the active-like and inactive-like conformations in both apo and antagonist-bound A2AAR, explaining the receptor's constitutive activity. For the agonist-bound A2AAR, we detected faster (390 ± 80 µs) ligand efficacy-dependent dynamics. Our work establishes a general smFRET platform for GPCR investigations that can potentially be used for drug screening and/or mechanism-of-action studies.

Structural insights into the effects of glycerol on ligand binding to cytochrome P450.

New antitubercular drugs are vital due to the spread of resistant strains. Carbethoxyhexyl imidazole (CHImi) inhibits cytochrome P450 CYP124, which is a steroid-metabolizing enzyme that is important for the survival of Mycobacterium tuberculosis in macrophages. The available crystal structure of the CYP124-CHImi complex reveals two glycerol molecules in the active site. A 1.15 Šresolution crystal structure of the glycerol-free CYP124-CHimi complex reported here shows multiple conformations of CHImi and the CYP124 active site which were previously restricted by glycerol. Complementary molecular dynamics simulations show coherence of the ligand and enzyme conformations. Spectrophotometric titration confirmed the influence of glycerol on CHImi binding: the affinity decreases more than tenfold in glycerol-containing buffer. In addition, it also showed that glycerol has a similar effect on other azole and triazole CYP124 ligands. Together, these data show that glycerol may compromise structural-functional studies and impede rational drug-design campaigns.

Structural basis for receptor selectivity and inverse agonism in S1P5 receptors.

The bioactive lysophospholipid sphingosine-1-phosphate (S1P) acts via five different subtypes of S1P receptors (S1PRs) - S1P1-5. S1P5 is predominantly expressed in nervous and immune systems, regulating the egress of natural killer cells from lymph nodes and playing a role in immune and neurodegenerative disorders, as well as carcinogenesis. Several S1PR therapeutic drugs have been developed to treat these diseases; however, they lack receptor subtype selectivity, which leads to side effects. In this article, we describe a 2.2 Å resolution room temperature crystal structure of the human S1P5 receptor in complex with a selective inverse agonist determined by serial femtosecond crystallography (SFX) at the Pohang Accelerator Laboratory X-Ray Free Electron Laser (PAL-XFEL) and analyze its structure-activity relationship data. The structure demonstrates a unique ligand-binding mode, involving an allosteric sub-pocket, which clarifies the receptor subtype selectivity and provides a template for structure-based drug design. Together with previously published S1PR structures in complex with antagonists and agonists, our structure with S1P5-inverse agonist sheds light on the activation mechanism and reveals structural determinants of the inverse agonism in the S1PR family.

High-pressure crystallography shows noble gas intervention into protein-lipid interaction and suggests a model for anaesthetic action.

In this work we examine how small hydrophobic molecules such as inert gases interact with membrane proteins (MPs) at a molecular level. High pressure atmospheres of argon and krypton were used to produce noble gas derivatives of crystals of three well studied MPs (two different proton pumps and a sodium light-driven ion pump). The structures obtained using X-ray crystallography showed that the vast majority of argon and krypton binding sites were located on the outer hydrophobic surface of the MPs - a surface usually accommodating hydrophobic chains of annular lipids (which are known structural and functional determinants for MPs). In conformity with these results, supplementary in silico molecular dynamics (MD) analysis predicted even greater numbers of argon and krypton binding positions on MP surface within the bilayer. These results indicate a potential importance of such interactions, particularly as related to the phenomenon of noble gas-induced anaesthesia.

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.

Structure and dynamics of the SARS-CoV-2 envelope protein monomer.

Coronaviruses, especially severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), present an ongoing threat to human wellbeing. Consequently, elucidation of molecular determinants of their function and interaction with the host is an important task. Whereas some of the coronaviral proteins are extensively characterized, others remain understudied. Here, we use molecular dynamics simulations to analyze the structure and dynamics of the SARS-CoV-2 envelope (E) protein (a viroporin) in the monomeric form. The protein consists of the hydrophobic α-helical transmembrane domain (TMD) and amphiphilic α-helices H2 and H3, connected by flexible linkers. We show that TMD has a preferable orientation in the membrane, while H2 and H3 reside at the membrane surface. Orientation of H2 is strongly influenced by palmitoylation of cysteines Cys40, Cys43, and Cys44. Glycosylation of Asn66 affects the orientation of H3. We also observe that the monomeric E protein both generates and senses the membrane curvature, preferably localizing with the C-terminus at the convex regions of the membrane; the protein in the pentameric form displays these properties as well. Localization to curved regions may be favorable for assembly of the E protein oligomers, whereas induction of curvature may facilitate the budding of the viral particles. The presented results may be helpful for a better understanding of the function of the coronaviral E protein and viroporins in general, and for overcoming the ongoing SARS-CoV-2 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.

Ambiguities in and completeness of SAS data analysis of membrane proteins: the case of the sensory rhodopsin II-transducer complex.

Membrane proteins (MPs) play vital roles in the function of cells and are also major drug targets. Structural information on proteins is vital for understanding their mechanism of function and is critical for the development of drugs. However, obtaining high-resolution structures of membrane proteins, in particular, under native conditions is still a great challenge. In such cases, the low-resolution methods small-angle X-ray and neutron scattering (SAXS and SANS) might provide valuable structural information. However, in some cases small-angle scattering (SAS) provides ambiguous ab initio structural information if complementary measurements are not performed and/or a priori information on the protein is not taken into account. Understanding the nature of the limitations may help to overcome these problems. One of the main problems of SAS data analysis of solubilized membrane proteins is the contribution of the detergent belt surrounding the MP. Here, a comprehensive analysis of how the detergent belt contributes to the SAS data of a membrane-protein complex of sensory rhodopsin II with its cognate transducer from Natronomonas pharaonis (NpSRII-NpHtrII) was performed. The influence of the polydispersity of NpSRII-NpHtrII oligomerization is the second problem that is addressed here. It is shown that inhomogeneity in the scattering length density of the detergent belt surrounding a membrane part of the complex and oligomerization polydispersity significantly impacts on SAXS and SANS profiles, and therefore on 3D ab initio structures. It is described how both problems can be taken into account to improve the quality of SAS data treatment. Since SAS data for MPs are usually obtained from solubilized proteins, and their detergent belt and, to a certain extent, oligomerization polydispersity are sufficiently common phenomena, the approaches proposed in this work might be used in SAS studies of different MPs.

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.

Nitrate- and Nitrite-Sensing Histidine Kinases: Function, Structure, and Natural Diversity.

Under anaerobic conditions, bacteria may utilize nitrates and nitrites as electron acceptors. Sensitivity to nitrous compounds is achieved via several mechanisms, some of which rely on sensor histidine kinases (HKs). The best studied nitrate- and nitrite-sensing HKs (NSHKs) are NarQ and NarX from Escherichia coli. Here, we review the function of NSHKs, analyze their natural diversity, and describe the available structural information. In particular, we show that around 6000 different NSHK sequences forming several distinct clusters may now be found in genomic databases, comprising mostly the genes from Beta- and Gammaproteobacteria as well as from Bacteroidetes and Chloroflexi, including those from anaerobic ammonia oxidation (annamox) communities. We show that the architecture of NSHKs is mostly conserved, although proteins from Bacteroidetes lack the HAMP and GAF-like domains yet sometimes have PAS. We reconcile the variation of NSHK sequences with atomistic models and pinpoint the structural elements important for signal transduction from the sensor domain to the catalytic module over the transmembrane and cytoplasmic regions spanning more than 200 Å.

Molecular model of a sensor of two-component signaling system.

Two-component systems (TCS) are widespread signaling systems present in all domains of life. TCS typically consist of a signal receptor/transducer and a response regulator. The receptors (histidine kinases, chemoreceptors and photoreceptors) are often embedded in the membrane and have a similar modular structure. Chemoreceptors were shown to function in highly ordered arrays, with trimers of dimers being the smallest functional unit. However, much less is known about photoreceptors. Here, we use small-angle scattering (SAS) to show that detergent-solubilized sensory rhodopsin II in complex with its cognate transducer forms dimers at low salt concentration, which associate into trimers of dimers at higher buffer molarities. We then fit an atomistic model of the whole complex into the SAS data. The obtained results suggest that the trimer of dimers is "tripod"-shaped and that the contacts between the dimers occur only through their cytoplasmic regions, whereas the transmembrane regions remain unconnected.

Lipid Dynamics in Diisobutylene-Maleic Acid (DIBMA) Lipid Particles in Presence of Sensory Rhodopsin II.

Amphiphilic diisobutylene/maleic acid (DIBMA) copolymers extract lipid-encased membrane proteins from lipid bilayers in a detergent-free manner, yielding nanosized, discoidal DIBMA lipid particles (DIBMALPs). Depending on the DIBMA/lipid ratio, the size of DIBMALPs can be broadly varied which makes them suitable for the incorporation of proteins of different sizes. Here, we examine the influence of the DIBMALP sizes and the presence of protein on the dynamics of encased lipids. As shown by a set of biophysical methods, the stability of DIBMALPs remains unaffected at different DIBMA/lipid ratios. Coarse-grained molecular dynamics simulations confirm the formation of viable DIBMALPs with an overall size of up to 35 nm. Electron paramagnetic resonance spectroscopy of nitroxides located at the 5th, 12th or 16th carbon atom positions in phosphatidylcholine-based spin labels reveals that the dynamics of enclosed lipids are not altered by the DIBMALP size. The presence of the membrane protein sensory rhodopsin II from Natronomonas pharaonis (NpSRII) results in a slight increase in the lipid dynamics compared to empty DIBMALPs. The light-induced photocycle shows full functionality of DIBMALPs-embedded NpSRII and a significant effect of the protein-to-lipid ratio during preparation on the NpSRII dynamics. This study indicates a possible expansion of the applicability of the DIBMALP technology on studies of membrane protein-protein interaction and oligomerization in a constraining environment.

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.

Sensor Histidine Kinase NarQ Activates via Helical Rotation, Diagonal Scissoring, and Eventually Piston-Like Shifts.

Membrane-embedded sensor histidine kinases (HKs) and chemoreceptors are used ubiquitously by bacteria and archaea to percept the environment, and are often crucial for their survival and pathogenicity. The proteins can transmit the signal from the sensor domain to the catalytic kinase domain reliably over the span of several hundreds of angstroms, and regulate the activity of the cognate response regulator proteins, with which they form two-component signaling systems (TCSs). Several mechanisms of transmembrane signal transduction in TCS receptors have been proposed, dubbed (swinging) piston, helical rotation, and diagonal scissoring. Yet, despite decades of studies, there is no consensus on whether these mechanisms are common for all TCS receptors. Here, we extend our previous work on Escherichia coli nitrate/nitrite sensor kinase NarQ. We determined a crystallographic structure of the sensor-TM-HAMP fragment of the R50S mutant, which, unexpectedly, was found in a ligand-bound-like conformation, despite an inability to bind nitrate. Subsequently, we reanalyzed the structures of the ligand-free and ligand-bound NarQ and NarX sensor domains, and conducted extensive molecular dynamics simulations of ligand-free and ligand-bound wild type and mutated NarQ. Based on the data, we show that binding of nitrate to NarQ causes, first and foremost, helical rotation and diagonal scissoring of the α-helices at the core of the sensor domain. These conformational changes are accompanied by a subtle piston-like motion, which is amplified by a switch in the secondary structure of the linker between the sensor and TM domains. We conclude that helical rotation, diagonal scissoring, and piston are simply different degrees of freedom in coiled-coil proteins and are not mutually exclusive in NarQ, and likely in other nitrate sensors and TCS proteins as well.

Molecular mechanism of light-driven sodium pumping.

The light-driven sodium-pumping rhodopsin KR2 from Krokinobacter eikastus is the only non-proton cation active transporter with demonstrated potential for optogenetics. However, the existing structural data on KR2 correspond exclusively to its ground state, and show no sodium inside the protein, which hampers the understanding of sodium-pumping mechanism. Here we present crystal structure of the O-intermediate of the physiologically relevant pentameric form of KR2 at the resolution of 2.1 Å, revealing a sodium ion near the retinal Schiff base, coordinated by N112 and D116 of the characteristic NDQ triad. We also obtained crystal structures of D116N and H30A variants, conducted metadynamics simulations and measured pumping activities of putative pathway mutants to demonstrate that sodium release likely proceeds alongside Q78 towards the structural sodium ion bound between KR2 protomers. Our findings highlight the importance of pentameric assembly for sodium pump function, and may be used for rational engineering of enhanced optogenetic tools.

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.