Mutant analysis of Kcng4b reveals how the different functional states of the voltage-gated potassium channel regulate ear development.

The voltage gated (Kv) slow-inactivating delayed rectifier channel regulates the development of hollow organs of the zebrafish. The functional channel consists of the tetramer of electrically active Kcnb1 (Kv2.1) subunits and Kcng4b (Kv6.4) modulatory or electrically silent subunits. The two mutations in zebrafish kcng4b gene - kcng4b-C1 and kcng4b-C2 (Gasanov et al., 2021) - have been studied during ear development using electrophysiology, developmental biology and in silico structural modelling. kcng4b-C1 mutation causes a C-terminal truncation characterized by mild Kcng4b loss-of-function (LOF) manifested by failure of kinocilia to extend and formation of ectopic otoliths. In contrast, the kcng4b-C2-/- mutation causes the C-terminal domain to elongate and the ectopic seventh transmembrane (TM) domain to form, converting the intracellular C-terminus to an extracellular one. Kcng4b-C2 acts as a Kcng4b gain-of-function (GOF) allele. Otoliths fail to develop and kinocilia are reduced in kcng4b-C2-/-. These results show that different mutations of the silent subunit Kcng4 can affect the activity of the Kv channel and cause a wide range of developmental defects.

Molecular stripping underpins derepression of a toxin-antitoxin system.

Transcription factors control gene expression; among these, transcriptional repressors must liberate the promoter for derepression to occur. Toxin-antitoxin (TA) modules are bacterial elements that autoregulate their transcription by binding the promoter in a T:A ratio-dependent manner, known as conditional cooperativity. The molecular basis of how excess toxin triggers derepression has remained elusive, largely because monitoring the rearrangement of promoter-repressor complexes, which underpin derepression, is challenging. Here, we dissect the autoregulation of the Salmonella enterica tacAT3 module. Using a combination of assays targeting DNA binding and promoter activity, as well as structural characterization, we determine the essential TA and DNA elements required to control transcription, and we reconstitute a repression-to-derepression path. We demonstrate that excess toxin triggers molecular stripping of the repressor complex off the DNA through multiple allosteric changes causing DNA distortion and ultimately leading to derepression. Thus, our work provides important insight into the mechanisms underlying conditional cooperativity.

Dual-Activity Fluoroquinolone-Transportan 10 Conjugates Offer Alternative Leukemia Therapy during Hematopoietic Cell Transplantation.

Hematopoietic cell transplantation (HCT) is often considered a last resort leukemia treatment, fraught with limited success due to microbial infections, a leading cause of mortality in leukemia patients. To address this critical issue, we explored a novel approach by synthesizing antileukemic agents containing antibacterial substances. This innovative strategy involves conjugating fluoroquinolone antibiotics, such as ciprofloxacin (CIP) or levofloxacin (LVX), with the cell-penetrating peptide transportan 10 (TP10). Here, we demonstrate that the resultant compounds display promising biologic activities in preclinical studies. These novel conjugates not only exhibit potent antimicrobial effects but are also selective against leukemia cells. The cytotoxic mechanism involves rapid disruption of cell membrane asymmetry leading to membrane damage. Importantly, these conjugates penetrated mammalian cells, accumulating within the nuclear membrane without significant effect on cellular architecture or mitochondrial function. Molecular simulations elucidated the aggregation tendencies of TP10 conjugates within lipid bilayers, resulting in membrane disruption and permeabilization. Moreover, mass spectrometry analysis confirmed efficient reduction of disulfide bonds within TP10 conjugates, facilitating release and activation of the fluoroquinolone derivatives. Intriguingly, these compounds inhibited human topoisomerases, setting them apart from traditional fluoroquinolones. Remarkably, TP10 conjugates generated lower intracellular levels of reactive oxygen species compared with CIP and LVX. The combination of antibacterial and antileukemic properties, coupled with selective cytostatic effects and minimal toxicity toward healthy cells, positions TP10 derivatives as promising candidates for innovative therapeutic approaches in the context of antileukemic HCT. This study highlights their potential in search of more effective leukemia treatments. SIGNIFICANCE STATEMENT: Fluoroquinolones are commonly used antibiotics, while transportan 10 (TP10) is a cell-penetrating peptide (CPP) with anticancer properties. In HCT, microbial infections are the primary cause of illness and death. Combining TP10 with fluoroquinolones enhanced their effects on different cell types. The dual pharmacological action of these conjugates offers a promising proof-of-concept solution for leukemic patients undergoing HCT. Strategically designed therapeutics, incorporating CPPs with antibacterial properties, have the potential to reduce microbial infections in the treatment of malignancies.

Splice variants of mitofusin 2 shape the endoplasmic reticulum and tether it to mitochondria.

In eukaryotic cells, different organelles interact at membrane contact sites stabilized by tethers. Mitochondrial mitofusin 2 (MFN2) acts as a membrane tether that interacts with an unknown partner on the endoplasmic reticulum (ER). In this work, we identified the MFN2 splice variant ERMIT2 as the ER tethering partner of MFN2. Splicing of MFN2 produced ERMIT2 and ERMIN2, two ER-specific variants. ERMIN2 regulated ER morphology, whereas ERMIT2 localized at the ER-mitochondria interface and interacted with mitochondrial mitofusins to tether ER and mitochondria. This tethering allowed efficient mitochondrial calcium ion uptake and phospholipid transfer. Expression of ERMIT2 ameliorated the ER stress, inflammation, and fibrosis typical of liver-specific Mfn2 knockout mice. Thus, ER-specific MFN2 variants display entirely extramitochondrial MFN2 functions involved in interorganellar tethering and liver metabolic activities.

Molecular mechanism and energetics of coupling between substrate binding and product release in the F1-ATPase catalytic cycle.

F1-ATPase is a motor protein that couples the rotation of its rotary [Formula: see text] subunit with ATP synthesis or hydrolysis. Single-molecule experiments indicate that nucleotide binding and release events occur almost simultaneously during the synthesis cycle, allowing the energy gain due to spontaneous binding of ADP to one catalytic [Formula: see text] subunit to be directly harnessed for driving the release of ATP from another rather than being dissipated as heat. Here, we examine the unknown mechanism of this coupling that is critical for an exceptionally high mechanochemical efficiency of F1-ATPase by means of all-atom free-energy simulations. We find that nondissipative and kinetically fast progression of the motor in the synthesis direction requires a concerted conformational change involving the closure of the ADP-binding [Formula: see text] subunit followed by the gradual opening of the ATP-releasing [Formula: see text] subunit over the course of the 30 to 40° rotary substep of the [Formula: see text] subunit. This rotary substep, preceding the ATP-dependent metastable state, allows for the recovery of a large portion of the ADP binding energy in the conformation of ATP-bound [Formula: see text] that gradually adopts the low-affinity conformation, captured also by the recent cryo-EM structure of this elusive state. The release of ATP from this nearly open conformation leads to its further opening, which enables the progression of the motor to the next catalytic metastable state. Our simulations explain this energy conversion mechanism in terms of intersubunit and ligand-protein interactions.

How acidic amino acid residues facilitate DNA target site selection.

Despite the negative charge of the DNA backbone, acidic residues (Asp/Glu) commonly participate in the base readout, with a strong preference for cytosine. In fact, in the solved DNA/protein structures, cytosine is recognized almost exclusively by Asp/Glu through a direct hydrogen bond, while at the same time, adenine, regardless of its amino group, shows no propensity for Asp/Glu. Here, we analyzed the contribution of Asp/Glu to sequence-specific DNA binding using classical and ab initio simulations of selected transcription factors and found that it is governed by a fine balance between the repulsion from backbone phosphates and attractive interactions with cytosine. Specifically, Asp/Glu lower the affinity for noncytosine sites and thus act as negative selectors preventing off-target binding. At cytosine-containing sites, the favorable contribution does not merely rely on the formation of a single H-bond but usually requires the presence of positive potential generated by multiple cytosines, consistently with the observed excess of cytosine in the target sites. Finally, we show that the preference of Asp/Glu for cytosine over adenine is a result of the repulsion from the adenine imidazole ring and a tendency of purine-purine dinucleotides to adopt the BII conformation.

Omicron mutations increase interdomain interactions and reduce epitope exposure in the SARS-CoV-2 spike.

Omicron BA.1 is a highly infectious variant of SARS-CoV-2 that carries more than thirty mutations on the spike protein in comparison to the Wuhan wild type (WT). Some of the Omicron mutations, located on the receptor-binding domain (RBD), are exposed to the surrounding solvent and are known to help evade immunity. However, the impact of buried mutations on the RBD conformations and on the mechanics of the spike opening is less evident. Here, we use all-atom molecular dynamics (MD) simulations with metadynamics to characterize the thermodynamic RBD-opening ensemble, identifying significant differences between WT and Omicron. Specifically, the Omicron mutations S371L, S373P, and S375F make more RBD interdomain contacts during the spike's opening. Moreover, Omicron takes longer to reach the transition state than WT. It stabilizes up-state conformations with fewer RBD epitopes exposed to the solvent, potentially favoring immune or antibody evasion.

Correction: The structural role of SARS-CoV-2 genetic background in the emergence and success of spike mutations: The case of the spike A222V mutation.

[This corrects the article DOI: 10.1371/journal.ppat.1010631.].

DNAffinity: a machine-learning approach to predict DNA binding affinities of transcription factors.

We present a physics-based machine learning approach to predict in vitro transcription factor binding affinities from structural and mechanical DNA properties directly derived from atomistic molecular dynamics simulations. The method is able to predict affinities obtained with techniques as different as uPBM, gcPBM and HT-SELEX with an excellent performance, much better than existing algorithms. Due to its nature, the method can be extended to epigenetic variants, mismatches, mutations, or any non-coding nucleobases. When complemented with chromatin structure information, our in vitro trained method provides also good estimates of in vivo binding sites in yeast.

Pre-exascale HPC approaches for molecular dynamics simulations. Covid-19 research: A use case.

Exascale computing has been a dream for ages and is close to becoming a reality that will impact how molecular simulations are being performed, as well as the quantity and quality of the information derived for them. We review how the biomolecular simulations field is anticipating these new architectures, making emphasis on recent work from groups in the BioExcel Center of Excellence for High Performance Computing. We exemplified the power of these simulation strategies with the work done by the HPC simulation community to fight Covid-19 pandemics. This article is categorized under:Data Science > Computer Algorithms and ProgrammingData Science > Databases and Expert SystemsMolecular and Statistical Mechanics > Molecular Dynamics and Monte-Carlo Methods.

The structural role of SARS-CoV-2 genetic background in the emergence and success of spike mutations: The case of the spike A222V mutation.

The S:A222V point mutation, within the G clade, was characteristic of the 20E (EU1) SARS-CoV-2 variant identified in Spain in early summer 2020. This mutation has since reappeared in the Delta subvariant AY.4.2, raising questions about its specific effect on viral infection. We report combined serological, functional, structural and computational studies characterizing the impact of this mutation. Our results reveal that S:A222V promotes an increased RBD opening and slightly increases ACE2 binding as compared to the parent S:D614G clade. Finally, S:A222V does not reduce sera neutralization capacity, suggesting it does not affect vaccine effectiveness.

Mechanism of reaction of RNA-dependent RNA polymerase from SARS-CoV-2.

We combine molecular dynamics, statistical mechanics, and hybrid quantum mechanics/molecular mechanics simulations to describe mechanistically the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA-dependent RNA polymerase (RdRp). Our study analyzes the binding mode of both natural triphosphate substrates as well as remdesivir triphosphate (the active form of drug), which is bound preferentially over ATP by RdRp while being poorly recognized by human RNA polymerase II (RNA Pol II). A comparison of incorporation rates between natural and antiviral nucleotides shows that remdesivir is incorporated more slowly into the nascent RNA compared with ATP, leading to an RNA duplex that is structurally very similar to an unmodified one, arguing against the hypothesis that remdesivir is a competitive inhibitor of ATP. We characterize the entire mechanism of reaction, finding that viral RdRp is highly processive and displays a higher catalytic rate of incorporation than human RNA Pol II. Overall, our study provides the first detailed explanation of the replication mechanism of RdRp.

Molecular basis of Arginine and Lysine DNA sequence-dependent thermo-stability modulation.

We have used a variety of theoretical and experimental techniques to study the role of four basic amino acids-Arginine, Lysine, Ornithine and L-2,4-Diaminobutyric acid-on the structure, flexibility and sequence-dependent stability of DNA. We found that the presence of organic ions stabilizes the duplexes and significantly reduces the difference in stability between AT- and GC-rich duplexes with respect to the control conditions. This suggests that these amino acids, ingredients of the primordial soup during abiogenesis, could have helped to equalize the stability of AT- and GC-rich DNA oligomers, facilitating a general non-catalysed self-replication of DNA. Experiments and simulations demonstrate that organic ions have an effect that goes beyond the general electrostatic screening, involving specific interactions along the grooves of the double helix. We conclude that organic ions, largely ignored in the DNA world, should be reconsidered as crucial structural elements far from mimics of small inorganic cations.

Molecular mechanism of proton-coupled ligand translocation by the bacterial efflux pump EmrE.

The current surge in bacterial multi-drug resistance (MDR) is one of the largest challenges to public health, threatening to render ineffective many therapies we rely on for treatment of serious infections. Understanding different factors that contribute to MDR is hence crucial from the global "one health" perspective. In this contribution, we focus on the prototypical broad-selectivity proton-coupled antiporter EmrE, one of the smallest known ligand transporters that confers resistance to aromatic cations in a number of clinically relevant species. As an asymmetric homodimer undergoing an "alternating access" protomer-swap conformational change, it serves as a model for the mechanistic understanding of more complex drug transporters. Here, we present a free energy and solvent accessibility analysis that indicates the presence of two complementary ligand translocation pathways that remain operative in a broad range of conditions. Our simulations show a previously undescribed desolvated apo state and anticorrelated accessibility in the ligand-bound state, explaining on a structural level why EmrE does not disrupt the pH gradient through futile proton transfer. By comparing the behavior of a number of model charged and/or aromatic ligands, we also explain the origin of selectivity of EmrE towards a broad class of aromatic cations. Finally, we explore unbiased pathways of ligand entry and exit to identify correlated structural changes implicated in ligand binding and release, as well as characterize key intermediates of occupancy changes.

Effective Drug Concentration and Selectivity Depends on Fraction of Primitive Cells.

Poor efficiency of chemotherapeutics in the eradication of Cancer Stem Cells (CSCs) has been driving the search for more active and specific compounds. In this work, we show how cell density-dependent stage culture profiles can be used in drug development workflows to achieve more robust drug activity (IC50 and EC50) results. Using flow cytometry and light microscopy, we characterized the cytological stage profiles of the HL-60-, A-549-, and HEK-293-derived sublines with a focus on their primitive cell content. We then used a range of cytotoxic substances-C-123, bortezomib, idarubicin, C-1305, doxorubicin, DMSO, and ethanol-to highlight typical density-related issues accompanying drug activity determination. We also showed that drug EC50 and selectivity indices normalized to primitive cell content are more accurate activity measurements. We tested our approach by calculating the corrected selectivity index of a novel chemotherapeutic candidate, C-123. Overall, our study highlights the usefulness of accounting for primitive cell fractions in the assessment of drug efficiency.

Light-Modulated Sunscreen Mechanism in the Retina of the Human Eye.

The functioning of the human eye in the extreme range of light intensity requires a combination of the high sensitivity of photoreceptors with their photostability. Here, we identify a regulatory mechanism based on dynamic modulation of light absorption by xanthophylls in the retina, realized by reorientation of pigment molecules induced by trans-cis photoisomerization. We explore this photochemically switchable system using chromatographic analysis coupled with microimaging based on fluorescence lifetime and Raman scattering, showing it at work in both isolated human retina and model lipid membranes. The molecular mechanism underlying xanthophyll reorientation is explained in terms of hydrophobic mismatch using molecular dynamics simulations. Overall, we show that xanthophylls in the human retina act as "molecular blinds", opening and closing on a submillisecond timescale to dynamically control the intensity of light reaching the photoreceptors, thus enabling vision at a very low light intensity and protecting the retina from photodegradation when suddenly exposed to strong light.

Defining a novel domain that provides an essential contribution to site-specific interaction of Rep protein with DNA.

An essential feature of replication initiation proteins is their ability to bind to DNA. In this work, we describe a new domain that contributes to a replication initiator sequence-specific interaction with DNA. Applying biochemical assays and structure prediction methods coupled with DNA-protein crosslinking, mass spectrometry, and construction and analysis of mutant proteins, we identified that the replication initiator of the broad host range plasmid RK2, in addition to two winged helix domains, contains a third DNA-binding domain. The phylogenetic analysis revealed that the composition of this unique domain is typical within the described TrfA-like protein family. Both in vitro and in vivo experiments involving the constructed TrfA mutant proteins showed that the newly identified domain is essential for the formation of the protein complex with DNA, contributes to the avidity for interaction with DNA, and the replication activity of the initiator. The analysis of mutant proteins, each containing a single substitution, showed that each of the three domains composing TrfA is essential for the formation of the protein complex with DNA. Furthermore, the new domain, along with the winged helix domains, contributes to the sequence specificity of replication initiator interaction within the plasmid replication origin.

Self-assembly, stability and conductance of amphotericin B channels: bridging the gap between structure and function.

Amphotericin B (AmB), one of the most powerful but also toxic drugs used to treat systemic mycoses, is believed to selectively permeabilize fungal cell membranes to ions in a sterol-dependent manner. Unfortunately, the structure of the biologically active AmB channels has long eluded researchers, obstructing the design of safer alternatives. Here, we investigate the structural and thermodynamic aspects of channel formation, stability, and selective ion conduction. We combine fluorescence lifetime imaging and molecular simulations to trace the process of channel assembly until the formation of stable, roughly octameric double-length channels (DLCs). This stoichiometry is confirmed by matching the predicted channel conductances with the past results of patch-clamp measurements. We then use free energy calculations to explain the effect of sterols on DLC stability and discuss the observed cation selectivity in structural terms, addressing several long-standing controversies in the context of their physiological relevance. Simulations of ion permeation indicate that only solvated ions pass through DLCs, revealing surprising solvation patterns in the channel lumen. We conclude our investigation by inspecting the role of the tail hydroxyl in the assembly of functional channels, pointing at possible origins of the cholesterol-ergosterol selectivity.

Molywood: streamlining the design and rendering of molecular movies.

MOTIVATION: High-quality dynamic visuals are needed at all levels of science communication, from the conference hall to the classroom. As scientific journals embrace new article formats, many key concepts-particularly, in structural biology-are also more easily conveyed as videos than still frames. Notwithstanding, the design and rendering of a complex molecular movie remain an arduous task. Here, we introduce Molywood, a robust and intuitive tool that builds on the capabilities of Visual Molecular Dynamics (VMD) to automate all stages of movie rendering. RESULTS: Molywood is a Python-based script that uses an integrated workflow to give maximal flexibility in movie design. It implements the basic concepts of actions, layers, grids and concurrency and requires no programming experience to run. AVAILABILITY AND IMPLEMENTATION: The script is freely available on GitLab (gitlab.com/KomBioMol/molywood) and PyPI (through pip), and features an extended documentation, tutorial and gallery hosted on mmb.irbbarcelona.org/molywood.

A shear stress micromodel of urinary tract infection by the Escherichia coli producing Dr adhesin.

In this study, we established a dynamic micromodel of urinary tract infection to analyze the impact of UT-segment-specific urinary outflow on the persistence of E. coli colonization. We found that the adherence of Dr+ E. coli to bladder T24 transitional cells and type IV collagen is maximal at lowest shear stress and is reduced by any increase in flow velocity. The analyzed adherence was effective in the whole spectrum of physiological shear stress and was almost irreversible over the entire range of generated shear force. Once Dr+ E. coli bound to host cells or collagen, they did not detach even in the presence of elevated shear stress or of chloramphenicol, a competitive inhibitor of binding. Investigating the role of epithelial surface architecture, we showed that the presence of budding cells-a model microarchitectural obstacle-promotes colonization of the urinary tract by E. coli. We report a previously undescribed phenomenon of epithelial cell "rolling-shedding" colonization, in which the detached epithelial cells reattach to the underlying cell line through a layer of adherent Dr+ E. coli. This rolling-shedding colonization progressed continuously due to "refilling" induced by the flow-perturbing obstacle. The shear stress of fluid containing free-floating bacteria fueled the rolling, while providing an uninterrupted supply of new bacteria to be trapped by the rolling cell. The progressive rolling allows for transfer of briefly attached bacteria onto the underlying monolayer in a repeating cascading event.