Critical review of general guidelines for membrane proteins model building and analysis

Molecular Modeling of Transporters: From Low Resolution Cryo-Electron Microscopy Map to Conformational Exploration. The Example of TSPO

This chapter describes a protocol to establish a three-dimensional (3D) model of a protein and to explore its conformational landscape. It combines predictions from up-to-date bioinformatics methods with low-resolution experimental data. It also proposes to examine rapidly the dynamics of the protein using molecular dynamics simulations with a coarse-grained force field. Tools for analyzing these trajectories are suggested as well as those for constructing all-atoms models. Thus, starting from a protein sequence and using free software, the user can get important conformational information, which might improve the knowledge about the protein function.

Identification of the intracellular gate for a member of the equilibrative nucleoside transporter (ENT) family

Equilibrative nucleoside transporters of the SLC29 family play important roles in many physiological and pharmacological processes, including import of drugs for treatment of cancer, AIDS, cardiovascular, and parasitic diseases. However, no crystal structure is available for any member of this family. In previous studies we generated a computational model of the Leishmania donovani nucleoside transporter 1.1 (LdNT1.1) that captured this permease in the outward-closed conformation, and we identified the extracellular gate. In the present study we have modeled the inward-closed conformation of LdNT1.1 using the crystal structure of the Escherichia coli fucose transporter FucP and have identified four transmembrane helices whose ends close to form a predicted intracellular gate. We have tested this prediction by site-directed mutagenesis of relevant helix residues and by cross-linking of introduced cysteine pairs. The results are consistent with the predictions of the computational model and suggest that a similarly constituted gate operates in other members of the equilibrative nucleoside transporter family.

Membrane Protein Simulations Using AMBER Force Field and Berger Lipid Parameters

AMBER force fields are among the most commonly used in molecular dynamics (MD) simulations of proteins. Unfortunately, they lack a specific set of lipid parameters, thus limiting its use in membrane protein simulations. In order to overcome this limitation we assessed whether the widely used united-atom lipid parameters described by Berger and co-workers could be used in conjunction with AMBER force fields in simulations of membrane proteins. Thus, free energies of solvation in water and in cyclohexane, and free energies of water to cyclohexane transfer, were computed by thermodynamic integration procedures for neutral amino acid side-chains employing AMBER99, AMBER03, and OPLS-AA amino acid force fields. In addition, MD simulations of three membrane proteins in a POPC lipid bilayer, the β2 adrenergic G protein-coupled receptor, Aquaporin-1, and the outer membrane protein Omp32, were performed with the aim of comparing the AMBER99SB/Berger combination of force fields with the OPLS-AA/Berger combination. We have shown that AMBER99SB and Berger force fields are compatible, they provide reliable free energy estimations relative to experimental values, and their combination properly describes both membrane and protein structural properties. We then suggest that the AMBER99SB/Berger combination is a reliable choice for the simulation of membrane proteins, which links the easiness of ligand parametrization and the ability to reproduce secondary structure of AMBER99SB force field with the largely validated Berger lipid parameters.