Statistical analyses and computational prediction of helical kinks in membrane proteins

Photoreception and signaling in bacterial phytochrome revealed by single-particle cryo-EM

Phytochromes are red-light photoreceptors discovered in plants with homologs in bacteria and fungi that regulate a variety of physiological responses. They display a reversible photocycle between two distinct states: a red-light-absorbing Pr state and a far-red light-absorbing Pfr state. The photoconversion regulates the activity of an enzymatic domain, usually a histidine kinase (HK). The molecular mechanism that explains how light controls the HK activity is not understood because structures of unmodified bacterial phytochromes with HK activity are missing. Here, we report three cryo-electron microscopy structures of a wild-type bacterial phytochrome with HK activity determined as Pr and Pfr homodimers and as a Pr/Pfr heterodimer with individual subunits in distinct states. We propose that the Pr/Pfr heterodimer is a physiologically relevant signal transduction intermediate. Our results offer insight into the molecular mechanism that controls the enzymatic activity of the HK as part of a bacterial two-component system that perceives and transduces light signals.

Signal Transduction in an Enzymatic Photoreceptor Revealed by Cryo-Electron Microscopy

Phytochromes are essential photoreceptor proteins in plants with homologs in bacteria and fungi that regulate a variety of important environmental responses. They display a reversible photocycle between two distinct states, the red-light absorbing Pr and the far-red light absorbing Pfr, each with its own structure. The reversible Pr to Pfr photoconversion requires covalently bound bilin chromophore and regulates the activity of a C-terminal enzymatic domain, which is usually a histidine kinase (HK). In plants, phytochromes translocate to nucleus where the C-terminal effector domain interacts with protein interaction factors (PIFs) to induce gene expression. In bacteria, the HK phosphorylates a response-regulator (RR) protein triggering downstream gene expression through a two-component signaling pathway. Although plant and bacterial phytochromes share similar structural composition, they have contrasting activity in the presence of light with most BphPs being active in the dark. The molecular mechanism that explains bacterial and plant phytochrome signaling has not been well understood due to limited structures of full-length phytochromes with enzymatic domain resolved at or near atomic resolution in both Pr and Pfr states. Here, we report the first Cryo-EM structures of a wild-type bacterial phytochrome with a HK enzymatic domain, determined in both Pr and Pfr states, between 3.75 and 4.13 Å resolution, respectively. Furthermore, we capture a distinct Pr/Pfr heterodimer of the same protein as potential signal transduction intermediate at 3.75 Å resolution. Our three Cryo-EM structures of the distinct signaling states of BphPs are further reinforced by Cryo-EM structures of the truncated PCM of the same protein determined for the Pr/Pfr heterodimer as well as Pfr state. These structures provide insight into the different light-signaling mechanisms that could explain how bacteria and plants see the light.

Structural basis for the Pr-Pfr long-range signaling mechanism of a full-length bacterial phytochrome at the atomic level

Phytochromes constitute a widespread photoreceptor family that typically interconverts between two photostates called Pr (red light–absorbing) and Pfr (far-red light–absorbing). The lack of full-length structures solved at the (near-)atomic level in both pure Pr and Pfr states leaves gaps in the structural mechanisms involved in the signal transmission pathways during the photoconversion. Here, we present the crystallographic structures of three versions from the plant pathogen Xanthomonas campestris virulence regulator XccBphP bacteriophytochrome, including two full-length proteins, in the Pr and Pfr states. The structures show a reorganization of the interaction networks within and around the chromophore-binding pocket, an α-helix/β-sheet tongue transition, and specific domain reorientations, along with interchanging kinks and breaks at the helical spine as a result of the photoswitching, which subsequently affect the quaternary assembly. These structural findings, combined with multidisciplinary studies, allow us to describe the signaling mechanism of a full-length bacterial phytochrome at the atomic level.

Identification of residues critical for topology inversion of the transmembrane protein TM4SF20 through regulated alternative translocation

Adopting a proper topology is crucial for transmembrane proteins to perform their functions. We previously reported that ceramide regulates a transmembrane protein called TM4SF20 (transmembrane 4 L six family member 20) through topological inversion by altering the direction through which the protein is translocated across membranes during translation. This regulatory mechanism, denoted regulated alternative translocation (RAT), depends on a GXXXN motif present in the first transmembrane helix of TM4SF20. Here, using site-directed mutagenesis, we show that Asn-26 in the motif is crucial for RAT of TM4SF20, as it cannot be replaced even by Gln. In contrast, Gly-22 in the motif could be substituted by other small residues such as Ala and Ser without affecting RAT of TM4SF20. We further demonstrate that the GXXXN motif alone is insufficient to induce RAT of a transmembrane protein because TM4SF4, a relative of TM4SF20 that also contains the motif in the first transmembrane helix, did not undergo RAT. Using TM4SF40-TM4SF20 chimeras, we identified Pro-29 of TM4SF20 as another important element required for RAT of the protein. Substituting Pro-29 alone did not affect RAT of TM4SF20, whereas replacing Pro-29 together with either Leu-25 or Val-17 of TM4SF20 with the corresponding residues of TM4SF4 abolished RAT of TM4SF20. Because Val-17, Gly-22, Leu-25, Asn-26, and Pro-29 are predicted to reside along the same surface of the transmembrane helix, our results suggest that interactions with other proteins mediated by this surface during translocation may be critical for RAT of TM4SF20.

Exploring conformational states and helical packings in the P2X receptor transmembrane domain by molecular dynamics simulation

The P2X receptor is a trimeric transmembrane protein that acts as an ATP-gated ion channel. Its transmembrane domain (TMD) contains only six helices and three of them, the M2 helices, line the ion conduction pathway. Here, using molecular dynamics simulation, I identify four conformational states of the TMD that are associated with four types of packing between M2 helices. Packing in the extracellular half of the M2 helix produces closed conformations, while packing in the intracellular half produces both open and closed conformations. State transition is observed and supports a mechanism where iris-like twisting of the M2 helices switches the location of helical packing between the extracellular and the intracellular halves of the helices. In addition, this twisting motion alters the position and orientation of residue side-chains relative to the pore and therefore influences the pore geometry and possibly ion permeation. Helical packing, on the other hand, may restrict the twisting motion and generate discrete conformational states.

Clustering and visualizing similarity networks of membrane proteins

We proposed a fast and unsupervised clustering method, minimum span clustering (MSC), for analyzing the sequence-structure-function relationship of biological networks, and demonstrated its validity in clustering the sequence/structure similarity networks (SSN) of 682 membrane protein (MP) chains. The MSC clustering of MPs based on their sequence information was found to be consistent with their tertiary structures and functions. For the largest seven clusters predicted by MSC, the consistency in chain function within the same cluster is found to be 100%. From analyzing the edge distribution of SSN for MPs, we found a characteristic threshold distance for the boundary between clusters, over which SSN of MPs could be properly clustered by an unsupervised sparsification of the network distance matrix. The clustering results of MPs from both MSC and the unsupervised sparsification methods are consistent with each other, and have high intracluster similarity and low intercluster similarity in sequence, structure, and function. Our study showed a strong sequence-structure-function relationship of MPs. We discussed evidence of convergent evolution of MPs and suggested applications in finding structural similarities and predicting biological functions of MP chains based on their sequence information.

Understanding intramembrane proteolysis: from protein dynamics to reaction kinetics

Intramembrane proteolysis - cleavage of proteins within the plane of a membrane - is a widespread phenomenon that can contribute to the functional activation of substrates and is involved in several diseases. Although different families of intramembrane proteases have been discovered and characterized, we currently do not know how these enzymes discriminate between substrates and non-substrates, how site-specific cleavage is achieved, or which factors determine the rate of proteolysis. Focusing on γ-secretase and rhomboid proteases, we argue that answers to these questions may emerge from connecting experimental readouts, such as reaction kinetics and the determination of cleavage sites, to the structures and the conformational dynamics of substrates and enzymes.

Computational prediction of kink properties of helices in membrane proteins

We have combined molecular dynamics simulations and fold identification procedures to investigate the structure of 696 kinked and 120 unkinked transmembrane (TM) helices in the PDBTM database. Our main aim of this study is to understand the formation of helical kinks by simulating their quasi-equilibrium heating processes, which might be relevant to the prediction of their structural features. The simulated structural features of these TM helices, including the position and the angle of helical kinks, were analyzed and compared with statistical data from PDBTM. From quasi-equilibrium heating processes of TM helices with four very different relaxation time constants, we found that these processes gave comparable predictions of the structural features of TM helices. Overall, 95 % of our best kink position predictions have an error of no more than two residues and 75 % of our best angle predictions have an error of less than 15°. Various structure assessments have been carried out to assess our predicted models of TM helices in PDBTM. Our results show that, in 696 predicted kinked helices, 70 % have a RMSD less than 2 Å, 71 % have a TM-score greater than 0.5, 69 % have a MaxSub score greater than 0.8, 60 % have a GDT-TS score greater than 85, and 58 % have a GDT-HA score greater than 70. For unkinked helices, our predicted models are also highly consistent with their crystal structure. These results provide strong supports for our assumption that kink formation of TM helices in quasi-equilibrium heating processes is relevant to predicting the structure of TM helices.

Minireview: Toward the establishment of a link between melatonin and glucose homeostasis: association of melatonin MT2 receptor variants with type 2 diabetes

The existence of interindividual variations in G protein-coupled receptor sequences has been recognized early on. Recent advances in large-scale exon sequencing techniques are expected to dramatically increase the number of variants identified in G protein-coupled receptors, giving rise to new challenges regarding their functional characterization. The current minireview will illustrate these challenges based on the MTNR1B gene, which encodes the melatonin MT2 receptor, for which exon sequencing revealed 40 rare nonsynonymous variants in the general population and in type 2 diabetes (T2D) cohorts. Functional characterization of these MT2 mutants revealed 14 mutants with loss of Gi protein activation that associate with increased risk of T2D development. This repertoire of disease-associated mutants is a rich source for structure-activity studies and will help to define the still poorly understood role of melatonin in glucose homeostasis and T2D development in humans. Defining the functional defects in carriers of rare MT2 mutations will help to provide personalized therapies to these patients in the future.