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.

Chemical Investigation on the Mechanism and Kinetics of the Atmospheric Degradation Reaction of Trichlorofluoroethene by OH⋠and Its Subsequent Fate in the Presence of O2 /NOx.

The M06-2X/6-311++G(d,p) level of theory was used to examine the degradation of Trichlorofluoroethene (TCFE) initiated by OH⋅ radicals. Additionally, the coupled-cluster single-double with triple perturbative [CCSD(T)] method was employed to refine the single-point energies using the complete basis set extrapolation approach. The results indicated that OH-addition is the dominant pathway. OH⋅ adds to both the C1 and C2 carbons, resulting in the formation of the C(OH)Cl2 -⋅CClF and ⋅CCl2 -C(OH)ClF species. The associated barrier heights were determined to be 1.11 and -0.99 kcal mol-1 , respectively. Furthermore, the energetic and thermodynamic parameters show that pathway 1 exhibits greater exothermicity and exergonicity compared to pathway 2, with differences of 8.11 and 8.21 kcal mol-1 , correspondingly. The primary pathway involves OH addition to the C2 position, with a rate constant of 6.2×10-13 cm3 molecule-1 sec-1 at 298 K. This analysis served to estimate the atmospheric lifetime, along with the photochemical ozone creation potential (POCP) and ozone depletion potential (ODP). It yielded an atmospheric lifetime of 8.49 days, an ODP of 4.8×10-4 , and a POCP value of 2.99, respectively. Radiative forcing efficiencies were also estimated at the M06-2X/6-311++G(d,p) level. Global warming potentials (GWPs) were calculated for 20, 100, and 500 years, resulting in values of 9.61, 2.61, and 0.74, respectively. TCFE is not expected to make a significant contribution to the radiative forcing of climate change. The results obtained from the time-dependent density functional theory (TDDFT) indicated that TCFE and its energized adducts are unable to photolysis under sunlight in the UV and visible spectrum. Secondary reactions involve the [TCFE-OH-O2 ]⋅ peroxy radical, leading subsequently to the [TCFE-OH-O]⋅ alkoxy radical. It was found that the alkoxy radical resulting from the peroxy radical can lead to the formation of phosgene (COCl2 ) and carbonyl chloride fluoride (CClFO), with phosgene being the primary product.

Screening Stability, Thermochemistry, and Chemical Kinetics of 3-Hydroxybutanoic Acid as a Bifunctional Biodiesel Additive.

The thermo-kinetic aspects of 3-hydroxybutyric acid (3-HBA) pyrolysis in the gas phase were investigated using density functional theory (DFT), specifically the M06-2X theoretical level in conjunction with the cc-pVTZ basis set. The obtained data were compared with benchmark CBS-QB3 results. The degradation mechanism was divided into 16 pathways, comprising 6 complex fissions and 10 barrierless reactions. Energy profiles were calculated and supplemented with computations of rate coefficients and branching ratios over the temperature range of 600-1700 K at a pressure of 1 bar using transition state theory (TST) and Rice-Ramsperger-Kassel-Marcus (RRKM) methods. Thermodynamics results indicated the presence of six stable conformers within a 4 kcal mol-1 energy range. The estimated chemical kinetics results suggested that TST and RRKM approaches are comparable, providing confidence in our calculations. The branching ratio analysis reveals that the dehydration reaction pathway leading to the formation of H2O and CH3CH═CHCO2H dominates entirely at T ≤ 650 K. At these temperatures, there is a minor contribution from the simple homolytic bond fission reaction, yielding related radicals [CH3•CHOH + •CH2CO2H]. However, at T ≥ 700 K, this reaction becomes the primary decomposition route. At T = 1700 K, there is a minor involvement of a reaction pathway resulting in the formation of CH3CH(OH)•CH2 + •CHO(OH) with an approximate contribution of 16%, and a reaction leading to [•CH3 + •CH2OHCH2CO2H] with around 9%.

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.

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.

Two-step mechanism of J-domain action in driving Hsp70 function.

J-domain proteins (JDPs), obligatory Hsp70 cochaperones, play critical roles in protein homeostasis. They promote key allosteric transitions that stabilize Hsp70 interaction with substrate polypeptides upon hydrolysis of its bound ATP. Although a recent crystal structure revealed the physical mode of interaction between a J-domain and an Hsp70, the structural and dynamic consequences of J-domain action once bound and how Hsp70s discriminate among its multiple JDP partners remain enigmatic. We combined free energy simulations, biochemical assays and evolutionary analyses to address these issues. Our results indicate that the invariant aspartate of the J-domain perturbs a conserved intramolecular Hsp70 network of contacts that crosses domains. This perturbation leads to destabilization of the domain-domain interface-thereby promoting the allosteric transition that triggers ATP hydrolysis. While this mechanistic step is driven by conserved residues, evolutionarily variable residues are key to initial JDP/Hsp70 recognition-via electrostatic interactions between oppositely charged surfaces. We speculate that these variable residues allow an Hsp70 to discriminate amongst JDP partners, as many of them have coevolved. Together, our data points to a two-step mode of J-domain action, a recognition stage followed by a mechanistic stage.

Why do G-quadruplexes dimerize through the 5'-ends? Driving forces for G4 DNA dimerization examined in atomic detail.

G-quadruplexes (G4) are secondary structures formed by guanine-rich nucleic acid sequences and shown to exist in living cells where they participate in regulation of gene expression and chromosome maintenance. G-quadruplexes with solvent-exposed guanine tetrads show the tendency to associate together through cofacial stacking, which may be important for packaging of G4-forming sequences and allows for the design of higher-order G4 DNA structures. To understand the molecular driving forces for G4 association, here, we study the binding interaction between two parallel-stranded G-quadruplexes using all-atom molecular dynamics simulations. The predicted dimerization free energies show that direct binding through the 5'-G-tetrads is the most preferred of all possible end-to-end stacking orientations, consistently with all available experimental data. Decomposition of dimerization enthalpies in combination with simulations at varying ionic strength further indicate that the observed orientational preferences arise from a fine balance between the electrostatic repulsion of the sugar-phosphate backbones and favorable counterion binding at the dimeric interface. We also demonstrate how these molecular-scale findings can be used to devise means of controlling G4 dimerization equilibrium, e.g., by altering salt concentration and using G4-targeted ligands.

The Product of Matrix Metalloproteinase Cleavage of Doxorubicin Conjugate for Anticancer Drug Delivery: Calorimetric, Spectroscopic, and Molecular Dynamics Studies on Peptide-Doxorubicin Binding to DNA.

Matrix metalloproteinases (MMPs) are extracellular matrix degradation factors, promoting cancer progression. Hence, they could provide an enzyme-assisted delivery of doxorubicin (DOX) in cancer treatment. In the current study, the intercalation process of DOX and tetrapeptide-DOX, the product of the MMPs' cleavage of carrier-linked DOX, into dsDNA was investigated using stationary and time-resolved fluorescence spectroscopy, UV-Vis spectrophotometry and isothermal titration calorimetry (ITC). The molecular dynamics (MD) simulations on the same tetrapeptide-DOX…DNA and DOX…DNA systems were also performed. The undertaken studies indicate that DOX and tetrapeptide-DOX can effectively bond with dsDNA through the intercalation mode; however, tetrapeptide-DOX forms less stable complexes than free DOX. Moreover, the obtained results demonstrate that the differences in DNA affinity of both forms of DOX can be attributed to different intercalation modes. Tetrapeptide-DOX shows a preference to intercalate into DNA through the major groove, whereas DOX does it through the minor one. In summary, we can conclude that the tetrapeptide-DOX intercalation to DNA is significant and that even the lack of non-specific proteases releasing DOX from the tetrapeptide conjugate, the presence of which is suggested by the literature for the efficient release of DOX, should not prevent the cytostatic action of the anthracycline.

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.

Effect of osmolytes on the thermal stability of proteins: replica exchange simulations of Trp-cage in urea and betaine solutions.

Although osmolytes are known to modulate the folding equilibrium, the molecular mechanism of their effect on thermal denaturation of proteins is still poorly understood. Here, we simulated the thermal denaturation of a small model protein (Trp-cage) in the presence of denaturing (urea) and stabilizing (betaine) osmolytes, using the all-atom replica exchange molecular dynamics simulations. We found that urea destabilizes Trp-cage by enthalpically-driven association with the protein, acting synergistically with temperature to induce unfolding. In contrast, betaine is sterically excluded from the protein surface thereby exerting entropic depletion forces that contribute to the stabilization of the native state. In fact, we find that while at low temperatures betaine slightly increases the folding free energy of Trp-cage by promoting another near-native conformation, it protects the protein against temperature-induced denaturation. This, in turn, can be attributed to enhanced exclusion of betaine at higher temperatures that arises from less attractive interactions with the protein surface.

Molecular dynamics simulations reveal the balance of forces governing the formation of a guanine tetrad-a common structural unit of G-quadruplex DNA.

G-quadruplexes (G4) are nucleic acid conformations of guanine-rich sequences, in which guanines are arranged in the square-planar G-tetrads, stacked on one another. G4 motifs formin vivoand are implicated in regulation of such processes as gene expression and chromosome maintenance. The structure and stability of various G4 topologies were determined experimentally; however, the driving forces for their formation are not fully understood at the molecular level. Here, we used all-atom molecular dynamics to probe the microscopic origin of the G4 motif stability. By computing the free energy profiles governing the dissociation of the 3'-terminal G-tetrad in the telomeric parallel-stranded G4, we examined the thermodynamic and kinetic stability of a single G-tetrad, as a common structural unit of G4 DNA. Our results indicate that the energetics of guanine association alone does not explain the overall stability of the G-tetrad and that interactions involving sugar-phosphate backbone, in particular, the constrained minimization of the phosphate-phosphate repulsion energy, are crucial in providing the observed enthalpic stabilization. This enthalpic gain is largely compensated by the unfavorable entropy change due to guanine association and optimization of the backbone topology.

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.

Iron-Sulfur Cluster Biogenesis Chaperones: Evidence for Emergence of Mutational Robustness of a Highly Specific Protein-Protein Interaction.

Biogenesis of iron-sulfur clusters (FeS) is a highly conserved process involving Hsp70 and J-protein chaperones. However, Hsp70 specialization differs among species. In most eukaryotes, including Schizosaccharomyces pombe, FeS biogenesis involves interaction between the J-protein Jac1 and the multifunctional Hsp70 Ssc1. But, in Saccharomyces cerevisiae and closely related species, Jac1 interacts with the specialized Hsp70 Ssq1, which emerged through duplication of SSC1. As little is known about how gene duplicates affect the robustness of their protein interaction partners, we analyzed the functional and evolutionary consequences of Ssq1 specialization on the ubiquitous J-protein cochaperone Jac1, by comparing S. cerevisiae and S. pombe. Although deletion of JAC1 is lethal in both species, alanine substitutions within the conserved His-Pro-Asp (HPD) motif, which is critical for Jac1:Hsp70 interaction, have species-specific effects. They are lethal in S. pombe, but not in S. cerevisiae. These in vivo differences correlated with in vitro biochemical measurements. Charged residues present in the J-domain of S. cerevisiae Jac1, but absent in S. pombe Jac1, are important for tolerance of S. cerevisiae Jac1 to HPD alterations. Moreover, Jac1 orthologs from species that encode Ssq1 have a higher sequence divergence. The simplest interpretation of our results is that Ssq1's coevolution with Jac1 resulted in expansion of their binding interface, thus increasing the efficiency of their interaction. Such an expansion could in turn compensate for negative effects of HPD substitutions. Thus, our results support the idea that the robustness of Jac1 emerged as consequence of its highly efficient and specific interaction with Ssq1.

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.