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.

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.

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.].

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.

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.

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.

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.

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.

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.

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.

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.

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.

How proteins bind to DNA: target discrimination and dynamic sequence search by the telomeric protein TRF1.

Target search as performed by DNA-binding proteins is a complex process, in which multiple factors contribute to both thermodynamic discrimination of the target sequence from overwhelmingly abundant off-target sites and kinetic acceleration of dynamic sequence interrogation. TRF1, the protein that binds to telomeric tandem repeats, faces an intriguing variant of the search problem where target sites are clustered within short fragments of chromosomal DNA. In this study, we use extensive (>0.5 ms in total) MD simulations to study the dynamical aspects of sequence-specific binding of TRF1 at both telomeric and non-cognate DNA. For the first time, we describe the spontaneous formation of a sequence-specific native protein-DNA complex in atomistic detail, and study the mechanism by which proteins avoid off-target binding while retaining high affinity for target sites. Our calculated free energy landscapes reproduce the thermodynamics of sequence-specific binding, while statistical approaches allow for a comprehensive description of intermediate stages of complex formation.

Role of the disulfide bond in stabilizing and folding of the fimbrial protein DraE from uropathogenic Escherichia coli.

Dr fimbriae are homopolymeric adhesive organelles of uropathogenic Escherichia coli composed of DraE subunits, responsible for the attachment to host cells. These structures are characterized by enormously high stability resulting from the structural properties of an Ig-like fold of DraE. One feature of DraE and other fimbrial subunits that makes them peculiar among Ig-like domain-containing proteins is a conserved disulfide bond that joins their A and B strands. Here, we investigated how this disulfide bond affects the stability and folding/unfolding pathway of DraE. We found that the disulfide bond stabilizes self-complemented DraE (DraE-sc) by ∼50 kJ mol-1 in an exclusively thermodynamic manner, i.e. by lowering the free energy of the native state and with almost no effect on the free energy of the transition state. This finding was confirmed by experimentally determined folding and unfolding rate constants of DraE-sc and a disulfide bond-lacking DraE-sc variant. Although the folding of both proteins exhibited similar kinetics, the unfolding rate constant changed upon deletion of the disulfide bond by 10 orders of magnitude, from ∼10-17 s-1 to 10-7 s-1 Molecular simulations revealed that unfolding of the disulfide bond-lacking variant is initiated by strands A or G and that disulfide bond-mediated joining of strand A to the core strand B cooperatively stabilizes the whole protein. We also show that the disulfide bond in DraE is recognized by the DraB chaperone, indicating a mechanism that precludes the incorporation of less stable, non-oxidized DraE forms into the fimbriae.

Mechanism of Binding of Antifungal Antibiotic Amphotericin B to Lipid Membranes: An Insight from Combined Single-Membrane Imaging, Microspectroscopy, and Molecular Dynamics.

Amphotericin B is a lifesaving polyene antibiotic used in the treatment of systemic mycoses. Unfortunately, the pharmacological applicability of this drug is limited because of its severe toxic side effects. At the same time, the lack of a well-defined mechanism of selectivity hampers the efforts to rationally design safer derivatives. As the drug primarily targets the biomembranes of both fungi and humans, new insights into the binding of amphotericin B to lipid membranes can be helpful in unveiling the molecular mechanisms underlying both its pharmacological activity and toxicity. We use fluorescence-lifetime-imaging microscopy combined with fluorescence-emission spectroscopy in the microscale to study the interaction of amphotericin B with single lipid bilayers, using model systems based on giant unilamellar liposomes formed with three lipids: dipalmitoylphosphatidylcholine (DPPC), dimirystoylphosphatidylcholine (DMPC), and 1-palmitoyl-2-oleoylphosphatidylcholine (POPC). The results show that amphotericin B introduced into the water phase as a DMSO solution binds to the membrane as dimers and small-molecular aggregates that we identify as tetramers and trimers. Fluorescence-detected linear-dichroism measurements revealed high orientational freedom of all the molecular-organization forms with respect to the membrane plane, which suggests that the drug partially binds to the membrane surface. The presence of sterols in the lipid phase (cholesterol but particularly ergosterol at 30 mol %) promotes the penetration of drug molecules into the lipid membrane, as concluded on the basis of the decreased orientation angle of amphotericin B molecules with respect to the axis normal to the membrane plane. Moreover, ergosterol facilitates the association of amphotericin B dimers into aggregated structures that can play a role in membrane destabilization or permeabilization. The presence of cholesterol inhibits the formation of small aggregates in the lipid phase of liposomes, making this system a promising candidate for a low-toxicity antibiotic-delivery system. Our conclusions are supported with molecular simulations that reveal the conformational properties of AmB oligomers in both aqueous solution and lipid bilayers of different compositions.

Mechanochemical Energy Transduction during the Main Rotary Step in the Synthesis Cycle of F1-ATPase.

F1-ATPase is a highly efficient molecular motor that can synthesize ATP driven by a mechanical torque. Its ability to function reversibly in either direction requires tight mechanochemical coupling between the catalytic domain and the rotating central shaft, as well as temporal control of substrate binding and product release. Despite great efforts and significant progress, the molecular details of this synchronized and fine-tuned energy conversion mechanism are not fully understood. Here, we use extensive molecular dynamics simulations to reconcile recent single-molecule experiments with structural data and provide a consistent thermodynamic, kinetic and mechanistic description of the main rotary substep in the synthetic cycle of mammalian ATP synthase. The calculated free energy profiles capture a discrete pattern in the rotation of the central γ-shaft, with a metastable intermediate located-consistently with recent experimental findings-at 70° relative to the X-ray position. We identify this rotary step as the ATP-dependent substep, and find that the associated free energy input supports the mechanism involving concurrent nucleotide binding and release. During the main substep, our simulations show no significant opening of the ATP-bound ß subunit; instead, we observe that mechanical energy is transmitted to its nucleotide binding site, thus lowering the affinity for ATP. Simultaneously, the empty subunit assumes a conformation that enables the enzyme to harness the free energy of ADP binding to drive ATP release. Finally, we show that ligand exchange is regulated by a checkpoint mechanism, an apparent prerequisite for high efficiency in protein nanomotors.

Molecular basis of the osmolyte effect on protein stability: a lesson from the mechanical unfolding of lysozyme.

Osmolytes are a class of small organic molecules that shift the protein folding equilibrium. For this reason, they are accumulated by organisms under environmental stress and find applications in biotechnology where proteins need to be stabilized or dissolved. However, despite years of research, debate continues over the exact mechanisms underpinning the stabilizing and denaturing effect of osmolytes. Here, we simulated the mechanical denaturation of lysozyme in different solvent conditions to study the molecular mechanism by which two biologically relevant osmolytes, denaturing (urea) and stabilizing (betaine), affect the folding equilibrium. We found that urea interacts favorably with all types of residues via both hydrogen bonds and dispersion forces, and therefore accumulates in a diffuse solvation shell around the protein. This not only provides an enthalpic stabilization of the unfolded state, but also weakens the hydrophobic effect, as hydrophobic forces promote the association of urea with nonpolar residues, facilitating the unfolding. In contrast, we observed that betaine is excluded from the protein backbone and nonpolar side chains, but is accumulated near the basic residues, yielding a nonuniform distribution of betaine molecules at the protein surface. Spatially resolved solvent-protein interaction energies further suggested that betaine behaves in a ligand- rather than solvent-like manner and its exclusion from the protein surface arises mostly from the scarcity of favorable binding sites. Finally, we found that, in the presence of betaine, the reduced ability of water molecules to solvate the protein results in an additional enthalpic contribution to the betaine-induced stabilization.

Membrane Sterols Modulate the Binding Mode of Amphotericin B without Affecting Its Affinity for a Lipid Bilayer.

Membrane-active antibiotics are known to selectively target certain pathogens based on cell membrane properties, such as fluidity, lipid ordering, and phase behavior. These are in turn modulated by the composition of a lipid bilayer and in particular by the presence and type of membrane sterols. Amphotericin B (AmB), the golden standard of antifungal treatment, exhibits higher activity toward ergosterol-rich fungal membranes, which permits its use against systemic mycoses; however, the selectivity for fungal membranes is far from satisfactory leading to severe side effects. Despite decades of research, no consensus has emerged on the origin of AmB specificity for fungal cells and its actual mode of action at the molecular level. Previously, it has been proposed that the specific action of AmB is related to differences in its affinity for membranes of different composition. In this work, we investigate this relationship by employing molecular dynamics simulations to compare the free energy of insertion of AmB into three types of membranes: a pure DMPC bilayer and DMPC bilayers containing 30% of cholesterol or ergosterol. We analyze the orientation of AmB molecules within the bilayer in order to unambiguously establish their membrane binding mode and relate the orientational freedom to the sterol-dependent tightness of lipid packing. Our results strongly indicate that the membrane insertion of AmB proceeds virtually to completion independent of membrane type, and hence the higher toxicity against fungal membranes may rather result from differences in subsequent oligomerization in the membrane and assembly of monomers into functional transmembrane pores. In particular, the latter could be facilitated by sterol-induced ordering of AmB molecules along the membrane normal, revealed by our free energy profiles. Moreover--in contrast to certain claims--we find no stable binding mode corresponding to the horizontal adsorption of AmB on the membrane surface.

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.