Helix perturbations in membrane proteins assist in inter-helical interactions and optimal helix positioning in the bilayer

Transmembrane (TM) helices in integral membrane proteins are primarily α-helical in structure. Here we analyze 1134 TM helices in 90 high resolution membrane proteins and find that apart from the widely prevalent α-helices, TM regions also contain stretches of 3₁₀ (3 to 8 residues) and π-helices (5 to 19 residues) with distinct sequence signatures. The various helix perturbations in TM regions comprise of helices with kinked geometry, as well as those with an interspersed 3₁₀/π-helical fragment and show high occurrence in a few membrane proteins. Proline is frequently present at sites of these perturbations, but it is neither a necessary nor a sufficient requirement. Helix perturbations are also conserved within a family of membrane proteins despite low sequence identity in the perturbed region. Furthermore, a perturbation influences the geometry of the TM helix, mediates inter-helical interactions within and across protein chains and avoids hydrophobic mismatch of the helix termini with the bilayer. An analysis of π-helices in the TM regions of the heme copper oxidase superfamily shows that interspersed π-helices can vary in length from 6 to 19 amino acids or be entirely absent, depending upon the protein function. The results presented here would be helpful for prediction of 3₁₀ and π-helices in TM regions and can assist the computational design of membrane proteins.

Examining the Conservation of Kinks in Alpha Helices

Kinks are a structural feature of alpha-helices and many are known to have functional roles. Kinks have previously tended to be defined in a binary fashion. In this paper we have deliberately moved towards defining them on a continuum, which given the unimodal distribution of kink angles is a better description. From this perspective, we examine the conservation of kinks in proteins. We find that kink angles are not generally a conserved property of homologs, pointing either to their not being functionally critical or to their function being related to conformational flexibility. In the latter case, the different structures of homologs are providing snapshots of different conformations. Sequence identity between homologous helices is informative in terms of kink conservation, but almost equally so is the sequence identity of residues in spatial proximity to the kink. In the specific case of proline, which is known to be prevalent in kinked helices, loss of a proline from a kinked helix often also results in the loss of a kink or reduction in its kink angle. We carried out a study of the seven transmembrane helices in the GPCR family and found that changes in kinks could be related both to subfamilies of GPCRs and also, in a particular subfamily, to the binding of agonists or antagonists. These results suggest conformational change upon receptor activation within the GPCR family. We also found correlation between kink angles in different helices, and the possibility of concerted motion could be investigated further by applying our method to molecular dynamics simulations. These observations reinforce the belief that helix kinks are key, functional, flexible points in structures.

Enhancing the Reliability of GPCR Models by Accounting for Flexibility of Their Pro-Containing Helices: the Case of the Human mAChR1 Receptor

To better investigate the GPCR structures, we have recently proposed to explore their flexibility by simulating the bending of their Pro-containing TM helices so generating a set of models (the so-called chimeras) which exhaustively combine the two conformations (bent and straight) of these helices. The primary objective of the study is to investigate whether such an approach can be exploited to enhance the reliability of the GPCR models generated by distant templates. The study was focused on the human mAChR1 receptor for which a presumably reliable model was generated using the congener mAChR3 as the template along with a second less reliable model based on the distant β2-AR template. The second model was then utilized to produce the chimeras by combining the conformations of its Pro-containing helices (i.e., TM4, TM5, TM6 and TM7 with 16 modeled chimeras). The reliability of such chimeras was assessed by virtual screening campaigns as evaluated using a novel skewness metric where they surpassed the predictive power of the more reliable mAChR1 model. Finally, the virtual screening campaigns emphasize the opportunity of synergistically combining the scores of more chimeras using a specially developed tool which generates highly predictive consensus functions by maximizing the corresponding enrichment factors.

Helix kinks are equally prevalent in soluble and membrane proteins

Helix kinks are a common feature of α-helical membrane proteins, but are thought to be rare in soluble proteins. In this study we find that kinks are a feature of long α-helices in both soluble and membrane proteins, rather than just transmembrane α-helices. The apparent rarity of kinks in soluble proteins is due to the relative infrequency of long helices (≥20 residues) in these proteins. We compare length-matched sets of soluble and membrane helices, and find that the frequency of kinks, the role of Proline, the patterns of other amino acid around kinks (allowing for the expected differences in amino acid distributions between the two types of protein), and the effects of hydrogen bonds are the same for the two types of helices. In both types of protein, helices that contain Proline in the second and subsequent turns are very frequently kinked. However, there are a sizeable proportion of kinked helices that do not contain a Proline in either their sequence or sequence homolog. Moreover, we observe that in soluble proteins, kinked helices have a structural preference in that they typically point into the solvent.

Toward an understanding of agonist binding to human Orexin-1 and Orexin-2 receptors with G-protein-coupled receptor modeling and site-directed mutagenesis

The class A G-protein-coupled receptors (GPCRs) Orexin-1 (OX1) and Orexin-2 (OX2) are located predominantly in the brain and are linked to a range of different physiological functions, including the control of feeding, energy metabolism, modulation of neuro-endocrine function, and regulation of the sleep-wake cycle. The natural agonists for OX1 and OX2 are two neuropeptides, Orexin-A and Orexin-B, which have activity at both receptors. Site-directed mutagenesis (SDM) has been reported on both the receptors and the peptides and has provided important insight into key features responsible for agonist activity. However, the structural interpretation of how these data are linked together is still lacking. In this work, we produced and used SDM data, homology modeling followed by MD simulation, and ensemble-flexible docking to generate binding poses of the Orexin peptides in the OX receptors to rationalize the SDM data. We also developed a protein pairwise similarity comparing method (ProS) and a GPCR-likeness assessment score (GLAS) to explore the structural data generated within a molecular dynamics simulation and to help distinguish between different GPCR substates. The results demonstrate how these newly developed methods of structural assessment for GPCRs can be used to provide a working model of neuropeptide-Orexin receptor interaction.

Fighting obesity with a sugar-based library: discovery of novel MCH-1R antagonists by a new computational-VAST approach for exploration of GPCR binding sites

Obesity is an increasingly common disease. While antagonism of the melanin-concentrating hormone-1 receptor (MCH-1R) has been widely reported as a promising therapeutic avenue for obesity treatment, no MCH-1R antagonists have reached the market. Discovery and optimization of new chemical matter targeting MCH-1R is hindered by reduced HTS success rates and a lack of structural information about the MCH-1R binding site. X-ray crystallography and NMR, the major experimental sources of structural information, are very slow processes for membrane proteins and are not currently feasible for every GPCR or GPCR-ligand complex. This situation significantly limits the ability of these methods to impact the drug discovery process for GPCR targets in "real-time", and hence, there is an urgent need for other practical and cost-efficient alternatives. We present here a conceptually pioneering approach that integrates GPCR modeling with design, synthesis, and screening of a diverse library of sugar-based compounds from the VAST technology (versatile assembly on stable templates) to provide structural insights on the MCH-1R binding site. This approach creates a cost-efficient new avenue for structure-based drug discovery (SBDD) against GPCR targets. In our work, a primary VAST hit was used to construct a high-quality MCH-1R model. Following model validation, a structure-based virtual screen yielded a 14% hit rate and 10 novel chemotypes of potent MCH-1R antagonists, including EOAI3367472 (IC50 = 131 nM) and EOAI3367474 (IC50 = 213 nM).