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.
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.
Thirteen tubulin protofilaments, made of αβ-tubulin heterodimers, interact laterally to produce cytoskeletal microtubules. Microtubules exhibit the striking property of dynamic instability, ...manifested in their intermittent growth and shrinkage at both ends. This behavior is key to many cellular processes, such as cell division, migration, maintenance of cell shape, etc. Although assembly and disassembly of microtubules is known to be linked to hydrolysis of a guanosine triphosphate molecule in the pocket of β-tubulin, detailed mechanistic understanding of corresponding conformational changes is still lacking. Here we take advantage of the recent generation of in-microtubule structures of tubulin to examine the properties of protofilaments, which serve as important microtubule assembly and disassembly intermediates. We find that initially straight tubulin protofilaments, relax to similar non-radially curved and slightly twisted conformations. Our analysis further suggests that guanosine triphosphate hydrolysis primarily affects the flexibility and conformation of the inter-dimer interface, without a strong impact on the shape or flexibility of αβ-heterodimer. Inter-dimer interfaces are significantly more flexible compared to intra-dimer interfaces. We argue that such a difference in flexibility could be key for distinct stability of the plus and minus microtubule ends. The higher flexibility of the inter-dimer interface may have implications for development of pulling force by curving tubulin protofilaments during microtubule disassembly, a process of major importance for chromosome motions in mitosis.
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.
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.
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.
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.
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.
Ambiguities in and the completeness of SAS data analysis of membrane proteins (MPs) are considered here in the case of the sensory rhodopsin II–transducer complex solubilized in a detergent. The contribution of the detergent belt surrounding MPs is one of the main problems in SAS data analysis of MPs; the second problem addressed here is the influence of oligomerization polydispersity (which is a sufficiently common phenomenon for MPs) on the quality of SAS structural analysis.
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.
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.