Structural properties of POPC monolayers under lateral compression: computer simulations analysis.

1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), a lipid comprising a saturated and an unsaturated acyl chain, belongs to the class of glycerophosphatidylcholines, major lipids in eukaryotic cell membranes. To get insight into the structural properties of this lipid within monolayers as membrane models, we performed molecular dynamics (MD) simulations of POPC monolayers under compression at the air/water interface. MD simulations were carried out at 300 K and at different surface pressures using the all-atom general Amber force field (GAFF). A good agreement was found between the simulated data and experimental isotherms. At surface pressures greater than 15 mN/m, two orientations of the head groups clearly appear: one nearly parallel to the monolayer interface and another one pointing toward the water. On the basis of the analysis of headgroup dihedral angles, we propose that the conformational variations around the bonds connecting the phosphorus atom to the adjacent oxygens are involved in these two orientations of the headgroup. The glycerol group orientation is characterized by a large distribution centered around 50° with respect to the monolayer normal. The acyl chains are predominantly in trans configuration from 7.5 to 43 mN/m surface pressures. Moreover, the calculated order parameter profiles of both chains suggest an independent behavior of the saturated and unsaturated chains that could be correlated with the formation of chain-type clusters observed along the simulated trajectories.

Membrane interface composition drives the structure and the tilt of the single transmembrane helix protein PMP1: MD studies.

PMP1, a regulatory subunit of the yeast plasma membrane H(+)-ATPase, is a single transmembrane helix protein. Its cytoplasmic C-terminus possesses several positively charged residues and interacts with phosphatidylserine lipids as shown through both (1)H- and (2)H-NMR experiments. We used all-atom molecular dynamics simulations to obtain atomic-scale data on the effects of membrane interface lipid composition on PMP1 structure and tilt. PMP1 was embedded in two hydrated bilayers, differing in the composition of the interfacial region. The neutral bilayer is composed of POPC (1-palmitoyl-2-oleoyl-3-glycero-phosphatidylcholine) lipids and the negatively charged bilayer is composed of POPC and anionic POPS (1-palmitoyl-2-oleoyl-3-glycero-phosphatidylserine) lipids. Our results were consistent with NMR data obtained previously, such as a lipid sn-2 chain lying on the W28 aromatic ring and in the groove formed on one side of the PMP1 helix. In pure POPC, the transmembrane helix is two residues longer than the initial structure and the helix tilt remains constant at 6 ± 3°. By contrast, in mixed POPC-POPS, the initial helical structure of PMP1 is stable throughout the simulation time even though the C-terminal residues interact strongly with POPS headgroups, leading to a significant increase of the helix tilt within the membrane to 20 ± 5°.

Cell uptake of a biosensor detected by hyperpolarized 129Xe NMR: the transferrin case.

For detection of biological events in vitro, sensors using hyperpolarized (129)Xe NMR can become a powerful tool, provided the approach can bridge the gap in sensitivity. Here we propose constructs based on the non-selective grafting of cryptophane precursors on holo-transferrin. This biological system was chosen because there are many receptors on the cell surface, and endocytosis further increases this density. The study of these biosensors with K562 cell suspensions via fluorescence microscopy and (129)Xe NMR indicates a strong interaction, as well as interesting features such as the capacity of xenon to enter the cryptophane even when the biosensor is endocytosed, while keeping a high level of polarization. Despite a lack of specificity for transferrin receptors, undoubtedly due to the hydrophobic character of the cryptophane moiety that attracts the biosensor into the cell membrane, these biosensors allow the first in-cell probing of biological events using hyperpolarized xenon.

Alteration of the tertiary structure of the major bee venom allergen Api m 1 by multiple mutations is concomitant with low IgE reactivity.

We have engineered a recombinant form of the major bee venom allergen (Api m 1) with the final goal of reducing its IgE reactivity. This molecule (Api mut) contains 24 mutations and one deletion of 10 amino acids. The successive introduction of these sequence modifications led to a progressive loss of specific IgE and IgG reactivity and did not reveal any immunodominant epitopes. However, Api mut exhibited a clear loss of reactivity for Api m 1-specific IgE and IgG. Injection of Api mut into mice induced specific antibody production. This humoral response was as high as that induced by the Api m 1 but the cross-reactivity of the antibodies was weak. As inferred by far UV circular dichroism, this mutant was correctly folded. However, near UV circular dichroism and denaturation curves of Api mut showed that it exhibits a dynamic tertiary structure and that it is a highly flexible molecule. Finally, as all the sequence modifications have been introduced outside the human and murine T cell epitope regions, we investigated its T cell properties in mice. We showed that Api mut-specific T lymphocytes induced in vivo were stimulated in vitro by both proteins. These data provide new insights in the design of hypoallergenic molecules.

Most of the structural elements of the globular domain of murine prion protein form fibrils with predominant beta-sheet structure.

The conversion of the cellular prion protein into the beta-sheet-rich scrapie prion protein is thought to be the key step in the pathogenesis of prion diseases. To gain insight into this structural conversion, we analyzed the intrinsic structural propensity of the amino acid sequence of the murine prion C-terminal domain. For that purpose, this globular domain was dissected into its secondary structural elements and the structural propensity of the protein fragments was determined. Our results show that all these fragments, excepted that strictly encompassing helix 1, have a very high propensity to form structured aggregates with a dominant content of beta-sheet structures.

New force field for calcium binding sites in annexin-membrane complexes.

For accurate classical molecular dynamics (MD) simulations of the calcium mediated bound complexes of annexin and membrane we have developed new force-field parameters correctly describing the interaction of the Ca ion with its environment. We have used quantum chemical calculations to investigate the potential energy surface experienced by the Ca ion within the three different binding sites found in domain 1 of annexin V (ANX V/1). Based on these calculations we were able to quantify the charge polarization of atoms within the binding sites, and to determine the geometry and force constants of harmonic restraints between the Ca ion and its coordinating oxygen atoms. Harmonic restraints were introduced to compensate for the deviations between the quantum mechanical potential energy surface and that of the classical force field. Our analysis has shown that using the refined force field for the Ca binding sites enables long-time MD simulations that conserve the initial structure of ANX V/1 significantly better than MD simulations using the standard force field.

Molecular dynamics simulations of E. coli MsbA transmembrane domain: formation of a semipore structure.

The human P-glycoprotein (MDR1/P-gp) is an ATP-binding cassette (ABC) transporter involved in cellular response to chemical stress and failures of anticancer chemotherapy. In the absence of a high-resolution structure for P-gp, we were interested in the closest P-gp homolog for which a crystal structure is available: the bacterial ABC transporter MsbA. Here we present the molecular dynamics simulations performed on the transmembrane domain of the open-state MsbA in a bilayer composed of palmitoyl oleoyl phosphatidylethanolamine lipids. The system studied contained more than 90,000 atoms and was simulated for 50 ns. This simulation shows that the open-state structure of MsbA can be stable in a membrane environment and provides invaluable insights into the structural relationships between the protein and its surrounding lipids. This study reveals the formation of a semipore-like structure stabilized by two key phospholipids which interact with the hinge region of the protein during the entire simulation. Multiple sequence alignments of ABC transporters reveal that one of the residues involved in the interaction with these two phospholipids are under a strong selection pressure specifically applied on the bacterial homologs of MsbA. Hence, comparison of molecular dynamics simulation and phylogenetic data appears as a powerful approach to investigate the functional relevance of molecular events occurring during simulations.

Protein unfolding transitions in an intrinsically unstable annexin domain: molecular dynamics simulation and comparison with nuclear magnetic resonance data.

Unfolding transitions of an intrinsically unstable annexin domain and the unfolded state structure have been examined using multiple approximately 10-ns molecular dynamics simulations. Three main basins are observed in the configurational space: native-like state, compact partially unfolded or intermediate compact state, and the unfolded state. In the native-like state fluctuations are observed that are nonproductive for unfolding. During these fluctuations, after an initial loss of approximately 20% of the core residue native contacts, the core of the protein transiently completely refolds to the native state. The transition from the native-like basin to the partially unfolded compact state involves approximately 75% loss of native contacts but little change in the radius of gyration or core hydration properties. The intermediate state adopts for part of the time in one of the trajectories a novel highly compact salt-bridge stabilized structure that can be identified as a conformational trap. The intermediate-to-unfolded state transition is characterized by a large increase in the radius of gyration. After an initial relaxation the unfolded state recovers a native-like topology of the domain. The simulated unfolded state ensemble reproduces in detail experimental nuclear magnetic resonance data and leads to a convincing complete picture of the unfolded domain.

Deciphering the role of individual acyl chains in the interaction network between phosphatidylserines and a single-spanning membrane protein.

PMP1 is a small single-spanning membrane protein functioning as a regulatory subunit of the yeast plasma membrane H(+)-ATPase. This protein forms a unique helix and exhibits a positively charged cytoplasmic domain that is able to specifically segregate phosphatidylserines (PSs). A marked groove formed at the helix surface is thought to play a major role in the related lipid-protein interaction network. Mutational analysis and (1)H NMR experiments were therefore performed on a synthetic PMP1 fragment using DPC-d(38) micelles as a membrane-like environment, in the presence of small amounts of POPS. A mutation designed for altering the helix groove was shown to disfavor the POPS binding specificity as much as that affecting the electrostatic interaction network. From POPS titration experiments monitored by a full set of one- and two-dimensional NOESY spectra, the association between the phospholipids and the PMP1 peptide has been followed. Our data reveal that the clustering of POPS molecules is promoted from a stabilized framework obtained by coupling the PMP1 helix groove to a POPS sn-2 chain. To our knowledge, the NOE-based titration plots displayed in this report constitute the first NMR data that directly distinguish the role of the sn-1 and sn-2 acyl chains in a lipid-protein interaction. The results are discussed while taking into account our accurate knowledge of the yeast plasma membrane composition and its ability to form functional lipid rafts.

A new consensus sequence for phosphatidylserine recognition by annexins.

Annexins are abundant and ubiquitous proteins that bind, by their four structurally identical domain cores, to phosphatidylserine-containing membranes in the presence of Ca2+. Using molecular simulation and mutagenesis, we have identified a new phosphatidylserine-binding site in annexin V domain 1 and established its structure. The residues involved in this site constitute a consensus sequence highly conserved in all annexins. Remarkably, this consensus sequence is exclusively found in domains 1 or 2, sometimes in both, but never in domains 3 and 4. Such a pattern actually delineates three classes of annexins, shedding new light on the role played by the four-domain core of annexins that could encode specific information discriminating the different annexins that compete within a given cell for membrane binding. Our findings thus provide new strategies for understanding the regulation of the cellular functions of annexins.

HLA-DP4, the most frequent HLA II molecule, defines a new supertype of peptide-binding specificity.

Among HLA-DP specificities, HLA-DP4 specificity involves at least two molecules, HLA-DPA1*0103/DPB1*0401 (DP401) and HLA-DPA1*0103/DPB1*0402 (DP402), which differ from each other by only three residues. Together, they are present worldwide at an allelic frequency of 20-60% and are the most abundant human HLA II alleles. Strikingly, the peptide-binding specificities of these molecules have never been investigated. Hence, in this study, we report the peptide-binding motifs of both molecules. We first set up a binding assay specific for the immunopurified HLA-DP4 molecules. Using multiple sets of synthetic peptides, we successfully defined the amino acid preferences of the anchor residues. With these assays, we were also able to identify new peptide ligands from allergens and viral and tumor Ags. DP401 and DP402 exhibit very similar patterns of recognition in agreement with molecular modeling of the complexes. Pockets P1 and P6 accommodate the main anchor residues and interestingly contain only two polymorphic residues, beta86 and beta11, respectively. Both positions are almost dimorphic and thus produce a limited number of pocket combinations. Taken together, our results support the existence of three main binding supertypes among HLA-DP molecules and should significantly contribute to the identification of universal epitopes to be used in peptide-based vaccines for cancer, as well as for allergic or infectious diseases.