ProMod3-A versatile homology modelling toolbox.

Computational methods for protein structure modelling are routinely used to complement experimental structure determination, thus they help to address a broad spectrum of scientific questions in biomedical research. The most accurate methods today are based on homology modelling, i.e. detecting a homologue to the desired target sequence that can be used as a template for modelling. Here we present a versatile open source homology modelling toolbox as foundation for flexible and computationally efficient modelling workflows. ProMod3 is a fully scriptable software platform that can perform all steps required to generate a protein model by homology. Its modular design aims at fast prototyping of novel algorithms and implementing flexible modelling pipelines. Common modelling tasks, such as loop modelling, sidechain modelling or generating a full protein model by homology, are provided as production ready pipelines, forming the starting point for own developments and enhancements. ProMod3 is the central software component of the widely used SWISS-MODEL web-server.

Conformational dynamics and role of the acidic pocket in ASIC pH-dependent gating.

Acid-sensing ion channels (ASICs) are proton-activated Na+ channels expressed in the nervous system, where they are involved in learning, fear behaviors, neurodegeneration, and pain sensation. In this work, we study the role in pH sensing of two regions of the ectodomain enriched in acidic residues: the acidic pocket, which faces the outside of the protein and is the binding site of several animal toxins, and the palm, a central channel domain. Using voltage clamp fluorometry, we find that the acidic pocket undergoes conformational changes during both activation and desensitization. Concurrently, we find that, although proton sensing in the acidic pocket is not required for channel function, it does contribute to both activation and desensitization. Furthermore, protonation-mimicking mutations of acidic residues in the palm induce a dramatic acceleration of desensitization followed by the appearance of a sustained current. In summary, this work describes the roles of potential pH sensors in two extracellular domains, and it proposes a model of acidification-induced conformational changes occurring in the acidic pocket of ASIC1a.

Implementation of a methodology for determining elastic properties of lipid assemblies from molecular dynamics simulations.

BACKGROUND: The importance of the material properties of membranes for diverse cellular processes is well established. Notably, the elastic properties of the membrane, which depend on its composition, can directly influence membrane reshaping and fusion processes as well as the organisation and function of membrane proteins. Determining these properties is therefore key for a mechanistic understanding of how the cell functions. RESULTS: We have developed a method to determine the bending rigidity and tilt modulus, for lipidic assemblies of arbitrary lipid composition and shape, from molecular dynamics simulations. The method extracts the elastic moduli from the distributions of microscopic tilts and splays of the lipid components. We present here an open source implementation of the method as a set of Python modules using the computational framework OpenStructure. These modules offer diverse algorithms typically used in the calculatation the elastic moduli, including routines to align MD trajectories of complex lipidic systems, to determine the water/lipid interface, to calculate lipid tilts and splays, as well as to fit the corresponding distributions to extract the elastic properties. We detail the implementation of the method and give several examples of how to use the modules in specific cases. CONCLUSIONS: The method presented here is, to our knowledge, the only available computational approach allowing to quantify the elastic properties of lipidic assemblies of arbitrary shape and composition (including lipid mixtures). The implementation as python modules offers flexibility, which has already allowed the method to be applied to diverse lipid assembly types, ranging from bilayers in the liquid ordered and disordered phases to a study of the inverted-hexagonal phase, and with different force-fields (both all-atom and coarse grained representations). The modules are freely available through GitHub at https://github.com/njohner/ost_pymodules/ while OpenStructure can be obtained at http://www.openstructure.org .

Computational modeling of the N-terminus of the human dopamine transporter and its interaction with PIP2 -containing membranes.

The dopamine transporter (DAT) is a transmembrane protein belonging to the family of neurotransmitter:sodium symporters (NSS). Members of the NSS are responsible for the clearance of neurotransmitters from the synaptic cleft, and for their translocation back into the presynaptic nerve terminal. The DAT contains long intracellular N- and C-terminal domains that are strongly implicated in the transporter function. The N-terminus (N-term), in particular, regulates the reverse transport (efflux) of the substrate through DAT. Currently, the molecular mechanisms of the efflux remain elusive in large part due to lack of structural information on the N-terminal segment. Here we report a computational model of the N-term of the human DAT (hDAT), obtained through an ab initio structure prediction, in combination with extensive atomistic molecular dynamics (MD) simulations in the context of a lipid membrane. Our analysis reveals that whereas the N-term is a highly dynamic domain, it contains secondary structure elements that remain stable in the long MD trajectories of interactions with the bilayer (totaling >2.2 µs). Combining MD simulations with continuum mean-field modeling we found that the N-term engages with lipid membranes through electrostatic interactions with the charged lipids PIP2 (phosphatidylinositol 4,5-Biphosphate) or PS (phosphatidylserine) that are present in these bilayers. We identify specific motifs along the N-term implicated in such interactions and show that differential modes of N-term/membrane association result in differential positioning of the structured segments on the membrane surface. These results will inform future structure-based studies that will elucidate the mechanistic role of the N-term in DAT function.

Protein and lipid interactions driving molecular mechanisms of in meso crystallization.

The recent advances in the in meso crystallization technique for the structural characterization of G-protein coupled receptor (GPCR) proteins have established the usefulness of the lipidic-cubic phases (LCPs) in the field of crystallography of membrane proteins. It is surprising that despite the success of the approach, the molecular mechanisms of the in meso method are still not well understood. Therefore, the approach must rely on extensive screening for a suitable protein construct, for host and additive lipids, and for the appropriate precipitants and temperature. To shed light on the in meso crystallization mechanisms, we used extensive coarse-grained molecular dynamics simulations to study, in molecular detail, LCPs under different conditions (compositions and temperatures relevant to crystallogenesis) and their interactions with different types of GPCR constructs. The results presented show how the modulation of the lattice constant of the LCP (triggered by the addition of precipitant during the in meso assay), or of the host lipid type, can destabilize monomeric proteins in the bilayer of the LCP and thus drive their aggregation into the stacked lamellae, where the residual hydrophobic mismatch between the protein and the membrane can drive the formation of lateral contacts leading to nucleation and crystal growth. Moreover, we demonstrate how particular protein designs (such as transmembrane proteins engineered to contain large polar regions) can promote protein stacking interactions in the third, out-of-plane, dimension. The insights provided by the new aspects of the specific molecular mechanisms responsible for protein-protein interactions inside the cubic phase presented here should be helpful in guiding the rational design of future in meso trials with successful outcomes.

How the dynamic properties and functional mechanisms of GPCRs are modulated by their coupling to the membrane environment.

Experimental observations of the dependence of function and organization of G protein-coupled receptors (GPCRs) on their lipid environment have stimulated new quantitative studies of the coupling between the proteins and the membrane. It is important to develop such a quantitative understanding at the molecular level because the effects of the coupling are seen to be physiologically and clinically significant. Here we review findings that offer insight into how membrane-GPCR coupling is connected to the structural characteristics of the GPCR, from sequence to 3D structural detail, and how this coupling is involved in the actions of ligands on the receptor. The application of a recently developed computational approach designed for quantitative evaluation of membrane remodeling and the energetics of membrane-protein interactions brings to light the importance of the radial asymmetry of the membrane-facing surface of GPCRs in their interaction with the surrounding membrane. As the radial asymmetry creates adjacencies of hydrophobic and polar residues at specific sites of the GPCR, the ability of membrane remodeling to achieve complete hydrophobic matching is limited, and the residual mismatch carries a significant energy cost. The adjacencies are shown to be affected by ligand-induced conformational changes. Thus, functionally important organization of GPCRs in the cell membrane can depend both on ligand-determined properties and on the lipid composition of various membrane regions with different remodeling capacities. That this functionally important reorganization can be driven by oligomerization patterns that reduce the energy cost of the residual mismatch, suggests a new perspective on GPCR dimerization and ligand-GPCR interactions. The relation between the modulatory effects on GPCRs from the binding of specific cell-membrane components, e.g., cholesterol, and those produced by the non-local energetics of hydrophobic mismatch are discussed in this context.

Why GPCRs behave differently in cubic and lamellar lipidic mesophases.

Recent successes in the crystallographic determination of structures of transmembrane proteins in the G protein-coupled receptor (GPCR) family have established the lipidic cubic phase (LCP) environment as the medium of choice for growing structure-grade crystals by the method termed "in meso". The understanding of in meso crystallogenesis is currently at a descriptive level. To enable an eventual quantitative, energy-based description of the nucleation and crystallization mechanism, we have examined the properties of the lipidic cubic phase system and the dynamics of the GPCR rhodopsin reconstituted into the LCP with coarse-grained molecular dynamics simulations with the Martini force-field. Quantifying the differences in the hydrophobic/hydrophilic exposure of the GPCR to lipids in the cubic and lamellar phases, we found that the highly curved geometry of the cubic phase provides more efficient shielding of the protein from unfavorable hydrophobic exposure, which leads to a lesser hydrophobic mismatch and less unfavorable hydrophobic-hydrophilic interactions between the protein and lipid-water interface in the LCP, compared to the lamellar phase. Since hydrophobic mismatch is considered a driving force for oligomerization, the differences in exposure mismatch energies between the LCP and the lamellar structures suggest that the latter provide a more favorable setting in which GPCRs can oligomerize as a prelude to nucleation and crystal growth. These new findings lay the foundation for future investigations of in meso crystallization mechanisms related to the transition from the LCP to the lamellar phase and studies aimed at an improved rational approach for generating structure-quality crystals of membrane proteins.

Calcium regulates acid-sensing ion channel 3 activation by competing with protons in the channel pore and at an allosteric binding site.

The extracellular Ca2+ concentration changes locally under certain physiological and pathological conditions. Such variations affect the function of ion channels of the nervous system and consequently also neuronal signalling. We investigated here the mechanisms by which Ca2+ controls the activity of acid-sensing ion channel (ASIC) 3. ASICs are neuronal, H+-gated Na+ channels involved in several physiological and pathological processes, including the expression of fear, learning, pain sensation and neurodegeneration after ischaemic stroke. It was previously shown that Ca2+ negatively modulates the ASIC pH dependence. While protons are default activators of ASIC3, this channel can also be activated at pH7.4 by the removal of the extracellular Ca2+. Two previous studies concluded that low pH opens ASIC3 by displacing Ca2+ ions that block the channel pore at physiological pH. We show here that an acidic residue, distant from the pore, together with pore residues, controls the modulation of ASIC3 by Ca2+. Our study identifies a new regulatory site in ASIC3 and demonstrates that ASIC3 activation involves an allosteric mechanism together with Ca2+ unbinding from the channel pore. We provide a molecular analysis of a regulatory mechanism found in many ion channels.

Optimal percolation of disordered segregated composites.

We evaluate the percolation threshold values for a realistic model of continuum segregated systems, where random spherical inclusions forbid the percolating objects, modeled by hardcore spherical particles surrounded by penetrable shells, to occupy large regions inside the composite. We find that the percolation threshold is generally a nonmonotonous function of segregation, and that an optimal (i.e., minimum) critical concentration exists well before maximum segregation is reached. We interpret this feature as originating from a competition between reduced available volume effects and enhanced concentrations needed to ensure percolation in the highly segregated regime. The relevance with existing segregated materials is discussed.

Percolative properties of hard oblate ellipsoids of revolution with a soft shell.

We present an in-depth analysis of the geometrical percolation behavior in the continuum of random assemblies of hard oblate ellipsoids of revolution. Simulations were carried out by considering a broad range of aspect ratios, from spheres up to aspect-ratio-100 platelike objects, and with various limiting two-particle interaction distances, from 0.05 times the major axis up to 4.0 times the major axis. We confirm the widely reported trend of a consistent lowering of the hard particle critical volume fraction with increase of the aspect ratio. Moreover, by assimilating the limiting interaction distance to a shell of constant thickness surrounding the ellipsoids, we propose a simple relation based on the total excluded volume of these objects which allows us to estimate the critical concentration from a quantity that is quasi-invariant over a large spectrum of limiting interaction distances. Excluded volume and volume quantities are derived explicitly.

An Amphipathic Helix Directs Cellular Membrane Curvature Sensing and Function of the BAR Domain Protein PICK1.

BAR domains are dimeric protein modules that sense, induce, and stabilize lipid membrane curvature. Here, we show that membrane curvature sensing (MCS) directs cellular localization and function of the BAR domain protein PICK1. In PICK1, and the homologous proteins ICA69 and arfaptin2, we identify an amphipathic helix N-terminal to the BAR domain that mediates MCS. Mutational disruption of the helix in PICK1 impaired MCS without affecting membrane binding per se. In insulin-producing INS-1E cells, super-resolution microscopy revealed that disruption of the helix selectively compromised PICK1 density on insulin granules of high curvature during their maturation. This was accompanied by reduced hormone storage in the INS-1E cells. In Drosophila, disruption of the helix compromised growth regulation. By demonstrating size-dependent binding on insulin granules, our finding highlights the function of MCS for BAR domain proteins in a biological context distinct from their function, e.g., at the plasma membrane during endocytosis.

Structure of Dimeric and Tetrameric Complexes of the BAR Domain Protein PICK1 Determined by Small-Angle X-Ray Scattering.

PICK1 is a neuronal scaffolding protein containing a PDZ domain and an auto-inhibited BAR domain. BAR domains are membrane-sculpting protein modules generating membrane curvature and promoting membrane fission. Previous data suggest that BAR domains are organized in lattice-like arrangements when stabilizing membranes but little is known about structural organization of BAR domains in solution. Through a small-angle X-ray scattering (SAXS) analysis, we determine the structure of dimeric and tetrameric complexes of PICK1 in solution. SAXS and biochemical data reveal a strong propensity of PICK1 to form higher-order structures, and SAXS analysis suggests an offset, parallel mode of BAR-BAR oligomerization. Furthermore, unlike accessory domains in other BAR domain proteins, the positioning of the PDZ domains is flexible, enabling PICK1 to perform long-range, dynamic scaffolding of membrane-associated proteins. Together with functional data, these structural findings are compatible with a model in which oligomerization governs auto-inhibition of BAR domain function.

Curvature and Lipid Packing Modulate the Elastic Properties of Lipid Assemblies: Comparing HII and Lamellar Phases.

Accumulating evidence indicates that membrane reshaping and fusion processes, as well as regulation of membrane protein function, depend on lipid composition. Although it is widely accepted that cell membranes are under considerable stress and frustration and can be locally highly curved, experimental approaches to determine the material properties of lipids usually rely on their study in a relaxed environment or in flat bilayers. Here, we propose a computational method to determine the elastic properties of lipid assemblies of arbitrarily shaped interfaces and apply it to lipidic mixtures in the inverted hexagonal and lamellar phases. We find that the bending rigidity critically depends on the geometry of the system and correlates with the changes in lipid chain order imposed by the specific environment. Our results are relevant for resolving local lipid properties of deformed, stressed, or frustrated membranes that notably emerge around integral membrane proteins or during different membrane remodeling processes.

Synergistic substrate binding determines the stoichiometry of transport of a prokaryotic H(+)/Cl(-) exchanger.

Active exchangers dissipate the gradient of one substrate to accumulate nutrients, export xenobiotics and maintain cellular homeostasis. Mechanistic studies have suggested that two fundamental properties are shared by all exchangers: substrate binding is antagonistic, and coupling is maintained by preventing shuttling of the empty transporter. The CLC H(+)/Cl(-) exchangers control the homeostasis of cellular compartments in most living organisms, but their transport mechanism remains unclear. We show that substrate binding to CLC-ec1 is synergistic rather than antagonistic: chloride binding induces protonation of a crucial glutamate. The simultaneous binding of H(+) and Cl(-) gives rise to a fully loaded state that is incompatible with conventional transport mechanisms. Mutations in the Cl(-) transport pathway identically alter the stoichiometries of H(+)/Cl(-) exchange and binding. We propose that the thermodynamics of synergistic substrate binding, rather than the kinetics of conformational changes and ion binding, determine the stoichiometry of transport.