Mono and biallelic variants in HCN2 cause severe neurodevelopmental disorders.

Hyperpolarization activated Cyclic Nucleotide (HCN) gated channels are crucial for various neurophysiological functions, including learning and sensory functions, and their dysfunction are responsible for brain disorders, such as epilepsy. To date, HCN2 variants have only been associated with mild epilepsy and recently, one monoallelic missense variant has been linked to developmental and epileptic encephalopathy. Here, we expand the phenotypic spectrum of HCN2- related disorders by describing twenty-one additional individuals from fifteen unrelated families carrying HCN2 variants. Seventeen individuals had developmental delay/intellectual disability (DD/ID), two had borderline DD/ID, and one had borderline DD. Ten individuals had epilepsy with DD/ID, with median age of onset of 10 months, and one had epilepsy with normal development. Molecular diagnosis identified thirteen different pathogenic HCN2 variants, including eleven missense variants affecting highly conserved amino acids, one frameshift variant, and one in-frame deletion. Seven variants were monoallelic of which five occurred de novo, one was not maternally inherited, one was inherited from a father with mild learning disabilities, and one was of unknown inheritance. The remaining six variants were biallelic, with four homozygous and two compound heterozygous variants. Functional studies using two-electrode voltage-clamp recordings in Xenopus laevis oocytes were performed on three monoallelic variants, p.(Arg324His), p.(Ala363Val), and p.(Met374Leu), and three biallelic variants, p.(Leu377His), p.(Pro493Leu) and p.(Gly587Asp). The p.(Arg324His) variant induced a strong increase of HCN2 conductance, while p.(Ala363Val) and p.(Met374Leu) displayed dominant negative effects, leading to a partial loss of HCN2 channel function. By confocal imaging, we found that the p.(Leu377His), p.(Pro493Leu) and p.(Gly587Asp) pathogenic variants impaired membrane trafficking, resulting in a complete loss of HCN2 elicited currents in Xenopus oocytes. Structural 3D-analysis in depolarized and hyperpolarized states of HCN2 channels, revealed that the pathogenic variants p.(His205Gln), p.(Ser409Leu), p.(Arg324Cys), p.(Asn369Ser) and p.(Gly460Asp) modify molecular interactions altering HCN2 function. Taken together, our data broadens the clinical spectrum associated with HCN2 variants, and disclose that HCN2 is involved in developmental encephalopathy with or without epilepsy.

Homology Modeling in the Twilight Zone: Improved Accuracy by Sequence Space Analysis.

The analysis of the relationship between sequence and structure similarities during the evolution of a protein family has revealed a limit of sequence divergence for which structural conservation can be confidently assumed and homology modeling is reliable. Below this limit, the twilight zone corresponds to sequence divergence for which homology modeling becomes increasingly difficult and requires specific methods. Either with conventional threading methods or with recent deep learning methods, such as AlphaFold, the challenge relies on the identification of a template that shares not only a common ancestor (homology) but also a conserved structure with the query. As both homology and structural conservation are transitive properties, mining of sequence databases followed by multidimensional scaling (MDS) of the query sequence space can reveal intermediary sequences to infer homology and structural conservation between the query and the template. Here, as a case study, we studied the plethodontid receptivity factor isoform 1 (PRF1) from Plethodon jordani, a member of a pheromone protein family present only in lungless salamanders and weakly related to cytokines of the IL6 family. A variety of conventional threading methods led to the cytokine CNTF as a template. Sequence mining, followed by phylogenetic and MDS analysis, provided missing links between PRF1 and CNTF and allowed reliable homology modeling. In addition, we compared automated models obtained from web servers to a customized model to show how modeling can be improved by expert information.

Evolutionary information helps understand distinctive features of the angiotensin II receptors AT1 and AT2 in amniota.

In vertebrates, the octopeptide angiotensin II (AngII) is an important in vivo regulator of the cardiovascular system. It acts mainly through two G protein-coupled receptors, AT1 and AT2. To better understand distinctive features of these receptors, we carried out a phylogenetic analysis that revealed a mirror evolution of AT1 and AT2, each one split into two clades, separating fish from terrestrial receptors. It also revealed that hallmark mutations occurred at, or near, the sodium binding site in both AT1 and AT2. Electrostatics computations and molecular dynamics simulations support maintained sodium binding to human AT1 with slow ingress from the extracellular side and an electrostatic component of the binding free energy around -3kT, to be compared to around -2kT for human AT2 and the δ opioid receptor. Comparison of the sodium binding modes in wild type and mutated AT1 and AT2 from humans and eels indicates that the allosteric control by sodium in both AT1 and AT2 evolved during the transition from fish to amniota. The unusual S7.46N mutation in AT1 is mirrored by a L3.36M mutation in AT2. In the presence of sodium, the N7.46 pattern in amniota AT1 stabilizes the inward orientation of N3.35 in the apo receptor, which should contribute to efficient N3.35 driven biased signaling. The M3.36 pattern in amniota AT2 favours the outward orientation of N3.35 and the receptor promiscuity. Both mutations have physiological consequences for the regulation of the renin-angiotensin system.

Homology Modeling of Class A G-Protein-Coupled Receptors in the Age of the Structure Boom.

With 700 members, G protein-coupled receptors (GPCRs) of the rhodopsin family (class A) form the largest membrane receptor family in humans and are the target of about 30% of presently available pharmaceutical drugs. The recent boom in GPCR structures led to the structural resolution of 57 unique receptors in different states (39 receptors in inactive state only, 2 receptors in active state only and 16 receptors in different activation states). In spite of these tremendous advances, most computational studies on GPCRs, including molecular dynamics simulations, virtual screening and drug design, rely on GPCR models obtained by homology modeling. In this protocol, we detail the different steps of homology modeling with the MODELLER software, from template selection to model evaluation. The present structure boom provides closely related templates for most receptors. If, in these templates, some of the loops are not resolved, in most cases, the numerous available structures enable to find loop templates with similar length for equivalent loops. However, simultaneously, the large number of putative templates leads to model ambiguities that may require additional information based on multiple sequence alignments or molecular dynamics simulations to be resolved. Using the modeling of the human bradykinin receptor B1 as a case study, we show how several templates are managed by MODELLER, and how the choice of template(s) and of template fragments can improve the quality of the models. We also give examples of how additional information and tools help the user to resolve ambiguities in GPCR modeling.

Bios2cor: an R package integrating dynamic and evolutionary correlations to identify functionally important residues in proteins.

SUMMARY: Both dynamic correlations in protein sidechain motions during molecular dynamics (MD) simulations and evolutionary correlations in multiple sequence alignments (MSAs) of homologous proteins may reveal functionally important residues. We developed the R package Bios2cor that provides a unique framework to investigate and, possibly, integrate both analyses. Bios2cor starts with an MSA or an MD trajectory and computes correlation/covariation scores between positions in the MSA or between sidechain dihedral angles or rotamers in the MD trajectory. In addition, Bios2cor provides a variety of tools for the analysis, the visualization and the interpretation of the data. AVAILABILITY AND IMPLEMENTATION: The R package Bios2cor is available from the Comprehensive R Archive Network, at https://CRAN.R-project.org/package=Bios2cor.

Deciphering collaborative sidechain motions in proteins during molecular dynamics simulations.

The dynamic structure of proteins is essential for their functions and may include large conformational transitions which can be studied by molecular dynamics (MD) simulations. However, details of these transitions are difficult to automatically track. To facilitate their analysis, we developed two scores of correlation between sidechain dihedral angles. The CIRCULAR and OMES scores are computed from, respectively, dihedral angle values and rotamer distributions. As a case study, we applied our methods to an activation-like transition of the chemokine receptor CXCR4, observed during accelerated MD simulations. The principal component analysis of the correlation matrices was consistent with the networking structure of the top ranking pairs. Both scores identify a set of residues whose "collaborative" sidechain rotamerization immediately preceded or accompanied the conformational transition of CXCR4. Detailed analysis of the sequential order of these rotamerizations suggests that an allosteric mechanism, involving the outward motion of an asparagine residue in transmembrane helix 3, might be a prerequisite to the large scale conformational transition of CXCR4. This case study provides the proof-of-concept that the correlation methods developed here are valuable exploratory techniques to help decipher complex reactional pathways.

Evolution of chemokine receptors is driven by mutations in the sodium binding site.

Chemokines and their receptors (members of the GPCR super-family) are involved in a wide variety of physiological processes and diseases; thus, understanding the specificity of the chemokine receptor family could help develop new receptor specific drugs. Here, we explore the evolutionary mechanisms that led to the emergence of the chemokine receptors. Based on GPCR hierarchical classification, we analyzed nested GPCR sets with an eigen decomposition approach of the sequence covariation matrix and determined three key residues whose mutation was crucial for the emergence of the chemokine receptors and their subsequent divergence into homeostatic and inflammatory receptors. These residues are part of the allosteric sodium binding site. Their structural and functional roles were investigated by molecular dynamics simulations of CXCR4 and CCR5 as prototypes of homeostatic and inflammatory chemokine receptors, respectively. This study indicates that the three mutations crucial for the evolution of the chemokine receptors dramatically altered the sodium binding mode. In CXCR4, the sodium ion is tightly bound by four protein atoms and one water molecule. In CCR5, the sodium ion is mobile within the binding pocket and moves between different sites involving from one to three protein atoms and two to five water molecules. Analysis of chemokine receptor evolution reveals that a highly constrained sodium binding site characterized most ancient receptors, and that the constraints were subsequently loosened during the divergence of this receptor family. We discuss the implications of these findings for the evolution of the chemokine receptor functions and mechanisms of action.

Nationwide French Study of RET Variants Detected from 2003 to 2013 Suggests a Possible Influence of Polymorphisms as Modifiers.

BACKGROUND: The presence of single nucleotide polymorphisms (SNPs) in the REarranged during Transfection (RET) gene has been investigated with regard to their potential role in the development or progression of medullary thyroid cancer or pheochromocytomas (PHEO) in patients with the multiple endocrine neoplasia type 2 (MEN2) syndrome. The aim of this study was to evaluate the spectrum of RET variants in France between 2003 and 2013, and to evaluate the impact of SNPs on the MEN2 A phenotype. METHODS: In this retrospective cohort study, RET variants were screened in 5109 index cases, and RET pathogenic variants were screened in 2214 relatives. Exons 5, 8, 10, 11, 13, 14, 15, and 16 were characterized by Sanger sequencing. RET pathogenic variants, RET variants with unknown functional significance (VUS), and four RET SNP variants-G691S (rs1799939), L769L (rs1800861), S836S (rs1800862), and S904S (rs1800863)-were characterized and are reported in index cases. In silico analysis and classification following the recommendation of the American College of Medical Genetics and Genomics was performed for RET VUS. Each patient's age at the time of diagnosis, sex, and the endocrine neoplasias present at molecular diagnosis were recorded. RESULTS: Twenty-six single VUS in RET without any well-defined risk profiles were found in 33 patients. Nine of these were considered probably pathogenic, 11 of uncertain significance, and six as probably benign. Three double pathogenic variants found in three patients were classified as pathogenic. A study of the entire cohort showed that patients carrying pathogenic variants or VUS in RET together with PHEO were diagnosed earlier than the others. The presence of the G691S SNP, or a combination of SNPs, increased the risk of developing PHEO but did not modify the date of the diagnosis. No association was found between SNPs and medullary thyroid cancer or hyperparathyroidism. CONCLUSIONS: The findings propose a classification of 15 of the 26 VUS in RET without any well-defined risk profiles and suggest that the G691S SNP, or a combination of SNPs, may be associated with the development of PHEO.

The G Protein-Coupled Receptor UT of the Neuropeptide Urotensin II Displays Structural and Functional Chemokine Features.

The urotensinergic system was previously considered as being linked to numerous physiopathological states, including atherosclerosis, heart failure, hypertension, pre-eclampsia, diabetes, renal disease, as well as brain vascular lesions. Thus, it turns out that the actions of the urotensin II (UII)/G protein-coupled receptor UT system in animal models are currently not predictive enough in regard to their effects in human clinical trials and that UII analogs, established to target UT, were not as beneficial as expected in pathological situations. Thus, many questions remain regarding the overall signaling profiles of UT leading to complex involvement in cardiovascular and inflammatory responses as well as cancer. We address the potential UT chemotactic structural and functional definition under an evolutionary angle, by the existence of a common conserved structural feature among chemokine receptorsopioïdergic receptors and UT, i.e., a specific proline position in the transmembrane domain-2 TM2 (P2.58) likely responsible for a kink helical structure that would play a key role in chemokine functions. Even if the last decade was devoted to the elucidation of the cardiovascular control by the urotensinergic system, we also attempt here to discuss the role of UII on inflammation and migration, likely providing a peptide chemokine status for UII. Indeed, our recent work established that activation of UT by a gradient concentration of UII recruits Gαi/o and Gα13 couplings in a spatiotemporal way, controlling key signaling events leading to chemotaxis. We think that this new vision of the urotensinergic system should help considering UT as a chemotactic therapeutic target in pathological situations involving cell chemoattraction.

Microvascular vasodilator properties of the angiotensin II type 2 receptor in a mouse model of type 1 diabetes.

Diabetes Mellitus is associated with severe cardiovascular disorders involving the renin-angiotensin system, mainly through activation of the angiotensin II type 1 receptor (AT1R). Although the type 2 receptor (AT2R) opposes the effects of AT1R, with vasodilator and anti-trophic properties, its role in diabetes is debatable. Thus we investigated AT2R-mediated dilatation in a model of type 1 diabetes induced by streptozotocin in 5-month-old male mice lacking AT2R (AT2R-/y). Glucose tolerance was reduced and markers of inflammation and oxidative stress (cyclooxygenase-2, gp91phox p22phox and p67phox) were increased in AT2R-/y mice compared to wild-type (WT) animals. Streptozotocin-induced hyperglycaemia was higher in AT2R-/y than in WT mice. Arterial gp91phox and MnSOD expression levels in addition to blood 8-isoprostane and creatinine were further increased in diabetic AT2R-/y mice compared to diabetic WT mice. AT2R-dependent dilatation in both isolated mesenteric resistance arteries and perfused kidneys was greater in diabetic mice than in non-diabetic animals. Thus, in type 1 diabetes, AT2R may reduce glycaemia and display anti-oxidant and/or anti-inflammatory properties in association with greater vasodilatation in mesenteric arteries and in the renal vasculature, a major target of diabetes. Therefore AT2R might represent a new therapeutic target in diabetes.

The angiotensin II type 2 receptor activates flow-mediated outward remodelling through T cells-dependent interleukin-17 production.

AIMS: The angiotensin II type 1 receptor (AT1R) through the activation of immune cells plays a key role in arterial inward remodelling and reduced blood flow in cardiovascular disorders. On the other side, flow (shear stress)-mediated outward remodelling (FMR), involved in collateral arteries growth in ischaemic diseases, allows revascularization. We hypothesized that the type 2 receptor (AT2R), described as opposing the effects of AT1R, could be involved in FMR. METHODS AND RESULTS: We studied FMR using a model of ligation of feed arteries supplying collateral pathways in the mouse mesenteric arterial bed in vivo. Seven days after ligation, diameter increased by 30% in high flow (HF) arteries compared with normal flow vessels. FMR was absent in mice lacking AT2R. At Day 2, T lymphocytes expressing AT2R were present preferentially around HF arteries. FMR did not occur in athymic (nude) mice lacking T cells and in mice treated with anti-CD3ε antibodies. AT2R activation induced interleukin-17 production by memory T cells. Treatment of nude mice or AT2R-deficient mice with interleukin-17 restored diameter enlargement in HF arteries. Interleukin-17 increased NO-dependent relaxation and matrix metalloproteinases activity, both important in FMR. Remodelling of feeding arteries in the skin flap model of ischaemia was also absent in AT2R-deficient mice and in anti-interleukin-17-treated mice. Finally, remodelling, absent in 12-month-old mice, was restored by a treatment with the AT2R non-peptidic agonist C21. CONCLUSION: AT2R-dependent interleukin-17 production by T lymphocyte is necessary for collateral artery growth and could represent a new therapeutic target in ischaemic disorders.

In Vitro Effects of the Endocrine Disruptor p,p'-DDT on Human Follitropin Receptor.

BACKGROUND: 1-chloro-4-[2,2,2-trichloro-1-(4-chlorophenyl)ethyl]benzene (p,p'-DDT) is a persistent environmental endocrine disruptor (ED). Several studies have shown an association between p,p'-DDT exposure and reproductive abnormalities. OBJECTIVES: To investigate the putative effects of p,p'-DDT on the human follitropin receptor (FSHR) function. METHODS AND RESULTS: We used Chinese hamster ovary (CHO) cells stably expressing human FSHR to investigate the impact of p,p'-DDT on FSHR activity and its interaction with the receptor. At a concentration of 5 µM p,p'-DDT increased the maximum response of the FSHR to follitropin by 32 ± 7.45%. However, 5 µM p,p'-DDT decreased the basal activity and did not influence the maximal response of the closely related LH/hCG receptor to human chorionic gonadotropin (hCG). The potentiating effect of p,p'-DDT was specific for the FSHR. Moreover, in cells that did not express FSHR, p,p'-DDT had no effect on cAMP response. Thus, the potentiating effect of p,p'-DDT was dependent on the FSHR. In addition, p,p'-DDT increased the sensitivity of FSHR to hCG and to a low molecular weight agonist of the FSHR, 3-((5methyl)-2-(4-benzyloxy-phenyl)-5-{[2-[3-ethoxy-4-methoxy-phenyl)-ethylcarbamoyl]-methyl}-4-oxo-thiazolidin-3-yl)-benzamide (16a). Basal activity in response to p,p'-DDT and potentiation of the FSHR response to FSH by p,p'-DDT varied among FSHR mutants with altered transmembrane domains (TMDs), consistent with an effect of p,p'-DDT via TMD binding. This finding was corroborated by the results of simultaneously docking p,p'-DDT and 16a into the FSHR transmembrane bundle. CONCLUSION: p,p'-DDT acted as a positive allosteric modulator of the FSHR in our experimental model. These findings suggest that G protein-coupled receptors are additional targets of endocrine disruptors. CITATION: Munier M, Grouleff J, Gourdin L, Fauchard M, Chantreau V, Henrion D, Coutant R, Schiøtt B, Chabbert M, Rodien P. 2016. In vitro effects of the endocrine disruptor p,p'-DDT on human follitropin receptor. Environ Health Perspect 124:991-999; http://dx.doi.org/10.1289/ehp.1510006.

Molecular Insights into the Transmembrane Domain of the Thyrotropin Receptor.

The thyrotropin receptor (TSHR) is a G protein-coupled receptor (GPCR) that is member of the leucine-rich repeat subfamily (LGR). In the absence of crystal structure, the success of rational design of ligands targeting the receptor internal cavity depends on the quality of the TSHR models built. In this subfamily, transmembrane helices (TM) 2 and 5 are characterized by the absence of proline compared to most receptors, raising the question of the structural conformation of these helices. To gain insight into the structural properties of these helices, we carried out bioinformatics and experimental studies. Evolutionary analysis of the LGR family revealed a deletion in TM5 but provided no information on TM2. Wild type residues at positions 2.58, 2.59 or 2.60 in TM2 and/or at position 5.50 in TM5 were substituted to proline. Depending on the position of the proline substitution, different effects were observed on membrane expression, glycosylation, constitutive cAMP activity and responses to thyrotropin. Only proline substitution at position 2.59 maintained complex glycosylation and high membrane expression, supporting occurrence of a bulged TM2. The TSHR transmembrane domain was modeled by homology with the orexin 2 receptor, using a protocol that forced the deletion of one residue in the TM5 bulge of the template. The stability of the model was assessed by molecular dynamics simulations. TM5 straightened during the equilibration phase and was stable for the remainder of the simulations. Our data support a structural model of the TSHR transmembrane domain with a bulged TM2 and a straight TM5 that is specific of glycoprotein hormone receptors.

Comparative analysis of sequence covariation methods to mine evolutionary hubs: examples from selected GPCR families.

Covariation between positions in a multiple sequence alignment may reflect structural, functional, and/or phylogenetic constraints and can be analyzed by a wide variety of methods. We explored several of these methods for their ability to identify covarying positions related to the divergence of a protein family at different hierarchical levels. Specifically, we compared seven methods on a model system composed of three nested sets of G-protein-coupled receptors (GPCRs) in which a divergence event occurred. The covariation methods analyzed were based on: χ2 test, mutual information, substitution matrices, and perturbation methods. We first analyzed the dependence of the covariation scores on residue conservation (measured by sequence entropy), and then we analyzed the networking structure of the top pairs. Two methods out of seven--OMES (Observed minus Expected Squared) and ELSC (Explicit Likelihood of Subset Covariation)--favored pairs with intermediate entropy and a networking structure with a central residue involved in several high-scoring pairs. This networking structure was observed for the three sequence sets. In each case, the central residue corresponded to a residue known to be crucial for the evolution of the GPCR family and the subfamily specificity. These central residues can be viewed as evolutionary hubs, in relation with an epistasis-based mechanism of functional divergence within a protein family.

Structural evolution of G-protein-coupled receptors: a sequence space approach.

Class A G-protein-coupled receptors (GPCRs) provide a fascinating example of evolutionary success. In this review, we discuss how metric multidimensional scaling (MDS), a multivariate analysis method, complements traditional tree-based phylogenetic methods and helps decipher the mechanisms that drove the evolution of class A GPCRs. MDS provides low-dimensional representations of a distance matrix. Applied to a multiple sequence alignment, MDS represents the sequences in a Euclidean space as points whose interdistances are as close as possible to the distances in the alignment (the so-called sequence space). We detail how to perform the MDS analysis of a multiple sequence alignment and how to analyze and interpret the resulting sequence space. We also show that the projection of supplementary data (a property of the MDS method) can be used to straightforwardly monitor the evolutionary drift of specific subfamilies. The sequence space of class A GPCRs reveals the key role of mutations at the level of the TM2 and TM5 proline residues in the evolution of class A GPCRs.

Bios2mds: an R package for comparing orthologous protein families by metric multidimensional scaling.

BACKGROUND: The distance matrix computed from multiple alignments of homologous sequences is widely used by distance-based phylogenetic methods to provide information on the evolution of protein families. This matrix can also be visualized in a low dimensional space by metric multidimensional scaling (MDS). Applied to protein families, MDS provides information complementary to the information derived from tree-based methods. Moreover, MDS gives a unique opportunity to compare orthologous sequence sets because it can add supplementary elements to a reference space. RESULTS: The R package bios2mds (from BIOlogical Sequences to MultiDimensional Scaling) has been designed to analyze multiple sequence alignments by MDS. Bios2mds starts with a sequence alignment, builds a matrix of distances between the aligned sequences, and represents this matrix by MDS to visualize a sequence space. This package also offers the possibility of performing K-means clustering in the MDS derived sequence space. Most importantly, bios2mds includes a function that projects supplementary elements (a.k.a. "out of sample" elements) onto the space defined by reference or "active" elements. Orthologous sequence sets can thus be compared in a straightforward way. The data analysis and visualization tools have been specifically designed for an easy monitoring of the evolutionary drift of protein sub-families. CONCLUSIONS: The bios2mds package provides the tools for a complete integrated pipeline aimed at the MDS analysis of multiple sets of orthologous sequences in the R statistical environment. In addition, as the analysis can be carried out from user provided matrices, the projection function can be widely used on any kind of data.

Multidimensional scaling reveals the main evolutionary pathways of class A G-protein-coupled receptors.

Class A G-protein-coupled receptors (GPCRs) constitute the largest family of transmembrane receptors in the human genome. Understanding the mechanisms which drove the evolution of such a large family would help understand the specificity of each GPCR sub-family with applications to drug design. To gain evolutionary information on class A GPCRs, we explored their sequence space by metric multidimensional scaling analysis (MDS). Three-dimensional mapping of human sequences shows a non-uniform distribution of GPCRs, organized in clusters that lay along four privileged directions. To interpret these directions, we projected supplementary sequences from different species onto the human space used as a reference. With this technique, we can easily monitor the evolutionary drift of several GPCR sub-families from cnidarians to humans. Results support a model of radiative evolution of class A GPCRs from a central node formed by peptide receptors. The privileged directions obtained from the MDS analysis are interpretable in terms of three main evolutionary pathways related to specific sequence determinants. The first pathway was initiated by a deletion in transmembrane helix 2 (TM2) and led to three sub-families by divergent evolution. The second pathway corresponds to the differentiation of the amine receptors. The third pathway corresponds to parallel evolution of several sub-families in relation with a covarion process involving proline residues in TM2 and TM5. As exemplified with GPCRs, the MDS projection technique is an important tool to compare orthologous sequence sets and to help decipher the mutational events that drove the evolution of protein families.

Structural determinants stabilizing helical distortions related to proline.

The distortions related to proline in a alpha-helix can be accommodated by different structural elements. To obtain an exhaustive view of these distortions, we data mined a non-redundant subset of the Protein Data Bank in search of proline residues included either within contiguous helices or within structural motifs in which two helices are joined by a few linker residues with backbone dihedral angles in the generous alpha region. The distortions correspond to "typical" and "non-typical" proline distortions, with relative ratio of 0.65 and 0.35, respectively. Analysis of "non-typical" proline distortions indicates that most linkers have one (75%) or two residues (20%) and that proline is preferentially located at the second or third position of the second helix (95%). The dihedral angles of the linker residues are located in two areas of the generous alpha region. Structures with linker(s) in the alpha1 area, which is characterized by very negative phi values, possess i to i-5 H-bonds and correspond to pi bulges. Structures with linker(s) in the alpha2 area, which links the alpha and beta regions, possess i to i-3 H-bonds and correspond to tight turns. Further classification of bulges and turns as a function of the linker length and proline position yields five canonical structures, representing 85% of "non-typical" proline distortions. These structures are characterized by distinct H-bonding patterns and structural determinants and correspond to different classes of pi bulges and tight turns. This hierarchical approach provides a straightforward and robust classification of proline-related helical distortions.

An indel in transmembrane helix 2 helps to trace the molecular evolution of class A G-protein-coupled receptors.

Class A G-protein-coupled receptors (GPCRs) constitute a large family of transmembrane receptors. Helical distortions play a major role in the overall fold of these receptors. Most are related to conserved proline residues. However, in transmembrane helix 2, the proline pattern is not conserved, and when present, proline may be located at position 2.58, 2.59, or 2.60. Sequence analysis, three-dimensional data mining, and molecular modeling were undertaken to investigate the origin of this unusual pattern. Taken together, the data strongly support the assumption that an indel led to two structural motifs for helix 2: a bulged structure in P2.59 and P2.60 receptors and a "typical" proline kink in P2.58 receptors. The proline pattern of helix 2 can be used as an evolutionary marker and helps to trace the molecular evolution of class A GPCRs. Two indel events yielding functional receptors occurred independently. One indel arose very early in GPCR evolution, in a bilaterian ancestor, before the protostome-deuterostome divergence. This indel led to the split between the P2.58 somatostatin/opioid receptors and other peptide receptors with the P2.59 pattern. A second indel also occurred in insect opsins and corresponds to a deletion. Subfamilies with proline at position 2.59 or no proline expanded earlier, whereas P2.60 receptors remained marginal throughout evolution. P2.58 receptors underwent rapid expansion in vertebrates with the development of the chemokine and purinergic receptor subfamilies from somatostatin/opioid-related ancestors.

Comprehensive analysis of the helix-X-helix motif in soluble proteins.

Alpha-helices are the most common secondary structures found in globular proteins. In this report, we analyze the stereochemical and sequence properties of helix-X-helix (HXH) motifs in which two alpha-helices are linked by a single residue, in search of characteristic structures and sequence signals. The analysis is carried out on a database of 837 nonredundant HXH motifs. The kinks are characterized by the bend angle between the axes of the N-terminal and C-terminal helices and the wobble angle corresponding to the rotation of C-terminal helix axis on the plane perpendicular to the N-terminal one. The phi-psi dihedral angles of the linker residue are clustered in six distinct areas of the Ramachandran plot: two areas are located in the additional allowed alpha region (alpha(1) and alpha(2)), two areas are in the additional allowed beta region (beta(1) and beta(2)) and two areas have positive phi values (alpha(L) and beta(M)). Each phi/psi region corresponds to characteristic bend and wobble angles and amino acid distributions. Bend angles can vary from 0 degrees to 160 degrees. Most wobble angles correspond to a counter-clockwise rotation of the C-terminal helix. Proline residues are rigorously excluded from the linker position X but have a high propensity at position X+1 of the beta(1) and beta(2) motifs (12 and 7, respectively) and at position X+3 of the alpha(1) motifs (9). Glycine linkers are located either in the alpha(L) region (20%) or in the beta(M) region (80%). This latter conformation is characterized by a marked bend angle (124 degrees +/- 18 degrees) and a clockwise wobble. Among other amino acids, Asn is remarkable for its high propensity (>3) at the linker position of the alpha(2), beta(1), and beta(2) motifs. Stabilization of HXH motifs by H-bonds between polar side chains of the linker and polar groups of the backbone is determined. A method based on position-specific scoring matrices is developed for conformational prediction. The accuracy of the predictions reaches 80% when the method is applied to proline-induced kinks or to kinks with bend angles in the 50 degrees-100 degrees range.