Bh3 induced conformational changes in Bcl‐Xl revealed by crystal structure and comparative analysis

Apoptosis or programmed cell death is a regulatory process in cells in response to stimuli perturbing physiological conditions. The Bcl‐2 family of proteins plays an important role in regulating homeostasis during apoptosis. In the process, the molecular interactions among the three members of this family, the pro‐apoptotic, anti‐apoptotic and BH3‐only proteins at the mitochondrial outer membrane define the fate of a cell. Here, we report the crystal structures of the human anti‐apoptotic protein Bcl‐X_L in complex with BH3‐only BID^BH3 and BIM^BH3 peptides determined at 2.0 Å and 1.5 Å resolution, respectively. The BH3 peptides bind to the canonical hydrophobic pocket in Bcl‐X_L and adopt an alpha helical conformation in the bound form. Despite a similar structural fold, a comparison with other BH3 complexes revealed structural differences due to their sequence variations. In the Bcl‐X_L‐BID^BH3 complex we observed a large pocket, in comparison with other BH3 complexes, lined by residues from helices α1, α2, α3, and α5 located adjacent to the canonical hydrophobic pocket. These results suggest that there are differences in the mode of interactions by the BH3 peptides that may translate into functional differences in apoptotic regulation. Proteins 2015; 83:1262–1272. © 2015 Wiley Periodicals, Inc.     

Distance restraints from crosslinking mass spectrometry: Mining a molecular dynamics simulation database to evaluate lysine–lysine distances

Integrative structural biology attempts to model the structures of protein complexes that are challenging or intractable by classical structural methods (due to size, dynamics, or heterogeneity) by combining computational structural modeling with data from experimental methods. One such experimental method is chemical crosslinking mass spectrometry (XL‐MS), in which protein complexes are crosslinked and characterized using liquid chromatography‐mass spectrometry to pinpoint specific amino acid residues in close structural proximity. The commonly used lysine‐reactive N‐hydroxysuccinimide ester reagents disuccinimidylsuberate (DSS) and bis(sulfosuccinimidyl)suberate (BS³) have a linker arm that is 11.4 Å long when fully extended, allowing C_α (alpha carbon of protein backbone) atoms of crosslinked lysine residues to be up to ～24 Å apart. However, XL‐MS studies on proteins of known structure frequently report crosslinks that exceed this distance. Typically, a tolerance of ～3 Å is added to the theoretical maximum to account for this observation, with limited justification for the chosen value. We used the Dynameomics database, a repository of high‐quality molecular dynamics simulations of 807 proteins representative of diverse protein folds, to investigate the relationship between lysine–lysine distances in experimental starting structures and in simulation ensembles. We conclude that for DSS/BS³, a distance constraint of 26–30 Å between C_α atoms is appropriate. This analysis provides a theoretical basis for the widespread practice of adding a tolerance to the crosslinker length when comparing XL‐MS results to structures or in modeling. We also discuss the comparison of XL‐MS results to MD simulations and known structures as a means to test and validate experimental XL‐MS methods.     

Predicting the disruption by of a protein-ligand interaction

The uranyl cation (UO₂²⁺) can be suspected to interfere with the binding of essential metal cations to proteins, underlying some mechanisms of toxicity. A dedicated computational screen was used to identify UO₂²⁺ binding sites within a set of nonredundant protein structures. The list of potential targets was compared to data from a small molecules interaction database to pinpoint specific examples where UO₂²⁺ should be able to bind in the vicinity of an essential cation, and would be likely to affect the function of the corresponding protein. The C‐reactive protein appeared as an interesting hit since its structure involves critical calcium ions in the binding of phosphorylcholine. Biochemical experiments confirmed the predicted binding site for UO₂²⁺ and it was demonstrated by surface plasmon resonance assays that UO₂²⁺ binding to CRP prevents the calcium‐mediated binding of phosphorylcholine. Strikingly, the apparent affinity of UO₂²⁺ for native CRP was almost 100‐fold higher than that of Ca²⁺. This result exemplifies in the case of CRP the capability of our computational tool to predict effective binding sites for UO₂²⁺ in proteins and is a first evidence of calcium substitution by the uranyl cation in a native protein.     

Synthetic prions with novel strain-specified properties

Prions are infectious proteins that possess multiple self-propagating structures. The information for strains and structural specific barriers appears to be contained exclusively in the folding of the pathological isoform, PrP^Sc. Many recent studies determined that de novo prion strains could be generated in vitro from the structural conversion of recombinant (rec) prion protein (PrP) into amyloidal structures. Our aim was to elucidate the conformational diversity of pathological recPrP amyloids and their biological activities, as well as to gain novel insights in characterizing molecular events involved in mammalian prion conversion and propagation. To this end we generated infectious materials that possess different conformational structures. Our methodology for the prion conversion of recPrP required only purified rec full-length mouse (Mo) PrP and common chemicals. Neither infected brain extracts nor amplified PrP^Sc were used. Following two different in vitro protocols recMoPrP converted to amyloid fibrils without any seeding factor. Mouse hypothalamic GT1 and neuroblastoma N2a cell lines were infected with these amyloid preparations as fast screening methodology to characterize the infectious materials. Remarkably, a large number of amyloid preparations were able to induce the conformational change of endogenous PrP^C to harbor several distinctive proteinase-resistant PrP forms. One such preparation was characterized in vivo habouring a synthetic prion with novel strain specified neuropathological and biochemical properties.     

Conformational plasticity of the calcium‐binding pocket in the Burkholderia glumae lipase: Remodeling induced by mutation of calcium coordinating residues

Most bacterial lipases bind one or more Ca²⁺ atoms at different locations and are a suitable case of study for investigating structural effects related to calcium binding, depletion, or mutation of calcium‐binding sites. Generally Ca²⁺ in microbial lipases can play a crucial role in the stabilization of the whole three‐dimensional structure by mediating long‐range effects. It has been recently demonstrated that calcium binding influences thermal stability of Burkholderia glumae lipase (BGL) through the restriction of conformational plasticity of specific regions. Moreover, calcium depletion results in a highly cooperative protein unfolding, eliciting protein aggregation. To further shed light on molecular mechanisms and structural features connected to calcium binding in microbial lipases, we present a molecular dynamics investigation, based on multiple‐replica approach at different temperatures, of BGL mutants targeting the calcium‐binding site. It turns out that additional acidic residues, which are conserved in other microbial lipases, help in overcoming effects induced by mutation of D241 Ca²⁺‐coordinating residue, upon rearrangements induced in the calcium binding site. © 2010 Wiley Periodicals, Inc. Biopolymers 95: 117–126, 2011.     

Differential stability of the bovine prion protein upon urea unfolding

Prion diseases, or transmissible spongiform encephalopathies, are a group of infectious neurological diseases associated with the structural conversion of an endogenous protein (PrP) in the central nervous system. There are two major forms of this protein: the native and noninfectious cellular form, PrP^C; and the misfolded, infectious, and proteinase K‐resistant form, PrP^Sc. The C‐terminal domain of PrP^C is mainly α‐helical in structure, whereas PrP^Sc in known to aggregate into an assembly of β‐sheets, forming amyloid fibrils. To identify the regions of PrP^C potentially involved in the initial steps of the conversion to the infectious conformation, we have used high‐resolution NMR spectroscopy to characterize the stability and structure of bovine recombinant PrP^C (residues 121 to 230) during unfolding with the denaturant urea. Analysis of the 800 MHz ¹H NMR spectra reveals region‐specific information about the structural changes occurring upon unfolding. Our data suggest that the dissociation of the native β‐sheet of PrP^C is a primary step in the urea‐induced unfolding process, while strong hydrophobic interactions between helices α1 and α3, and between α2 and α3, stabilize these regions even at very high concentrations of urea.     

Protein surface roughness accounts for binding free energy of Plasmepsin II‐ligand complexes

The calculation of absolute binding affinities for protein‐inhibitor complexes remains as one of the main challenges in computational structure‐based ligand design. The present work explored the calculations of surface fractal dimension (as a measure of surface roughness) and the relationship with experimental binding free energies of Plasmepsin II complexes. Plasmepsin II is an attractive target for novel therapeutic compounds to treat malaria. However, the structural flexibility of this enzyme is a drawback when searching for specific inhibitors. Concerning that, we performed separate explicitly solvated molecular dynamics simulations using the available high‐resolution crystal structures of different Plasmepsin II complexes. Molecular dynamics simulations allowed a better approximation to systems dynamics and, therefore, a more reliable estimation of surface roughness. This constitutes a novel approximation in order to obtain more realistic values of fractal dimension, because previous works considered only x‐ray structures. Binding site fractal dimension was calculated considering the ensemble of structures generated at different simulation times. A linear relationship between binding site fractal dimension and experimental binding free energies of the complexes was observed within 20 ns. Previous studies of the subject did not uncover this relationship. Regression model, coined FD model, was built to estimate binding free energies from binding site fractal dimension values. Leave‐one‐out cross‐validation showed that our model reproduced accurately the absolute binding free energies for our training set (R² = 0.76; <|error|> =0.55 kcal/mol; SD_error = 0.19 kcal/mol). The fact that such a simple model may be applied raises some questions that are addressed in the article.     

New Vistas on the Recruiting of Ferrous Iron and Dioxygen by Ferritins: A Case Study of the Escherichia coli 24‐mer Ferritin by All‐Atom Molecular Dynamics in Aqueous Medium

It is shown here that Fe²⁺ and O₂ ligands are displaced from the ferroxidase center of the C1 four‐helix bundle of E. coli 24‐mer ferritin under molecular dynamics (MD) aided by a randomly oriented external force applied to the ligand. Under these conditions, ligand egress toward the external aqueous medium occurs preferentially from the same four‐helix bundle, in the case of O₂, or other bundle, in the case of Fe². Viewing ligand egress from the protein as the microscopic reverse of ligand influx into the protein under unbiased MD, these findings challenge current views that preferential gates for recruitment of Fe²⁺ are 3‐fold channels with human ferritin, or the short path from the ferroxidase center to H93 with bacterial ferritins.     

Ornithokinin (avian bradykinin) from the skin of the Chinese bamboo odorous frog,Odorrana versabilis

One of the most widespread and abundant families of pharmacologically active peptides in amphibian defensive skin secretions is the bradykinins and related peptides. Despite retaining certain primary structural attributes that assign them to this peptide family, bradykinins and related peptides are unique among amphibian skin peptides in that they exhibit a wide range of primary structural variations, post‐translational modifications and/or N‐terminal or C‐terminal extensions. Initially it was believed that their high degree of primary structural heterogeneity was reflective of random gene mutations within species, but latterly, there is an increasing body of evidence that the spectrum of structural modifications found within this peptide family is reflective of the vertebrate predator spectrum of individual species. Here we report the discovery of ornithokinin (avian bradykinin – Thr⁶, Leu⁸‐bradykinin) in the skin secretion of the Chinese bamboo odorous frog, Odorrana versabilis. Molecular cloning of its biosynthetic precursor‐encoding cDNA from a skin secretion‐derived cDNA library revealed a deduced open‐reading frame of 86 amino acid residues, encoding a single copy of ornithokinin towards its C‐terminus. The domain architecture of this ornithokinin precursor protein was consistent with that of a typical amphibian skin peptide and quite different to that of the ornithokininogen from chicken plasma. Ornithokinin was reported to induce hypotension in the chicken and to contract the chicken oviduct but to have no obvious effect on the rat uterus. However, in this study, synthetic ornithokinin was found to contract the rat ileum (EC₅₀ = 539 nM) and to increase contraction frequency in the rat uterus (EC₅₀ = 1.87 μM). Copyright © 2014 European Peptide Society and John Wiley & Sons, Ltd.     

A DNA‐promoted amyloid proteinopathy in Escherichia coli

Protein amyloids arise from the conformational conversion and assembly of a soluble protein into fibrilar aggregates with a crossed β‐sheet backbone. Amyloid aggregates are able to replicate by acting as a template for the structural transformation and accretion of further protein molecules. In physicochemical terms, amyloids arguably constitute the simplest self‐replicative macromolecular assemblies. Similarly to the mammalian proteins PrP and α‐synuclein, the winged‐helix dimerization (WH1) domain of the bacterial, plasmid‐encoded protein RepA can assemble into amyloid fibres upon binding to DNA in vitro. Here we report that a hyper‐amyloidogenic functional variant (A31V) of RepA, fused to a red fluorescent protein, causes an amyloid proteinopathy in Escherichia coli with the following features: (i) in the presence of multiple copies of the specific DNA sequence opsp, WH1(A31V) accumulates as cytoplasmatic inclusions segregated from the nucleoid; (ii) such aggregates are amyloid in nature; (iii) bacteria carrying the amyloid inclusions age, exhibiting a fivefold expanded generation time; (iv) before cytokinesis, small inclusions are assembled de novo and transferred to the daughter cells, in which transmission failures cure amyloidosis; and (v) in the absence of inducer DNA, purified cellular WH1(A31V) inclusions seed amyloid fibre growth in vitro from the soluble protein. RepA‐WH1 is a suitable bacterial model system for amyloid proteinopathies.     

An insight of early PrP-E200K aggregation by combined molecular dynamics/fragment molecular orbital approaches

Unveiling the events leading to the formation of prion particles is a nowadays challenge in the field of neurochemistry. Pathogenic mutants of prion protein (PrP) are characterized by both an intrinsic tendency to aggregation and scrapie conversion propensity. However, the question about a possible correlation between these two events lasts still unanswered. Here, a multilayered computational workflow was employed to investigate structure, stability, and molecular interaction properties of a dimer of PrP^C-E200K, a well-known mutant of the PrP that represents a reduced model of early aggregates of this protein. Based on the combination of molecular dynamics and quantum mechanical approaches, this study provided for an in depth insight of PrP^C-E200K dimer in terms of residue-residue interactions. Assembly hypotheses for the early aggregation of PrP^C-E200K are paved and compared with PrP^Sc models reported in the literature to find a structural link between early and late (scrapie) aggregates of this protein.     

The CBL–CIPK signaling module in plants: a mechanistic perspective

In a given environment, plants are constantly exposed to multitudes of stimuli. These stimuli are sensed and transduced to generate a diverse array of responses by several signal transduction pathways. Calcium (Ca²⁺) signaling is one such important pathway involved in transducing a large number of stimuli or signals in both animals and plants. Ca²⁺ engages a plethora of decoders to mediate signaling in plants. Among these groups of decoders, the sensor responder complex of calcineurin B‐like protein (CBL) and CBL‐interacting protein kinases (CIPKs) play a very significant role in transducing these signals. The signal transduction mechanism in most cases is phosphorylation events, but some structural role for the pair has also come to light recently. In this review, we discuss the structural nature of the sensor‐responder duo; their mechanism of substrate phosphorylation and also their structural role in modulating targets. Moreover, the mechanism of complex formation and mechanistic role of protein phosphatases with CBL–CIPK module has been mentioned. A comparison of CBL–CIPK with other decoders of Ca²⁺ signaling in plants also signifies the relatedness and diversity in signaling pathways. Further an attempt has been made to compare this aspect of Ca²⁺ signaling pathways in different plant species to develop a holistic understanding of conservation of stimulus–response‐coupling mediated by this Ca²⁺–CBL–CIPK module.     

Virtual analysis of structurally diverse synthetic analogs as inhibitors of snake venom secretory phospholipase A2

Due to the toxic pathophysiological role of snake venom phospholipase A₂ (PLA₂), its compelling limitations to anti‐venom therapy in humans and the need for alternative therapy foster considerable pharmacological interest towards search of PLA₂ specific inhibitors. In this study, an integrated approach involving homology modeling, molecular dynamics and molecular docking studies on VRV‐PL‐V (Vipera russellii venom phospholipase A₂ fraction—V) belonging to Group II‐B secretory PLA₂ from Daboia russelli pulchella is carried out in order to study the structure‐based inhibitor design. The accuracy of the model was validated using multiple computational approaches. The molecular docking study of this protein was undertaken using different classes of experimentally proven, structurally diverse synthetic inhibitors of secretory PLA₂ whose selection is based on IC₅₀ value that ranges from 25 μM to 100 μM. Estimation of protein–ligand contacts by docking analysis sheds light on the importance of His 47 and Asp 48 within the VRV‐PL‐V binding pocket as key residue for hydrogen bond interaction with ligands. Our virtual analysis revealed that compounds with different scaffold binds to the same active site region. ADME analysis was also further performed to filter and identify the best potential specific inhibitor against VRV‐PL‐V. Additionally, the e‐pharmacophore was generated for the best potential specific inhibitor against VRV‐PL‐V and reported here. The present study should therefore play a guiding role in the experimental design of VRV‐PL‐V inhibitors that may provide better therapeutic molecular models for PLA₂ recognition and anti‐ophidian activity. Copyright © 2015 John Wiley & Sons, Ltd.     

The role of conformational selection in the molecular recognition of the wild type and mutants XPA67‐80 peptides by ERCC1

Molecular recognition is a fundamental step in the coordination of biomolecular pathways. Understanding how recognition and binding occur between highly flexible protein domains is a complex task. The conformational selection theory provides an elegant rationalization of the recognition mechanism, especially valid in cases when unstructured protein regions are involved. The recognition of a poorly structured peptide, namely XPA_67‐80, by its target receptor ERCC1, falls in this challenging study category. The microsecond molecular dynamics (MD) simulations, discussed in this work, show that the conformational propensity of the wild type XPA_67‐80 peptide in solution supports conformational selection as the key mechanism driving its molecular recognition by ERCC1. Moreover, all the mutations of the XPA_67‐80 peptide studied here cause a significant increase of its conformational disorder, relative to the wild type. Comparison to experimental data suggests that the loss of the recognized structural motifs at the microscopic time scale can contribute to the critical decrease in binding observed for one of the mutants, further substantiating the key role of conformational selection in recognition. Ultimately, because of the high sequence identity and analogy in binding, it is conceivable that the conclusions of this study on the XPA_67‐80 peptide also apply to the ERCC1‐binding domain of the XPA protein. Proteins 2015; 83:1341–1351. © 2015 Wiley Periodicals, Inc.     

Molecular mechanism of β-sheet self-organization at water-hydrophobic interfaces

The capacity to form β‐sheet structure and to self‐organize into amyloid aggregates is a property shared by many proteins. Severe neurodegenerative pathologies such as Alzheimer's disease are thought to involve the interaction of amyloidogenic protein oligomers with neuronal membranes. To understand the experimentally observed catalysis of amyloid formation by lipid membranes and other water‐hydrophobic interfaces, we examine the physico‐chemical basis of peptide adsorption and aggregation in a model membrane using atomistic molecular simulations. Blocked octapeptides with simple, repetitive sequences, (Gly‐Ala)₄, and (Gly‐Val)₄, are used as models of β‐sheet‐forming polypeptide chains found in the core of amyloid fibrils. In the presence of an n‐octane phase mimicking the core of lipid membranes, the peptides spontaneously partition at the octane‐water interface. The adsorption of nonpolar sidechains displaces the peptides' conformational equilibrium from a heterogeneous ensemble characterized by a high degree of structural disorder toward a more ordered ensemble favoring β‐hairpins and elongated β‐strands. At the interface, peptides spontaneously aggregate and rapidly evolve β‐sheet structure on a 10 to 100 ns time scale, while aqueous aggregates remain amorphous. Catalysis of β‐sheet formation results from the combination of the hydrophobic effect and of reduced conformational entropy of the polypeptide chain. While the former drives interfacial partition and displaces the conformational equilibrium of monomeric peptides, the planar interface further facilitates β‐sheet organization by increasing peptide concentration and reducing the dimensionality of self‐assembly from three to two. These findings suggest a general mechanism for the formation of β‐sheets on the surface of globular proteins and for amyloid self‐organization at hydrophobic interfaces. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Copper (II) promotes the formation of soluble neurotoxic PrP oligomers in acidic environment

Prion diseases are classically considered to be “amyloid diseases” caused by the deposition of amyloid fibrils in the brain. Recent studies identified soluble oligomers of PrP (prion protein) in damaged neuronal tissue. However, the details of PrP oligomerization are still unclear. In this study, we demonstrate that Cu²⁺ plays a vital role in the formation of soluble PrP oligomers. A Cu²⁺‐triggered structural conversion of PrP (90–231) to a β‐sheet isoform in pH 5.0 buffer was revealed by circular dichroism spectra and fluorescence measurement. Soluble oligomers were isolated by size exclusion chromatography from experimental solutions, allowing atomic force microscopy to reveal their morphology. 3‐(4,5‐Dimethylthiazol‐2‐yl)‐2,5‐diphenyltetrazolium bromide (MTT) and flow cytometry assays demonstrated that oligomeric PrP induced significant damage in and apoptosis of neuroblastoma SK‐N‐SH cells. These results indicate that in an acidic environment, Cu²⁺ promotes the formation of neurotoxic soluble PrP oligomers. J. Cell. Biochem. 111: 627–633, 2010. © 2010 Wiley‐Liss, Inc.     

Equilibrium transitions between side‐chain conformations in leucine and isoleucine

Despite recent improvements in computational methods for protein design, we still lack a quantitative, predictive understanding of the intrinsic probabilities for amino acids to adopt particular side‐chain conformations. Surprisingly, this question has remained unsettled for many years, in part because of inconsistent results from different experimental approaches. To explicitly determine the relative populations of different side‐chain dihedral angles, we performed all‐atom hard‐sphere Langevin Dynamics simulations of leucine (Leu) and isoleucine (Ile) dipeptide mimetics with stereo‐chemical constraints and repulsive‐only steric interactions between non‐bonded atoms. We determine the relative populations of the different χ₁ and χ₂ dihedral angle combinations as a function of the backbone dihedral angles ? and ψ. We also propose, and test, a mechanism for inter‐conversion between the different side‐chain conformations. Specifically, we discover that some of the transitions between side‐chain dihedral angle combinations are very frequent, whereas others are orders of magnitude less frequent, because they require rare coordinated motions to avoid steric clashes. For example, to transition between different values of χ₂, the Leu side‐chain bond angles κ₁ and κ₂ must increase, whereas to transition in χ₁, the Ile bond angles λ₁ and λ₂ must increase. These results emphasize the importance of computational approaches in stimulating further experimental studies of the conformations of side‐chains in proteins. Moreover, our studies emphasize the power of simple steric models to inform our understanding of protein structure, dynamics, and design. Proteins 2015; 83:1488–1499. © 2015 Wiley Periodicals, Inc.     

Hierarchical domain‐motion analysis of conformational changes in sarcoplasmic reticulum Ca2+‐ATPase

Sarco(endo)plasmic reticulum Ca²⁺‐ATPase transports two Ca²⁺ per ATP‐hydrolyzed across biological membranes against a large concentration gradient by undergoing large conformational changes. Structural studies with X‐ray crystallography revealed functional roles of coupled motions between the cytoplasmic domains and the transmembrane helices in individual reaction steps. Here, we employed “Motion Tree (MT),” a tree diagram that describes a conformational change between two structures, and applied it to representative Ca²⁺‐ATPase structures. MT provides information of coupled rigid‐body motions of the ATPase in individual reaction steps. Fourteen rigid structural units, “common rigid domains (CRDs)” are identified from seven MTs throughout the whole enzymatic reaction cycle. CRDs likely act as not only the structural units, but also the functional units. Some of the functional importance has been newly revealed by the analysis. Stability of each CRD is examined on the morphing trajectories that cover seven conformational transitions. We confirmed that the large conformational changes are realized by the motions only in the flexible regions that connect CRDs. The Ca²⁺‐ATPase efficiently utilizes its intrinsic flexibility and rigidity to response different switches like ligand binding/dissociation or ATP hydrolysis. The analysis detects functional motions without extensive biological knowledge of experts, suggesting its general applicability to domain movements in other membrane proteins to deepen the understanding of protein structure and function. Proteins 2015; 83:746–756. © 2015 Wiley Periodicals, Inc.     

Unfolding of beta‐lactoglobulin on the surface of polystyrene nanoparticles: Experimental and computational approaches

Structural changes ensuing from the non‐covalent absorption of bovine beta‐lactoglobulin (BLG) on the surface of polystyrene nanoparticles were investigated by using spectroscopic approaches, by assessing the reactivity of specific residues, and by limited proteolysis/mass spectrometry. Also, the immunoreactivity of absorbed and free BLG was compared. All these approaches indicated substantial rearrangements of the protein structure in the absorbed state, in spite of the reported structural rigidity of BLG. Changes made evident by experimental measurements were confirmed by computational approaches. These indicate that adsorption‐related changes are most marked in the area between the main C‐terminal alpha helix and the beta‐barrel, and lead to full exposure of the thiol on Cys₁₂₁, consistent with experimental measurements. In the computational model of bound BLG, both Trp₆₁ and Trp₁₉ also move away from their neighboring quenchers and become solvent‐exposed, as indicated by fluorescence measurement. Upon binding, the beta‐barrel also loosens, with a substantial increase in immunoreactivity and with noticeable changes in the trypsinolytic pattern. The possible general significance of the structural changes reported here for non‐covalently adsorbed BLG is discussed with respect to recognition events involving surface‐bound proteins, as are aspects related to the carrier function(s) of BLG, and to its use as a common ingredient in many food systems. Proteins 2014; 82:1272–1282. © 2013 Wiley Periodicals, Inc.     

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 PIP₂ (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. Proteins 2015; 83:952–969. © 2015 Wiley Periodicals, Inc.