Explicit Lipid Research Articles

Low serum total cholesterol levels predict inferior prognosis of patients with POEMS syndrome

Low serum cholesterol levels are associated with increased tumor morbidity and mortality. However, the relationship between serum lipid profile and POEMS syndrome (polyneuropathy, organomegaly, endocrinopathy, M-protein, skin changes) is still unclear. The aim of our study was to clarify the importance of the serum lipid profile in predicting the severity and prognosis of patients with POEMS syndrome. Forty-three patients with newly diagnosed POEMS syndrome admitted to the Department of Hematology of Jiangsu Provincial People's Hospital between August 2013 and February 2023 were selected. They had explicit serum lipid profiles. There were 27 males and 16 females with a median age of 54 years (range, 28–77 years). Survival curves were plotted using the Kaplan–Meier method, and comparisons between the two groups were performed using the log-rank test. The Cox proportional-hazards model examined risk factors associated with the prognosis of POEMS syndrome. Receiver-operator characteristic (ROC) curves assessed the predictive accuracy. 23 (53.5%) patients had low total cholesterol (TC) levels. Low levels of TC were concerned with unfavorable progression-free survival (PFS) (p = 0.007) and overall survival (OS) (p = 0.004), and at the same time, the low circulating TC concentration was an independent risk factor for PFS (p = 0.020) and OS (p = 0.011). Low TC values could improve the risk stratification, especially in high-risk patients. In conclusion, low serum TC levels may predict inferior prognosis in patients with POEMS syndrome; in future clinical application, low TC may be a reliable indicator of prognosis.

Deciphering insights into the binding mechanism and plasticity of Telacebec with M. tuberculosis cytochrome bcc-aa3 supercomplex through an unbiased molecular dynamics simulation, free-energy analysis, and DFT study

The cytochrome bcc-aa3 supercomplex, a key component in the electron transport chain pathway involved in bacterial energy production and homeostasis, is a clinically validated target for tuberculosis (TB), leading to Telacebec (Q203). Telacebec is a potent candidate drug under Phase II clinical development for the treatment of drug-sensitive and drug-resistant TB. Recently, the cryo-electron microscopy structure of this supercomplex from Mycobacterium tuberculosis (Mtb) complexed with Q203 was resolved at 6.9 Å resolution (PDB ID: 7E1W). To understand the binding site (QP site) flexibility and Q203’s stability at the QP site of the Mtb cytochrome bcc complex, we conducted molecular dynamics (MD) simulation and free energy analysis on this complex in an explicit hydrated lipid bilayer environment for 500 ns. Through this study, the persistence of a range of direct and indirect interactions was observed over the course of the simulation. The significance of the interactions with His375, Tyr161, Ala178, Ala179, Ile183, His355, Leu356, and Thr313 is underlined. Electrostatic energy was the primary source of the net binding free energy, regardless of the important interacting residues. The overall binding free energy for Q203 was −112.84 ± 7.73 kcal/mol, of which the electrostatic and lipophilic energy contributions were −116.31 ± 1.14 and −21.32 ± 2.35 kcal/mol, respectively. Meanwhile, DFT calculations were utilized to elucidate Q203’s molecular properties. Overall, this study deciphers key insights into the cytochrome bcc-aa3 supercomplex with Q203 on the ground of molecular mechanics and quantum mechanics that may facilitate structure-based drug design and optimization for the discovery of the next-generation antitubercular drug(s). Communicated by Ramaswamy H. Sarma

Structural Modeling of Cytokine-Receptor-JAK2 Signaling Complexes Using AlphaFold Multimer.

Homodimeric class 1 cytokine receptors include the erythropoietin (EPOR), thrombopoietin (TPOR), granulocyte colony-stimulating factor 3 (CSF3R), growth hormone (GHR), and prolactin receptors (PRLR). These cell-surface single-pass transmembrane (TM) glycoproteins regulate cell growth, proliferation, and differentiation and induce oncogenesis. An active TM signaling complex consists of a receptor homodimer, one or two ligands bound to the receptor extracellular domains, and two molecules of Janus Kinase 2 (JAK2) constitutively associated with the receptor intracellular domains. Although crystal structures of soluble extracellular domains with ligands have been obtained for all of the receptors except TPOR, little is known about the structure and dynamics of the complete TM complexes that activate the downstream JAK-STAT signaling pathway. Three-dimensional models of five human receptor complexes with cytokines and JAK2 were generated here by using AlphaFold Multimer. Given the large size of the complexes (from 3220 to 4074 residues), the modeling required a stepwise assembly from smaller parts, with selection and validation of the models through comparisons with published experimental data. The modeling of active and inactive complexes supports a general activation mechanism that involves ligand binding to a monomeric receptor followed by receptor dimerization and rotational movement of the receptor TM α-helices, causing proximity, dimerization, and activation of associated JAK2 subunits. The binding mode of two eltrombopag molecules to the TM α-helices of the active TPOR dimer was proposed. The models also help elucidate the molecular basis of oncogenic mutations that may involve a noncanonical activation route. Models equilibrated in explicit lipids of the plasma membrane are publicly available.

OPRLM: Orientations of proteins in realistic lipid membranes.

A new web tool and a database, OPRLM (Orientations of Proteins in Realistic Lipid Membranes), has been created based on the OPM database and the broadly used GHARMM-GUI toolset. The database provides a complete set of experimental 3D structures of integral and peripheral membrane proteins from the Protein Data Bank that are positioned in planar and curved lipid bilayers using implicit and explicit membrane representations. Membrane systems vary from single-component artificial lipid bilayers to multi-component bilayers with complex lipid composition corresponding to 18 biomembranes with either asymmetric or symmetric distribution of lipid types between two leaflets: 14 “average” eukaryotic cellular and organelle membranes (plasma membranes of mammals, plants, and fungi, endoplasmic reticulum and apparatus Golgi membranes of mammals and fungi, mammalian membranes of endosomes and lysosomes, mitochondrial outer and inner membranes, plant vacuole membranes, thylakoid membranes of plants and cyanobacteria), the outer and the inner membranes of Gram-negative bacteria (E. coli), and “average” cell membranes of Gram-positive bacteria and archaebacteria. All these membrane systems are accessible for public use in CHARMM-GUI Archive (https://www.charmm-gui.org/docs/archive/biomembrane). The web tool allows users of any level of expertise an easy setup of protein-membrane systems for atomic molecular dynamics (MD) simulations in pre-assembled realistic lipid membranes formed by mixtures of explicit lipids. OPRLM is publicly available at https://oprlm.org/.

Understanding the conformational changes in the influenza B M2 ion channel at various protonation states

The characterization of influenza (A/B M2) ion channels is very important as they are potential binding sites for the drugs. We report the all-atom molecular dynamics study of the influenza B M2 ion channel in the presence of explicit solvent and lipid bilayers using the high resolution solid-state NMR structures. The importance of the various protonation states of histidine in the activation of the ion channel is discussed. The conformational changes at the closed and the open structures clearly show that the increase in tilt angle is necessary for the activation of the ion channel. Additionally, the free energy surfaces of the eight systems show the importance of the protonation state of the histidine residues in the activation of the influenza B M2 ion channel. The protonation of the histidine residues increases the tilt angle and the intra-helix distance which is evident from the superimposition of the structures corresponding to the maxima and the minima in the free energy landscape. The findings imply differences in the singly protonated and double protonated conformational states of BM2 ion channel and provide insights to help further studies of these ion channels as the drug targets for the influenza virus.

Assessing hERG1 Blockade from Bayesian Machine-Learning-Optimized Site Identification by Ligand Competitive Saturation Simulations.

Drug-induced cardiotoxicity is a potentially lethal and yet one of the most common side effects with the drugs in clinical use. Most of the drug-induced cardiotoxicity is associated with an off-target pharmacological blockade of K+ currents carried out by the cardiac Human-Ether-a-go-go-Related (hERG1) potassium channel. There is a compulsory preclinical stage safety assessment for the hERG1 blockade for all classes of drugs, which adds substantially to the cost of drug development. The availability of a high-resolution cryogenic electron microscopy (cryo-EM) structure for the channel in its open/depolarized state solved in 2017 enabled the application of molecular modeling for rapid assessment of drug blockade by molecular docking and simulation techniques. More importantly, if successful, in silico methods may allow a path to lead-compound salvaging by mapping out key block determinants. Here, we report the blind application of the site identification by the ligand competitive saturation (SILCS) protocol to map out druggable/regulatory hotspots in the hERG1 channel available for blockers and activators. The SILCS simulations use small solutes representative of common functional groups to sample the chemical space for the entire protein and its environment using all-atom simulations. The resulting chemical maps, FragMaps, explicitly account for receptor flexibility, protein-fragment interactions, and fragment desolvation penalty allowing for rapid ranking of potential ligands as blockers or nonblockers of hERG1. To illustrate the power of the approach, SILCS was applied to a test set of 55 blockers with diverse chemical scaffolds and pIC50 values measured under uniform conditions. The original SILCS model was based on the all-atom modeling of the hERG1 channel in an explicit lipid bilayer and was further augmented with a Bayesian-optimization/machine-learning (BML) stage employing an independent literature-derived training set of 163 molecules. BML approach was used to determine weighting factors for the FragMaps contributions to the scoring function. pIC50 predictions from the combined SILCS/BML approach to the 55 blockers showed a Pearson correlation (PC) coefficient of >0.535 relative to the experimental data. SILCS/BML model was shown to yield substantially improved performance as compared to commonly used rigid and flexible molecular docking methods for a well-established cohort of hERG1 blockers, where no correlation with experimental data was recorded. SILCS/BML results also suggest that a proper weighting of protonation states of common blockers present at physiological pH is essential for accurate predictions of blocker potency. The precalculated and optimized SILCS FragMaps can now be used for the rapid screening of small molecules for their cardiotoxic potential as well as for exploring alternative binding pockets in the hERG1 channel with applications to the rational design of activators.

Atomistic Structure and Dynamics of the Ca2+-ATPase Bound to Phosphorylated Phospholamban.

Sarcoplasmic reticulum Ca2+-ATPase (SERCA) and phospholamban (PLB) are essential components of the cardiac Ca2+ transport machinery. PLB phosphorylation at residue Ser16 (pSer16) enhances SERCA activity in the heart via an unknown structural mechanism. Here, we report a fully atomistic model of SERCA bound to phosphorylated PLB and study its structural dynamics on the microsecond time scale using all-atom molecular dynamics simulations in an explicit lipid bilayer and water environment. The unstructured N-terminal phosphorylation domain of PLB samples different orientations and covers a broad area of the cytosolic domain of SERCA but forms a stable complex mediated by pSer16 interactions with a binding site formed by SERCA residues Arg324/Lys328. PLB phosphorylation does not affect the interaction between the transmembrane regions of the two proteins; however, pSer16 stabilizes a disordered structure of the N-terminal phosphorylation domain that releases key inhibitory contacts between SERCA and PLB. We found that PLB phosphorylation is sufficient to guide the structural transitions of the cytosolic headpiece that are required to produce a competent structure of SERCA. We conclude that PLB phosphorylation serves as an allosteric molecular switch that releases inhibitory contacts and strings together the catalytic elements required for SERCA activation. This atomistic model represents a vivid atomic-resolution visualization of SERCA bound to phosphorylated PLB and provides previously inaccessible insights into the structural mechanism by which PLB phosphorylation releases SERCA inhibition in the heart.

An implicit lipid model for efficient reaction-diffusion simulations of protein binding to surfaces of arbitrary topology.

Localization of proteins to a membrane is an essential step in a broad range of biological processes such as signaling, virion formation, and clathrin-mediated endocytosis. The strength and specificity of proteins binding to a membrane depend on the lipid composition. Single-particle reaction-diffusion methods offer a powerful tool for capturing lipid-specific binding to membrane surfaces by treating lipids explicitly as individual diffusible binding sites. However, modeling lipid particle populations is expensive. Here, we present an algorithm for reversible binding of proteins to continuum surfaces with implicit lipids, providing dramatic speed-ups to many body simulations. Our algorithm can be readily integrated into most reaction-diffusion software packages. We characterize changes to kinetics that emerge from explicit vs implicit lipids as well as surface adsorption models, showing excellent agreement between our method and the full explicit lipid model. Compared to models of surface adsorption, which couple together binding affinity and lipid concentration, our implicit lipid model decouples them to provide more flexibility for controlling surface binding properties and lipid inhomogeneity, thus reproducing binding kinetics and equilibria. Crucially, we demonstrate our method's application to membranes of arbitrary curvature and topology, modeled via a subdivision limit surface, again showing excellent agreement with explicit lipid simulations. Unlike adsorption models, our method retains the ability to bind lipids after proteins are localized to the surface (through, e.g., a protein-protein interaction), which can greatly increase the stability of multiprotein complexes on the surface. Our method will enable efficient cell-scale simulations involving proteins localizing to realistic membrane models, which is a critical step for predictive modeling and quantification of in vitro and in vivo dynamics.

Familial L723P Mutation Can Shift the Distribution between the Alternative APP Transmembrane Domain Cleavage Cascades by Local Unfolding of the Ε-Cleavage Site Suggesting a Straightforward Mechanism of Alzheimer's Disease Pathogenesis.

Alzheimer's disease is an age-related pathology associated with accumulation of amyloid-β peptides, products of enzymatic cleavage of amyloid-β precursor protein (APP) by secretases. Several familial mutations causing early onset of the disease have been identified in the APP transmembrane (TM) domain. The mutations influence production of amyloid-β, but the molecular mechanisms of this effect are unclear. The "Australian" (L723P) mutation located in the C-termini of APP TM domain is associated with autosomal-dominant, early onset Alzheimer's disease. Herein, we describe the impact of familial L723P mutation on the structural-dynamic behavior of APP TM domain studied by high-resolution NMR in membrane-mimicking micelles and augmented by molecular dynamics simulations in explicit lipid bilayer. We found L723P mutation to cause local unfolding of the C-terminal turn of the APP TM domain helix and increase its accessibility to water required for cleavage of the protein backbone by γ-secretase in the ε-site, thus switching between alternative ("pathogenic" and "non-pathogenic") cleavage cascades. These findings suggest a straightforward mechanism of the pathogenesis associated with this mutation, and are of generic import for understanding the molecular-level events associated with APP sequential proteolysis resulting in accumulation of the pathogenic forms of amyloid-β. Moreover, age-related onset of Alzheimer's disease can be explained by a similar mechanism, where the effect of mutation is emulated by the impact of local environmental factors, such as oxidative stress and/or membrane lipid composition. Knowledge of the mechanisms regulating generation of amyloidogenic peptides of different lengths is essential for development of novel treatment strategies of the Alzheimer's disease.

Kinetic and thermodynamic insights into sodium ion translocation through the μ-opioid receptor from molecular dynamics and machine learning analysis

The differential modulation of agonist and antagonist binding to opioid receptors (ORs) by sodium (Na+) has been known for decades. To shed light on the molecular determinants, thermodynamics, and kinetics of Na+ translocation through the μ-OR (MOR), we used a multi-ensemble Markov model framework combining equilibrium and non-equilibrium atomistic molecular dynamics simulations of Na+ binding to MOR active or inactive crystal structures embedded in an explicit lipid bilayer. We identify an energetically favorable, continuous ion pathway through the MOR active conformation only, and provide, for the first time: i) estimates of the energy differences and required timescales of Na+ translocation in inactive and active MORs, ii) estimates of Na+-induced changes to agonist binding validated by radioligand measurements, and iii) testable hypotheses of molecular determinants and correlated motions involved in this translocation, which are likely to play a key role in MOR signaling.

Structure refinement of membrane proteins via molecular dynamics simulations.

A refinement protocol based on physics-based techniques established for water soluble proteins is tested for membrane protein structures. Initial structures were generated by homology modeling and sampled via molecular dynamics simulations in explicit lipid bilayer and aqueous solvent systems. Snapshots from the simulations were selected based on scoring with either knowledge-based or implicit membrane-based scoring functions and averaged to obtain refined models. The protocol resulted in consistent and significant refinement of the membrane protein structures similar to the performance of refinement methods for soluble proteins. Refinement success was similar between sampling in the presence of lipid bilayers and aqueous solvent but the presence of lipid bilayers may benefit the improvement of lipid-facing residues. Scoring with knowledge-based functions (DFIRE and RWplus) was found to be as good as scoring using implicit membrane-based scoring functions suggesting that differences in internal packing is more important than orientations relative to the membrane during the refinement of membrane protein homology models.

Structural Conservation and Effects of Alterations in T Cell Receptor Transmembrane Interfaces

T cell receptors (TCRs) are octameric assemblies of type-I membrane proteins in which a receptor heterodimer (αβ, δγ, or pre-Tαβ) is associated with three dimeric signaling modules (CD3δε, CD3γε, and ζζ) at the T cell or pre-T cell surface. In the human αβTCR, the α and β transmembrane (TM) domains form a specific structure that acts as a hub for assembly with the signaling modules inside the lipid bilayer. Conservation of key polar contacts across the C-terminal half of this TM interface suggests that the structure is a common feature of all TCR types. In this study, using molecular dynamics simulations in explicit lipid bilayers, we show that human δγ and pre-Tαβ TM domains also adopt stable αβ-like interfaces, yet each displays unique features that modulate the stability of the interaction and are related to sequences that are conserved within TCR types, but are distinct from the αβ sequences. We also performed simulations probing effects of previously reported mutations in the human αβ TM interface, and observed that the most disruptive mutations caused substantial departures from the wild-type TM structure and increased dynamics. These simulations show a strong correlation between structural instability, increased conformational variation, and the severity of structural defects in whole-TCR complexes measured in our previous biochemical assays. These results thus support the view that the stability of the core TM structure is a key determinant of TCR structural integrity and suggest that the interface has been evolutionarily optimized for different forms of TCRs.

The Flow of Information from Nucleus to Golgi Is Contingent upon Nuclear Membrane Synthesis and Protraction of the Ceramide-Containing Membrane to Endoplasmic Reticulum

The nucleus-initiated augmentation of ER membrane is reflected in a coordinated synthesis and intercalation of the explicit proteins and lipids required for the replacement, repair and function of the cell and its organelles. The direct connection between nucleus and the membranes containing labeled sphingosine (SphN) and ceramide (Cer) was affirmed by determining synthetic activity of serine palmitoyltransferase (SPT). The SPT and the newly synthesized serine-labeled lipid products were identified in the Outer- and Inner-Nuclear Membrane (ONM, INM) and ER. The pulse-chase experiments disclosed that the incorporation of radiolabeled lipids into both nuclear membranes declined upon their simultaneous increase in Endoplasmic Reticulum (ER). These results, and prior findings regarding metabolic transfer of nuclear membrane phosphoinositides to the outer leaflet of ER [Slomiany and Slomiany, Health, 2011, 3, 187-199], allowed us to reason that INM and ONM are not distinct entities, but uninterrupted continuum facing nucleosol and then cytosol when protracted into segment known as ER. Consequently, the identification of SPT and its products in the inner leaflet of nuclear and ER microsomes lent credence to the luminal presence of Cer in Golgi, luminal synthesis of glycosphingolipids (GSphLs), sphingomyelin (SM), and their delivery to the outer leaflet of apical and basolateral cell membrane, respectively. The findings presented in this communication provide further support to our concept that the factual intercalation of proteins and lipids into the cell membranes can only take place during their simultaneous synthesis that is guided by the nuclear and cytosolic processes enacted in nuclear-ER membrane continuum. At the nuclear stage, the signal-specific genes expression promotes active synthesis and intercalation of lipids into the organelles’ customized membrane that is protracted and articulated in ER in form of transport vesicles.

Death Receptor 5 Activation Is Energetically Coupled to Opening of the Transmembrane Domain Dimer

The precise mechanism by which binding of tumor necrosis factor ligands to the extracellular domain of their corresponding receptors transmits signals across the plasma membrane has remained elusive. Recent studies have proposed that activation of several tumor necrosis factor receptors, including Death Receptor 5, involves a scissorlike opening of the disulfide-linked transmembrane (TM) dimer. Using time-resolved fluorescence resonance energy transfer, we provide, to our knowledge, the first direct biophysical evidence that Death Receptor 5 TM-dimers open in response to ligand binding. Then, to probe the importance of the closed-to-open TM domain transition in the overall energetics of receptor activation, we designed point-mutants (alanine to phenylalanine) in the predicted, tightly packed TM domain dimer interface. We hypothesized that the bulky residues should destabilize the closed conformation and eliminate the ∼3 kcal/mol energy barrier to TM domain opening and the ∼2 kcal/mol energy difference between the closed and open states, thus oversensitizing the receptor. To test this, we used all-atom molecular dynamics simulations of the isolated TM domain in explicit lipid bilayers coupled to thermodynamic potential of mean force calculations. We showed that single point mutants at the interface altered the energy landscape as predicted, but were not enough to completely eliminate the barrier to opening. However, the computational model did predict that a double mutation at i, i+4 positions at the center of the TM domain dimer eliminates the barrier and stabilizes the open conformation relative to the closed. We tested these mutants in cells with time-resolved fluorescence resonance energy transfer and death assays, and show remarkable agreement with the calculations. The single mutants had a small effect on TM domain separation and cell death, whereas the double mutant significantly increased the TM domain separation and more than doubled the sensitivity of cells to ligand stimulation.

All atom NMDA receptor transmembrane domain model development and simulations in lipid bilayers and water.

N-methyl-d-aspartate receptors (NMDARs) are members of the ionotropic glutamate receptor family that mediate excitatory synaptic transmission in the central nervous system. The channels of NMDARs are permeable to Ca2+ but blocked by Mg2+, distinctive properties that underlie essential brain processes such as induction of synaptic plasticity. However, due to limited structural information about the NMDAR transmembrane ion channel forming domain, the mechanism of divalent cation permeation and block is understood poorly. In this paper we developed an atomistic model of the transmembrane domain (TMD) of NMDARs composed of GluN1 and GluN2A subunits (GluN1/2A receptors). The model was generated using (a) a homology model based on the structure of the NaK channel and a partially resolved structure of an AMPA receptor (AMPAR), and (b) a partially resolved X-ray structure of GluN1/2B NMDARs. Refinement and extensive Molecular Dynamics (MD) simulations of the NMDAR TMD model were performed in explicit lipid bilayer membrane and water. Targeted MD with simulated annealing was introduced to promote structure refinement. Putative positions of the Mg2+ and Ca2+ ions in the ion channel divalent cation binding site are proposed. Differences in the structural and dynamic behavior of the channel protein in the presence of Mg2+ or Ca2+ are analyzed. NMDAR protein conformational flexibility was similar with no ion bound to the divalent cation binding site and with Ca2+ bound, whereas Mg2+ binding reduced protein fluctuations. While bound at the binding site both ions retained their preferred ligand coordination numbers: 6 for Mg2+, and 7–8 for Ca2+. Four asparagine side chain oxygens, a back-bone oxygen, and a water molecule participated in binding a Mg2+ ion. The Ca2+ ion first coordination shell ligands typically included four to five side-chain oxygen atoms of the binding site asparagine residues, two water molecules and zero to two backbone oxygens of the GluN2B subunits. These results demonstrate the importance of high-resolution channel structures for elucidation of mechanisms of NMDAR permeation and block.

Temperature-Sensitive Gating of TRPV1 Channel as Probed by Atomistic Simulations

Heat-activated transient receptor potential channel TRPV1 is one of the most studied eukaryotic proteins involved in temperature sensation. Upon heating, it exhibits rapid reversible pore gating, which depolarizes neurons and generates action potentials. Underlying molecular details of such effects in the pore region of TRPV1 is of a crucial importance to control temperature responses of the organism. Despite the spatial structure of the channel in both open (O) and closed (C) states is known, microscopic nature of channel gating and mechanism of thermal sensitivity are still poorly understood. In this work, we used unrestrained atomistic molecular dynamics simulations of TRPV1 (without N- and C-terminal cytoplasmic domains) embedded into explicit lipid bilayer in its O- and C-states [1]. We found that the pore domain with its neighboring loops undergoes large temperature-dependent conformational transitions in an asymmetric way, when fragments only one monomer move with large amplitude, freeing the pore upon heating. Such an asymmetrical gating looks rather biologically relevant because it is faster and more reliable than traditionally proposed “iris-like” symmetric scheme of channel opening. Analysis of structural, dynamic, and hydrophobic organization of the pore domain revealed entropy growth upon TRPV1 gating, which is in line with current concepts of thermal sensitivity.Acknowledgements: Authors thank Russian Science Foundation (14-14-00871), Russian Foundation for Basic Research (16-04-00578), RAS MCB Program. Access to computational facilities of the Supercomputer Center “Polytechnical” at the St. Petersburg Polytecnic University is gratefully acknowledged.Reference:[1] A.O. Chugunov et al., Scientific Reports 6 (2016) 33112.

Role of trans to cis Transition in Viral Fusion Pore Dilation

Fusion of viral and host membranes is a key step during infection by membrane-enclosed viruses. The fusion pore plays a critical role, and must dilate to release the viral genome. Previous studies of fusion mediated by influenza A hemagglutinin (HA) revealed ∼2-5 nm pores that flickered before dilating to >10 nm. The mechanisms are unknown.Here we studied HA-mediated fusion pore dynamics using a novel single-pore assay, combined with computational simulations accessing extraordinarily long ms-s timescales. We measured pores between HA-expressing fibroblasts and bilayer nanodiscs. From pore currents we infer pore size with millisecond time resolution. Unlike previous in vitro studies, use of nanodiscs limited the membrane contact areas and maximum pore sizes, better mimicking the initial phases of virus-endosome fusion. With wild type (WT) HA, fusion pores flickered about a mean pore size ∼1 nm. By contrast, fusion pores formed by GPI-anchored HA nucleated at half the WT rate and were significantly larger.We developed radically coarse-grained, explicit lipid molecular dynamics simulations of the fusion pore reconstituted with post-fusion, trans HA hairpins. With WT HA, fusion pores were small, similar to experiment. Over time hairpins gradually converted from trans to cis, but contrary to a common view, cis hairpins accumulated on the “viral” membrane, not the pore waist, due to the low mobility HA transmembrane domains. With GPI-HA the anchoring lipids were far more mobile and the trans-cis transition much accelerated. Once most hairpins had converted to cis, because apposing membranes were released the fusion pore dilated significantly.Our results suggest pore dilation requires the trans-cis transition. We hypothesize that this transition is accelerated in GPI-HA by the more mobile lipid anchor, explaining the larger observed pores.

CHARMM-GUI MDFF/xMDFF Utilizer for Molecular Dynamics Flexible Fitting Simulations in Various Environments.

X-ray crystallography and cryo-electron microscopy are two popular methods for the structure determination of biological molecules. Atomic structures are derived through the fitting and refinement of an initial model into electron density maps constructed by both experiments. Two computational approaches, MDFF and xMDFF, have been developed to facilitate this process by integrating the experimental data with molecular dynamics simulation. However, the setup of an MDFF/xMDFF simulation requires knowledge of both experimental and computational methods, which is not straightforward for nonexpert users. In addition, sometimes it is desirable to include realistic environments, such as explicit solvent and lipid bilayers during the simulation, which poses another challenge even for expert users. To alleviate these difficulties, we have developed MDFF/xMDFF Utilizer in CHARMM-GUI that helps users to set up an MDFF/xMDFF simulation. The capability of MDFF/xMDFF Utilizer is greatly enhanced by integration with other CHARMM-GUI modules, including protein structure manipulation, a diverse set of lipid types, and all-atom CHARMM and coarse-grained PACE force fields. With this integration, various simulation environments are available for MDFF Utilizer (vacuum, implicit/explicit solvent, and bilayers) and xMDFF Utilizer (vacuum and solution). In this work, three examples are shown to demonstrate the usage of MDFF/xMDFF Utilizer.

Temperature-sensitive gating of TRPV1 channel as probed by atomistic simulations of its trans- and juxtamembrane domains

Heat-activated transient receptor potential channel TRPV1 is one of the most studied eukaryotic proteins involved in temperature sensation. Upon heating, it exhibits rapid reversible pore gating, which depolarizes neurons and generates action potentials. Underlying molecular details of such effects in the pore region of TRPV1 is of a crucial importance to control temperature responses of the organism. Despite the spatial structure of the channel in both open (O) and closed (C) states is known, microscopic nature of channel gating and mechanism of thermal sensitivity are still poorly understood. In this work, we used unrestrained atomistic molecular dynamics simulations of TRPV1 (without N- and C-terminal cytoplasmic domains) embedded into explicit lipid bilayer in its O- and C-states. We found that the pore domain with its neighboring loops undergoes large temperature-dependent conformational transitions in an asymmetric way, when fragments of only one monomer move with large amplitude, freeing the pore upon heating. Such an asymmetrical gating looks rather biologically relevant because it is faster and more reliable than traditionally proposed “iris-like” symmetric scheme of channel opening. Analysis of structural, dynamic, and hydrophobic organization of the pore domain revealed entropy growth upon TRPV1 gating, which is in line with current concepts of thermal sensitivity.

Mutational Analysis on Membrane Associated Transporter Protein (MATP) and Their Structural Consequences in Oculocutaeous Albinism Type 4 (OCA4)-A Molecular Dynamics Approach.

Oculocutaneous albinism type IV (OCA4) is an autosomal recessive inherited disorder which is characterized by reduced biosynthesis of melanin pigmentation in skin, hair, and eyes and caused by the genetic mutations in the membrane-associated transporter protein (MATP) encoded by SLC45A2 gene. The MATP protein consists of 530 amino acids which contains 12 putative transmembrane domains and plays an important role in pigmentation and probably functions as a membrane transporter in melanosomes. We scrutinized the most OCA4 disease-associated mutation and their structural consequences on SLC45A2 gene. To understand the atomic arrangement in 3D space, the native and mutant structures were modeled. Further the structural behavior of native and mutant MATP protein was investigated by molecular dynamics simulation (MDS) approach in explicit lipid and water background. We found Y317C as the most deleterious and disease-associated SNP on SLC45A2 gene. In MDS, mutations in MATP protein showed loss of stability and became more flexible, which alter its structural conformation and function. This phenomenon has indicated a significant role in inducing OCA4. Our study explored the understanding of molecular mechanism of MATP protein upon mutation at atomic level and further helps in the field of pharmacogenomics to develop a personalized medicine for OCA4 disorder. J. Cell. Biochem. 117: 2608-2619, 2016. © 2016 Wiley Periodicals, Inc.