Membrane Protein Simulations Using AMBER Force Field and Berger Lipid Parameters

Synthesis and Characterization of a Bidirectional Photoswitchable Antagonist Toolbox for Real-Time GPCR Photopharmacology

Noninvasive methods to modulate G protein-coupled receptors (GPCRs) with temporal and spatial precision are in great demand. Photopharmacology uses photons to control in situ the biological properties of photoswitchable small-molecule ligands, which bodes well for chemical biological precision approaches. Integrating the light-switchable configurational properties of an azobenzene into the ligand core, we developed a bidirectional antagonist toolbox for an archetypical family A GPCR, the histamine H₃ receptor (H₃R). From 16 newly synthesized photoswitchable compounds, VUF14738 (28) and VUF14862 (33) were selected as they swiftly and reversibly photoisomerize and show over 10-fold increased or decreased H₃R binding affinities, respectively, upon illumination at 360 nm. Both ligands combine long thermal half-lives with fast and high photochemical trans-/cis conversion, allowing their use in real-time electrophysiology experiments with oocytes to confirm dynamic photomodulation of H₃R activation in repeated second-scale cycles. VUF14738 and VUF14862 are robust and fatigue-resistant photoswitchable GPCR antagonists suitable for spatiotemporal studies of H₃R signaling.

Cross-communication between Gi and Gs in a G-protein-coupled receptor heterotetramer guided by a receptor C-terminal domain

Background

G-protein-coupled receptor (GPCR) heteromeric complexes have distinct properties from homomeric GPCRs, giving rise to new receptor functionalities. Adenosine receptors (A₁R or A_2AR) can form A₁R-A_2AR heteromers (A₁-A_2AHet), and their activation leads to canonical G-protein-dependent (adenylate cyclase mediated) and -independent (β-arrestin mediated) signaling. Adenosine has different affinities for A₁R and A_2AR, allowing the heteromeric receptor to detect its concentration by integrating the downstream G_i- and G_s-dependent signals. cAMP accumulation and β-arrestin recruitment assays have shown that, within the complex, activation of A_2AR impedes signaling via A₁R.

Results

We examined the mechanism by which A₁-A_2AHet integrates G_i- and G_s-dependent signals. A₁R blockade by A_2AR in the A₁-A_2AHet is not observed in the absence of A_2AR activation by agonists, in the absence of the C-terminal domain of A_2AR, or in the presence of synthetic peptides that disrupt the heteromer interface of A₁-A_2AHet, indicating that signaling mediated by A₁R and A_2AR is controlled by both G_i and G_s proteins.

Conclusions

We identified a new mechanism of signal transduction that implies a cross-communication between G_i and G_s proteins guided by the C-terminal tail of the A_2AR. This mechanism provides the molecular basis for the operation of the A₁-A_2AHet as an adenosine concentration-sensing device that modulates the signals originating at both A₁R and A_2AR.

Electronic supplementary material

The online version of this article (10.1186/s12915-018-0491-x) contains supplementary material, which is available to authorized users.

Development of novel fluorescent histamine H1-receptor antagonists to study ligand-binding kinetics in living cells

The histamine H1-receptor (H1R) is an important mediator of allergy and inflammation. H1R antagonists have particular clinical utility in allergic rhinitis and urticaria. Here we have developed six novel fluorescent probes for this receptor that are very effective for high resolution confocal imaging, alongside bioluminescence resonance energy transfer approaches to monitor H1R ligand binding kinetics in living cells. The latter technology exploits the opportunities provided by the recently described bright bioluminescent protein NanoLuc when it is fused to the N-terminus of a receptor. Two different pharmacophores (mepyramine or the fragment VUF13816) were used to generate fluorescent H1R antagonists conjugated via peptide linkers to the fluorophore BODIPY630/650. Kinetic properties of the probes showed wide variation, with the VUF13816 analogues having much longer H1R residence times relative to their mepyramine-based counterparts. The kinetics of these fluorescent ligands could also be monitored in membrane preparations providing new opportunities for future drug discovery applications.

Distinct modulation of inactivation by a residue in the pore domain of voltage-gated Na+ channels: mechanistic insights from recent crystal structures

Inactivation of voltage-gated Na⁺ channels (VGSC) is essential for the regulation of cellular excitability. The molecular rearrangement underlying inactivation is thought to involve the intracellular linker between domains III and IV serving as inactivation lid, the receptor for the lid (domain III S4-S5 linker) and the pore-lining S6 segements. To better understand the role of the domain IV S6 segment in inactivation we performed a cysteine scanning mutagenesis of this region in rNav 1.4 channels and screened the constructs for perturbations in the voltage-dependence of steady state inactivation. This screen was performed in the background of wild-type channels and in channels carrying the mutation K1237E, which profoundly alters both permeation and gating-properties. Of all tested constructs the mutation I1581C was unique in that the mutation-induced gating changes were strongly influenced by the mutational background. This suggests that I1581 is involved in specific short-range interactions during inactivation. In recently published crystal structures VGSCs the respective amino acids homologous to I1581 appear to control a bend of the S6 segment which is critical to the gating process. Furthermore, I1581 may be involved in the transmission of the movement of the DIII voltage-sensor to the domain IV S6 segment.

Identifying human diamine sensors for death related putrescine and cadaverine molecules

Pungent chemical compounds originating from decaying tissue are strong drivers of animal behavior. Two of the best-characterized death smell components are putrescine (PUT) and cadaverine (CAD), foul-smelling molecules produced by decarboxylation of amino acids during decomposition. These volatile polyamines act as ‘necromones’, triggering avoidance or attractive responses, which are fundamental for the survival of a wide range of species. The few studies that have attempted to identify the cognate receptors for these molecules have suggested the involvement of the seven-helix trace amine-associated receptors (TAARs), localized in the olfactory epithelium. However, very little is known about the precise chemosensory receptors that sense these compounds in the majority of organisms and the molecular basis of their interactions. In this work, we have used computational strategies to characterize the binding between PUT and CAD with the TAAR6 and TAAR8 human receptors. Sequence analysis, homology modeling, docking and molecular dynamics studies suggest a tandem of negatively charged aspartates in the binding pocket of these receptors which are likely to be involved in the recognition of these small biogenic diamines.

The distinctive dead smell comes largely from molecules like cadaverine and putrescine that are produced during decomposition of organic tissues. These volatile compounds act as powerful chemical signals important for the survival of a wide range of species. Previous studies have identified the trace amine-associated receptor 13c (or TAAR13c) in zebrafish as the cognate receptor of cadaverine in bony fishes. In this work, we employed computational strategies to disclose the human TAAR6 and TAAR8 receptors as sensors of the putrescine and cadaverine molecules. Our results indicate that several negatively charged residues in the ligand binding pocket of these receptors constitute the molecular basis for recognition of these necromones in humans.

Intracellular Transfer of Na+ in an Active-State G-Protein-Coupled Receptor

Playing a central role in cell signaling, G-protein-coupled receptors (GPCRs) are the largest superfamily of membrane proteins and form the majority of drug targets in humans. How extracellular agonist binding triggers the activation of GPCRs and associated intracellular effector proteins remains, however, poorly understood. Structural studies have revealed that inactive class A GPCRs harbor a conserved binding site for Na+ ions in the center of their transmembrane domain, accessible from the extracellular space. Here, we show that the opening of a conserved hydrated channel in the activated state receptors allows the Na+ ion to egress from its binding site into the cytosol. Coupled with protonation changes, this ion movement occurs without significant energy barriers, and can be driven by physiological transmembrane ion and voltage gradients. We propose that Na+ ion exchange with the cytosol is a key step in GPCR activation. Further, we hypothesize that this transition locks receptors in long-lived active-state conformations.

Why the Flavin Adenine Dinucleotide (FAD) Cofactor Needs To Be Covalently Linked to Complex II of the Electron-Transport Chain for the Conversion of FADH2 into FAD

A covalently bound flavin cofactor is predominant in the succinate‐ubiquinone oxidoreductase (SQR; Complex II), an essential component of aerobic electron transport, and in the menaquinol‐fumarate oxidoreductase (QFR), the anaerobic counterpart, although it is only present in approximately 10 % of the known flavoenzymes. This work investigates the role of this 8α‐N3‐histidyl linkage between the flavin adenine dinucleotide (FAD) cofactor and the respiratory Complex II. After parameterization with DFT calculations, classical molecular‐dynamics simulations and quantum‐mechanics calculations for Complex II:FAD and Complex II:FADH₂, with and without the covalent bond, were performed. It was observed that the covalent bond is essential for the active‐center arrangement of the FADH₂/FAD cofactor. Removal of this bond causes a displacement of the isoalloxazine group, which influences interactions with the protein, flavin solvation, and possible proton‐transfer pathways. Specifically, for the noncovalently bound FADH₂ cofactor, the N1 atom moves away from the His‐A365 and His‐A254 residues and the N5 atom moves away from the glutamine‐62A residue. Both of the histidine and glutamine residues interact with a chain of water molecules that cross the enzyme, which is most likely involved in proton transfer. Breaking this chain of water molecules could thereby compromise proton transfer across the two active sites of Complex II.

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.

Docking-based Screening of Ficus religiosa Phytochemicals as Inhibitors of Human Histamine H2 Receptor

Background:

Ficus religiosa L. is generally known as Peepal and belongs to family Moraceae. The tree is a source of many compounds having high medicinal value. In gastrointestinal tract, histamine H2 receptors have key role in histamine-stimulated gastric acid secretion. Their over stimulation causes its excessive production which is responsible for gastric ulcer.

Objective:

This study aims to screen the range of phytochemicals present in F. religiosa for binding with human histamine H2 and identify therapeutics for a gastric ulcer from the plant.

Materials and Methods:

In this work, a 3D-structure of human histamine H2 receptor was modeled by using homology modeling and the predicted model was validated using PROCHECK. Docking studies were also performed to assess binding affinities between modeled receptor and 34 compounds. Molecular dynamics simulations were done to identify most stable receptor-ligand complexes. Absorption, distribution, metabolism, excretion, and screening was done to evaluate pharmacokinetic properties of compounds.

Results:

The results suggest that seven ligands, namely, germacrene, bergaptol, lanosterol, Ergost-5-en-3beta-ol, α-amyrin acetate, bergapten, and γ-cadinene showed better binding affinities.

Conclusion:

Among seven phytochemicals, lanosterol and α-amyrin acetate were found to have greater stability during simulation studies. These two compounds may be a suitable therapeutic agent against histamine H2 receptor.

SUMMARY

This study was performed to screen antiulcer compounds from F. religiosa. Molecular modeling, molecular docking and MD simulation studies were performed with selected phytochemicals from F. religiosa. The analysis suggests that Lanosterol and α-amyrin may be a suitable therapeutic agent against histamine H2 receptor. This study facilitates initiation of the herbal drug discovery process for the antiulcer activity.

Abbreviations used: ADMET: Absorption, distribution, metabolism, excretion and toxicity, DOPE: Discrete Optimized Potential Energy, OPLS: Optimized potential for liquid simulations, RMSD: Root-mean-square deviation, HOA: Human oral absorption, MW: Molecular weight, SP: Standard-precision, XP: Extra-precision, GPCRs: G protein-coupled receptors, SASA: Solvent accessible surface area, Rg: Radius of gyration, NHB: Number of hydrogen bond

Alamandine reverses hyperhomocysteinemia-induced vascular dysfunction via PKA-dependent mechanisms

INTRODUCTION

Hyperhomocysteinemia (HHcy) impairs nitric oxide endothelium-dependent vasodilation, consequently leading to atherosclerosis, a risk factor for cardiovascular disease. Novel treatments for HHcy are necessary.

AIM

We tested the hypothesis that alamandine, a vasoactive peptide of the renin-angiotensin system (RAS), could reverse HHcy-induced vascular dysfunction through the MrgD receptor and that this is mediated by the protein kinase A (PKA) pathway. Furthermore, we sought to determine a putative binding model of alamandine to the MrgD receptor through docking and molecular dynamics simulations.

METHOD

The abdominal aorta was excised from New Zealand white rabbits (n = 15) and incubated with 3 mmol/L Hcy (to mimic HHcy) to induce vascular dysfunction in vitro. Vascular function was assessed by vasodilatory responses to cumulative doses of acetylcholine.

RESULT

Vasodilation was significantly impaired in HHcy-incubated aortic rings while alamandine reversed this effect (control, 74.2 ± 5.0%; Hcy, 30.3 ± 9.8%; alamandine + Hcy, 59.7 ± 4.8%, P < .0001). KT5720 (PKA inhibitor) significantly inhibited the ability of alamandine to attenuate the impaired vasodilation caused by HHcy (KT5720 + Hcy + alamandine, 27.1 ± 24.1, P < .01). Following immunohistochemistry analysis, the MrgD receptor was highly expressed within the media and endothelial layer of aortic rings in HHcy compared to control (media: 0.23 ± 0.003 vs control 0.16 ± 0.01, P < .05 and endothelium: 0.68 ± 0.07 vs control 0.13 ± 0.02, P < .01, in PA/I (A.U) units). Computational studies also propose certain interactions of alamandine within the MrgD transmembrane domain.

CONCLUSION

This study shows that alamandine is effective in reversing HHcy-induced vascular dysfunction, possibly through the PKA signaling pathway via MrgD. Our results indicate a therapeutic potential of alamandine in reversing the detrimental effects of HHcy.

A-type ECG and EGCG dimers inhibit 3T3-L1 differentiation by binding to cholesterol in lipid rafts

The present study aimed to explore the underlying mechanisms of epicatechin-3-gallate-(4β→8, 2β→O→7)-epicatechin-3-gallate (A-type ECG dimer) and epigallocatechin-3-gallate-(4β→8, 2β→O→7)-epigallocatechin-3-gallate (A-type EGCG dimer) involved in their strong inhibitory effects on 3T3-L1 preadipocytes differentiation. In the synthetic "lipid raft-like" liposome, A-type ECG and EGCG dimers incorporated into the liposome with high affinity and decreased the fluidity of the liposome. In 3T3-L1 preadipocytes, A-type ECG and EGCG dimers possibly bonded to lipid rafts cholesterol and disrupted the integrity of lipid rafts, thus exerting their notable inhibitory effects on 3T3-L1 preadipocytes differentiation by suppressing mitotic clonal expansion process and mRNA levels of PPARγ, C/EBPα and SREBP1C. A highly positive correlation between the cholesterol binding capacity of the two dimers and their inhibitory effect on 3T3-L1 preadipocytes differentiation (R²=0.9328) was observed. Molecular dynamics simulation further verified that A-type ECG and EGCG dimers could bond to cholesterol via hydrogen bonding. The results of this study suggested that the disruption of A-type ECG and EGCG dimers on membrane lipid rafts by targeting cholesterol in the lipid rafts was involved in the underlying mechanisms of their strong inhibitory effects on 3T3-L1 preadipocytes differentiation. This broadens the understanding of the molecular mechanisms of polyphenols on modulating and controlling of metabolic dysregulation, particularly adipocyte differentiation, which is a significant risk factor associated with the development of cardiovascular disease.

Decoupled side chain and backbone dynamics for proton translocation - M2 of influenza A

The 97 amino acid bitopic membrane protein M2 of influenza A forms a tetrameric bundle in which two of the monomers are covalently linked via a cysteine bridge. In its tetrameric assembly the protein conducts protons across the viral envelope and within intracellular compartments during the infectivity cycle of the virus. A key residue in the translocation of the protons is His-37 which forms a planar tetrad in the configuration of the bundle accepting and translocating the incoming protons from the N terminal side, exterior of the virus, to the C terminal side, inside the virus. With experimentally available data from NMR spectroscopy of the transmembrane domains of the tetrameric M2 bundle classical MD simulations are conducted with the protein bundle in different protonation stages in respect to His-37. A full correlation analysis (FCA) of the data sets with the His-37 tetrad either in a fully four times unprotonated or protonated state, assumed to mimic high and low pH in vivo, respectively, in both cases reveal asymmetric backbone dynamics. His-37 side chain rotation dynamics is increased at full protonation of the tetrad compared to the dynamics in the fully unprotonated state. The data suggest that proton translocation can be achieved by decoupled side chain or backbone dynamics. Graphical abstract Visualization of the tetrameric bundle of the transmembrane domains of M2 of influenza A after 200 ns of MD simulations (upper left). The four histidine residues 37 are either not protonated as in M2⁰ or fully protonated is in M2⁴⁺. The asymmetric dynamics of the backbones are shown after a full correlation analysis (FCA) in blue (lower left). The rotational dynamics of the χ2 dihedral angles of the histidines in M2⁰ (upper right) are less than those in M2⁴⁺ (lower right).

Temperature dependence of protein-water interactions in a gated yeast aquaporin

Regulation of aquaporins is a key process of living organisms to counteract sudden osmotic changes. Aqy1, which is a water transporting aquaporin of the yeast Pichia pastoris, is suggested to be gated by chemo-mechanical stimuli as a protective regulatory-response against rapid freezing. Here, we tested the influence of temperature by determining the X-ray structure of Aqy1 at room temperature (RT) at 1.3 Å resolution, and by exploring the structural dynamics of Aqy1 during freezing through molecular dynamics simulations. At ambient temperature and in a lipid bilayer, Aqy1 adopts a closed conformation that is globally better described by the RT than by the low-temperature (LT) crystal structure. Locally, for the blocking-residue Tyr31 and the water molecules inside the pore, both LT and RT data sets are consistent with the positions observed in the simulations at room-temperature. Moreover, as the temperature was lowered, Tyr31 adopted a conformation that more effectively blocked the channel, and its motion was accompanied by a temperature-driven rearrangement of the water molecules inside the channel. We therefore speculate that temperature drives Aqy1 from a loosely- to a tightly-blocked state. This analysis provides high-resolution structural evidence of the influence of temperature on membrane-transport channels.

Identification of potent cholecystokinin-B receptor antagonists: synthesis, molecular modeling and anti-cancer activity against pancreatic cancer cells

Advanced malignant stages of pancreatic cancer have poor prognosis and very few treatment strategies are available. Pancreatic cancer is known to possess unique growth-related receptors that when activated, stimulate tumour proliferation. Gastrin and its related peptide cholecystokinin (CCK) are also significantly involved in the growth of this cancer type as well as other malignancies through activation of the cholecystokinin-B receptor (CCK-BR). New treatment strategies with CCK-BR antagonists are being suggested that suppress the growth promoting effects of gastrin. In this paper, we report the development of two series of quinazolinone derivatives incorporating hydrazinecarbothioamide (compounds 3a-g) and the hydrazino group (compounds 4a-e) as linkers for developing CCK-BR antagonists. The affinities of the compounds were determined using docking into the CCK-BR homology modeled structure. The compounds were tested for in vitro CCK-BR binding and gastric acid secretion in an isolated lumen-perfused mouse stomach assay. The compounds exhibited CCK-BR binding activity (IC₅₀) in the range of 0.2-975 nM and showed good gastric acid secretion inhibitory activity. Molecular modeling of the compounds was done and pharmacophore mapping results showed good prediction of in vitro activity which correlated well with the experimental antagonistic activity. The compounds were further tested for their cytotoxicity on CCK-BR expressing pancreatic cancer cells. The results of the study provided two potent CCK-BR antagonists which also possess good to moderate growth inhibitory activities against pancreatic cancer cells.

A Molecular Basis for Selective Antagonist Destabilization of Dopamine D3 Receptor Quaternary Organization

The dopamine D3 receptor (D3R) is a molecular target for both first-generation and several recently-developed antipsychotic agents. Following stable expression of this mEGFP-tagged receptor, Spatial Intensity Distribution Analysis indicated that a substantial proportion of the receptor was present within dimeric/oligomeric complexes and that increased expression levels of the receptor favored a greater dimer to monomer ratio. Addition of the antipsychotics, spiperone or haloperidol, resulted in re-organization of D3R quaternary structure to promote monomerization. This action was dependent on ligand concentration and reversed upon drug washout. By contrast, a number of other antagonists with high affinity at the D3R, did not alter the dimer/monomer ratio. Molecular dynamics simulations following docking of each of the ligands into a model of the D3R derived from the available atomic level structure, and comparisons to the receptor in the absence of ligand, were undertaken. They showed that, in contrast to the other antagonists, spiperone and haloperidol respectively increased the atomic distance between reference α carbon atoms of transmembrane domains IV and V and I and II, both of which provide key interfaces for D3R dimerization. These results offer a molecular explanation for the distinctive ability of spiperone and haloperidol to disrupt D3R dimerization.

Ligand chain length drives activation of lipid G protein-coupled receptors

Sphingosine-1-phosphate (S1P) is a lipid mediator that can activate five cell membrane G protein-coupled receptors (GPCRs) which carry a variety of essential functions and are promising drug targets. S1P is composed of a polar zwitterionic head-group and a hydrophobic alkyl chain. This implies an activation mechanism of its cognate receptor that must be significantly different from what is known for prototypical GPCRs (ie receptor to small hydrophilic ligands). Here we aim to identify the structural features responsible for S1P agonism by combining molecular dynamics simulations and functional assays using S1P analogs of different alkyl chain lengths. We propose that high affinity binding involves polar interactions between the lipid head-group and receptor side chains while activation is due to hydrophobic interactions between the lipid tail and residues in a distinct binding site. We observe that ligand efficacy is directly related to alkyl chain length but also varies with receptor subtypes in correlation with the size of this binding pocket. Integrating experimental and computational data, we propose an activation mechanism for the S1P receptors involving agonist-induced conformational events that are conserved throughout class A GPCRs.

Adsorption of the natural protein surfactant Rsn-2 onto liquid interfaces

To stabilize foams, droplets and films at liquid interfaces a range of protein biosurfactants have evolved in nature. Compared to synthetic surfactants, these combine surface activity with biocompatibility and low solution aggregation. One recently studied example is Rsn-2, a component of the foam nest of the frog Engystomops pustulosus, which has been predicted to undergo a clamshell-like opening transition at the air-water interface. Using atomistic molecular dynamics simulations and surface tension measurements we study the adsorption of Rsn-2 onto air-water and cyclohexane-water interfaces. The protein adsorbs readily at both interfaces, with adsorption mediated by the hydrophobic N-terminus. At the cyclohexane-water interface the clamshell opens, due to the favourable interaction between hydrophobic residues and cyclohexane molecules and the penetration of cyclohexane molecules into the protein core. Simulations of deletion mutants showed that removal of the N-terminus inhibits interfacial adsorption, which is consistent with the surface tension measurements. Deletion of the hydrophilic C-terminus also affects adsorption, suggesting that this plays a role in orienting the protein at the interface. The characterisation of the interfacial behaviour gives insight into the factors that control the interfacial adsorption of proteins, which may inform new applications of this and similar proteins in areas including drug delivery and food technology and may also be used in the design of synthetic molecules showing similar changes in conformation at interfaces.

Potential of Mean Force Calculations of Solute Permeation Across UT-B and AQP1: A Comparison between Molecular Dynamics and 3D-RISM

Membrane channels facilitate the efficient and selective flux of various solutes across biological membranes. A common approach to investigate the selectivity of a channel has been the calculation of potentials of mean force (PMFs) for solute permeation across the pore. PMFs have been frequently computed from molecular dynamics (MD) simulations, yet the three-dimensional reference interaction site model (3D-RISM) has been suggested as a computationally efficient alternative to MD. Whether the two methods yield comparable PMFs for solute permeation has remained unclear. In this study, we calculated potentials of mean force for water, ammonia, urea, molecular oxygen, and methanol across the urea transporter B (UT-B) and aquaporin-1 (AQP1), using 3D-RISM, as well as using MD simulations and umbrella sampling. To allow direct comparison between the PMFs from 3D-RISM and MD, we ensure that all PMFs refer to a well-defined reference area in the bulk or, equivalently, to a well-defined density of channels in the membrane. For PMFs of water permeation, we found reasonable agreement between the two methods, with differences of ≲3 kJ mol-1. In contrast, we found stark discrepancies for the PMFs for all other solutes. Additional calculations confirm that discrepancies between MD and 3D-RISM are not explained by the choice for the closure relation, the definition the reaction coordinate (center of mass-based versus atomic site-based), details of the molecule force field, or fluctuations of the protein. Comparison of the PMFs suggests that 3D-RISM may underestimate effects from hydrophobic solute-channel interactions, thereby, for instance, missing the urea binding sites in UT-B. Furthermore, we speculate that the orientational averages inherent to 3D-RISM might lead to discrepancies in the narrow channel lumen. These findings suggest that current 3D-RISM solvers provide reasonable estimates for the PMF for water permeation, but that they are not suitable to study the selectivity of membrane channels with respect to uncharged nonwater solutes.

A new serotonin 5-HT6 receptor antagonist with procognitive activity - Importance of a halogen bond interaction to stabilize the binding

Serotonin 5-HT6 receptor has been proposed as a promising therapeutic target for cognition enhancement though the development of new antagonists is still needed to validate these molecules as a drug class for the treatment of Alzheimer's disease and other pathologies associated with memory deficiency. As part of our efforts to target the 5-HT6 receptor, new benzimidazole-based compounds have been designed and synthesized. Site-directed mutagenesis and homology models show the importance of a halogen bond interaction between a chlorine atom of the new class of 5-HT6 receptor antagonists identified herein and a backbone carbonyl group in transmembrane domain 4. In vitro pharmacological characterization of 5-HT6 receptor antagonist 7 indicates high affinity and selectivity over a panel of receptors including 5-HT2B subtype and hERG channel, which suggests no major cardiac issues. Compound 7 exhibited in vivo procognitive activity (1 mg/kg, ip) in the novel object recognition task as a model of memory deficit.

Dynamics of the EAG1 K+ channel selectivity filter assessed by molecular dynamics simulations

EAG1 channels belong to the KCNH family of voltage gated potassium channels. They are expressed in several brain regions and increased expression is linked to certain cancer types. Recent cryo-EM structure determination finally revealed the structure of these channels in atomic detail, allowing computational investigations. In this study, we performed molecular dynamics simulations to investigate the ion binding sites and the dynamical behavior of the selectivity filter. Our simulations suggest that sites S2 and S4 form stable ion binding sites, while ions placed at sites S1 and S3 rapidly switched to sites S2 and S4. Further, ions tended to dissociate away from S0 within less than 20 ns, due to increased filter flexibility. This was followed by water influx from the extracellular side, leading to a widening of the filter in this region, and likely non-conductive filter configurations. Simulations with the inactivation-enhancing mutant Y464A or Na⁺ ions lead to trapped water molecules behind the SF, suggesting that these simulations captured early conformational changes linked to C-type inactivation.

Molecular Insight into Affinities of Gallated and Nongallated Proanthocyanidins Dimers to Lipid Bilayers

Experimental studies have proved the beneficial effects of proanthocyanidins (Pas) relating to interaction with the cell membrane. But the detailed mechanisms and structure-function relationship was unclear. In present study, molecular dynamics (MD) simulations were used to study the interactions of four PA dimers with a lipid bilayer composed of 1:1 mixed 1-palmitoyl-2-oleoyl phosphatidylcholine (POPC) and 1-palmitoyl-2-oleoyl phosphatidylethanolamine (POPE). The results showed that the gallated PA dimers had much higher affinities to the bilayer with lower binding free energies compared with nongallated PA dimers. The gallated PA dimers penetrated deeper into the bilayer and formed more hydrogen bonds (H-bonds) with bilayer oxygen atoms, especially the deeper oxygen atoms of the lipids simultaneously, thus inducing stronger lateral expansion of the membrane and lipid tails disorder. The present results provided molecular insights into the interactions between PA dimers and bio-membranes and agreed with our experimental results well. These molecular interactions helped to elucidate the structure-function relationship of the PA dimers and provided a foundation for a better understanding of the underlying mechanisms of the bioactivities of PA oligomers.

Discovery of 8-Trifluoromethyl-3-cyclopropylmethyl-7-[(4-(2,4-difluorophenyl)-1-piperazinyl)methyl]-1,2,4-triazolo[4,3-a]pyridine (JNJ-46356479), a Selective and Orally Bioavailable mGlu2 Receptor Positive Allosteric Modulator (PAM)

Positive allosteric modulators of the metabotropic glutamate 2 receptor have generated great interest in the past decade. There is mounting evidence of their potential as therapeutic agents in the treatment of multiple central nervous system disorders. We have previously reported substantial efforts leading to potent and selective mGlu2 PAMs. However, finding compounds with the optimal combination of in vitro potency and good druglike properties has remained elusive, in part because of the hydrophobic nature of the allosteric binding site. Herein, we report on the lead optimization process to overcome the poor solubility inherent to the advanced lead 6. Initial prototypes already showed significant improvements in solubility while retaining good functional activity but displayed new liabilities associated with metabolism and hERG inhibition. Subsequent subtle modifications efficiently addressed those issues leading to the identification of compound 27 (JNJ-46356479). This new lead represents a more balanced profile that offers a significant improvement on the druglike attributes compared to previously reported leads.

Structural basis of control of inward rectifier Kir2 channel gating by bulk anionic phospholipids

Phospholipids are required to bind to two distinct sites on the inward rectifier potassium channel for maximal efficacy. Lee et al. show that a membrane-associating tryptophan residue in the second site can mimic the effect of phospholipid binding and cause a conformational change to reveal the primary binding site.

Inward rectifier potassium (Kir) channel activity is controlled by plasma membrane lipids. Phosphatidylinositol-4,5-bisphosphate (PIP₂) binding to a primary site is required for opening of classic inward rectifier Kir2.1 and Kir2.2 channels, but interaction of bulk anionic phospholipid (PL⁻) with a distinct second site is required for high PIP₂ sensitivity. Here we show that introduction of a lipid-partitioning tryptophan at the second site (K62W) generates high PIP₂ sensitivity, even in the absence of PL⁻. Furthermore, high-resolution x-ray crystal structures of Kir2.2[K62W], with or without added PIP₂ (2.8- and 2.0-Å resolution, respectively), reveal tight tethering of the C-terminal domain (CTD) to the transmembrane domain (TMD) in each condition. Our results suggest a refined model for phospholipid gating in which PL⁻ binding at the second site pulls the CTD toward the membrane, inducing the formation of the high-affinity primary PIP₂ site and explaining the positive allostery between PL⁻ binding and PIP₂ sensitivity.

Release of Entropic Spring Reveals Conformational Coupling Mechanism in the ABC Transporter BtuCD-F

Substrate translocation by ATP-binding cassette (ABC) transporters involves coupling of ATP binding and hydrolysis in the nucleotide-binding domains (NBDs) to conformational changes in the transmembrane domains. We used molecular dynamics simulations to investigate the atomic-level mechanism of conformational coupling in the ABC transporter BtuCD-F, which imports vitamin B12 across the inner membrane of Escherichia coli. Our simulations show how an engineered disulfide bond across the NBD dimer interface reduces conformational fluctuations and hence configurational entropy. As a result, the disulfide bond is under substantial mechanical stress. Releasing this entropic spring, as is the case in the wild-type transporter, combined with analyzing the pairwise forces between individual residues, unravels the coupling mechanism. The identified pathways along which force is propagated from the NBDs via the coupling helix to the transmembrane domains are composed of highly conserved residues, underlining their functional relevance. This study not only reveals the details of conformational coupling in BtuCD-F, it also provides a promising approach to other long-range conformational couplings, e.g., in ABC exporters or other ATP-driven molecular machines.

Structural Mechanisms of Voltage Sensing in G Protein-Coupled Receptors

G-protein-coupled receptors (GPCRs) form the largest superfamily of membrane proteins and one-third of all drug targets in humans. A number of recent studies have reported evidence for substantial voltage regulation of GPCRs. However, the structural basis of GPCR voltage sensing has remained enigmatic. Here, we present atomistic simulations on the δ-opioid and M2 muscarinic receptors, which suggest a structural and mechanistic explanation for the observed voltage-induced functional effects. The simulations reveal that the position of an internal Na(+) ion, recently detected to bind to a highly conserved aqueous pocket in receptor crystal structures, strongly responds to voltage changes. The movements give rise to gating charges in excellent agreement with previous experimental recordings. Furthermore, free energy calculations show that these rearrangements of Na(+) can be induced by physiological membrane voltages. Due to its role in receptor function and signal bias, the repositioning of Na(+) has important general implications for signal transduction in GPCRs.

Asymmetric dynamics of ion channel forming proteins - Hepatitis C virus (HCV) p7 bundles

Protein p7 of hepatitis C virus (HCV) is a short 63 amino acid membrane protein which homo-oligomerises in the lipid membrane to form ion and proton conducting bundles. Two different genotypes (GTs) of p7, 1a and 5a, are used to simulate hexameric bundles of the protein embedded in a fully hydrated lipid bilayer during 400 ns molecular dynamics (MD) simulations. Whilst the bundle of GT 1a is based on a fully computational derived structure, the bundle of GT 5a is based on NMR spectroscopic data. Results of a full correlation analysis (FCA) reveal that albeit structural differences both bundles screen local minima during the simulation. The collective motion of the protein domains is asymmetric. No 'breathing-mode'-like dynamics is observed. The presence of divalent ions, such as Ca-ions affects the dynamics of especially solvent exposed parts of the protein, but leaves the asymmetric domain motion unaffected.

Quaternary structure of a G-protein-coupled receptor heterotetramer in complex with Gi and Gs

Background

G-protein-coupled receptors (GPCRs), in the form of monomers or homodimers that bind heterotrimeric G proteins, are fundamental in the transfer of extracellular stimuli to intracellular signaling pathways. Different GPCRs may also interact to form heteromers that are novel signaling units. Despite the exponential growth in the number of solved GPCR crystal structures, the structural properties of heteromers remain unknown.

Results

We used single-particle tracking experiments in cells expressing functional adenosine A₁-A_2A receptors fused to fluorescent proteins to show the loss of Brownian movement of the A₁ receptor in the presence of the A_2A receptor, and a preponderance of cell surface 2:2 receptor heteromers (dimer of dimers). Using computer modeling, aided by bioluminescence resonance energy transfer assays to monitor receptor homomerization and heteromerization and G-protein coupling, we predict the interacting interfaces and propose a quaternary structure of the GPCR tetramer in complex with two G proteins.

Conclusions

The combination of results points to a molecular architecture formed by a rhombus-shaped heterotetramer, which is bound to two different interacting heterotrimeric G proteins (G_i and G_s). These novel results constitute an important advance in understanding the molecular intricacies involved in GPCR function.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-016-0247-4) contains supplementary material, which is available to authorized users.

Markovian and Non-Markovian Modeling of Membrane Dynamics with Milestoning

We exploit atomically detailed simulations and the milestoning theory to extract coarse grained models of membrane kinetics and thermodynamics. Non-Markovian and Markovian theories for the phosphate group displacements are used to coarsely represent membrane motions. The construction of the two models makes it possible to examine their consistency and accuracy. The equilibrium and fluctuations of the phosphate groups along the normal to the membrane plane are estimated, and milestoning equations are constructed and solved. An optimal Markovian model is constructed that reproduces exactly the equilibrium and mean first passage time (MFPT) of the non-Markovian model. The equilibrium solution of both models is favorably compared to distributions obtained from straightforward molecular dynamics simulations. The picture for the kinetics is complex. Multiple local relaxation times of the mass density are illustrated emphasizing the non-Markovian characteristics of the process. In Markovian modeling, only a single relaxation time is assumed for a state. Mapping of particle dynamics to the dynamics of a field density offers a new way of coarse graining complex systems as membranes that may bridge between atomically detailed models and phenomenological descriptions of macroscopic membranes.

Voltage-Gated Sodium Channels: Mechanistic Insights From Atomistic Molecular Dynamics Simulations

The permeation of ions and other molecules across biological membranes is an inherent requirement of all cellular organisms. Ion channels, in particular, are responsible for the conduction of charged species, hence modulating the propagation of electrical signals. Despite the universal physiological implications of this property, the molecular functioning of ion channels remains ambiguous. The combination of atomistic structural data with computational methodologies, such as molecular dynamics (MD) simulations, is now considered routine to investigate structure-function relationships in biological systems. A fuller understanding of conduction, selectivity, and gating, therefore, is steadily emerging due to the applicability of these techniques to ion channels. However, because their structure is known at atomic resolution, studies have consistently been biased toward K(+) channels, thus the molecular determinants of ionic selectivity, activation, and drug blockage in Na(+) channels are often overlooked. The recent increase of available crystallographic data has eminently encouraged the investigation of voltage-gated sodium (NaV) channels via computational methods. Here, we present an overview of simulation studies that have contributed to our understanding of key principles that underlie ionic conduction and selectivity in Na(+) channels, in comparison to the K(+) channel analogs.

Efficient preparation and analysis of membrane and membrane protein systems

Molecular dynamics (MD) simulations have become a highly important technique to consider lipid membrane systems, and quite often they provide considerable added value to laboratory experiments. Rapid development of both software and hardware has enabled the increase of time and size scales reachable by MD simulations to match those attainable by several accurate experimental techniques. However, until recently, the quality and maturity of software tools available for building membrane models for simulations as well as analyzing the results of these simulations have seriously lagged behind. Here, we discuss the recent developments of such tools from the end-users' point of view. In particular, we review the software that can be employed to build lipid bilayers and other related structures with or without embedded membrane proteins to be employed in MD simulations. Additionally, we provide a brief critical insight into force fields and MD packages commonly used for membrane and membrane protein simulations. Finally, we list analysis tools that can be used to study the properties of membrane and membrane protein systems. In all these points we comment on the respective compatibility of the covered tools. We also share our opinion on the current state of the available software. We briefly discuss the most commonly employed tools and platforms on which new software can be built. We conclude the review by providing a few ideas and guidelines on how the development of tools can be further boosted to catch up with the rapid pace at which the field of membrane simulation progresses. This includes improving the compatibility between software tools and promoting the openness of the codes on which these applications rely. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Molecular mechanism of positive allosteric modulation of the metabotropic glutamate receptor 2 by JNJ-46281222

BACKGROUND AND PURPOSE

Allosteric modulation of the mGlu2 receptor is a potential strategy for treatment of various neurological and psychiatric disorders. Here, we describe the in vitro characterization of the mGlu2 positive allosteric modulator (PAM) JNJ-46281222 and its radiolabelled counterpart [(3) H]-JNJ-46281222. Using this novel tool, we also describe the allosteric effect of orthosteric glutamate binding and the presence of a bound G protein on PAM binding and use computational approaches to further investigate the binding mode.

EXPERIMENTAL APPROACH

We have used radioligand binding studies, functional assays, site-directed mutagenesis, homology modelling and molecular dynamics to study the binding of JNJ-46281222.

KEY RESULTS

JNJ-46281222 is an mGlu2 -selective, highly potent PAM with nanomolar affinity (KD = 1.7 nM). Binding of [(3) H]-JNJ-46281222 was increased by the presence of glutamate and greatly reduced by the presence of GTP, indicating the preference for a G protein bound state of the receptor for PAM binding. Its allosteric binding site was visualized and analysed by a computational docking and molecular dynamics study. The simulations revealed amino acid movements in regions expected to be important for activation. The binding mode was supported by [(3) H]-JNJ-46281222 binding experiments on mutant receptors.

CONCLUSION AND IMPLICATIONS

Our results obtained with JNJ-46281222 in unlabelled and tritiated form further contribute to our understanding of mGlu2 allosteric modulation. The computational simulations and mutagenesis provide a plausible binding mode with indications of how the ligand permits allosteric activation. This study is therefore of interest for mGlu2 and class C receptor drug discovery.

Probing the Huntingtin 1-17 membrane anchor on a phospholipid bilayer by using all-atom simulations

Mislocalization and aggregation of the huntingtin protein are related to Huntington's disease. Its first exon-more specifically the first 17 amino acids (Htt17)-is crucial for the physiological and pathological functions of huntingtin. It regulates huntingtin's activity through posttranslational modifications and serves as an anchor to membrane-containing organelles of the cell. Recently, structure and orientation of the Htt17 membrane anchor were determined using a combined solution and solid-state NMR approach. This prompted us to refine this model by investigating the dynamics and thermodynamics of this membrane anchor on a POPC bilayer using all-atom, explicit solvent molecular dynamics and Hamiltonian replica exchange. Our simulations are combined with various experimental measurements to generate a high-resolution atomistic model for the huntingtin Htt17 membrane anchor on a POPC bilayer. More precisely, we observe that the single α-helix structure is more stable in the phospholipid membrane than the NMR model obtained in the presence of dodecylphosphocholine detergent micelles. The resulting Htt17 monomer has its hydrophobic plane oriented parallel to the bilayer surface. Our results further unveil the key residues interacting with the membrane in terms of hydrogen bonds, salt-bridges, and nonpolar contributions. We also observe that Htt17 equilibrates at a well-defined insertion depth and that it perturbs the physical properties-order parameter, thickness, and area per lipid-of the bilayer in a manner that could favor its dimerization. Overall, our observations reinforce and refine the NMR measurements on the Htt17 membrane anchor segment of huntingtin that is of fundamental importance to its biological functions.

How transmembrane peptides insert and orientate in biomembranes: a combined experimental and simulation study

After the synthesis of transmembrane peptides/proteins (TMPs), their insertion into a lipid bilayer is a fundamental biophysical process.

The molecular mechanism of secondary sodium symporters elucidated through the lens of the computational microscope

Transport of molecules across cellular membranes is a key biological process for normal cell function. In this review we describe current state-of-the-art knowledge on molecular mechanism of secondary active transporters obtained by molecular simulations studies.

Assembly of the MHC I peptide-loading complex determined by a conserved ionic lock-switch

Salt bridges in lipid bilayers play a decisive role in the dynamic assembly and downstream signaling of the natural killer and T-cell receptors. Here, we describe the identification of an inter-subunit salt bridge in the membrane within yet another key component of the immune system, the peptide-loading complex (PLC). The PLC regulates cell surface presentation of self-antigens and antigenic peptides via molecules of the major histocompatibility complex class I. We demonstrate that a single salt bridge in the membrane between the transporter associated with antigen processing TAP and the MHC I-specific chaperone tapasin is essential for the assembly of the PLC and for efficient MHC I antigen presentation. Molecular modeling and all-atom molecular dynamics simulations suggest an ionic lock-switch mechanism for the binding of TAP to tapasin, in which an unfavorable uncompensated charge in the ER-membrane is prevented through complex formation. Our findings not only deepen the understanding of the interaction network within the PLC, but also provide evidence for a general interaction principle of dynamic multiprotein membrane complexes in immunity.

Effects of Lys to Glu mutations in GsMTx4 on membrane binding, peptide orientation, and self-association propensity, as analyzed by molecular dynamics simulations

GsMTx4, a gating modifier peptide acting on cationic mechanosensitive channels, has a positive charge (+5e) due to six Lys residues. The peptide does not have a stereospecific binding site on the channel but acts from the boundary lipids within a Debye length of the pore probably by changing local stress. To gain insight into how these Lys residues interact with membranes, we performed molecular dynamics simulations of Lys to Glu mutants in parallel with our experimental work. In silico, K15E had higher affinity for 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine bilayers than wild-type (WT) peptide or any other mutant tested, and showed deeper penetration than WT, a finding consistent with the experimental data. Experimentally, the inhibitory activities of K15E and K25E were most compromised, whereas K8E and K28E inhibitory activities remained similar to WT peptide. Binding of WT in an interfacial mode did not influence membrane thickness. With interfacial binding, the direction of the dipole moments of K15E and K25E was predicted to differ from WT, whereas those of K8E and K28E oriented similarly to that of WT. These results support a model in which binding of GsMTx4 to the membrane acts like an immersible wedge that serves as a membrane expansion buffer reducing local stress and thus inhibiting channel activity. In simulations, membrane-bound WT attracted other WT peptides to form aggregates. This may account for the positive cooperativity observed in the ion channel experiments. The Lys residues seem to fine-tune the depth of membrane binding, the tilt angle, and the dipole moments.

Biomembranes in atomistic and coarse-grained simulations

The architecture of biological membranes is tightly coupled to the localization, organization, and function of membrane proteins. The organelle-specific distribution of lipids allows for the formation of functional microdomains (also called rafts) that facilitate the segregation and aggregation of membrane proteins and thus shape their function. Molecular dynamics simulations enable to directly access the formation, structure, and dynamics of membrane microdomains at the molecular scale and the specific interactions among lipids and proteins on timescales from picoseconds to microseconds. This review focuses on the latest developments of biomembrane force fields for both atomistic and coarse-grained molecular dynamics (MD) simulations, and the different levels of coarsening of biomolecular structures. It also briefly introduces scale-bridging methods applicable to biomembrane studies, and highlights selected recent applications.

Modulation of the Interaction between a Peptide Ligand and a G Protein-Coupled Receptor by Halogen Atoms

Systematic halogenation of two native opioid peptides has shown that halogen atoms can modulate peptide-receptor interactions in different manners. First, halogens may produce a steric hindrance that reduces the binding of the peptide to the receptor. Second, chlorine, bromine, or iodine may improve peptide binding if their positive σ-hole forms a halogen bond interaction with negatively charged atoms of the protein. Lastly, the negative electrostatic potential of fluorine can interact with positively charged atoms of the protein to improve peptide binding.

Structure and Dynamics of Phospholipid Nanodiscs from All-Atom and Coarse-Grained Simulations

We investigated structural and dynamical properties of nanodiscs comprising dimyristoylphosphatidylcholine (DMPC) lipids and major scaffold protein MSP1Δ(1-22) from human apolipoprotein A-1 using combined all-atom and coarse-grained (CG) molecular dynamics (MD) simulations. The computational efficiency of the Martini-CG force field enables the spontaneous self-assembly of lipids and scaffold proteins into stable nanodisc structures on time scales up to tens of microseconds. Subsequent all-atom and CG-MD simulations reveal that the lipids in the nanodisc have lower configurational entropy and higher acyl tail order than in a lamellar bilayer phase. These altered average properties arise from rather differential behavior of lipids, depending on their location in the nanodisc. Since the scaffold proteins exert constrictive forces from the outer rim of the disc toward its center, lipids at the center of the nanodisc are highly ordered, whereas annular lipids that are in contact with the MSP proteins are remarkably disordered due to perturbed packing. Although specific differences between all-atom and CG simulations are also evident, the results obtained at both levels of resolution are in overall good agreement with each other and provide atomic level interpretations of recent experiments. Thus, the present study highlights the applicability of multiscale simulation approaches for nanodisc systems and opens the way for future applications, including the study of nanodisc-embedded membrane proteins.