Search

Vortioxetine (VTX) is a recently approved antidepressant that targets a variety of serotonin receptors. Here, we investigate the drug's molecular mechanism of operation at the serotonin 5-HT 3 receptor (5-HT 3 R), which features two properties: VTX acts differently on rodent and human 5-HT 3 R, and VTX appears to suppress any subsequent response to agonists. Using a combination of cryo-EM, electrophysiology, voltage-clamp fluorometry and molecular dynamics, we show that VTX stabilizes a resting inhibited state of the mouse 5-HT 3 R and an agonist-bound-like state of human 5-HT 3 R, in line with the functional profile of the drug. We report four human 5-HT 3 R structures and show that the human receptor transmembrane domain is intrinsically fragile. We also explain the lack of recovery after VTX administration via a membrane partition mechanism., (© 2024. The Author(s), under exclusive licence to Springer Nature America, Inc.)

All-atom simulations are a powerful approach to study the structure and dynamics of biological membranes. However, sampling the atomic configurations of inhomogeneous membranes can be challenging due to the slow lateral diffusion of the various constituents. To address this issue, a hybrid nonequilibrium molecular dynamics Monte Carlo (neMD/MC) simulation method is proposed in which randomly chosen lipid molecules are swapped to generate configurations that are subsequently accepted or rejected according to a Metropolis criterion based on the alchemical work for the attempted swap calculated via a short trajectory. A dual-topology framework constraining the common atoms of the exchanging molecules yields a good acceptance probability using switching trajectories as short as 10 ps. The performance of the hybrid neMD/MC algorithm and its ability to sample the distribution of lipids near a transmembrane helix carrying a net charge are illustrated for a binary mixture of charged and zwitterionic lipids.

The ability of the Torpedo nicotinic acetylcholine receptor (nAChR) to undergo agonist-induced conformational transitions requires the presence of cholesterol and/or anionic lipids. Here we use recently solved structures along with multiscale molecular dynamics simulations to examine lipid binding to the nAChR in bilayers that have defined effects on nAChR function. We examine how phosphatidic acid and cholesterol, lipids that support conformational transitions, individually compete for binding with phosphatidylcholine, a lipid that does not. We also examine how the two lipids work synergistically to stabilize an agonist-responsive nAChR. We identify rapidly exchanging lipid binding sites, including both phospholipid sites with a high affinity for phosphatidic acid and promiscuous cholesterol binding sites in the grooves between adjacent transmembrane α-helices. A high affinity cholesterol site is confirmed in the inner leaflet framed by a key tryptophan residue on the MX α-helix. Our data provide insight into the dynamic nature of lipid-nAChR interactions and set the stage for a detailed understanding of the mechanisms by which lipids facilitate nAChR function at the neuromuscular junction., (© 2024. The Author(s).)

Glutathione transferases (GSTs) constitute a widespread superfamily of enzymes notably involved in detoxification processes and/or in specialized metabolism. In the cyanobacterium Synechocsytis sp. PCC 6803, SynGSTC1, a chi-class GST (GSTC), is thought to participate in the detoxification process of methylglyoxal, a toxic by-product of cellular metabolism. A comparative genomic analysis showed that GSTCs were present in all orders of cyanobacteria with the exception of the basal order Gloeobacterales. These enzymes were also detected in some marine and freshwater noncyanobacterial bacteria, probably as a result of horizontal gene transfer events. GSTCs were shorter of about 30 residues compared to most cytosolic GSTs and had a well-conserved SRAS motif in the active site ( 10 SRAS 13 in SynGSTC1). The crystal structure of SynGSTC1 in complex with glutathione adopted the canonical GST fold with a very open active site because the α4 and α5 helices were exceptionally short. A transferred multipolar electron-density analysis allowed a fine description of the solved structure. Unexpectedly, Ser10 did not have an electrostatic influence on glutathione as usually observed in serinyl-GSTs. The S10A variant was only slightly less efficient than the wild-type and molecular dynamics simulations suggested that S10 was a stabilizer of the protein backbone rather than an anchor site for glutathione.

Atomic-level information is essential to explain the formation of specific protein complexes in terms of structure and dynamics. The set of Dpr and DIP proteins, which play a key role in the neuromorphogenesis in the nervous system of Drosophila melanogaster , offer a rich paradigm to learn about protein-protein recognition. Many members of the DIP subfamily cross-react with several members of the Dpr family and vice versa. While there exists a total of 231 possible Dpr-DIP heterodimer complexes from the 21 Dpr and 11 DIP proteins, only 57 "cognate" pairs have been detected by surface plasmon resonance (SPR) experiments, suggesting that the remaining 174 pairs have low or unreliable binding affinity. Our goal is to assess the performance of computational approaches to characterize the global set of interactions between Dpr and DIP proteins and identify the specificity of binding between each DIP with their corresponding Dpr binding partners. In addition, we aim to characterize how mutations influence the specificity of the binding interaction. In this work, a wide range of knowledge-based and physics-based approaches are utilized, including mutual information, linear discriminant analysis, homology modeling, molecular dynamics simulations, Poisson-Boltzmann continuum electrostatics calculations, and alchemical free energy perturbation to decipher the origin of binding specificity of the Dpr-DIP complexes examined. Ultimately, the results show that those two broad strategies are complementary, with different strengths and limitations. Biological inter-relations are more clearly revealed through knowledge-based approaches combining evolutionary and structural features, the molecular determinants controlling binding specificity can be predicted accurately with physics-based approaches based on atomic models.

Perforin-2 (PFN2, MPEG1) is a pore-forming protein that acts as a first line of defense in the mammalian immune system, rapidly killing engulfed microbes within the phagolysosome in macrophages. PFN2 self-assembles into hexadecameric pre-pore rings that transition upon acidification into pores damaging target cell membranes. Here, using high-speed atomic force microscopy (HS-AFM) imaging and line-scanning and molecular dynamics simulation, we elucidate PFN2 pre-pore to pore transition pathways and dynamics. Upon acidification, the pre-pore rings (pre-pore-I) display frequent, 1.8 s -1 , ring-opening dynamics that eventually, 0.2 s -1 , initiate transition into an intermediate, short-lived, ~75 ms, pre-pore-II state, inducing a clockwise pre-pore-I to pre-pore-II propagation. Concomitantly, the first pre-pore-II subunit, undergoes a major conformational change to the pore state that propagates also clockwise at a rate ~15 s -1 . Thus, the pre-pore to pore transition is a clockwise hand-over-hand mechanism that is accomplished within ~1.3 s. Our findings suggest a clockwise mechanism of membrane insertion that with variations may be general for the MACPF/CDC superfamily., (© 2022. The Author(s).)

Deinococcus radiodurans is a spherical bacterium well-known for its outstanding resistance to DNA-damaging agents. Exposure to such agents leads to drastic changes in the transcriptome of D. radiodurans. In particular, four Deinococcus-specific genes, known as DNA Damage Response genes, are strongly up-regulated and have been shown to contribute to the resistance phenotype of D. radiodurans. One of these, DdrC, is expressed shortly after exposure to γ-radiation and is rapidly recruited to the nucleoid. In vitro, DdrC has been shown to compact circular DNA, circularize linear DNA, anneal complementary DNA strands and protect DNA from nucleases. To shed light on the possible functions of DdrC in D. radiodurans, we determined the crystal structure of the domain-swapped DdrC dimer at a resolution of 2.5 Å and further characterized its DNA binding and compaction properties. Notably, we show that DdrC bears two asymmetric DNA binding sites located on either side of the dimer and can modulate the topology and level of compaction of circular DNA. These findings suggest that DdrC may be a DNA damage-induced nucleoid-associated protein that enhances nucleoid compaction to limit the dispersion of the fragmented genome and facilitate DNA repair after exposure to severe DNA damaging conditions., (© The Author(s) 2022. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Fast synaptic communication requires receptors that respond to the presence of neurotransmitter by opening an ion channel across the post-synaptic membrane. The muscle-type nicotinic acetylcholine receptor from the electric fish, Torpedo, is the prototypic ligand-gated ion channel, yet the structural changes underlying channel activation remain undefined. Here we use cryo-EM to solve apo and agonist-bound structures of the Torpedo nicotinic receptor embedded in a lipid nanodisc. Using both a direct biochemical assay to define the conformational landscape and molecular dynamics simulations to assay flux through the pore, we correlate structures with functional states and elucidate the motions that lead to pore activation of a heteromeric nicotinic receptor. We highlight an underappreciated role for the complementary subunit in channel gating, establish the structural basis for the differential agonist affinities of α/δ versus α /γ sites, and explain why nicotine is less potent at muscle nicotinic receptors compared to neuronal ones., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 Elsevier Inc. All rights reserved.)

Designing a reliable computational methodology to calculate protein:ligand standard binding free energies is extremely challenging. The large change in configurational enthalpy and entropy that accompanies the association of ligand and protein is notoriously difficult to capture in naive brute-force simulations. Addressing this issue, the present protocol rests upon a rigorous statistical mechanical framework for the determination of protein:ligand binding affinities together with the comprehensive Binding Free-Energy Estimator 2 (BFEE2) application software. With the knowledge of the bound state, available from experiments or docking, application of the BFEE2 protocol with a reliable force field supplies in a matter of days standard binding free energies within chemical accuracy, for a broad range of protein:ligand complexes. Limiting undesirable human intervention, BFEE2 assists the end user in preparing all the necessary input files and performing the post-treatment of the simulations towards the final estimate of the binding affinity., (© 2022. The Author(s), under exclusive licence to Springer Nature Limited.)

The mitochondrial ADP/ATP carrier (AAC) performs the first and last step in oxidative phosphorylation by exchanging ADP and ATP across the mitochondrial inner membrane. Its optimal function has been shown to be dependent on cardiolipins (CLs), unique phospholipids located almost exclusively in the mitochondrial membrane. In addition, AAC exhibits an enthralling threefold pseudosymmetry, a unique feature of members of the SLC25 family. Recently, its conformation poised for binding of ATP was solved by x-ray crystallography referred to as the matrix state. Binding of the substrate leads to conformational changes that export of ATP to the mitochondrial intermembrane space. In this contribution, we investigate the influence of CLs on the structure, substrate-binding properties, and structural symmetry of the matrix state, employing microsecond-scale molecular dynamics simulations. Our findings demonstrate that CLs play a minor stabilizing role on the AAC structure. The interdomain salt bridges and hydrogen bonds forming the cytoplasmic network and tyrosine braces, which ensure the integrity of the global AAC scaffold, highly benefit from the presence of CLs. Under these conditions, the carrier is found to be organized in a more compact structure in its interior, as revealed by analyses of the electrostatic potential, measure of the AAC cavity aperture, and the substrate-binding assays. Introducing a convenient structure-based symmetry metric, we quantified the structural threefold pseudosymmetry of AAC, not only for the crystallographic structure, but also for conformational states of the carrier explored in the molecular dynamics simulations. Our results suggest that CLs moderately contribute to preserve the pseudosymmetric structure of AAC., (Copyright © 2021 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Inaccurately perceived as niche drugs, antiemetics are key elements of cancer treatment alleviating the most dreaded side effect of chemotherapy. Serotonin 5-HT3 receptor antagonists are the most commonly prescribed class of drugs to control chemotherapy-induced nausea and vomiting. These antagonists have been clinically successful drugs since the 1980s, yet our understanding of how they operate at the molecular level has been hampered by the difficulty of obtaining structures of drug-receptor complexes. Here, we report the cryoelectron microscopy structure of the palonosetron-bound 5-HT3 receptor. We investigate the binding of palonosetron, granisetron, dolasetron, ondansetron, and cilansetron using molecular dynamics, covering the whole set of antagonists used in clinical practice. The structural and computational results yield detailed atomic insight into the binding modes of the drugs. In light of our data, we establish a comprehensive framework underlying the inhibition mechanism by the -setron drug family., Competing Interests: Declaration of Interests The authors declare no competing interests., (Copyright © 2020 Elsevier Ltd. All rights reserved.)

The mitochondrial respiratory chain, formed by five protein complexes, utilizes energy from catabolic processes to synthesize ATP. Complex I, the first and the largest protein complex of the chain, harvests electrons from NADH to reduce quinone, while pumping protons across the mitochondrial membrane. Detailed knowledge of the working principle of such coupled charge-transfer processes remains, however, fragmentary due to bottlenecks in understanding redox-driven conformational transitions and their interplay with the hydrated proton pathways. Complex I from Thermus thermophilus encases 16 subunits with nine iron-sulfur clusters, reduced by electrons from NADH. Here, employing the latest crystal structure of T. thermophilus complex I, we have used microsecond-scale molecular dynamics simulations to study the chemo-mechanical coupling between redox changes of the iron-sulfur clusters and conformational transitions across complex I. First, we identify the redox switches within complex I, which allosterically couple the dynamics of the quinone binding pocket to the site of NADH reduction. Second, our free-energy calculations reveal that the affinity of the quinone, specifically menaquinone, for the binding-site is higher than that of its reduced, menaquinol form-a design essential for menaquinol release. Remarkably, the barriers to diffusive menaquinone dynamics are lesser than that of the more ubiquitous ubiquinone, and the naphthoquinone headgroup of the former furnishes stronger binding interactions with the pocket, favoring menaquinone for charge transport in T. thermophilus . Our computations are consistent with experimentally validated mutations and hierarchize the key residues into three functional classes, identifying new mutation targets. Third, long-range hydrogen-bond networks connecting the quinone-binding site to the transmembrane subunits are found to be responsible for proton pumping. Put together, the simulations reveal the molecular design principles linking redox reactions to quinone turnover to proton translocation in complex I.

Promoting drug delivery across the biological membrane is a common strategy to improve bioavailability. Inspired by the observation that carbonated alcoholic beverages can increase the absorption rate of ethanol, we speculate that carbon dioxide (CO 2 ) molecules could also enhance membrane permeability to drugs. In the present work, we have investigated the effect of CO 2 on the permeability of a model membrane formed by 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine lipids to three drug-like molecules, namely, ethanol, 2',3'-dideoxyadenosine, and trimethoprim. The free-energy and fractional-diffusivity profiles underlying membrane translocation were obtained from μs-timescale simulations and combined in the framework of the fractional solubility-diffusion model. We find that addition of CO 2 in the lipid environment results in an increase of the membrane permeability to the three substrates. Further analysis of the permeation events reveals that CO 2 expands and loosens the membrane, which, in turn, facilitates permeation of the drug-like molecules. © 2019 Wiley Periodicals, Inc., (© 2019 Wiley Periodicals, Inc.)

DNA compaction is essential to ensure the packaging of the genetic material in living cells and also plays a key role in the epigenetic regulation of gene expression. In both humans and bacteria, DNA packaging is achieved by specific well-conserved proteins. Here, by means of all-atom molecular dynamics simulations, including the determination of relevant free-energy profiles, we rationalize the molecular bases for this remarkable process in bacteria, illustrating the crucial role played by positively charged amino acids of a small histone-like protein. We also present compelling evidence that this histone-like protein alone can induce strong bending of a DNA duplex around its core domain, a process that requires overcoming a major free-energy barrier.

In humans, vision is limited to a small fraction of the whole electromagnetic spectrum. One possible strategy for enhancing vision in deep-red or poor-light conditions consists of recruiting chlorophyll derivatives in the rod photoreceptor cells of the eye, as suggested in the case of some deep-sea fish. Here, we employ all-atom molecular simulations and high-level quantum chemistry calculations to rationalize how chlorin e6 (Ce6), widely used in photodynamic therapy although accompanied by enhanced visual sensitivity, mediates vision in the dark, shining light on a fascinating but largely unknown molecular mechanism. First, we identify persistent interaction sites between Ce6 and the extracellular loops of rhodopsin, the transmembrane photoreceptor protein responsible for the first steps in vision. Triggered by Ce6 deep-red light absorption, the retinal within rhodopsin can be isomerized thus starting the visual phototransduction cascade. Our data largely exclude previously hypothesized energy-transfer mechanisms while clearly lending credence to a retinal isomerization indirectly triggered by singlet oxygen, proposing an alternative mechanism to rationalize photosensitizer-mediated night vision.

Potent and broadly neutralizing antibodies (bnAbs) are the hallmark of HIV-1 protection by vaccination. The membrane-proximal external region (MPER) of the HIV-1 gp41 fusion protein is targeted by the most broadly reactive HIV-1 neutralizing antibodies. Here, we examine the structural and molecular mechansims of neutralization by anti-MPER bnAb, LN01, which was isolated from lymph-node-derived germinal center B cells of an elite controller and exhibits broad neutralization breadth. LN01 engages both MPER and the transmembrane (TM) region, which together form a continuous helix in complex with LN01. The tilted TM orientation allows LN01 to interact simultaneously with the peptidic component of the MPER epitope and membrane via two specific lipid binding sites of the antibody paratope. Although LN01 carries a high load of somatic mutations, most key residues interacting with the MPER epitope and lipids are germline encoded, lending support for the LN01 epitope as a candidate for lineage-based vaccine development., (Copyright © 2019 The Authors. Published by Elsevier Inc. All rights reserved.)

The DNA-binding of the natural benzophenanthridine alkaloid chelerythrine (CHE) has been assessed by combining molecular modeling and optical absorption spectroscopy. Specifically, both double-helical (B-DNA) and G-quadruplex sequences-representative of different topologies and possessing biological relevance, such as telomeric or regulatory sequences-have been considered. An original multiscale protocol, making use of molecular dynamics (MD) simulations and quantum mechanics/molecular mechanics (QM/MM) calculations, allowed us to compare the theoretical and experimental circular dichroism spectra of the different DNA topologies, readily providing atomic-level details of the CHE-DNA binding modes. The binding selectivity towards G-quadruplexes is confirmed by both experimental and theoretical determination of the binding free energies. Overall, our mixed computational and experimental approach is able to shed light on the interaction of small molecules with different DNA conformations. In particular, CHE may be seen as the building block of promising drug candidates specifically targeting G-quadruplexes for both antitumoral and antiviral purposes.

The design of conjugates displaying simultaneously high selectivity and high affinity for different subtypes of integrins is a current challenge. The arginine-glycine-aspartic acid amino acid sequence (RGD) is one of the most efficient short peptides targeting these receptors. We report herein the development of linear and cyclic fluoro-C-glycoside"RGD" conjugates, taking advantage of the robustness and hydrophilicity of C-glycosides. As attested by in vitro evaluation, the design of these C-glyco"RGD" with a flexible three-carbon triazolyl linker allows distinct profiles towards α IIb β 3 and α v β 3 integrins. Molecular-dynamics simulations confirm the suitability of cyclic C-glyco-c(RGDfC) to target α v β 3 integrin. These C-glyco"RGD" could become promising biological tools in particular for Positron Emission Tomography imaging., (Copyright © 2019 Elsevier Ltd. All rights reserved.)

Coordinated conformational transitions in oligomeric enzymatic complexes modulate function in response to substrates and play a crucial role in enzyme inhibition and activation. Caseinolytic protease (ClpP) is a tetradecameric complex, which has emerged as a drug target against multiple pathogenic bacteria. Activation of different ClpPs by inhibitors has been independently reported from drug development efforts, but no rationale for inhibitor-induced activation has been hitherto proposed. Using an integrated approach that includes x-ray crystallography, solid- and solution-state nuclear magnetic resonance, molecular dynamics simulations, and isothermal titration calorimetry, we show that the proteasome inhibitor bortezomib binds to the ClpP active-site serine, mimicking a peptide substrate, and induces a concerted allosteric activation of the complex. The bortezomib-activated conformation also exhibits a higher affinity for its cognate unfoldase ClpX. We propose a universal allosteric mechanism, where substrate binding to a single subunit locks ClpP into an active conformation optimized for chaperone association and protein processive degradation.

Anabaena Sensory Rhodopsin (ASR), a microbial photoactive protein featuring the retinal chromophore in two different conformations, exhibits a pH-dependent electronic absorption spectrum. Using the recently developed CpHMD-then-QM/MM multiscale protocol applied to ASR embedded in a membrane model, the pH-induced changes in its maximum absorption wavelength have been reproduced and analyzed. While the acidic tiny red-shift is essentially correlated with the deprotonation of an aspartic acid located on the ASR extracellular side, the larger blue-shift experimentally reported at pH values larger than 5 involves a cluster of titrating residues sitting on the cytoplasmic side. The ASR pH-dependent spectrum is the consequence of the competitive stabilization of retinal ground and excited states by the protein electrostatic potential.

In cells, many constituents are able to assemble resulting in large macromolecular machineries possessing very specific biological and physiological functions, e.g. ribosome, spliceosome and proteasome. Assembly of such entities is commonly mediated by transient protein factors. SPAG1 is a multidomain protein, known to participate in the assembly of both the inner and outer dynein arms. These arms are required for the function of sensitive and motile cells. Together with RUVBL1, RUVBL2 and PIH1D2, SPAG1 is a key element of R2SP, a protein complex assisting the quaternary assembly of specific protein clients in a tissue-specific manner and associating with heat shock proteins (HSPs) and regulators. In this study, we have investigated the role of TPR domains of SPAG1 in the recruitment of HSP chaperones by combining biochemical assays, ITC, NMR spectroscopy and molecular dynamics (MD) simulations. First, we propose that only two, out of the three TPR domains, are able to recruit the protein chaperones HSP70 and HSP90. We then focused on one of these TPR domains and elucidated its 3D structure using NMR spectroscopy. Relying on an NMR-driven docking approach and MD simulations, we deciphered its binding interface with the C-terminal tails of both HSP70 and HSP90. Finally, we addressed the biological function of SPAG1 and specifically demonstrated that a SPAG1 sub-fragment, containing a putative P-loop motif, cannot efficiently bind and hydrolyze GTP in vitro Our data challenge the interpretation of SPAG1 possessing GTPase activity. We propose instead that SPAG1 regulates nucleotide hydrolysis activity of the HSP and RUVBL1/2 partners., (© 2019 The Author(s). Published by Portland Press Limited on behalf of the Biochemical Society.)

The hepatitis C virus (HCV) viroporin p7 oligomerizes to form ion channels, which are required for the assembly and secretion of infectious viruses. The 63-amino acid p7 monomer has two putative transmembrane domains connected by a cytosolic loop, and has both N- and C- termini exposed to the endoplasmic reticulum (ER) lumen. NMR studies have indicated differences between p7 structures of distantly related HCV genotypes. A critical question is whether these differences arise from the high sequence variation between the different isolates and if so, how the divergent structures can support similar biological functions. Here, we present a side-by-side characterization of p7 derived from genotype 1b (isolate J4) in the detergent 6-cyclohexyl-1-hexylphosphocholine (Cyclofos-6) and p7 derived from genotype 5a (isolate EUH1480) in n-dodecylphosphocholine (DPC). The 5a isolate p7 in conditions previously associated with a disputed oligomeric form exhibits secondary structure, dynamics, and solvent accessibility broadly like those of the monomeric 1b isolate p7. The largest differences occur at the start of the second transmembrane domain, which is destabilized in the 5a isolate. The results show a broad consensus among the p7 variants that have been studied under a range of different conditions and indicate that distantly related HCVs preserve key features of structure and dynamics.

Endonuclease III (EndoIII) is a bifunctional DNA glycosylase that removes oxidized pyrimidines from DNA. The genome of Deinococcus radiodurans encodes for an unusually high number of DNA glycosylases, including three EndoIII enzymes (drEndoIII1-3). Here, we compare the properties of these enzymes to those of their well-studied homologues from E. coli and human. Our biochemical and mutational data, reinforced by MD simulations of EndoIII-DNA complexes, reveal that drEndoIII2 exhibits a broad substrate specificity and a catalytic efficiency surpassing that of its counterparts. In contrast, drEndoIII1 has much weaker and uncoupled DNA glycosylase and AP-lyase activities, a characteristic feature of eukaryotic DNA glycosylases, and was found to present a relatively robust activity on single-stranded DNA substrates. To our knowledge, this is the first report of such an activity for an EndoIII. In the case of drEndoIII3, no catalytic activity could be detected, but its ability to specifically recognize lesion-containing DNA using a largely rearranged substrate binding pocket suggests that it may play an alternative role in genome maintenance. Overall, these findings reveal that D. radiodurans possesses a unique set of DNA repair enzymes, including three non-redundant EndoIII variants with distinct properties and complementary activities, which together contribute to genome maintenance in this bacterium., (Copyright © 2019 The Authors. Published by Elsevier B.V. All rights reserved.)

B- to A-DNA transition is known to be sensitive to the macroscopic properties of the solution, such as salt and ethanol concentrations. Microenvironmental effects on DNA conformational transition have been broadly studied. Providing an intuitive picture of how DNA responds to environmental changes is, however, still needed. Analyzing the chemical equilibrium of B-to-A DNA transition at critical concentrations, employing explicit-solvent simulations, is envisioned to help understand such microenvironmental effects. In the present study, free-energy calculations characterizing the B- to A-DNA transition and the distribution of cations were carried out in solvents with different ethanol concentrations. With the addition of ethanol, the most stable structure of DNA changes from the B- to A-form, in agreement with previous experimental observation. In 60% ethanol, a chemical equilibrium is found, showing reversible transition between B- and A-DNA. Analysis of the microenvironment around DNA suggests that with the increase of ethanol concentration, the cations exhibit a significant tendency to move toward the backbone, and mobility of water molecules around the major groove and backbone decreases gradually, leading eventually to a B-to-A transition. The present results provide a free-energy view of DNA microenvironment and of the role of cation motion in the conformational transition.

Aromatic groups can engage in an interesting class of noncovalent interactions termed π-π interactions, which play a pivotal role in stabilizing a variety of molecular architectures, including nucleic acids, proteins, and supramolecular assemblies. When the aromatic compounds interact with each other in an aqueous environment, their association is facilitated by the hydrophobic effect-the trend of nonpolar solutes to aggregate in a polar solution. To develop an in-depth understanding of hydrophobic association, we investigate in the present work π-π interactions in water, employing as a paradigm the benzene dimer. Using DFT-CES, a mean-field QM/MM method recently developed by our group, we describe the benzene solute at a quantum-mechanical level. Full consideration of detailed solute-electron density enables an optimal description of the solute-solvent interactions, leading to an accurate prediction of hydration free energies. In π-stacking of benzene, we find an entropic stabilization associated with the shrinkage of the solvent-excluded volume, which agrees with the theory of hydrophobic effect at subnanoscales. However, at the equilibrium binding distance of the benzene dimer, we find that the entropic stabilization nearly cancels out due to the enthalpic cost required for dewetting the internal space. Such an enthalpy-entropy compensation leads the association free energy to be predominantly dictated by the solute-solute interaction enthalpy. The present work offers new insight into the mechanistic role of water and the primary thermodynamic driving force of hydrophobic association.

The serotonin 5-HT 3 receptor is a pentameric ligand-gated ion channel (pLGIC). It belongs to a large family of receptors that function as allosteric signal transducers across the plasma membrane 1,2 ; upon binding of neurotransmitter molecules to extracellular sites, the receptors undergo complex conformational transitions that result in transient opening of a pore permeable to ions. 5-HT 3 receptors are therapeutic targets for emesis and nausea, irritable bowel syndrome and depression 3 . In spite of several reported pLGIC structures 4-8 , no clear unifying view has emerged on the conformational transitions involved in channel gating. Here we report four cryo-electron microscopy structures of the full-length mouse 5-HT 3 receptor in complex with the anti-emetic drug tropisetron, with serotonin, and with serotonin and a positive allosteric modulator, at resolutions ranging from 3.2 Å to 4.5 Å. The tropisetron-bound structure resembles those obtained with an inhibitory nanobody 5 or without ligand 9 . The other structures include an 'open' state and two ligand-bound states. We present computational insights into the dynamics of the structures, their pore hydration and free-energy profiles, and characterize movements at the gate level and cation accessibility in the pore. Together, these data deepen our understanding of the gating mechanism of pLGICs and capture ligand binding in unprecedented detail.

In biological environments and in aqueous solution, DNA generally adopts the canonical B conformation. Recently, an azobenzene photoswitch containing a polyamine chain with three positive charges was shown to induce a reversible conformational transition between the A and B forms of DNA, the transition being triggered by trans-cis isomerization of the photoswitch upon non-covalent intercalation. It was proposed that, in its trans conformation, azobenzene stabilizes the A form of DNA. The structural details and the mechanism upon which trans-azobenzene induces the B-to-A DNA transition remain, however, unclear. In the present work, two possible intercalating modes of trans-azobenzene, from the minor groove and from the major groove, were investigated with all-atom molecular-dynamics simulations. Intercalation from the major groove was found to be the most probable binding mode due to favorable electrostatic and π-π stacking interactions. The free-energy profile associated with the B-to-A conformational transition reveals that intercalation from the major groove leads to a conformational change of DNA, showing a slight tendency to interconvert from B- to A-DNA, in agreement with the CD spectrum obtained from the experiment. However, the presence of only one interacting azobenzene is not sufficient to lead to a global conformational change to A-DNA. The present results are expected to serve in the design of DNA switches, which can induce reversible DNA conformational changes.

Membrane proteins perform a host of vital cellular functions. Deciphering the molecular mechanisms whereby they fulfill these functions requires detailed biophysical and structural investigations. Detergents have proven pivotal to extract the protein from its native surroundings. Yet, they provide a milieu that departs significantly from that of the biological membrane, to the extent that the structure, the dynamics, and the interactions of membrane proteins in detergents may considerably vary, as compared to the native environment. Understanding the impact of detergents on membrane proteins is, therefore, crucial to assess the biological relevance of results obtained in detergents. Here, we review the strengths and weaknesses of alkyl phosphocholines (or foscholines), the most widely used detergent in solution-NMR studies of membrane proteins. While this class of detergents is often successful for membrane protein solubilization, a growing list of examples points to destabilizing and denaturing properties, in particular for α-helical membrane proteins. Our comprehensive analysis stresses the importance of stringent controls when working with this class of detergents and when analyzing the structure and dynamics of membrane proteins in alkyl phosphocholine detergents.

The low bioavailability of most therapeutic compounds is often counterbalanced by association with molecular vectors capable of crossing cell membranes. Previous studies demonstrated that for vectors bearing titratable chemical groups, the translocation process might be accompanied by a change in the protonation state. For simple compounds e.g. a lysine analog, free energy calculations, using a single collective variable, namely the insertion depth, suggest that such a transition could only take place if the amino acid diffuses deep enough into the hydrophobic core of the membrane, a situation thermodynamically unfavorable. Here, we determined the 2D potential of mean force associated with the translocation of lysine across a model membrane using as reaction coordinates not only its location in the bilayer but also its hydration. Our results cogently demonstrate that the change in protonation can result from a small fluctuation in the latter, even at low insertion depth.

Characterizing the structure of membrane proteins (MPs) generally requires extraction from their native environment, most commonly with detergents. Yet, the physicochemical properties of detergent micelles and lipid bilayers differ markedly and could alter the structural organization of MPs, albeit without general rules. Dodecylphosphocholine (DPC) is the most widely used detergent for MP structure determination by NMR, but the physiological relevance of several prominent structures has been questioned, though indirectly, by other biophysical techniques, e.g., functional/thermostability assay (TSA) and molecular dynamics (MD) simulations. Here, we resolve unambiguously this controversy by probing the functional relevance of three different mitochondrial carriers (MCs) in DPC at the atomic level, using an exhaustive set of solution-NMR experiments, complemented by functional/TSA and MD data. Our results provide atomic-level insight into the structure, substrate interaction and dynamics of the detergent-membrane protein complexes and demonstrates cogently that, while high-resolution NMR signals can be obtained for MCs in DPC, they systematically correspond to nonfunctional states.

DNA is the target of chemical compounds (drugs, pollutants, photosensitizers, etc.), which bind through non-covalent interactions. Depending on their structure and their chemical properties, DNA binders can associate to the minor or to the major groove of double-stranded DNA. They can also intercalate between two adjacent base pairs, or even replace one or two base pairs within the DNA double helix. The subsequent biological effects are strongly dependent on the architecture of the binding motif. Discriminating between the different binding patterns is of paramount importance to predict and rationalize the effect of a given compound on DNA. The structural characterization of DNA complexes remains, however, cumbersome at the experimental level. In this contribution, we employed all-atom molecular dynamics simulations to determine the standard binding free energy of DNA with netropsin, a well-characterized antiviral and antimicrobial drug, which associates to the minor groove of double-stranded DNA. To overcome the sampling limitations of classical molecular dynamics simulations, which cannot capture the large change in configurational entropy that accompanies binding, we resort to a series of potentials of mean force calculations involving a set of geometrical restraints acting on collective variables., Competing Interests: The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

Biophysical investigation of membrane proteins generally requires their extraction from native sources using detergents, a step that can lead, possibly irreversibly, to protein denaturation. The propensity of dodecylphosphocholine (DPC), a detergent widely utilized in NMR studies of membrane proteins, to distort their structure has been the subject of much controversy. It has been recently proposed that the binding specificity of the yeast mitochondrial ADP/ATP carrier (yAAC3) toward cardiolipins is preserved in DPC, thereby suggesting that DPC is a suitable environment in which to study membrane proteins. In this communication, we used all-atom molecular dynamics simulations to investigate the specific binding of cardiolipins to yAAC3. Our data demonstrate that the interaction interface observed in a native-like environment differs markedly from that inferred from an NMR investigation in DPC, implying that in this detergent, the protein structure is distorted. We further investigated yAAC3 solubilized in DPC and in the milder dodecylmaltoside with thermal-shift assays. The loss of thermal transition observed in DPC confirms that the protein is no longer properly folded in this environment., (Copyright © 2017 Biophysical Society. All rights reserved.)

Resorting to state-of-the art molecular modeling and simulation techniques we provide full characterization of the photophysical properties of the naturally occurring hypericin chromophore, currently used in photodynamic therapy. In particular, we reveal the different photophysical pathways leading to intersystem-crossing and hence, triplet manifold population that is necessary for the subsequent production of singlet oxygen. In particular we identify an extended region of quasi-degeneracy between the first singlet excited state and three triplet state surfaces. This energetic factor allows the occurrence of intersystem-crossing even in the presence of a relatively small spin-orbit coupling. Furthermore, thanks to extended all-atom molecular dynamics simulations we provide insight into the interaction of hypericin with lipid bilayers. We demonstrate the formation of stable interactions with the membrane and, in particular, the penetration of hypericin into its hydrophobic core. This organization allows a spatial overlap between hypericin and the lipid oxidizable double bond pointing towards the production of singlet oxygen in close spatial proximity to its reactant, hence favoring photosensitization.

DNA photolesions constitute a particularly deleterious class of molecular defects responsible for the insurgence of a vast majority of skin malignant tumors. Dimerization of two adjacent thymines or cytosines mostly gives rise to cyclobutane pyrimidine dimers (CPD) and pyrimidine(6-4)pyrimidone 64-PP as the most common defects. We perform all-atom classical simulations, up to 2 μs, of CPD and 64-PP embedded in a 16-bp duplex, which reveal the constrasted behavior of the two lesions. In particular we evidence a very limited structural deformation induced by CPD while 64-PP is characterized by a complex structural polymorphism. Our simulations also allow to unify the contrasting experimental structural results obtained by nuclear magnetic resonance or Förster Resonant Energy Transfer method, showing that both low and high bent structures are indeed accessible. These contrasting behaviors can also explain repair resistance or the different replication obstruction, and hence the genotoxicity of these two photolesions., (© The Author(s) 2017. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Cation-π interactions of aromatic rings and positively charged groups are among the most important interactions in structural biology. The role and energetic characteristics of these interactions are well established. However, the occurrence of cation-π-cation interactions is an unexpected motif, which raises intriguing questions about its functional role in proteins. We present a statistical analysis of the occurrence, composition and geometrical preferences of cation-π-cation interactions identified in a set of non-redundant protein structures taken from the Protein Data Bank. Our results demonstrate that this structural motif is observed at a small, albeit non-negligible frequency in proteins, and suggest a preference to establish cation-π-cation motifs with Trp, followed by Tyr and Phe. Furthermore, we have found that cation-π-cation interactions tend to be highly conserved, which supports their structural or functional role. Finally, we have performed an energetic analysis of a representative subset of cation-π-cation complexes combining quantum-chemical and continuum solvation calculations. Our results point out that the protein environment can strongly screen the cation-cation repulsion, leading to an attractive interaction in 64% of the complexes analyzed. Together with the high degree of conservation observed, these results suggest a potential stabilizing role in the protein fold, as demonstrated recently for a miniature protein (Craven et al., J. Am. Chem. Soc. 2016, 138, 1543). From a computational point of view, the significant contribution of non-additive three-body terms challenges the suitability of standard additive force fields for describing cation-π-cation motifs in molecular simulations.

β-Galactosidase (β-Gal) is one of the most important enzymes used in milk processing for improving their nutritional quality and digestibility. Herein, β-Gal has been entrapped into a meso-macroporous material (average pore size 9 and 200 nm, respectively) prepared by a sol-gel method from a silica precursor and a dispersion of solid lipid nanoparticles in a micelle phase. The physisorption of the enzyme depends on the concentration of the feed solution and on the pore size of the support. The enzyme is preferentially adsorbed either in mesopores or in macropores, depending on its initial concentration. Moreover, this selective adsorption, arising from the oligomeric complexation of the enzyme (monomer/dimer/tetramer), has an effect on the catalytic activity of the material. Indeed, the enzyme encapsulated in macropores is more active than the enzyme immobilized in mesopores. Designed materials containing β-Gal are of particular interest for food applications and potentially extended to bioconversion, bioremediation, or biosensing when coupling the designed support with other enzymes.