Membrane specificity of the human cholesterol transfer protein STARD4

STARD4 regulates cholesterol homeostasis by transferring cholesterol between the plasma membrane and endoplasmic reticulum. The STARD4 structure features a helix-grip fold surrounding a large hydrophobic cavity holding the sterol. Its access is controlled by a gate formed by the Ω1 and Ω4 loops and the C-terminal α-helix. Little is known about the mechanisms by which STARD4 binds to membranes and extracts/releases cholesterol. All available structures of STARD4 are without a bound sterol and display the same closed conformation of the gate. The cholesterol transfer activity of the mouse STARD4 is enhanced in the presence of anionic lipids, and in particular of phosphatidylinositol biphosphates (PIP2) for which two binding sites were proposed on the mouse STARD4 surface. Yet only one of these sites is conserved in human STARD4. We here report the results of a liposome microarray-based assay and microseconds-long molecular dynamics simulations of human STARD4 with complex lipid bilayers mimicking the composition of the donor and acceptor membranes. We show that the binding of apo form of human STARD4 is sensitive to the presence of PIP2 through two specific binding sites, one of which was not identified on mouse STARD4. We report two novel conformations of the gate in holo-STARD4: a yet-unobserved close conformation and an open conformation of Ω4 shedding light on the opening/closure mechanism needed for cholesterol uptake/release. Overall, the modulation of human STARD4 membrane-binding by lipid composition, and by the presence of the cargo supports the capacity of human STARD4 to achieve directed transfer between specific organelle membranes.

Investigating Polypharmacology through Targeting Known Human Neutrophil Elastase Inhibitors to Proteinase 3

Using a combination of multisite λ-dynamics (MSλD) together with in vitro IC50 assays, we evaluated the polypharmacological potential of a scaffold currently in clinical trials for inhibition of human neutrophil elastase (HNE), targeting cardiopulmonary disease, for efficacious inhibition of Proteinase 3 (PR3), a related neutrophil serine proteinase. The affinities we observe suggest that the dihydropyrimidinone scaffold can serve as a suitable starting point for the establishment of polypharmacologically targeting both enzymes and enhancing the potential for treatments addressing diseases like chronic obstructive pulmonary disease.

Overlay databank unlocks data-driven analyses of biomolecules for all

Tools based on artificial intelligence (AI) are currently revolutionising many fields, yet their applications are often limited by the lack of suitable training data in programmatically accessible format. Here we propose an effective solution to make data scattered in various locations and formats accessible for data-driven and machine learning applications using the overlay databank format. To demonstrate the practical relevance of such approach, we present the NMRlipids Databank-a community-driven, open-for-all database featuring programmatic access to quality-evaluated atom-resolution molecular dynamics simulations of cellular membranes. Cellular membrane lipid composition is implicated in diseases and controls major biological functions, but membranes are difficult to study experimentally due to their intrinsic disorder and complex phase behaviour. While MD simulations have been useful in understanding membrane systems, they require significant computational resources and often suffer from inaccuracies in model parameters. Here, we demonstrate how programmable interface for flexible implementation of data-driven and machine learning applications, and rapid access to simulation data through a graphical user interface, unlock possibilities beyond current MD simulation and experimental studies to understand cellular membranes. The proposed overlay databank concept can be further applied to other biomolecules, as well as in other fields where similar barriers hinder the AI revolution.

Standard Binding Free Energy and Membrane Desorption Mechanism for a Phospholipase C

Peripheral membrane proteins (PMPs) bind temporarily to cellular membranes and play important roles in signaling, lipid metabolism, and membrane trafficking. Obtaining accurate membrane-PMP affinities using experimental techniques is more challenging than for protein-ligand affinities in an aqueous solution. At the theoretical level, calculation of the standard protein-membrane binding free energy using molecular dynamics simulations remains a daunting challenge owing to the size of the biological objects at play, the slow lipid diffusion, and the large variation in configurational entropy that accompanies the binding process. To overcome these challenges, we used a computational framework relying on a series of potential-of-mean-force (PMF) calculations including a set of geometrical restraints on collective variables. This methodology allowed us to determine the standard binding free energy of a PMP to a phospholipid bilayer using an all-atom force field. Bacillus thuringiensis phosphatidylinositol-specific phospholipase C (BtPI-PLC) was chosen due to its importance as a virulence factor and owing to the host of experimental affinity data available. We computed a standard binding free energy of -8.2 ± 1.4 kcal/mol in reasonable agreement with the reported experimental values (-6.6 ± 0.2 kcal/mol). In light of the 2.3-μs separation PMF calculation, we investigated the mechanism whereby BtPI-PLC disengages from interactions with the lipid bilayer during separation. We describe how a short amphipathic helix engages in transitory interactions to ease the passage of its hydrophobes through the interfacial region upon desorption from the bilayer.

Dissecting peripheral protein-membrane interfaces

Peripheral membrane proteins (PMPs) include a wide variety of proteins that have in common to bind transiently to the chemically complex interfacial region of membranes through their interfacial binding site (IBS). In contrast to protein-protein or protein-DNA/RNA interfaces, peripheral protein-membrane interfaces are poorly characterized. We collected a dataset of PMP domains representative of the variety of PMP functions: membrane-targeting domains (Annexin, C1, C2, discoidin C2, PH, PX), enzymes (PLA, PLC/D) and lipid-transfer proteins (START). The dataset contains 1328 experimental structures and 1194 AphaFold models. We mapped the amino acid composition and structural patterns of the IBS of each protein in this dataset, and evaluated which were more likely to be found at the IBS compared to the rest of the domains' accessible surface. In agreement with earlier work we find that about two thirds of the PMPs in the dataset have protruding hydrophobes (Leu, Ile, Phe, Tyr, Trp and Met) at their IBS. The three aromatic amino acids Trp, Tyr and Phe are a hallmark of PMPs IBS regardless of whether they protrude on loops or not. This is also the case for lysines but not arginines suggesting that, unlike for Arg-rich membrane-active peptides, the less membrane-disruptive lysine is preferred in PMPs. Another striking observation was the over-representation of glycines at the IBS of PMPs compared to the rest of their surface, possibly procuring IBS loops a much-needed flexibility to insert in-between membrane lipids. The analysis of the 9 superfamilies revealed amino acid distribution patterns in agreement with their known functions and membrane-binding mechanisms. Besides revealing novel amino acids patterns at protein-membrane interfaces, our work contributes a new PMP dataset and an analysis pipeline that can be further built upon for future studies of PMPs properties, or for developing PMPs prediction tools using for example, machine learning approaches.

Phosphatidylcholine Cation—Tyrosine π Complexes: Motifs for Membrane Binding by a Bacterial Phospholipase C

Phosphatidylinositol-specific phospholipase C (PI-PLC) enzymes are a virulence factor in many Gram-positive organisms. The specific activity of the Bacillus thuringiensis PI-PLC is significantly increased by adding phosphatidylcholine (PC) to vesicles composed of the substrate phosphatidylinositol, in part because the inclusion of PC reduces the apparent Kd for the vesicle binding by as much as 1000-fold when comparing PC-rich vesicles to PI vesicles. This review summarizes (i) the experimental work that localized a site on BtPI-PLC where PC is bound as a PC choline cation—Tyr-π complex and (ii) the computational work (including all-atom molecular dynamics simulations) that refined the original complex and found a second persistent PC cation—Tyr-π complex. Both complexes are critical for vesicle binding. These results have led to a model for PC functioning as an allosteric effector of the enzyme by altering the protein dynamics and stabilizing an ‘open’ active site conformation.

Specificity of Loxosceles α clade phospholipase D enzymes for choline-containing lipids: Role of a conserved aromatic cage

Spider venom GDPD-like phospholipases D (SicTox) have been identified to be one of the major toxins in recluse spider venom. They are divided into two major clades: the α clade and the β clade. Most α clade toxins present high activity against lipids with choline head groups such as sphingomyelin, while activities in β clade toxins vary and include preference for substrates containing ethanolamine headgroups (Sicarius terrosus, St_βIB1). A structural comparison of available structures of phospholipases D (PLDs) reveals a conserved aromatic cage in the α clade. To test the potential influence of the aromatic cage on membrane-lipid specificity we performed molecular dynamics (MD) simulations of the binding of several PLDs onto lipid bilayers containing choline headgroups; two SicTox from the α clade, Loxosceles intermedia αIA1 (Li_αIA) and Loxosceles laeta αIII1 (Ll_αIII1), and one from the β clade, St_βIB1. The simulation results reveal that the aromatic cage captures a choline-headgroup and suggest that the cage plays a major role in lipid specificity. We also simulated an engineered St_βIB1, where we introduced the aromatic cage, and this led to binding with choline-containing lipids. Moreover, a multiple sequence alignment revealed the conservation of the aromatic cage among the α clade PLDs. Here, we confirmed that the i-face of α and β clade PLDs is involved in their binding to choline and ethanolamine-containing bilayers, respectively. Furthermore, our results suggest a major role in choline lipid recognition of the aromatic cage of the α clade PLDs. The MD simulation results are supported by in vitro liposome binding assay experiments.

Envenomation following bites from recluse spiders (Loxosceles) causes loxoscelism, a necrotic tissue breakdown in mammals, and leads to skin degeneration and systemic reactions in the worst case. Recluse spiders belong to the Sicariidae family which also includes six-eyed sand spiders in the genera Sicarius and Hexopthalma. While sicariid spiders are found natively on all continents except Australia, treatments of loxoscelism are typically antibody based and available in some regions of the Americas. Sphingomyelinase D/phospholipase D enzymes are one of the major toxins in venom of sicariid spiders, and have been divided in two clades called α and β. The activity of α and β clades toxins differs; most α clade toxins present high activity against lipids with choline headgroups (-N (CH₃)₃⁺) such as sphingomyelin, while activities in β clade toxins vary and include preference for substrates containing ethanolamine headgroups (-NH₃⁺). When comparing the structures of two α clade toxins and one β clade toxin, we noticed the presence in the α clade toxins only of a cage consisting of three aromatic amino acids. In this work we used numerical molecular simulations to probe the role of this cage in the preference of α clade toxins for choline head groups over ethanolamine head groups.

PDBe-KB: collaboratively defining the biological context of structural data

The Protein Data Bank in Europe - Knowledge Base (PDBe-KB, https://pdbe-kb.org) is an open collaboration between world-leading specialist data resources contributing functional and biophysical annotations derived from or relevant to the Protein Data Bank (PDB). The goal of PDBe-KB is to place macromolecular structure data in their biological context by developing standardised data exchange formats and integrating functional annotations from the contributing partner resources into a knowledge graph that can provide valuable biological insights. Since we described PDBe-KB in 2019, there have been significant improvements in the variety of available annotation data sets and user functionality. Here, we provide an overview of the consortium, highlighting the addition of annotations such as predicted covalent binders, phosphorylation sites, effects of mutations on the protein structure and energetic local frustration. In addition, we describe a library of reusable web-based visualisation components and introduce new features such as a bulk download data service and a novel superposition service that generates clusters of superposed protein chains weekly for the whole PDB archive.

Phospholipids in Motion: High-Resolution 31P NMR Field Cycling Studies

Diverse phospholipid motions are key to membrane function but can be quite difficult to untangle and quantify. High-resolution field cycling ³¹P NMR spin-lattice relaxometry, where the sample is excited at high field, shuttled in the magnet bore for low-field relaxation, then shuttled back to high field for readout of the residual magnetization, provides data on phospholipid dynamics and structure. This information is encoded in the field dependence of the ³¹P spin-lattice relaxation rate (R₁). In the field range from 11.74 down to 0.003 T, three dipolar nuclear magnetic relaxation dispersions (NMRDs) and one due to ³¹P chemical shift anisotropy contribute to R₁ of phospholipids. Extraction of correlation times and maximum relaxation amplitudes for these NMRDs provides (1) lateral diffusion constants for different phospholipids in the same bilayer, (2) estimates of how additives alter the motion of the phospholipid about its long axis, and (3) an average ³¹P-¹H angle with respect to the bilayer normal, which reveals that polar headgroup motion is not restricted on a microsecond timescale. Relative motions within a phospholipid are also provided by comparing ³¹P NMRD profiles for specifically deuterated molecules as well as ¹³C and ¹H field dependence profiles to that of ³¹P. Although this work has dealt exclusively with phospholipids in small unilamellar vesicles, these same NMRDs can be measured for phospholipids in micelles and nanodisks, making this technique useful for monitoring lipid behavior in a variety of structures and assessing how additives alter specific lipid motions.

Martini 3: a general purpose force field for coarse-grained molecular dynamics

The coarse-grained Martini force field is widely used in biomolecular simulations. Here we present the refined model, Martini 3 ( http://cgmartini.nl ), with an improved interaction balance, new bead types and expanded ability to include specific interactions representing, for example, hydrogen bonding and electronic polarizability. The updated model allows more accurate predictions of molecular packing and interactions in general, which is exemplified with a vast and diverse set of applications, ranging from oil/water partitioning and miscibility data to complex molecular systems, involving protein-protein and protein-lipid interactions and material science applications as ionic liquids and aedamers.

Classification and phylogeny for the annotation of novel eukaryotic GNAT acetyltransferases

The enzymes of the GCN5-related N-acetyltransferase (GNAT) superfamily count more than 870 000 members through all kingdoms of life and share the same structural fold. GNAT enzymes transfer an acyl moiety from acyl coenzyme A to a wide range of substrates including aminoglycosides, serotonin, glucosamine-6-phosphate, protein N-termini and lysine residues of histones and other proteins. The GNAT subtype of protein N-terminal acetyltransferases (NATs) alone targets a majority of all eukaryotic proteins stressing the omnipresence of the GNAT enzymes. Despite the highly conserved GNAT fold, sequence similarity is quite low between members of this superfamily even when substrates are similar. Furthermore, this superfamily is phylogenetically not well characterized. Thus functional annotation based on sequence similarity is unreliable and strongly hampered for thousands of GNAT members that remain biochemically uncharacterized. Here we used sequence similarity networks to map the sequence space and propose a new classification for eukaryotic GNAT acetyltransferases. Using the new classification, we built a phylogenetic tree, representing the entire GNAT acetyltransferase superfamily. Our results show that protein NATs have evolved more than once on the GNAT acetylation scaffold. We use our classification to predict the function of uncharacterized sequences and verify by in vitro protein assays that two fungal genes encode NAT enzymes targeting specific protein N-terminal sequences, showing that even slight changes on the GNAT fold can lead to change in substrate specificity. In addition to providing a new map of the relationship between eukaryotic acetyltransferases the classification proposed constitutes a tool to improve functional annotation of GNAT acetyltransferases.

Enzymes of the GCN5-related N-acetyltransferase (GNAT) superfamily transfer an acetyl group from one molecule to another. This reaction is called acetylation and is one of the most common reactions inside the cell. The GNAT superfamily counts more than 870 000 members through all kingdoms of life. Despite sharing the same fold the GNAT superfamily is very diverse in terms of amino acid sequence and substrates. The eight N-terminal acetyltransferases (NatA, NatB, etc.. to NatH) are a GNAT subtype which acetylates the free amine group of polypeptide chains. This modification is called N-terminal acetylation and is one of the most abundant protein modifications in eukaryotic cells. This subtype is also characterized by a high sequence diversity even though they share the same substrate. In addition, the phylogeny of the superfamily is not characterized. This hampers functional annotation based on sequence similarity, and discovery of novel NATs. In this work we set out to solve the problem of the classification of eukaryotic GCN5-related acetyltransferases and report the first classification framework of the superfamily. This framework can be used as a tool for annotation of all GCN5-related acetyltransferases. As an example of what can be achieved we report in this paper the computational prediction and in vitro verification of the function of two previously uncharacterized N-terminal acetyltransferases. We also report the first acetyltransferase phylogenetic tree of the GCN5 superfamily. It indicates that N-terminal acetyltransferases do not constitute one homogeneous protein family, but that the ability to bind and acetylate protein N-termini had evolved more than once on the same acetylation scaffold. We also show that even small changes in key positions can lead to altered enzyme specificity.

Capturing Choline-Aromatics Cation-π Interactions in the MARTINI Force Field

Cation−π interactions play an important role in biomolecular recognition, including interactions between membrane phosphatidylcholine lipids and aromatic amino acids of peripheral proteins. While molecular mechanics coarse grain (CG) force fields are particularly well suited to simulate membrane proteins in general, they are not parameterized to explicitly reproduce cation−π interactions. We here propose a modification of the polarizable MARTINI coarse grain (CG) model enabling it to model membrane binding events of peripheral proteins whose aromatic amino acid interactions with choline headgroups are crucial for their membrane binding. For this purpose, we first collected and curated a dataset of eight peripheral proteins from different families. We find that the MARTINI CG model expectedly underestimates aromatics–choline interactions and is unable to reproduce membrane binding of the peripheral proteins in our dataset. Adjustments of the relevant interactions in the polarizable MARTINI force field yield significant improvements in the observed binding events. The orientation of each membrane-bound protein is comparable to reference data from all-atom simulations and experimental binding data. We also use negative controls to ensure that choline–aromatics interactions are not overestimated. We finally check that membrane properties, transmembrane proteins, and membrane translocation potential of mean force (PMF) of aromatic amino acid side-chain analogues are not affected by the new parameter set. This new version “MARTINI 2.3P” is a significant improvement over its predecessors and is suitable for modeling membrane proteins including peripheral membrane binding of peptides and proteins.

Dynamics-function relationship in the catalytic domains of N-terminal acetyltransferases

N-terminal acetyltransferases (NATs) belong to the superfamily of acetyltransferases. They are enzymes catalysing the transfer of an acetyl group from acetyl coenzyme A to the N-terminus of polypeptide chains. N-terminal acetylation is one of the most common protein modifications. To date, not much is known on the molecular basis for the exclusive substrate specificity of NATs. All NATs share a common fold called GNAT. A characteristic of NATs is the β6β7 hairpin loop covering the active site and forming with the α1α2 loop a narrow tunnel surrounding the catalytic site in which cofactor and polypeptide meet and exchange an acetyl group. We investigated the dynamics-function relationships of all available structures of NATs covering the three domains of Life. Using an elastic network model and normal mode analysis, we found a common dynamics pattern conserved through the GNAT fold; a rigid V-shaped groove formed by the β4 and β5 strands and splitting the fold in two dynamical subdomains. Loops α1α2, β3β4 and β6β7 all show clear displacements in the low frequency normal modes. We characterized the mobility of the loops and show that even limited conformational changes of the loops along the low-frequency modes are able to significantly change the size and shape of the ligand binding sites. Based on the fact that these movements are present in most low-frequency modes, and common to all NATs, we suggest that the α1α2 and β6β7 loops may regulate ligand uptake and the release of the acetylated polypeptide.

The Arabidopsis (ASHH2) CW domain binds monomethylated K4 of the histone H3 tail through conformational selection

Chromatin post-translational modifications are thought to be important for epigenetic effects on gene expression. Methylation of histone N-terminal tail lysine residues constitutes one of many such modifications, executed by families of histone lysine methyltransferase (HKMTase). One such protein is ASHH2 from the flowering plant Arabidopsis thaliana, equipped with the interaction domain, CW, and the HKMTase domain, SET. The CW domain of ASHH2 is a selective binder of monomethylation at lysine 4 on histone H3 (H3K4me1) and likely helps the enzyme dock correctly onto chromatin sites. The study of CW and related interaction domains has so far been emphasizing lock-key models, missing important aspects of histone-tail CW interactions. We here present an analysis of the ASHH2 CW-H3K4me1 complex using NMR and molecular dynamics, as well as mutation and affinity studies of flexible coils. β-augmentation and rearrangement of coils coincide with changes in the flexibility of the complex, in particular the η1, η3 and C-terminal coils, but also in the β1 and β2 strands and the C-terminal part of the ligand. Furthermore, we show that mutating residues with outlier dynamic behaviour affect the complex binding affinity despite these not being in direct contact with the ligand. Overall, the binding process is consistent with conformational selection. We propose that this binding mechanism presents an advantage when searching for the correct post-translational modification state among the highly modified and flexible histone tails, and also that the binding shifts the catalytic SET domain towards the nucleosome. DATABASES: Structural data are available in the PDB database under the accession code 6QXZ. Resonance assignments for CW42 in its apo- and holo-forms are available in the BMRB database under the accession code 27251.

Interfacial Aromatics Mediating Cation-π Interactions with Choline-Containing Lipids Can Contribute as Much to Peripheral Protein Affinity for Membranes as Aromatics Inserted below the Phosphates

Membrane-binding interfaces of peripheral proteins are restricted to a small part of their exposed surface, so the ability to engage in strong selective interactions with membrane lipids at various depths in the interface, both below and above the phosphates, is an advantage. Driven by their hydrophobicity, aromatic amino acids preferentially partition into membrane interfaces, often below the phosphates, yet enthalpically favorable interactions with the lipid headgroups, above the phosphate plane, are likely to further stabilize high interfacial positions. Using free-energy perturbation, we calculate the energetic cost of alanine substitution for 11 interfacial aromatic amino acids from 3 peripheral proteins. We show that the involvement in cation-π interactions with the headgroups (i) increases the ΔΔG_transfer as compared with insertion at the same depth without cation-π stabilization and (ii) can contribute at least as much as deeper insertion below the phosphates, highlighting the multiple roles of aromatics in peripheral membrane protein affinity.

Cation-π Interactions between Methylated Ammonium Groups and Tryptophan in the CHARMM36 Additive Force Field

Cation-π interactions between tryptophan and choline or trimethylated lysines are vital for many biological processes. The performance of the additive CHARMM36 force field against target quantum mechanical data is shown to reproduce QM equilibrium geometries but required modified Lennard-Jones potentials to accurately reproduce the QM interaction energies. The modified parameter set allows accurate modeling, including free energies, of cation-π indole-choline and indole-trimethylated lysines interactions relevant for protein-ligand, protein-membrane, and protein-protein interfaces.

d-Peptides as inhibitors of PR3-membrane interactions

Proteinase 3 (PR3) is a neutrophil serine protease present in cytoplasmic granules but also expressed at the neutrophil surface where it mediates proinflammatory effects. Studies of the underlying molecular mechanisms have been hampered by the lack of inhibitors of the PR3 membrane anchorage. Indeed while there exist inhibitors of the catalytic activity of PR3, its membrane interfacial binding site (IBS) is distinct from its catalytic site. The IBS has been characterized both by mutagenesis experiments and molecular modeling. Through docking and molecular dynamics simulations we have designed d-peptides targeting the PR3 IBS. We used surface plasmon resonance to evaluate their effect on the binding of PR3 to phospholipid bilayers. Next, we verified their ability of binding to PR3 via fluorescence spectroscopy and isothermal titration calorimetry. The designed peptides did not affect the catalytic activity of PR3. A few peptides bound to PR3 hydrophobic pockets and inhibited PR3 binding to lipids. While the (KFF)₃K d-peptide inconveniently showed a significant affinity for the lipids, another d-peptide (SAKEAFFKLLAS) did not and it inhibited the PR3-membrane binding site with IC₅₀ of about 40μM. Our work puts forward d-peptides as promising inhibitors of peripheral protein-membrane interactions, which remain high-hanging fruits in drug design.

Molecular determinants of the N-terminal acetyltransferase Naa60 anchoring to the Golgi membrane

Nα-Acetyltransferase 60 (Naa60 or NatF) was recently identified as an unconventional N-terminal acetyltransferase (NAT) because it localizes to organelles, in particular the Golgi apparatus, and has a preference for acetylating N termini of the transmembrane proteins. This knowledge challenged the prevailing view of N-terminal acetylation as a co-translational ribosome-associated process and suggested a new mechanistic functioning for the enzymes responsible for this increasingly recognized protein modification. Crystallography studies on Naa60 were unable to resolve the C-terminal tail of Naa60, which is responsible for the organellar localization. Here, we combined modeling, in vitro assays, and cellular localization studies to investigate the secondary structure and membrane interacting capacity of Naa60. The results show that Naa60 is a peripheral membrane protein. Two amphipathic helices within the Naa60 C terminus bind the membrane directly in a parallel position relative to the lipid bilayer via hydrophobic and electrostatic interactions. A peptide corresponding to the C terminus was unstructured in solution and only folded into an α-helical conformation in the presence of liposomes. Computational modeling and cellular mutational analysis revealed the hydrophobic face of two α-helices to be critical for membranous localization. Furthermore, we found a strong and specific binding preference of Naa60 toward membranes containing the phosphatidylinositol PI(4)P, thus possibly explaining the primary residency of Naa60 at the PI(4)P-rich Golgi. In conclusion, we have defined the mode of cytosolic Naa60 anchoring to the Golgi apparatus, most likely occurring post-translationally and specifically facilitating post-translational N-terminal acetylation of many transmembrane proteins.

A Role for Weak Electrostatic Interactions in Peripheral Membrane Protein Binding

Bacillus thuringiensis phosphatidylinositol-specific phospholipase C (BtPI-PLC) is a secreted virulence factor that binds specifically to phosphatidylcholine (PC) bilayers containing negatively charged phospholipids. BtPI-PLC carries a negative net charge and its interfacial binding site has no obvious cluster of basic residues. Continuum electrostatic calculations show that, as expected, nonspecific electrostatic interactions between BtPI-PLC and membranes vary as a function of the fraction of anionic lipids present in the bilayers. Yet they are strikingly weak, with a calculated ΔGel below 1 kcal/mol, largely due to a single lysine (K44). When K44 is mutated to alanine, the equilibrium dissociation constant for small unilamellar vesicles increases more than 50 times (∼2.4 kcal/mol), suggesting that interactions between K44 and lipids are not merely electrostatic. Comparisons of molecular-dynamics simulations performed using different lipid compositions reveal that the bilayer composition does not affect either hydrogen bonds or hydrophobic contacts between the protein interfacial binding site and bilayers. However, the occupancies of cation-π interactions between PC choline headgroups and protein tyrosines vary as a function of PC content. The overall contribution of basic residues to binding affinity is also context dependent and cannot be approximated by a rule-of-thumb value because these residues can contribute to both nonspecific electrostatic and short-range protein-lipid interactions. Additionally, statistics on the distribution of basic amino acids in a data set of membrane-binding domains reveal that weak electrostatics, as observed for BtPI-PLC, might be a less unusual mechanism for peripheral membrane binding than is generally thought.

Improving the Force Field Description of Tyrosine-Choline Cation-π Interactions: QM Investigation of Phenol-N(Me)4+ Interactions

Cation-π interactions between tyrosine amino acids and compounds containing N,N,N-trimethylethanolammonium (N(CH3)3) are involved in the recognition of histone tails by chromodomains and in the recognition of phosphatidylcholine (PC) phospholipids by membrane-binding proteins. Yet, the lack of explicit polarization or charge transfer effects in molecular mechanics force fields raises questions about the reliability of the representation of these interactions in biomolecular simulations. Here, we investigate the nature of phenol-tetramethylammonium (TMA) interactions using quantum mechanical (QM) calculations, which we also use to evaluate the accuracy of the additive CHARMM36 and Drude polarizable force fields in modeling tyrosine-choline interactions. We show that the potential energy surface (PES) obtained using SAPT2+/aug-cc-pVDZ compares well with the large basis-set CCSD(T) PES when TMA approaches the phenol ring perpendicularly. Furthermore, the SAPT energy decomposition reveals comparable contributions from electrostatics and dispersion in phenol-TMA interactions. We then compared the SAPT2+/aug-cc-pVDZ PES obtained along various approach directions to the corresponding PES obtained with CHARMM, and we show that the force field accurately reproduces the minimum distances while the interaction energies are underestimated. The use of the Drude polarizable force field significantly improves the interaction energies but decreases the agreement on distances at energy minima. The best agreement between force field and QM PES is obtained by modifying the Lennard-Jones terms for atom pairs involved in the phenol-TMA cation-π interactions. This is further shown to improve the correlation between the occupancy of tyrosine-choline cation-π interactions obtained from molecular dynamics simulations of a bilayer-bound bacterial phospholipase and experimental affinity data of the wild-type protein and selected mutants.

Similarity in Shape Dictates Signature Intrinsic Dynamics Despite No Functional Conservation in TIM Barrel Enzymes

The conservation of the intrinsic dynamics of proteins emerges as we attempt to understand the relationship between sequence, structure and functional conservation. We characterise the conservation of such dynamics in a case where the structure is conserved but function differs greatly. The triosephosphate isomerase barrel fold (TBF), renowned for its 8 β-strand-α-helix repeats that close to form a barrel, is one of the most diverse and abundant folds found in known protein structures. Proteins with this fold have diverse enzymatic functions spanning five of six Enzyme Commission classes, and we have picked five different superfamily candidates for our analysis using elastic network models. We find that the overall shape is a large determinant in the similarity of the intrinsic dynamics, regardless of function. In particular, the β-barrel core is highly rigid, while the α-helices that flank the β-strands have greater relative mobility, allowing for the many possibilities for placement of catalytic residues. We find that these elements correlate with each other via the loops that link them, as opposed to being directly correlated. We are also able to analyse the types of motions encoded by the normal mode vectors of the α-helices. We suggest that the global conservation of the intrinsic dynamics in the TBF contributes greatly to its success as an enzymatic scaffold both through evolution and enzyme design.

Proteinase 3 Is a Phosphatidylserine-binding Protein That Affects the Production and Function of Microvesicles

Proteinase 3 (PR3), the autoantigen in granulomatosis with polyangiitis, is expressed at the plasma membrane of resting neutrophils, and this membrane expression increases during both activation and apoptosis. Using surface plasmon resonance and protein-lipid overlay assays, this study demonstrates that PR3 is a phosphatidylserine-binding protein and this interaction is dependent on the hydrophobic patch responsible for membrane anchorage. Molecular simulations suggest that PR3 interacts with phosphatidylserine via a small number of amino acids, which engage in long lasting interactions with the lipid heads. As phosphatidylserine is a major component of microvesicles (MVs), this study also examined the consequences of this interaction on MV production and function. PR3-expressing cells produced significantly fewer MVs during both activation and apoptosis, and this reduction was dependent on the ability of PR3 to associate with the membrane as mutating the hydrophobic patch restored MV production. Functionally, activation-evoked MVs from PR3-expressing cells induced a significantly larger respiratory burst in human neutrophils compared with control MVs. Conversely, MVs generated during apoptosis inhibited the basal respiratory burst in human neutrophils, and those generated from PR3-expressing cells hampered this inhibition. Given that membrane expression of PR3 is increased in patients with granulomatosis with polyangiitis, MVs generated from neutrophils expressing membrane PR3 may potentiate oxidative damage of endothelial cells and promote the systemic inflammation observed in this disease.

Evolution of oligomeric state through allosteric pathways that mimic ligand binding

Evolution and design of protein complexes are almost always viewed through the lens of amino acid mutations at protein interfaces. We showed previously that residues not involved in the physical interaction between proteins make important contributions to oligomerization by acting indirectly or allosterically. In this work, we sought to investigate the mechanism by which allosteric mutations act, using the example of the PyrR family of pyrimidine operon attenuators. In this family, a perfectly sequence-conserved helix that forms a tetrameric interface is exposed as solvent-accessible surface in dimeric orthologs. This means that mutations must be acting from a distance to destabilize the interface. We identified 11 key mutations controlling oligomeric state, all distant from the interfaces and outside ligand-binding pockets. Finally, we show that the key mutations introduce conformational changes equivalent to the conformational shift between the free versus nucleotide-bound conformations of the proteins.

WEBnm@ v2.0: Web server and services for comparing protein flexibility

Background

Normal mode analysis (NMA) using elastic network models is a reliable and cost-effective computational method to characterise protein flexibility and by extension, their dynamics. Further insight into the dynamics–function relationship can be gained by comparing protein motions between protein homologs and functional classifications. This can be achieved by comparing normal modes obtained from sets of evolutionary related proteins.

Results

We have developed an automated tool for comparative NMA of a set of pre-aligned protein structures. The user can submit a sequence alignment in the FASTA format and the corresponding coordinate files in the Protein Data Bank (PDB) format. The computed normalised squared atomic fluctuations and atomic deformation energies of the submitted structures can be easily compared on graphs provided by the web user interface. The web server provides pairwise comparison of the dynamics of all proteins included in the submitted set using two measures: the Root Mean Squared Inner Product and the Bhattacharyya Coefficient. The Comparative Analysis has been implemented on our web server for NMA, WEBnm@, which also provides recently upgraded functionality for NMA of single protein structures. This includes new visualisations of protein motion, visualisation of inter-residue correlations and the analysis of conformational change using the overlap analysis. In addition, programmatic access to WEBnm@ is now available through a SOAP-based web service. Webnm@ is available at http://apps.cbu.uib.no/webnma.

Conclusion

WEBnm@ v2.0 is an online tool offering unique capability for comparative NMA on multiple protein structures. Along with a convenient web interface, powerful computing resources, and several methods for mode analyses, WEBnm@ facilitates the assessment of protein flexibility within protein families and superfamilies. These analyses can give a good view of how the structures move and how the flexibility is conserved over the different structures.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0427-6) contains supplementary material, which is available to authorized users.

Quantifying transient interactions between Bacillus phosphatidylinositol-specific phospholipase-C and phosphatidylcholine-rich vesicles

Bacillus thuringiensis secretes the virulence factor phosphatidylinositol-specific phospholipase C (BtPI-PLC), which specifically binds to phosphatidylcholine (PC) and cleaves GPI-anchored proteins off eukaryotic plasma membranes. To elucidate how BtPI-PLC searches for GPI-anchored proteins on the membrane surface, we measured residence times of single fluorescently labeled proteins on PC-rich small unilamellar vesicles (SUVs). BtPI-PLC interactions with the SUV surface are transient with a lifetime of 379 ± 49 ms. These data also suggest that BtPI-PLC does not directly sense curvature, but rather prefers to bind to the numerous lipid packing defects in SUVs. Despite this preference for defects, all-atom molecular dynamics simulations of BtPI-PLC interacting with PC-rich bilayers show that the protein is shallowly anchored with the deepest insertions ∼18 Å above the bilayer center. Membrane partitioning is mediated, on average, by 41 hydrophobic, 8 hydrogen-bonding, and 2 cation−π (between PC choline headgroups and Tyr residues) transient interactions with phospholipids. These results lead to a quantitative model for BtPI-PLC interactions with cell membranes where protein binding is mediated by lipid packing defects, possibly near GPI-anchored proteins, and the protein diffuses on the membrane for ∼100–380 ms, during which time it may cleave ∼10 GPI-anchored proteins before dissociating. This combination of short two-dimensional scoots followed by three-dimensional hops may be an efficient search strategy on two-dimensional surfaces with obstacles.

Biochemical and cellular analysis of Ogden syndrome reveals downstream Nt-acetylation defects

The X-linked lethal Ogden syndrome was the first reported human genetic disorder associated with a mutation in an N-terminal acetyltransferase (NAT) gene. The affected males harbor an Ser37Pro (S37P) mutation in the gene encoding Naa10, the catalytic subunit of NatA, the major human NAT involved in the co-translational acetylation of proteins. Structural models and molecular dynamics simulations of the human NatA and its S37P mutant highlight differences in regions involved in catalysis and at the interface between Naa10 and the auxiliary subunit hNaa15. Biochemical data further demonstrate a reduced catalytic capacity and an impaired interaction between hNaa10 S37P and Naa15 as well as Naa50 (NatE), another interactor of the NatA complex. N-Terminal acetylome analyses revealed a decreased acetylation of a subset of NatA and NatE substrates in Ogden syndrome cells, supporting the genetic findings and our hypothesis regarding reduced Nt-acetylation of a subset of NatA/NatE-type substrates as one etiology for Ogden syndrome. Furthermore, Ogden syndrome fibroblasts display abnormal cell migration and proliferation capacity, possibly linked to a perturbed retinoblastoma pathway. N-Terminal acetylation clearly plays a role in Ogden syndrome, thus revealing the in vivo importance of N-terminal acetylation in human physiology and disease.

Reversible ketomethylene-based inhibitors of human neutrophil proteinase 3

Neutrophil serine proteases, proteinase 3 (PR3) and human neutrophil elastase (HNE), are considered as targets for chronic inflammatory diseases. Despite sharing high sequence similarity, the two enzymes have different substrate specificities and functions. While a plethora of HNE inhibitors exist, PR3 specific inhibitors are still in their infancy. We have designed ketomethylene-based inhibitors for PR3 that show low micromolar IC50 values. Their synthesis was made possible by amending a previously reported synthesis of ketomethylene dipeptide isosteres to allow for the preparation of derivatives suitable for solid phase peptide synthesis. The best inhibitor (Abz-VADnV[Ψ](COCH2)ADYQ-EDDnp) was found to be selective for PR3 over HNE and to display a competitive and reversible inhibition mechanism. Molecular dynamics simulations show that the interactions between enzyme and ketomethylene-containing inhibitors are similar to those with the corresponding substrates. We also confirm that N- and C-terminal FRET groups are important for securing high inhibitory potency toward PR3.

Comparing the intrinsic dynamics of multiple protein structures using elastic network models

BACKGROUND

Elastic network models (ENMs) are based on the simple idea that a protein can be described as a set of particles connected by springs, which can then be used to describe its intrinsic flexibility using, for example, normal mode analysis. Since the introduction of the first ENM by Monique Tirion in 1996, several variants using coarser protein models have been proposed and their reliability for the description of protein intrinsic dynamics has been widely demonstrated. Lately an increasing number of studies have focused on the meaning of slow dynamics for protein function and its potential conservation through evolution. This leads naturally to comparisons of the intrinsic dynamics of multiple protein structures with varying levels of similarity.

SCOPE OF REVIEW

We describe computational strategies for calculating and comparing intrinsic dynamics of multiple proteins using elastic network models, as well as a selection of examples from the recent literature.

MAJOR CONCLUSIONS

The increasing interest for comparing dynamics across protein structures with various levels of similarity, has led to the establishment and validation of reliable computational strategies using ENMs. Comparing dynamics has been shown to be a viable way for gaining greater understanding for the mechanisms employed by proteins for their function. Choices of ENM parameters, structure alignment or similarity measures will likely influence the interpretation of the comparative analysis of protein motion.

GENERAL SIGNIFICANCE

Understanding the relation between protein function and dynamics is relevant to the fundamental understanding of protein structure-dynamics-function relationship. This article is part of a Special Issue entitled Recent developments of molecular dynamics.

Evaluation of Protein Elastic Network Models Based on an Analysis of Collective Motions

Elastic network models (ENMs) are valuable tools for investigating collective motions of proteins, and a rich variety of simple models have been proposed over the past decade. A good representation of the collective motions requires a good approximation of the covariances between the fluctuations of the individual atoms. Nevertheless, most studies have validated such models only by the magnitudes of the single-atom fluctuations they predict. In the present study, we have quantified the agreement between the covariance structure predicted by molecular dynamics (MD) simulations and those predicted by a representative selection of proposed coarse-grained ENMs. We then contrast this approach with the comparison to MD-predicted atomic fluctuations and comparison to crystallographic B-factors. While all the ENMs yield approximations to the MD-predicted covariance structure, we report large and consistent differences between proposed models. We also find that the ability of the ENMs to predict atomic fluctuations is correlated with their ability to capture the covariance structure. In contrast, we find that the models that agree best with B-factors model collective motions less reliably and recommend against using B-factors as a benchmark.

Visual cavity analysis in molecular simulations

Molecular surfaces provide a useful mean for analyzing interactions between biomolecules; such as identification and characterization of ligand binding sites to a host macromolecule. We present a novel technique, which extracts potential binding sites, represented by cavities, and characterize them by 3D graphs and by amino acids. The binding sites are extracted using an implicit function sampling and graph algorithms. We propose an advanced cavity exploration technique based on the graph parameters and associated amino acids. Additionally, we interactively visualize the graphs in the context of the molecular surface. We apply our method to the analysis of MD simulations of Proteinase 3, where we verify the previously described cavities and suggest a new potential cavity to be studied.

Cation-π interactions as lipid-specific anchors for phosphatidylinositol-specific phospholipase C

Amphitropic proteins, such as the virulence factor phosphatidylinositol-specific phospholipase C (PI-PLC) from Bacillus thuringiensis , often depend on lipid-specific recognition of target membranes. However, the recognition mechanisms for zwitterionic lipids, such as phosphatidylcholine, which is enriched in the outer leaflet of eukaryotic cells, are not well understood. A 500 ns long molecular dynamics simulation of PI-PLC at the surface of a lipid bilayer revealed a strikingly high number of interactions between tyrosines at the interfacial binding site and lipid choline groups with structures characteristic of cation-π interactions. Membrane affinities of PI-PLC tyrosine variants mostly tracked the simulation results, falling into two classes: (i) those with minor losses in affinity, Kd(mutant)/Kd(wild-type) ≤ 5 and (ii) those where the apparent Kd was 50-200 times higher than wild-type. Estimating ΔΔG for these Tyr/PC interactions from the apparent Kd values reveals that the free energy associated with class I is ~1 kcal/mol, comparable to the value predicted by the Wimley-White hydrophobicity scale. In contrast, removal of class II tyrosines has a higher energy cost: ~2.5 kcal/mol toward pure PC vesicles. These higher energies correlate well with the occupancy of the cation-π adducts throughout the MD simulation. Together, these results strongly indicate that PI-PLC interacts with PC headgroups via cation-π interactions with tyrosine residues and suggest that cation-π interactions at the interface may be a mechanism for specific lipid recognition by amphitropic and membrane proteins.

Does changing the predicted dynamics of a phospholipase C alter activity and membrane binding?

The enzymatic activity of secreted phosphatidylinositol-specific phospholipase C (PI-PLC) enzymes is associated with bacterial virulence. Although the PI-PLC active site has no obvious lid, molecular-dynamics simulations suggest that correlated loop motions may limit access to the active site, and two Pro residues, Pro(245) and Pro(254), are associated with these correlated motions. Whereas the region containing both Pro residues is quite variable among PI-PLCs, it shows high conservation in virulence-associated, secreted PI-PLCs that bind to the surface of cells. These regions of the protein are also associated with phosphatidylcholine binding, which enhances PI-PLC activity. In silico mutagenesis of Pro(245) disrupts correlated motions between the two halves of Bacillus thuringiensis PI-PLC, and Pro(245) variants show significantly reduced enzymatic activity in all assay systems. PC still enhanced activity, but not to the level of wild-type enzyme. Mutagenesis of Pro(254) appears to stiffen the PI-PLC structure, but experimental mutations had minor effects on activity and membrane binding. With the exception of P245Y, reduced activity was not associated with reduced membrane affinity. This combination of simulations and experiments suggests that correlated motions between the two halves of PI-PLC may be more important for enzymatic activity than for vesicle binding.