Structural determinants for activity of the antidepressant vortioxetine at human and rodent 5-HT3 receptors.

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-HT3 receptor (5-HT3R), which features two properties: VTX acts differently on rodent and human 5-HT3R, 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-HT3R and an agonist-bound-like state of human 5-HT3R, in line with the functional profile of the drug. We report four human 5-HT3R 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.

Configurational Sampling of All-Atom Solvated Membranes Using Hybrid Nonequilibrium Molecular Dynamics Monte Carlo Simulations.

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.

State-dependent binding of cholesterol and an anionic lipid to the muscle-type Torpedo nicotinic acetylcholine receptor.

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.

Biochemical and Structural Characterization of Chi-Class Glutathione Transferases: A Snapshot on the Glutathione Transferase Encoded by sll0067 Gene in the Cyanobacterium Synechocystis sp. Strain PCC 6803.

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 (10SRAS13 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.

Perforin-2 clockwise hand-over-hand pre-pore to pore transition mechanism.

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.

Structural and functional characterization of DdrC, a novel DNA damage-induced nucleoid associated protein involved in DNA compaction.

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.

Computational Assessment of Protein-Protein Binding Specificity within a Family of Synaptic Surface Receptors.

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.

Accurate determination of protein:ligand standard binding free energies from molecular dynamics simulations.

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.

Conformational transitions and ligand-binding to a muscle-type nicotinic acetylcholine receptor.

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.

Structure, substrate binding, and symmetry of the mitochondrial ADP/ATP carrier in its matrix-open state.

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.

Structure of HIV-1 gp41 with its membrane anchors targeted by neutralizing antibodies.

The HIV-1 gp120/gp41 trimer undergoes a series of conformational changes in order to catalyze gp41-induced fusion of viral and cellular membranes. Here, we present the crystal structure of gp41 locked in a fusion intermediate state by an MPER-specific neutralizing antibody. The structure illustrates the conformational plasticity of the six membrane anchors arranged asymmetrically with the fusion peptides and the transmembrane regions pointing into different directions. Hinge regions located adjacent to the fusion peptide and the transmembrane region facilitate the conformational flexibility that allows high-affinity binding of broadly neutralizing anti-MPER antibodies. Molecular dynamics simulation of the MPER Ab-stabilized gp41 conformation reveals a possible transition pathway into the final post-fusion conformation with the central fusion peptides forming a hydrophobic core with flanking transmembrane regions. This suggests that MPER-specific broadly neutralizing antibodies can block final steps of refolding of the fusion peptide and the transmembrane region, which is required for completing membrane fusion.

The Binding of Palonosetron and Other Antiemetic Drugs to the Serotonin 5-HT3 Receptor.

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.

Charge Transfer and Chemo-Mechanical Coupling in Respiratory Complex I.

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.

Modulation of membrane permeability by carbon dioxide.

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 (CO2 ) molecules could also enhance membrane permeability to drugs. In the present work, we have investigated the effect of CO2 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 CO2 in the lipid environment results in an increase of the membrane permeability to the three substrates. Further analysis of the permeation events reveals that CO2 expands and loosens the membrane, which, in turn, facilitates permeation of the drug-like molecules. © 2019 Wiley Periodicals, Inc.

Molecular Bases of DNA Packaging in Bacteria Revealed by All-Atom Molecular Dynamics Simulations: The Case of Histone-Like Proteins in Borrelia burgdorferi.

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.

Targeting G-quadruplexes with Organic Dyes: Chelerythrine-DNA Binding Elucidated by Combining Molecular Modeling and Optical Spectroscopy.

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.

Induced Night Vision by Singlet-Oxygen-Mediated Activation of Rhodopsin.

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.

Structural Basis for Broad HIV-1 Neutralization by the MPER-Specific Human Broadly Neutralizing Antibody LN01.

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.

Mechanism of the allosteric activation of the ClpP protease machinery by substrates and active-site inhibitors.

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.