Binding of Cholesterol to the N-Terminal Domain of the NPC1L1 Transporter: Analysis of the Epimerization-Related Binding Selectivity and Loop Mutations.

Cholesterol is a fat-like substance with a pivotal physiological relevance in humans, and its homeostasis is tightly regulated by various cellular processes, including the import in the small intestine and the reabsorption in the biliary ducts by the Niemann-Pick C1 Like 1 (NPC1L1) importer. NPC1L1 can mediate the absorption of a variety of sterols but strikingly exhibits a large sensitivity to cholesterol epimerization. This study examines the molecular basis of the epimerization-related selective binding of cholesterol by combining extended unbiased molecular dynamics simulations of the apo and holo species of the N-terminal domain of wild-type NPC1L1, in conjunction with relative binding free energy, umbrella sampling, and well-tempered metadynamics calculations. The analysis of the results discloses the existence of two distinct binding modes for cholesterol and epi-cholesterol. The former binds deeper in the cavity, forming key hydrogen-bond interactions with Q95, S56, and a water molecule. In contrast, epi-cholesterol is shifted ca. 3 Å to the mouth of the cavity and the transition to the Q95 site is prevented by an energetic barrier of 4.1 kcal·mol-1. Thus, the configuration of the hydroxyl group of cholesterol, together with the presence of a structural water molecule, is a key feature for effective absorption. Finally, whereas these findings may seemingly be challenged by single-point mutations that impair cholesterol transport but have a mild impact on the binding of cholesterol to the Q95 binding site, our results reveal that they have a drastic influence on the conformational landscape of the α8/ß7 loop in the apo species compared to the wild-type protein. Overall, the results give support to the functional role played by the α8/ß7 loop in regulating the access of ligands to NPC1L1, and hence to interpreting the impact of these mutations on diseases related to disruption of sterol absorption, paving the way to understanding certain physiological dysfunctions.

Disease related single point mutations alter the global dynamics of a tetratricopeptide (TPR) α-solenoid domain.

Tetratricopeptide repeat (TPR) proteins belong to the class of α-solenoid proteins, in which repetitive units of α-helical hairpin motifs stack to form superhelical, often highly flexible structures. TPR domains occur in a wide variety of proteins, and perform key functional roles including protein folding, protein trafficking, cell cycle control and post-translational modification. Here, we look at the TPR domain of the enzyme O-linked GlcNAc-transferase (OGT), which catalyses O-GlcNAcylation of a broad range of substrate proteins. A number of single-point mutations in the TPR domain of human OGT have been associated with the disease Intellectual Disability (ID). By extended steered and equilibrium atomistic simulations, we show that the OGT-TPR domain acts as an elastic nanospring, and that each of the ID-related local mutations substantially affect the global dynamics of the TPR domain. Since the nanospring character of the OGT-TPR domain is key to its function in binding and releasing OGT substrates, these changes of its biomechanics likely lead to defective substrate interaction. We find that neutral mutations in the human population, selected by analysis of the gnomAD database, do not incur these changes. Our findings may not only help to explain the ID phenotype of the mutants, but also aid the design of TPR proteins with tailored biomechanical properties.

Combining Virtual Reality Visualization with Ensemble Molecular Dynamics to Study Complex Protein Conformational Changes.

Molecular dynamics (MD) simulations are increasingly used to elucidate relationships between protein structure, dynamics, and their biological function. Currently, it is extremely challenging to perform MD simulations of large-scale structural rearrangements in proteins that occur on millisecond timescales or beyond, as this requires very significant computational resources, or the use of cumbersome "collective variable" enhanced sampling protocols. Here, we describe a framework that combines ensemble MD simulations and virtual reality visualization (eMD-VR) to enable users to interactively generate realistic descriptions of large amplitude, millisecond timescale protein conformational changes in proteins. Detailed tests demonstrate that eMD-VR substantially decreases the computational cost of folding simulations of a WW domain, without the need to define collective variables a priori. We further show that eMD-VR generated pathways can be combined with Markov state models to describe the thermodynamics and kinetics of large-scale loop motions in the enzyme cyclophilin A. Our results suggest eMD-VR is a powerful tool for exploring protein energy landscapes in bioengineering efforts.

On the Binding of Congo Red to Amyloid Fibrils.

Amyloids are characterized by their capacity to bind Congo red (CR), one of the most used amyloid-specific dyes. The structural features of CR binding were unknown for years, mainly because of the lack of amyloid structures solved at high resolution. In the last few years, solid-state NMR spectroscopy enabled the determination of the structural features of amyloids, such as the HET-s prion forming domain (HET-s PFD), which also has recently been used to determine the amyloid-CR interface at atomic resolution. Herein, we combine spectroscopic data with molecular docking, molecular dynamics, and excitonic quantum/molecular mechanics calculations to examine and rationalize CR binding to amyloids. In contrast to a previous assumption on the binding mode, our results suggest that CR binding to the HET-s PFD involves a cooperative process entailing the formation of a complex with 1:1 stoichiometry. This provides a molecular basis to explain the bathochromic shift in the maximal absorbance wavelength when CR is bound to amyloids.

Mechanism of the Pseudoirreversible Binding of Amantadine to the M2 Proton Channel.

The M2 proton channel of influenza A virus is an integral membrane protein involved in the acidification of the viral interior, a step necessary for the release of the viral genetic material and replication of new virions. The aim of this study is to explore the mechanism of drug (un)binding to the M2 channel in order to gain insight into the structural and energetic features relevant for the development of novel inhibitors. To this end, we have investigated the binding of amantadine (Amt) to the wild type (wt) M2 channel and its V27A variant using multiple independent molecular dynamics simulations, exploratory conventional metadynamics, and multiple-walkers well-tempered metadynamics calculations. The results allow us to propose a sequential mechanism for the (un)binding of Amt to the wt M2 channel, which involves the adoption of a transiently populated intermediate (up state) leading to the thermodynamically favored down binding mode in the channel pore. Furthermore, they suggest that chloride anions play a relevant role in stabilizing the down binding mode of Amt to the wt channel, giving rise to a kinetic trapping that explains the experimentally observed pseudoirreversible inhibition of the wt channel by Amt. We propose that this trapping mechanism underlies the inhibitory activity of potent M2 channel blockers, as supported by the experimental confirmation of the irreversible binding of a pyrrolidine analogue from electrophysiological current assays. Finally, the results reveal that the thermodynamics and kinetics of Amt (un)binding is very sensitive to the V27A mutation, providing a quantitative rationale to the drastic decrease in inhibitory potency against the V27A variant. Overall, these findings pave the way to explore the inhibitory activity of Amt-related analogues in mutated M2 channel variants, providing guidelines for the design of novel inhibitors against resistant virus strains.

Evolution of a multicomponent system: computational and mechanistic studies on the chemo- and stereoselectivity of a divergent process.

The evolution of a ternary molecular system (imine, diene and nitrile) is analyzed to disclose the pathways leading to a divergent synthetic outcome. The Lewis acid catalyzed reaction between cyclohexadiene, 2-phenyl-indol-3-one and acetonitrile yields the imino-Diels-Alder adduct as the major product together with minor amounts of the Mannich-Ritter-amidine product. The experimental and computational data show that the relative orientation of the initial reactants dictates the synthetic outcome. The exo approach between imine and diene leads to the Diels-Alder adduct in a concerted process, whereas the endo mode leads to a polarized intermediate, which is trapped by acetonitrile to yield the multicomponent adduct.

Exploration of forbidden Povarov processes as a source of unexpected reactivity: a multicomponent Mannich-Ritter transformation.

When a door closes, a window opens! The use of geometrically or electronically restricted imines for Povarov-type processes does not afford the anti-Bredt tetrahydroquinolines, but leads instead to highly functionalized structures through novel reaction pathways (see picture; LA=Lewis acid). The exploration of "forbidden" routes constitutes a valuable approach in the search for new multicomponent reactions.

Energetics of a protein disorder-order transition in small molecule recognition.

Many proteins recognise other proteins via mechanisms that involve the folding of intrinsically disordered regions upon complex formation. Here we investigate how the selectivity of a drug-like small molecule arises from its modulation of a protein disorder-to-order transition. Binding of the compound AM-7209 has been reported to confer order upon an intrinsically disordered 'lid' region of the oncoprotein MDM2. Calorimetric measurements revealed that truncation of the lid region of MDM2 increases the apparent dissociation constant of AM-7209 250-fold. By contrast, lid truncation has little effect on the binding of the ligand Nutlin-3a. Insights into these differential binding energetics were obtained via a complete thermodynamic analysis that featured adaptive absolute alchemical free energy of binding calculations with enhanced-sampling molecular dynamics simulations. The simulations reveal that in apo MDM2 the ordered lid state is energetically disfavoured. AM-7209, but not Nutlin-3a, shows a significant energetic preference for ordered lid conformations, thus shifting the balance towards ordering of the lid in the AM-7209/MDM2 complex. The methodology reported herein should facilitate broader targeting of intrinsically disordered regions in medicinal chemistry.

An antibiotic-resistance conferring mutation in a neisserial porin: Structure, ion flux, and ampicillin binding.

Gram-negative bacteria cause the majority of highly drug-resistant bacterial infections. To cross the outer membrane of the complex Gram-negative cell envelope, antibiotics permeate through porins, trimeric channel proteins that enable the exchange of small polar molecules. Mutations in porins contribute to the development of drug-resistant phenotypes. In this work, we show that a single point mutation in the porin PorB from Neisseria meningitidis, the causative agent of bacterial meningitis, can strongly affect the binding and permeation of beta-lactam antibiotics. Using X-ray crystallography, high-resolution electrophysiology, atomistic biomolecular simulation, and liposome swelling experiments, we demonstrate differences in drug binding affinity, ion selectivity and drug permeability of PorB. Our work further reveals distinct interactions between the transversal electric field in the porin eyelet and the zwitterionic drugs, which manifest themselves under applied electric fields in electrophysiology and are altered by the mutation. These observations may apply more broadly to drug-porin interactions in other channels. Our results improve the molecular understanding of porin-based drug-resistance in Gram-negative bacteria.

High-resolution experimental and computational electrophysiology reveals weak ß-lactam binding events in the porin PorB.

The permeation of most antibiotics through the outer membrane of Gram-negative bacteria occurs through porin channels. To design drugs with increased activity against Gram-negative bacteria in the face of the antibiotic resistance crisis, the strict constraints on the physicochemical properties of the permeants imposed by these channels must be better understood. Here we show that a combination of high-resolution electrophysiology, new noise-filtering analysis protocols and atomistic biomolecular simulations reveals weak binding events between the ß-lactam antibiotic ampicillin and the porin PorB from the pathogenic bacterium Neisseria meningitidis. In particular, an asymmetry often seen in the electrophysiological characteristics of ligand-bound channels is utilised to characterise the binding site and molecular interactions in detail, based on the principles of electro-osmotic flow through the channel. Our results provide a rationale for the determinants that govern the binding and permeation of zwitterionic antibiotics in porin channels.

The O-GlcNAc Transferase Intellectual Disability Mutation L254F Distorts the TPR Helix.

O-linked ß-N-acetyl-D-glucosamine (O-GlcNAc) transferase (OGT) regulates protein O-GlcNAcylation, an essential post-translational modification that is abundant in the brain. Recently, OGT mutations have been associated with intellectual disability, although it is not understood how they affect OGT structure and function. Using a multi-disciplinary approach we show that the L254F OGT mutation leads to conformational changes of the tetratricopeptide repeats and reduced activity, revealing the molecular mechanisms contributing to pathogenesis.

Modulation of the Neisseria gonorrhoeae drug efflux conduit MtrE.

Widespread antibiotic resistance, especially of Gram-negative bacteria, has become a severe concern for human health. Tripartite efflux pumps are one of the major contributors to resistance in Gram-negative pathogens, by efficiently expelling a broad spectrum of antibiotics from the organism. In Neisseria gonorrhoeae, one of the first bacteria for which pan-resistance has been reported, the most expressed efflux complex is MtrCDE. Here we present the electrophysiological characterisation of the outer membrane component MtrE and the membrane fusion protein MtrC, obtained by a combination of planar lipid bilayer recordings and in silico techniques. Our in vitro results show that MtrE can be regulated by periplasmic binding events and that the interaction between MtrE and MtrC is sufficient to stabilize this complex in an open state. In contrast to other efflux conduits, the open complex only displays a slight preference for cations. The maximum conductance we obtain in the in vitro recordings is comparable to that seen in our computational electrophysiology simulations conducted on the MtrE crystal structure, indicating that this state may reflect a physiologically relevant open conformation of MtrE. Our results suggest that the MtrC/E binding interface is an important modulator of MtrE function, which could potentially be targeted by new efflux inhibitors.

Structure of the MacAB-TolC ABC-type tripartite multidrug efflux pump.

The MacA-MacB-TolC assembly of Escherichia coli is a transmembrane machine that spans the cell envelope and actively extrudes substrates, including macrolide antibiotics and polypeptide virulence factors. These transport processes are energized by the ATPase MacB, a member of the ATP-binding cassette (ABC) superfamily. We present an electron cryo-microscopy structure of the ABC-type tripartite assembly at near-atomic resolution. A hexamer of the periplasmic protein MacA bridges between a TolC trimer in the outer membrane and a MacB dimer in the inner membrane, generating a quaternary structure with a central channel for substrate translocation. A gating ring found in MacA is proposed to act as a one-way valve in substrate transport. The MacB structure features an atypical transmembrane domain with a closely packed dimer interface and a periplasmic opening that is the likely portal for substrate entry from the periplasm, with subsequent displacement through an allosteric transport mechanism.

New polycyclic dual inhibitors of the wild type and the V27A mutant M2 channel of the influenza A virus with unexpected binding mode.

Two new polycyclic scaffolds were synthesized and evaluated as anti-influenza A compounds. The 5-azapentacyclo[6.4.0.0(2,10).0(3,7).0(9,11)]dodecane derivatives were only active against the wild-type M2 channel in the low-micromolar range. However, some of the 14-azaheptacyclo[8.6.1.0(2,5).0(3,11).0(4,9).0(6,17).0(12,16)]heptadecane derivatives were dual inhibitors of the wild-type and the V27A mutant M2 channels. The antiviral activity of these molecules was confirmed by cell culture assays. Their binding mode was analysed through molecular dynamics simulations, which showed the existence of distinct binding modes in the wild type M2 channel and its V27A variant.

Molecular basis of the selective binding of MDMA enantiomers to the alpha4beta2 nicotinic receptor subtype: synthesis, pharmacological evaluation and mechanistic studies.

The α4ß2 nicotinic acetylcholine receptor (nAChR) is a molecular target of 3,4-methylenedioxymethamphetamine (MDMA), a synthetic drug also known as ecstasy, and it modulates the MDMA-mediated reinforcing properties. However, the enantioselective preference of the α4ß2 nAChR subtype still remains unknown. Since the two enantiomers exhibit different pharmacological profiles and stereoselective metabolism, the aim of this study is to assess a possible difference in the interaction of the MDMA enantiomers with this nAChR subtype. To this end, we report a novel simple, yet highly efficient enantioselective synthesis of the MDMA enantiomers, in which the key step is the diastereoselective reduction of imides derived from optically pure tert-butylsulfinamide. The enantioselective binding to the receptor is examined using [(3)H]epibatidine in a radioligand assay. Even though the two enantiomers induced a concentration-dependent binding displacement, (S)-MDMA has an inhibition constant 13-fold higher than (R)-MDMA, which shows a Hill's coefficient not significantly different from unity, implying a competitive interaction. Furthermore, when NGF-differentiated PC12 cells were pretreated with the compounds, a significant increase in binding of [(3)H]epibatidine was found for (R)-MDMA, indicating up-regulation of heteromeric nAChR in the cell surface. Finally, docking and molecular dynamics studies have been used to identify the binding mode of the two enantiomers, which provides a structural basis to justify the differences in affinity from the differential interactions played by the substituents at the stereogenic centre of MDMA. The results provide a basis to explore the distinct psychostimulant profiles of the MDMA enantiomers mediated by the α4ß2 nAChR subtype.

Easily accessible polycyclic amines that inhibit the wild-type and amantadine-resistant mutants of the M2 channel of influenza A virus.

Amantadine inhibits the M2 proton channel of influenza A virus, yet most of the currently circulating strains of the virus carry mutations in the M2 protein that render the virus amantadine-resistant. While most of the research on novel amantadine analogues has revolved around the synthesis of novel adamantane derivatives, we have recently found that other polycyclic scaffolds effectively block the M2 proton channel, including amantadine-resistant mutant channels. In this work, we have synthesized and characterized a series of pyrrolidine derivatives designed as analogues of amantadine. Inhibition of the wild-type M2 channel and the A/M2-S31N, A/M2-V27A, and A/M2-L26F mutant forms of the channel were measured in Xenopus oocytes using two-electrode voltage clamp assays. Most of the novel compounds inhibited the wild-type ion channel in the low micromolar range. Of note, two of the compounds inhibited the amantadine-resistant A/M2-V27A and A/M2-L26F mutant ion channels with submicromolar and low micromolar IC50, respectively. None of the compounds was found to inhibit the S31N mutant ion channel.