Miłosz Wieczór

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.

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

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

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

Thermodynamics and kinetics of amphotericin B self-association in aqueous solution characterized in molecular detail.

Amphotericin B (AmB) is a potent but toxic drug commonly used to treat systemic mycoses. Its efficiency as a therapeutic agent depends on its ability to discriminate between mammalian and fungal cell membranes. The association of AmB monomers in an aqueous environment plays an important role in drug selectivity, as oligomers formed prior to membrane insertion – presumably dimers – are believed to act differently on fungal (ergosterol-rich) and mammalian (cholesterol-rich) membranes. In this work, we investigate the initial steps of AmB self-association by studying the structural, thermodynamic and spectral properties of AmB dimers in aqueous medium using molecular dynamics simulations. Our results show that in water, the hydrophobic aggregation of AmB monomers yields almost equiprobable populations of parallel and antiparallel dimers that rapidly interconvert into each other, and the dipole-dipole interaction between zwitterionic head groups plays a minor role in determining the drug’s tendency for self-aggregation. A simulation of circular dichroism (CD) spectra indicates that in experimental measurements, the signature CD spectrum of AmB aggregates should be attributed to higher-order oligomers rather than dimers. Finally, we suggest that oligomerization can impair the selectivity of AmB molecules for fungal membranes by increasing their hydrophobic drive for non-specific membrane insertion.

Correlation between the number of Pro-Ala repeats in the EmrA homologue of Acinetobacter baumannii and resistance to netilmicin, tobramycin, imipenem and ceftazidime.

Acinetobacter baumannii coccobacilli are dangerous to patients in intensive care units because of their multidrug resistance to antibiotics, developed mainly in the past decade. This study aimed to examine whether there is a significant correlation between the number of Pro-Ala repeats in the CAP01997 protein, the EmrA homologue of A. baumannii, and resistance to antibiotics. A total of 79 multidrug-resistant A. baumannii strains isolated from patients were analysed. Resistance to antibiotics was determined on Mueller-Hinton agar plates using the Kirby-Bauer disk diffusion method. The number of CCTGCA repeats encoding Pro-Ala repeats in CAP01997 was determined by PCR and capillary electrophoresis. The 3D models of CAP01997 containing Pro-Ala repeats were initially generated using RaptorX Structure Prediction server and were assembled with EasyModeller 4.0. The models were embedded in a model bacterial membrane based on structural information from homologous proteins and were refined using 100-ns molecular dynamics simulations. The results of this research show significant correlation between susceptibility to netilmicin, tobramycin and imipenem and the number of repeated Pro-Ala sequences in the CAP01997 protein, a homologue of the Escherichia coli transporter EmrA. Predicted structures suggest potential mechanisms that confer drug resistance by reshaping the cytoplasmic interface between CAP01997 protein and the critical component of the multidrug efflux pump homologous to EmrB. Based on these results, we can conclude that the CAP01997 protein, an EmrA homologue of A. baumannii, confers resistance to netilmicin, tobramycin and imipenem, depending on the number of Pro-Ala repeats.

Molecular recognition in complexes of TRF proteins with telomeric DNA.

Telomeres are specialized nucleoprotein assemblies that protect the ends of linear chromosomes. In humans and many other species, telomeres consist of tandem TTAGGG repeats bound by a protein complex known as shelterin that remodels telomeric DNA into a protective loop structure and regulates telomere homeostasis. Shelterin recognizes telomeric repeats through its two major components known as Telomere Repeat-Binding Factors, TRF1 and TRF2. These two homologous proteins are therefore essential for the formation and normal function of telomeres. Indeed, TRF1 and TRF2 are implicated in a plethora of different cellular functions and their depletion leads to telomere dysfunction with chromosomal fusions, followed by apoptotic cell death. More specifically, it was found that TRF1 acts as a negative regulator of telomere length, and TRF2 is involved in stabilizing the loop structure. Consequently, these proteins are of great interest, not only because of their key role in telomere maintenance and stability, but also as potential drug targets. In the current study, we investigated the molecular basis of telomeric sequence recognition by TRF1 and TRF2 and their DNA binding mechanism. We used molecular dynamics (MD) to calculate the free energy profiles for binding of TRFs to telomeric DNA. We found that the predicted binding free energies were in good agreement with experimental data. Further, different molecular determinants of binding, such as binding enthalpies and entropies, the hydrogen bonding pattern and changes in surface area, were analyzed to decompose and examine the overall binding free energies at the structural level. With this approach, we were able to draw conclusions regarding the consecutive stages of sequence-specific association, and propose a novel aspartate-dependent mechanism of sequence recognition. Finally, our work demonstrates the applicability of computational MD-based methods to studying protein-DNA interactions.

N-(4-Methyl­piperazin-4-ium-1-yl)dithio­carbamate sesquihydrate

In the crystal structure of the title compound, C6H13N3S2·1.5H2O, weak N—H⋯S interactions between the zwitterionic molecules are observed, leading to an extensively folded layered arrangement parallel to (100). There are three crystallographically independent water molecules in the asymmetric unit, which are disordered and only half occupied.

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.

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.