In order to associate with the cytoplasmic leaflet of the plasma membrane, many cytosolic signalling proteins possess a distinct lipid binding domain as part of their overall fold. Here, a multiscale simulation approach has been used to investigate three membrane-binding proteins involved in cellular processes such as growth and proliferation. The pleckstrin homology (PH) domain from the general receptor for phosphoinositides 1 (GRP1-PH) binds phosphatidylinositol (3,4,5)-trisphosphate (PI(3,4,5)P₃) with high affinity and specificity. To investigate how this peripheral protein is able to locate its target lipid in the complex membrane environment, Brownian dynamics (BD) simulations were employed to explore association pathways for GRP1-PH binding to PI(3,4,5)P₃ embedded in membranes with different surface charge densities and distributions. The results indicated that non-PI(3,4,5)P₃ lipids can act as decoys to disrupt PI(3,4,5)P₃ binding, but that at approximately physiological anionic lipid concentrations steering towards PI(3,4,5)P₃ is actually enhanced. Atomistic molecular dynamics (MD) simulations revealed substantial membrane penetration of membrane-bound GRP1-PH, evident when non-equilibrium, steered MD simulations were used to forcibly dissociate the protein from the membrane surface. Atomistic and coarse grained (CG) MD simulations of the phosphatase and tensin homologue deleted on chromosome ten (PTEN) tumour suppressor, which also binds PI(3,4,5)P₃, detected numerous non-specific protein-lipid contacts and anionic lipid clustering around PTEN that can be modulated by selective in silico mutagenesis. These results suggested a dual recognition model of membrane binding, with non-specific membrane interactions complementing the protein-ligand interaction. Molecular docking and MD simulations were used to characterise the lipid binding properties of kindlin-1 PH. Simulations demonstrated that a dynamic salt bridge was responsible for controlling the accessibility of the binding site. Electrostatics calculations applied to a variety of PH domains suggested that their molecular dipole moments are typically aligned with their ligand binding sites, which has implications for steering and ligand electrostatic funnelling.

Computational methods, including homology modelling, in-silico dockings, and molecular dynamics simulations have been used to study the functional dynamics and interactions of K⁺ channels. Molecular models were built of the inwardly rectifying K⁺ channel Kir2.2, the bacterial homolog K⁺ channel KirBac3.1, and the twin pore (K2P) K⁺ channels TREK-1 and TRESK. To investigate the electrostatic energy profile of K⁺ permeating through these homology models, continuum electrostatic calculations were performed. The primary mechanism of KirBac3.1 gating is believed to involve an opening at the helix bundle crossing (HBC). However, simulations of Kir channels have not yet revealed opening at the HBC. Here, in simulations of the new KirBac3.1-S129R X-ray crystal structure, in which the HBC was trapped open by the S129R mutation in the inner pore-lining helix (TM2), the HBC was found to exhibit considerable mobility. In a simulation of the new KirBac3.1-S129R-S205L double mutant structure, if the S129R and the S205L mutations were converted back to the wild-type serine, the HBC would close faster than in the simulations of the KirBac3.1-S129R single mutant structure. The double mutant structure KirBac3.1-S129R-S205L therefore likely represents a higher-energy state than the single mutant KirBac3.1-S129R structure, and these simulations indicate a staged pathway of gating in KirBac channels. Molecular modelling and MD simulations of the Kir2.2 channel structure demonstrated that the HBC would tend to open if the C-linker between the transmembrane and cytoplasmic domain was modelled helical. The electrostatic energy barrier for K⁺ permeation at the helix bundle crossing was found to be sensitive to subtle structural changes in the C-linker. Charge neutralization or charge reversal of the PIP2-binding residue R186 on the C-linker decreased the electrostatic barrier for K⁺ permeation through the HBC, suggesting an electrostatic contribution to the PIP2-dependent gating mechanism. Multi-scale simulations determined the PIP2 binding site in Kir2.2, in good agreement with crystallographic predictions. A TREK-1 homology model was built, based on the TRAAK structure. Two PIP2 binding sites were found in this TREK-1 model, at the C-terminal end, in line with existing functional data, and between transmembrane helices TM2 and TM3. The TM2-TM3 site is in reasonably good agreement with electron density attributed to an acyl tail in a recently deposited TREK-2 structure.

L'objectif de ce travail est de comprendre la structure et les interactions des liquides ioniques au contact de surfaces métalliques à l’échelle moléculaire en ayant recours aux méthodes de dynamique moléculaire. Il s’agit également d’étudier l’impact de ces caractéristiques microscopiques sur les propriétés tribologiques du système. Les liquides ioniques choisis en tant qu‘huiles lubrifiantes potentielles présentent des propriétés biodégradables et des caractéristiques tribologiques appropriées. Ils reposent sur des cations alkylammonium combinés avec des anions alkylsulfonate et bistriflamide. Notre étude est structurée en quatre parties. Elle commence par l’analyse des liquides ioniques purs puis, des liquides ioniques confinés entre deux surfaces de fer à l’équilibre et sous cisaillement, et enfin, en présence d’eau. Les propriétés structurales et dynamiques des liquides ioniques sont étudiées à travers la fonction de distribution radiale et les coefficients d’auto-diffusion. L’organisation des charges ainsi que la formation de micro-domaines en solution sont étudiées conjointement au comportement diffusif des espèces ioniques. Un champ de forces atomique, basé sur des méthodes quantiques, a été développé pour modéliser les interactions entre les liquides ioniques et la surface métallique. Des calculs DFT ont été réalisés sur des fragments de liquides ioniques en interaction avec un cluster de fer en fonction de la distance et de leur orientation. Une fonction modélisant des interactions site-site a été ajustée aux valeurs d’énergies fragment–cluster calculées par DFT afin d’obtenir les paramètres du champ de forces. Finalement, la polarisation du métal par les ions a été prise en compte en utilisant un modèle de dipôles induits afin de reproduire l’énergie d’interaction entre les charges et la surface conductrice. Avec ce modèle d’interaction, les simulations de dynamique moléculaire ont permis d’étudier la structure de l’interface entre une surface de fer plane et différents liquides ioniques. Cette analyse s’est concentrée sur l’étude du positionnement des différentes espèces au niveau de la surface, sur l’orientation des chaines alkyles et sur les profils de densité de charge. Des simulations de dynamique moléculaire hors-équilibre de liquides ioniques en interaction avec des surfaces de fer ont été réalisées en utilisant le champ de forces développé précédemment. Un protocole de simulation, basé sur une définition locale de la pression, a été développé pour prédire de manière quantitative le coefficient de friction en fonction de la valeur de la charge et du taux de cisaillement. La dépendance de la friction avec la charge, la vitesse de cisaillement, la topologie de la surface et la taille de la chaine alkyle du liquide ionique a été étudiée. La variation des forces de friction s’explique par l’arrangement spécifique des ions et l’orientation des groupements du liquide ionique à proximité de la surface. Finalement, l’effet de la présence d’eau en petite quantité dans une solution de liquide ionique a aussi été étudié à l’équilibre et hors-équilibre. Un potentiel a été construit pour décrire les interactions entre l’eau et une surface de fer en utilisant la même approche que celle décrite précédemment. Des résultats préliminaires concernant la structure de l’interface liquide-métal et la valeur du coefficient de friction ont été présentés et comparés avec ceux obtenus pour les liquides ioniques purs. / The aim of the present work is to understand at the molecular level the structure and interactions of ionic liquids at metallic surfaces, using molecular dynamics simulations, and to investigate the impact that these microscopic features can have in the tribological properties of the system. The chosen ionic liquids as potential lubricant oils present suitable ecotoxic and biodegradable properties and appropriate tribological characteristics. They are based in alkylammonium cations combined with alkylsulfonate and bistriflamide anions. Our study is divided in four parts, starting from the analyses of pure ionic liquids solutions and evolving to systems of ionic liquids confined between surfaces of iron, at the equilibrium, under shear and also in the presence of water. Structural and dynamic properties of ionic liquids are investigated in terms of the site-site radial distribution functions and the self-diffusion coefficients. The presence of charge-ordering and the formation of micro-domains in solution are discussed, as well as the diffusive behavior of the ionic species. An atomistic force field for ionic liquids interacting with a metal surface was built based on quantum methods. Density functional calculations of alkylammonium cations, alkylsulfonate and bistriflamide anions interacting with a cluster of iron atoms are performed, at a series of distances and orientations. A site-site potential function was then adjusted to the DFT interactions energies, to obtain the force field parameters. Finally, the polarization of the metal by the ions was taken into account using induced dipoles to reproduce the interaction energy between charges and a conductor surface. Using this interaction model, molecular dynamics simulations were performed to study the structure of the interfacial layer of several ionic liquids at a flat iron surface, including analyses of the positional and orientational ordering of the ions near the surface, and charge density profiles. Non-equilibrium molecular dynamics simulations of ionic liquids interacting with iron surfaces were carried out using the specific set of interaction parameters developed previously. A procedure was developed for a quantitative prediction of the friction coefficient at different loads and shear rates, based in a definition of pressure measured locally. The dependence of friction on the load, shear velocity, surface topology and length of alkyl side chains in the ionic liquid was investigated. The changes in the frictional forces were explained in terms of the specific arrangements and orientations of groups forming the ionic liquid at the vicinity of the surface. Finally, the effect of the presence of water in a small quantity in an ionic liquid solution is also studied at equilibrium and non-equilibrium. An interaction potential was build that describes the interaction between water and an iron surface, using the same approach described previously. Preliminary results are presented on the structure at the metal–liquid interface and friction coefficient, and compared with the pure ionic liquids.

L’albumine sérique humaine (HSA) est une protéine connue pour ses propriétés de transport exceptionnelles et son contenu élevé en ponts disulfure. L’étude de sa dynamique conformationnelle représente un défi important dans la compréhension de ses fonctions physiologiques. Le but de notre travail a été d’étudier cette dynamique conformationnelle et de comprendre le rôle des ponts disulfure dans le maintien de la structure native de la protéine. Notre analyse est basée sur des simulations de dynamique moléculaire couplées à des analyses par composantes principales. Outre la validation de la méthode de simulation les résultats fournissent de nouveaux éclairages sur les principaux effets de la réduction des ponts disulfure dans les albumines sériques. Les processus de dépliement/repliement protéique ont été détaillés. La prédiction de la structure réduite d’équilibre a également fait l’objet d’une attention particulière. Une étude détaillée de la dynamique conformationnelle globale de la protéine ainsi que celle des deux sites principaux de complexation a été effectuée. D’éventuels effets allostériques entre ces deux sites ont été recherchés. Les résultats théoriques obtenus ont été discutés avec les données expérimentales disponibles / Human serum albumin (HSA) is a protein known for its exceptional transport properties and its high content of disulfide bridges. The study of the conformational dynamics represents a major challenge in the comprehension of its physiological functions. The aim of our work was to study the conformational dynamics and to understand the roleof disulfide bonds in the stability of the native protein structure. Our analysis is based on simulations of molecular dynamics coupled with principal component analysis. Beyond the validation of the simulation method, the results provide new insights on the main effects of the disulfide bonds reduction in serum albumins. Protein unfolding/refolding processes were detailed. A special attention is paid to the prediction of the reduced structure at the equilibrium. A detailed study of the global protein conformational dynamics as well as the two main binding sites were performed. Possible allosteric effects between these two sites were researched. The theoretical results have been discussed with the available experimental data

No description available.

Οι νικοτινικοί υποδοχείς της ακετυλοχολίνης (nAChRs) ανήκουν στην υπερ-οικογένεια των ιοντικών καναλιών που ενεργοποιούνται από τη δέσμευση ενός προσδέτη (LGICs) και αποτελούνται από πέντε ομόλογες υπομονάδες. Κάθε μονομερής υπομονάδα αποτελείται από μία Ν-τελική εξωκυττάρια περιοχή (ΕΚΠ), από τέσσερεις διαμεμβρανικές α-έλικες και από μία κυτταροπλασματική περιοχή. Στην ΕΚΠ βρίσκεται η χαρακτηριστική Cys-θηλιά της υπερ-οικογένειας, καθώς και οι θέσεις πρόσδεσης αγωνιστών και ανταγωνιστών του υποδοχέα. Οι γνώσεις μας γύρω από τη δομή των nAChRs προέρχονται κυρίως από κρυσταλλογραφικές δομές ομολόγων πρωτεϊνών δέσμευσης της ACh (AChBP) μαλακίων, από μια δομή του nAChR από ιχθείς του γένους Torpedo που προέρχεται από ηλεκτρονική μικροσκοπία, από την κρυσταλλογραφική δομή της α1-ΕΚΠ ποντικού σε σύμπλοκο με α-μπουγκαροτοξίνη (α-Btx) και από κρυσταλλογραφικές δομές δύο προκαρυωτικών LGICs. Παρά τη μεγάλη πρόοδο που πραγματοποιήθηκε με τα παραπάνω επιτεύγματα, ακόμη δεν έχει επιλυθεί πειραματικά η δομή ανθρώπινου υποδοχέα. Επίσης λίγα είναι γνωστά για την επίδραση της γλυκοζυλίωσης των ΕΚΠ στη λειτουργία του nAChR. Χρησιμοποιώντας ως εκμαγείο την κρυσταλλογραφική δομή του συμπλόκου α1-ΕΚΠ ποντικού/α-Btx δημιουργήθηκαν υπολογιστικά μοντέλα της ανθρώπινης α1-ΕΚΠ προσδεμένης στις τοξίνες α-μπουγκαροτοξίνη (α-Btx), α-κομπρατοξίνη (α-Cbtx), α-κωνοτοξίνη (α-Ctx) ImI και α-κωνοτοξίνη GI. Στα σύμπλοκα με α-Btx και α-Cbtx προστέθηκε η υδατανθρακική αλυσίδα, συνδεδεμένη με το κατάλοιπο Asn141, που συγκρυσταλλώθηκε μαζί με την α1-ΕΚΠ ποντικού. Για να μελετηθεί η δυναμική συμπεριφορά της αλληλεπίδρασης υποδοχέα-τοξίνης καθώς και η συνεισφορά των σακχάρων σε αυτήν πραγματοποιήθηκαν προσομοιώσεις Μοριακής Δυναμικής σε υδατικό περιβάλλον. Με τη χρήση υπολογιστικών εργαλείων για τη μελέτη των συμπλόκων προσδιορίστηκαν σε ατομικό επίπεδο οι αλληλεπιδράσεις που καθοδηγούν την πρόσδεση τοξινών στην α1-ΕΚΠ. Βρέθηκε ότι η υδατανθρακική αλυσίδα συμμετέχει δυναμικά στη δέσμευση της τοξίνης στον υποδοχέα. Τα σάκχαρα συγκλίνουν προς την προσδεμένη τοξίνη στηριζόμενα στα κατάλοιπα Ser187 και Trp184 της α1 υπομονάδας. Η τοξίνη επίσης μετακινείται φέρνοντας τη θηλιά Ι σε επαφή με τα σάκχαρα. Αναγνωρίστηκαν σημαντικές αλληλεπιδράσεις των σακχάρων με τα τοξινικά κατάλοιπα Thr6, Ser9, και Th15 της α-Btx και Thr6 και Pro7 της α-Cbtx. Επίσης επιβεβαιώθηκε η ύπαρξη μιας υδρόφιλης κοιλότητας στο εσωτερικό του υδρόφοβου πυρήνα της α1-ΕΚΠ, η οποία πιθανόν εμπλέκεται στο άνοιγμα του ιοντικού καναλιού του nAChR. Τα αποτελέσματα αυτά παρέχουν σημαντικά δεδομένα για την κατανόηση της επίδρασης της υδατανθρακικής αλυσίδας στη λειτουργία του υποδοχέα, η οποία μπορεί να αξιοποιηθεί στην αντιμετώπιση των πολλών παθολογικών καταστάσεων στις οποίες εμπλέκονται οι nAChRs. / Nicotinic acetylcholine receptors (nAChRs) belong to the superfamily of ligand-gated ion channels (LGICs). LGICs form homo- or hetero-pentamers of related subunits, and each of them consists of a N-terminal extracellular ligand-binding domain (ECD), four transmembrane α-helixes and an intracellular region. The characteristic Cys-loop of the superfamily is found in the ECD of each subunit. The ECD also contains binding sites for agonists and competitive antagonists. Our knowledge regarding the nAChR structure mainly derives from the X-ray crystal structures of the molluscan ACh-binding proteins (AChBPs), the electron microscopy structure of the Torpedo nAChR, the X-ray crystal structure of the mouse nAChR α1-ECD bound to α-bungarotoxin (α-Btx), and the X-ray crystal structures of two prokaryotic LGICs. Despite the progress made by these achievements, the determination of any human nAChR structure has not yet been accomplished. Furthermore, the effect of glycosylation on nAChR function has not yet been explored. Based on the crystal structure of the extracellular domain of the mouse nAChR α1 subunit bound to α-Btx we have generated in silico models of the human nAChR α1-ECD bound to the toxins α-bungarotoxin (α-Btx), α-cobratoxin (α-Cbtx), α-conotoxin (α-Ctx) ImI and α-conotoxin GI. In the case of the α1-ECD/α-Btx and α-Cbtx complexes, a Asn141-linked carbohydrate chain was modeled, its coordinates taken from the crystal structure of the mouse α1-ECD. To gain further insight into the structural role of glycosylation molecular dynamics (MD) simulations were carried out in explicit solvent so as to compare the conformational dynamics of the binding interface between nAChR α1 and the two toxins. The use of computational methods allowed the monitoring of the interactions that govern toxin binding. The MD simulations revealed the strengthening of the receptor-toxin interaction in the presence of the carbohydrate chain. A shift in the position of the sugars towards the bound toxin was observed. Residues Ser187 and Trp184 of nAChR act as critical anchor points for the stabilization of the sugar chain in a close position to the toxin. Toxin Finger I shifts closer to the mannoses, forming important toxin-sugar interactions that implicate residues Thr6, Ser9, and Thr15 of α-Btx, as well as Thr6 and Pro7 of α-Cbtx. Additionally the MD simulations of the human α1 ECD–toxin complexes confirmed the possible accommodation of two water molecules into a hydration cavity inside the hydrophobic core of the subunit, which may contribute to the gating mechanism of the receptor. These findings provide additional structural data that are intended to inspire biophysical studies on the functional role of glycosylation in the gating mechanism of nAChR and also guide the development of novel therapeutic agents for the treatment of nAChR-associated diseases.

Simulation of processes in cellular membranes Abstract Many important processes in cells involve ions, e.g., fusion of synaptic vesi- cles with neuronal cell membranes is controlled by a divalent cation Ca2+ ; and the exchange of Na+ and K+ drives the the fast electrical signal transmis- sion in neurons. We have investigated model phospholipid membranes and their interactions with these biologically relevant ions. Using state-of-the-art molecular dynamics simulations, we accurately quantified their respective affinites towards neutral and negatively charged phospholipid bilayers. In order to achieve that, we developed a new model of phospholipids termed ECC-lipids, which accounts for the electronic polarization via the electronic continuum correction implemented as charge rescaling. Our simulations with this new force field reach for the first time a quantitative agreement with the experimental lipid electrometer concept for POPC as well as for POPS with all the studied cations. We have also examined the effects of transmembrane voltage on phospholipid bilayers. The electric field induced by the voltage exists exclusively in the hydrophobic region of the membrane, where it has an almost constant strength. This field affects the structure of nearby water molecules highlighting its importance in electroporation. 1

Recent experimental studies indicate that polymeric structures with a well-adjusted balance of amphiphilic parts may translocate through self-assembled phospholipid bilayers and enhance the passive trans-membrane transport of smaller molecules. Using a coarse grained lattice Monte Carlo model with explicit solvent we investigate self-assembled lipid bilayers interacting with a linear polymer chain under variation of the hydrophobicity of the chain. Here, we focus on the relationship between the chain's hydrophobicity and its translocation behavior through the membrane as well as induced membrane perturbations. We show, that there is an adsorption transition of the polymer at the bilayer interface, where effectively the solvent phase and the tail phase of the bilayer are equally repulsive for the polymer. Close to this adsorption threshold of the polymer both the translocation probability of the polymer as well as the permeability of the membrane with respect to solvent are enhanced significantly. The frequency of polymer translocation events can be understood quantitatively assuming a simple diffusion along a one-dimensional free energy profile, which is controlled by the effective lipophilicity of the chain and the tail-packing in the bilayer's core. / Dieser Beitrag ist mit Zustimmung des Rechteinhabers aufgrund einer (DFG-geförderten) Allianz- bzw. Nationallizenz frei zugänglich.

No description available.

The human immunodeficiency virus type 1 (HIV-1) protease (PR) is a critical drug target as it is responsible for virion maturation. Mutations within the active site (1°) of the PR directly interfere with inhibitor binding while mutations distal to the active site (2°) to restore enzymatic fitness. Increasing mutation number is not directly proportional to the severity of resistance, suggesting that resistance is not simply additive but that it is interdependent. The interdependency of both primary and secondary mutations to drive protease inhibitor (PI) resistance is grossly understudied. To structurally and dynamically characterize the direct role of secondary mutations in drug resistance, I selected a panel of single-site mutant protease crystal structures complexed with the PI darunavir (DRV). From these studies, I developed a network hypothesis that explains how mutations outside the active site are able to perpetuate changes to the active site of the protease to disrupt inhibitor binding. I then expanded the panel to include highly mutated multi-drug resistant variants. To elucidate the interdependency between primary and secondary mutations I used statistical and machine-learning techniques to determine which specific mutations underlie the perturbations of key inter-molecular interactions. From these studies, I have determined that mutations distal to the active site are able to perturb the global PR hydrogen bonding patterns, while primary and secondary mutations cooperatively perturb hydrophobic contacts between the PR and DRV. Discerning and exploiting the mechanisms that underlie drug resistance in viral targets could proactively ameliorate both current treatment and inhibitor design for HIV-1 targets.