Mireia Olivella

Serine and Threonine Residues Bend α-Helices in the χ1 = g− Conformation

The relationship between the Ser, Thr, and Cys side-chain conformation (chi(1) = g(-), t, g(+)) and the main-chain conformation (phi and psi angles) has been studied in a selection of protein structures that contain alpha-helices. The statistical results show that the g(-) conformation of both Ser and Thr residues decreases their phi angles and increases their psi angles relative to Ala, used as a control. The additional hydrogen bond formed between the O(gamma) atom of Ser and Thr and the i-3 or i-4 peptide carbonyl oxygen induces or stabilizes a bending angle in the helix 3-4 degrees larger than for Ala. This is of particular significance for membrane proteins. Incorporation of this small bending angle in the transmembrane alpha-helix at one side of the cell membrane results in a significant displacement of the residues located at the other side of the membrane. We hypothesize that local alterations of the rotamer configurations of these Ser and Thr residues may result in significant conformational changes across transmembrane helices, and thus participate in the molecular mechanisms underlying transmembrane signaling. This finding has provided the structural basis to understand the experimentally observed influence of Ser residues on the conformational equilibrium between inactive and active states of the receptor, in the neurotransmitter subfamily of G protein-coupled receptors.

Activation of CCR5 by chemokines involves an aromatic cluster between transmembrane helices 2 and 3

CCR5 is a G protein-coupled receptor responding to four natural agonists, the chemokines RANTES (regulated on activation normal T cell expressed and secreted), macrophage inflammatory protein (MIP)-1α, MIP-1β, and monocyte chemotactic protein (MCP)-2, and is the main co-receptor for the macrophage-tropic human immunodeficiency virus strains. We have previously identified a structural motif in the second transmembrane helix of CCR5, which plays a crucial role in the mechanism of receptor activation. We now report the specific role of aromatic residues in helices 2 and 3 of CCR5 in this mechanism. Using site-directed mutagenesis and molecular modeling in a combined approach, we demonstrate that a cluster of aromatic residues at the extracellular border of these two helices are involved in chemokine-induced activation. These aromatic residues are involved in interhelical interactions that are key for the conformation of the helices and govern the functional response to chemokines in a ligand-specific manner. We therefore suggest that transmembrane helices 2 and 3 contain important structural elements for the activation mechanism of chemokine receptors, and possibly other related receptors as well.

Ser and Thr Residues Modulate the Conformation of Pro-Kinked Transmembrane α-Helices

Functionally required conformational plasticity of transmembrane proteins implies that specific structural motifs have been integrated in transmembrane helices. Surveying a database of transmembrane helices and the large family of G-protein coupled receptors we identified a series of overrepresented motifs associating Pro with either Ser or Thr. Thus, we have studied the conformation of Pro-kinked transmembrane helices containing Ser or Thr residues, in both g+ and g- rotamers, by molecular dynamics simulations in a hydrophobic environment. Analysis of the simulations shows that Ser or Thr can significantly modulate the deformation of the Pro. A series of motifs, such as (S/T)P and (S/T)AP in the g+ rotamer and the TAP and PAA(S/T) motifs in the g- rotamer, induce an increase in bending angle of the helix compared to a standard Pro-kink, apparently due to the additional hydrogen bond formed between the side chain of Ser/Thr and the backbone carbonyl oxygen. In contrast, (S/T)AAP and PA(S/T) motifs, in both g+ and g-, and PAA(S/T) in g+ rotamers decrease the bending angle of the helix by either reducing the steric clash between the pyrrolidine ring of Pro and the helical backbone, or by adding a constrain in the form of a hydrogen bond in the curved-in face of the helix. Together with a number of available experimental data, our results strongly suggest that association of Ser and Thr with Pro is commonly used in transmembrane helices to accommodate the structural needs of specific functions.

Influence of the Environment in the Conformation of α-Helices Studied by Protein Database Search and Molecular Dynamics Simulations

The influence of the solvent on the main-chain conformation (phi and Psi dihedral angles) of alpha-helices has been studied by complementary approaches. A first approach consisted in surveying crystal structures of both soluble and membrane proteins. The residues of analysis were further classified as exposed to either the water (polar solvent) or the lipid (apolar solvent) environment or buried to the core of the protein (intermediate polarity). The statistical results show that the more polar the environment, the lower the value of phi(i) and the higher the value of Psi(i) are. The intrahelical hydrogen bond distance increases in water-exposed residues due to the additional hydrogen bond between the peptide carbonyl oxygen and the aqueous environment. A second approach involved nanosecond molecular dynamics simulations of poly-Ala alpha-helices in environments of different polarity: water to mimic hydrophilic environments that can form hydrogen bonds with the peptide carbonyl oxygen and methane to mimic hydrophobic environments without this hydrogen bond capabilities. These simulations reproduce similar effects in phi and Psi angles and intrahelical hydrogen bond distance and angle as observed in the protein survey analysis. The magnitude of the intrahelical hydrogen bond in the methane environment is stronger than in the water environment, suggesting that alpha-helices in membrane-embedded proteins are less flexible than in soluble proteins. There is a remarkable coincidence between the phi and Psi angles obtained in the analysis of residues exposed to the lipid in membrane proteins and the results from computer simulations in methane, which suggests that this simulation protocol properly mimic the lipidic cell membrane and reproduce several structural characteristics of membrane-embedded proteins. Finally, we have compared the phi and Psi torsional angles of Pro kinks in membrane protein crystal structures and in computer simulations.

Design and synthesis of S-(−)-2-[[4-(napht-1-yl)piperazin-1-yl]methyl]-1,4-dioxoperhydropyrrolo[1,2-a]pyrazine (CSP-2503) using computational simulation. A 5-HT1A receptor agonist

Based on a computational model for 5-HT(1A)R-ligand interaction and QSAR studies, we have designed and synthesized a new series of arylpiperazines 2-8 which exhibit high 5-HT(1A)R affinity and selectivity over alpha(1)-adrenergic receptors. Among them, compound CSP-2503 (4) has been pharmacologically characterized as a 5-HT(1A)R agonist at somatodendritic and postsynaptic sites, endowed with anxiolytic properties.

Influence of the g− conformation of Ser and Thr on the structure of transmembrane helices

In order to study the influence of Ser and Thr on the structure of transmembrane helices we have analyzed a database of helix stretches extracted from crystal structures of membrane proteins and an ensemble of model helices generated by molecular dynamics simulations. Both complementary analyses show that Ser and Thr in the g- conformation induce and/or stabilize a structural distortion in the helix backbone. Using quantum mechanical calculations, we have attributed this effect to the electrostatic repulsion between the side chain Ogamma atom of Ser and Thr and the backbone carbonyl oxygen at position i-3. In order to minimize the repulsive force between these negatively charged oxygens, there is a modest increase of the helix bend angle as well as a local opening of the helix turn preceding Ser/Thr. This small distortion can be amplified through the helix, resulting in a significant displacement of the residues located at the other side of the helix. The crystal structures of aquaporin Z and the beta(2)-adrenergic receptor are used to illustrate these effects. Ser/Thr-induced structural distortions can be implicated in processes as diverse as ligand recognition, protein function and protein folding.

Computational model of the complex between GR113808 and the 5-HT4 receptor guided by site-directed mutagenesis and the crystal structure of rhodopsin

A computational model of the transmembrane domain of the human 5-HT4 receptor complexed with the GR113808 antagonist was constructed from the crystal structure of rhodopsin and the putative residues of the ligand-binding site, experimentally determined by site-directed mutagenesis. The recognition mode of GR113808 consist of: (i) the ionic interaction between the protonated amine and Asp3.32; (ii) the hydrogen bond between the carbonylic oxygen and Ser5.43; (iii) the hydrogen bond between the ether oxygen and Asn6.55; (iv) the hydrogen bond between the C-H groups adjacent to the protonated piperidine nitrogen and the π electrons of Phe6.51; and (v) the π-σ aromatic-aromatic interaction between the indole ring and Phe6.52.This computational model offers structural indications about the role of Asp3.32, Ser5.43, Phe6.51, Phe6.52, and Asn6.55 in the experimental binding affinities. Asp3.32Asn mutation does not affect the binding of GR113808 because the loss of binding affinity from an ion pair to a charged hydrogen bond is compensated by the larger energetical penalty of Asp to disrupt its side chain environment in the ligand-free form, and the larger interaction between Phe6.51 and the piperidine ring of the ligand in the mutant receptor. In the Phe6.52Val mutant the indole ring of the ligand replaces the interaction with Phe6.52 by a similarly intense interaction with Tyr5.38, with no significant effect in the binding of GR113808. The mutation of Asn6.55 to Leu replaces the hydrogen bond of the ether oxygen of the ligand from Asn6.55 to Cys5.42, with a decrease of binding affinity that approximately equals the free energy difference between the SH⋯O and NH⋯O hydrogen bonds.Because these residues are also present in the other members of the neurotransmitter family of G protein-coupled receptors, these findings will also serve for our understanding of the binding of related ligands to their cognate receptors.

Analysis of the interactions of sulfur-containing amino acids in membrane proteins.

The interactions of Met and Cys with other amino acid side chains have received little attention, in contrast to aromatic–aromatic, aromatic–aliphatic or/and aliphatic–aliphatic interactions. Precisely, these are the only amino acids that contain a sulfur atom, which is highly polarizable and, thus, likely to participate in strong Van der Waals interactions. Analysis of the interactions present in membrane protein crystal structures, together with the characterization of their strength in small‐molecule model systems at the ab‐initio level, predicts that Met–Met interactions are stronger than Met–Cys ≈ Met–Phe ≈ Cys–Phe interactions, stronger than Phe–Phe ≈ Phe–Leu interactions, stronger than the Met–Leu interaction, and stronger than Leu–Leu ≈ Cys–Leu interactions. These results show that sulfur‐containing amino acids form stronger interactions than aromatic or aliphatic amino acids. Thus, these amino acids may provide additional driving forces for maintaining the 3D structure of membrane proteins and may provide functional specificity.

TMalphaDB and TMbetaDB: web servers to study the structural role of sequence motifs in α-helix and β-barrel domains of membrane proteins

BackgroundMembrane proteins represent over 25 % of human protein genes and account for more than 60 % of drug targets due to their accessibility from the extracellular environment. The increasing number of available crystal structures of these proteins in the Protein Data Bank permits an initial estimation of their structural properties.DescriptionWe have developed two web servers—TMalphaDB for α-helix bundles and TMbetaDB for β-barrels—to analyse the growing repertoire of available crystal structures of membrane proteins. TMalphaDB and TMbetaDB permit to search for these specific sequence motifs in a non-redundant structure database of transmembrane segments and quantify structural parameters such as ϕ and ψ backbone dihedral angles, χ1 side chain torsion angle, unit bend and unit twist.ConclusionsThe structural information offered by TMalphaDB and TMbetaDB permits to quantify structural distortions induced by specific sequence motifs, and to elucidate their role in the 3D structure. This specific structural information has direct implications in homology modeling of the growing sequences of membrane proteins lacking experimental structure. TMalphaDB and TMbetaDB are freely available at http://lmc.uab.cat/TMalphaDB and http://lmc.uab.cat/TMbetaDB.

GPCR-SAS: A web application for statistical analyses on G protein-coupled receptors sequences

G protein-coupled receptors (GPCRs) are one of the largest protein families in mammals. They mediate signal transduction across cell membranes and are important targets for the pharmaceutical industry. The G Protein-Coupled Receptors—Sequence Analysis and Statistics (GPCR-SAS) web application provides a set of tools to perform comparative analysis of sequence positions between receptors, based on a curated structural-informed multiple sequence alignment. The analysis tools include: (i) percentage of occurrence of an amino acid or motif and entropy at a position or range of positions, (ii) covariance of two positions, (iii) correlation between two amino acids in two positions (or two sequence motifs in two ranges of positions), and (iv) snake-plot representation for a specific receptor or for the consensus sequence of a group of selected receptors. The analysis of conservation of residues and motifs across transmembrane (TM) segments may guide the design of more selective ligands or help to rationalize activation mechanisms, among others. As an example, here we analyze the amino acids of the “transmission switch”, that initiates receptor activation following ligand binding. The tool is freely accessible at http://lmc.uab.cat/gpcrsas/.