Reengineering of a flavin-binding fluorescent protein using ProteinMPNN

Recent advances in machine learning techniques have led to development of a number of protein design and engineering approaches. One of them, ProteinMPNN, predicts an amino acid sequence that would fold and match user-defined backbone structure. Its performance was previously tested for proteins composed of standard amino acids, as well as for peptide- and protein-binding proteins. In this short report, we test whether ProteinMPNN can be used to reengineer a non-proteinaceous ligand-binding protein, flavin-based fluorescent protein CagFbFP. We fixed the native backbone conformation and the identity of 20 amino acids interacting with the chromophore (flavin mononucleotide, FMN) while letting ProteinMPNN predict the rest of the sequence. The software package suggested replacing 36-48 out of the remaining 86 amino acids so that the resulting sequences are 55%-66% identical to the original one. The three designs that we tested experimentally displayed different expression levels, yet all were able to bind FMN and displayed fluorescence, thermal stability, and other properties similar to those of CagFbFP. Our results demonstrate that ProteinMPNN can be used to generate diverging unnatural variants of fluorescent proteins, and, more generally, to reengineer proteins without losing their ligand-binding capabilities.

Overlay databank unlocks data-driven analyses of biomolecules for all

Tools based on artificial intelligence (AI) are currently revolutionising many fields, yet their applications are often limited by the lack of suitable training data in programmatically accessible format. Here we propose an effective solution to make data scattered in various locations and formats accessible for data-driven and machine learning applications using the overlay databank format. To demonstrate the practical relevance of such approach, we present the NMRlipids Databank-a community-driven, open-for-all database featuring programmatic access to quality-evaluated atom-resolution molecular dynamics simulations of cellular membranes. Cellular membrane lipid composition is implicated in diseases and controls major biological functions, but membranes are difficult to study experimentally due to their intrinsic disorder and complex phase behaviour. While MD simulations have been useful in understanding membrane systems, they require significant computational resources and often suffer from inaccuracies in model parameters. Here, we demonstrate how programmable interface for flexible implementation of data-driven and machine learning applications, and rapid access to simulation data through a graphical user interface, unlock possibilities beyond current MD simulation and experimental studies to understand cellular membranes. The proposed overlay databank concept can be further applied to other biomolecules, as well as in other fields where similar barriers hinder the AI revolution.

Two distinct mechanisms of flavoprotein spectral tuning revealed by low-temperature and time-dependent spectroscopy

Flavins such as flavin mononucleotide or flavin adenine dinucleotide are bound by diverse proteins, yet have very similar spectra when in the oxidized state. Recently, we developed new variants of flavin-binding protein CagFbFP exhibiting notable blue (Q148V) or red (I52V A85Q) shifts of fluorescence emission maxima. Here, we use time-resolved and low-temperature spectroscopy to show that whereas the chromophore environment is static in Q148V, an additional protein-flavin hydrogen bond is formed upon photoexcitation in the I52V A85Q variant. Consequently, in Q148V, excitation, emission, and phosphorescence spectra are shifted, whereas in I52V A85Q, excitation and low-temperature phosphorescence spectra are relatively unchanged, while emission spectrum is altered. We also determine the x-ray structures of the two variants to reveal the flavin environment and complement the spectroscopy data. Our findings illustrate two distinct color-tuning mechanisms of flavin-binding proteins and could be helpful for the engineering of new variants with improved optical properties.

Mechanisms of inward transmembrane proton translocation

Proton transport is indispensable for cell life. It is believed that molecular mechanisms of proton movement through different types of proton-conducting molecules have general universal features. However, elucidation of such mechanisms is a challenge. It requires true-atomic-resolution structures of all key proton-conducting states. Here we present a comprehensive function-structure study of a light-driven bacterial inward proton pump, xenorhodopsin, from Bacillus coahuilensis in all major proton-conducting states. The structures reveal that proton translocation is based on proton wires regulated by internal gates. The wires serve as both selectivity filters and translocation pathways for protons. The cumulative results suggest a general concept of proton translocation. We demonstrate the use of serial time-resolved crystallography at a synchrotron source with sub-millisecond resolution for rhodopsin studies, opening the door for principally new applications. The results might also be of interest for optogenetics since xenorhodopsins are the only alternative tools to fire neurons.

Development and Characterization of Flavin-Binding Fluorescent Proteins, Part I: Basic Characterization

Flavin-based fluorescent proteins (FbFPs) are small fluorescent proteins derived from light-oxygen-voltage (LOV) domains. The proteins bind ubiquitous endogenous flavins as chromophores and can be used as versatile in vivo reporter proteins under aerobic and anaerobic conditions. This chapter presents the methodology to identify LOV domain sequences in genomic databases; design new FbFPs; characterize their biochemical, spectroscopic, photophysical, and photochemical properties; and conduct basic fluorescence microscopy experiments.

Structural insights into light-driven anion pumping in cyanobacteria

Transmembrane ion transport is a key process in living cells. Active transport of ions is carried out by various ion transporters including microbial rhodopsins (MRs). MRs perform diverse functions such as active and passive ion transport, photo-sensing, and others. In particular, MRs can pump various monovalent ions like Na⁺, K⁺, Cl^-, I^-, NO₃^-. The only characterized MR proposed to pump sulfate in addition to halides belongs to the cyanobacterium Synechocystis sp. PCC 7509 and is named Synechocystis halorhodopsin (SyHR). The structural study of SyHR may help to understand what makes an MR pump divalent ions. Here we present the crystal structure of SyHR in the ground state, the structure of its sulfate-bound form as well as two photoreaction intermediates, the K and O states. These data reveal the molecular origin of the unique properties of the protein (exceptionally strong chloride binding and proposed pumping of divalent anions) and sheds light on the mechanism of anion release and uptake in cyanobacterial halorhodopsins. The unique properties of SyHR highlight its potential as an optogenetics tool and may help engineer different types of anion pumps with applications in optogenetics.

Structural basis for receptor selectivity and inverse agonism in S1P5 receptors

The bioactive lysophospholipid sphingosine-1-phosphate (S1P) acts via five different subtypes of S1P receptors (S1PRs) - S1P_1-5. S1P₅ is predominantly expressed in nervous and immune systems, regulating the egress of natural killer cells from lymph nodes and playing a role in immune and neurodegenerative disorders, as well as carcinogenesis. Several S1PR therapeutic drugs have been developed to treat these diseases; however, they lack receptor subtype selectivity, which leads to side effects. In this article, we describe a 2.2 Å resolution room temperature crystal structure of the human S1P₅ receptor in complex with a selective inverse agonist determined by serial femtosecond crystallography (SFX) at the Pohang Accelerator Laboratory X-Ray Free Electron Laser (PAL-XFEL) and analyze its structure-activity relationship data. The structure demonstrates a unique ligand-binding mode, involving an allosteric sub-pocket, which clarifies the receptor subtype selectivity and provides a template for structure-based drug design. Together with previously published S1PR structures in complex with antagonists and agonists, our structure with S1P₅-inverse agonist sheds light on the activation mechanism and reveals structural determinants of the inverse agonism in the S1PR family.

S1P₅ is a sphingosine-1-phosphate (S1P) receptor implicated in immune and neurodegenerative disorders. Here, authors report a crystal structure of the S1P₅ receptor in complex with a selective inverse agonist, revealing an allosteric subpocket and shedding light on inverse agonism in S1P receptors.

Structure and dynamics of the SARS-CoV-2 envelope protein monomer

Coronaviruses, especially severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), present an ongoing threat to human wellbeing. Consequently, elucidation of molecular determinants of their function and interaction with the host is an important task. Whereas some of the coronaviral proteins are extensively characterized, others remain understudied. Here, we use molecular dynamics simulations to analyze the structure and dynamics of the SARS-CoV-2 envelope (E) protein (a viroporin) in the monomeric form. The protein consists of the hydrophobic α-helical transmembrane domain (TMD) and amphiphilic α-helices H2 and H3, connected by flexible linkers. We show that TMD has a preferable orientation in the membrane, while H2 and H3 reside at the membrane surface. Orientation of H2 is strongly influenced by palmitoylation of cysteines Cys40, Cys43, and Cys44. Glycosylation of Asn66 affects the orientation of H3. We also observe that the monomeric E protein both generates and senses the membrane curvature, preferably localizing with the C-terminus at the convex regions of the membrane; the protein in the pentameric form displays these properties as well. Localization to curved regions may be favorable for assembly of the E protein oligomers, whereas induction of curvature may facilitate the budding of the viral particles. The presented results may be helpful for a better understanding of the function of the coronaviral E protein and viroporins in general, and for overcoming the ongoing SARS-CoV-2 pandemic.

Mechanisms of Formation, Structure, and Dynamics of Lipoprotein Discs Stabilized by Amphiphilic Copolymers: A Comprehensive Review

Amphiphilic copolymers consisting of alternating hydrophilic and hydrophobic units account for a major recent methodical breakthrough in the investigations of membrane proteins. Styrene–maleic acid (SMA), diisobutylene–maleic acid (DIBMA), and related copolymers have been shown to extract membrane proteins directly from lipid membranes without the need for classical detergents. Within the particular experimental setup, they form disc-shaped nanoparticles with a narrow size distribution, which serve as a suitable platform for diverse kinds of spectroscopy and other biophysical techniques that require relatively small, homogeneous, water-soluble particles of separate membrane proteins in their native lipid environment. In recent years, copolymer-encased nanolipoparticles have been proven as suitable protein carriers for various structural biology applications, including cryo-electron microscopy (cryo-EM), small-angle scattering, and conventional and single-molecule X-ray diffraction experiments. Here, we review the current understanding of how such nanolipoparticles are formed and organized at the molecular level with an emphasis on their chemical diversity and factors affecting their size and solubilization efficiency.

Microbial Rhodopsins

The first microbial rhodopsin, a light-driven proton pump bacteriorhodopsin from Halobacterium salinarum (HsBR), was discovered in 1971. Since then, this seven-α-helical protein, comprising a retinal molecule as a cofactor, became a major driver of groundbreaking developments in membrane protein research. However, until 1999 only a few archaeal rhodopsins, acting as light-driven proton and chloride pumps and also photosensors, were known. A new microbial rhodopsin era started in 2000 when the first bacterial rhodopsin, a proton pump, was discovered. Later it became clear that there are unexpectedly many rhodopsins, and they are present in all the domains of life and even in viruses. It turned out that they execute such a diversity of functions while being "nearly the same." The incredible evolution of the research area of rhodopsins and the scientific and technological potential of the proteins is described in the review with a focus on their function-structure relationships.

Extreme dependence of Chloroflexus aggregans LOV domain thermo- and photostability on the bound flavin species

Light-oxygen-voltage (LOV) domains are common photosensory modules that found many applications in fluorescence microscopy and optogenetics. Here, we show that the Chloroflexus aggregans LOV domain can bind different flavin species (lumichrome, LC; riboflavin, RF; flavin mononucleotide, FMN; flavin adenine dinucleotide, FAD) during heterologous expression and that its physicochemical properties depend strongly on the nature of the bound flavin. We show that whereas the dissociation constants for different chromophores are similar, the melting temperature of the protein reconstituted with single flavin species varies from ~ 60 °C for LC to ~ 81 °C for FMN, and photobleaching half-times vary almost 100-fold. These observations serve as a caution for future studies of LOV domains in non-native conditions yet raise the possibility of fine-tuning various properties of LOV-based fluorescent probes and optogenetic tools by manipulating the chromophore composition.

Crystal structure of human 14-3-3ζ complexed with the noncanonical phosphopeptide from proapoptotic BAD

Several signaling pathways control phosphorylation of the proapoptotic protein BAD and its phosphorylation-dependent association with 14-3-3 proteins in the cytoplasm. The stability of the 14-3-3/BAD complex determines the cell fate: unphosphorylated BAD escapes from 14-3-3, migrates to the mitochondria and initiates apoptosis. While the 14-3-3/BAD interaction represents a promising drug target, it lacks structural characterization. Among several phosphosites identified in vivo, Ser75 and Ser99 of human BAD match the consensus sequence RXXpSXP recognized by 14-3-3 and, therefore, represent canonical 14-3-3-binding sites. Yet, BAD contains other serines phosphorylatable in vivo, whose role is less understood. Here, we report a 2.36 Å crystal structure of 14-3-3ζ complexed with a BAD fragment which includes residues Ser74 and Ser75, both being substrates for protein kinases. While the BAD peptide is anchored to 14-3-3 by phosphoserine as expected, the BAD peptide was unexpectedly phosphorylated at Ser74 instead of Ser75, revealing noncanonical binding within the amphipathic groove and leading to a one-step positional shift and reorganization of the interface. This observation exemplifies plasticity of the amphipathic 14-3-3 groove in accommodating various peptides and suggests the redundancy of Ser74 and Ser75 phosphosites with respect to binding of BAD to 14-3-3.

Inverse Conformational Selection in Lipid-Protein Binding

Interest in lipid interactions with proteins and other biomolecules is emerging not only in fundamental biochemistry but also in the field of nanobiotechnology where lipids are commonly used, for example, in carriers of mRNA vaccines. The outward-facing components of cellular membranes and lipid nanoparticles, the lipid headgroups, regulate membrane interactions with approaching substances, such as proteins, drugs, RNA, or viruses. Because lipid headgroup conformational ensembles have not been experimentally determined in physiologically relevant conditions, an essential question about their interactions with other biomolecules remains unanswered: Do headgroups exchange between a few rigid structures, or fluctuate freely across a practically continuous spectrum of conformations? Here, we combine solid-state NMR experiments and molecular dynamics simulations from the NMRlipids Project to resolve the conformational ensembles of headgroups of four key lipid types in various biologically relevant conditions. We find that lipid headgroups sample a wide range of overlapping conformations in both neutral and charged cellular membranes, and that differences in the headgroup chemistry manifest only in probability distributions of conformations. Furthermore, the analysis of 894 protein-bound lipid structures from the Protein Data Bank suggests that lipids can bind to proteins in a wide range of conformations, which are not limited by the headgroup chemistry. We propose that lipids can select a suitable headgroup conformation from the wide range available to them to fit the various binding sites in proteins. The proposed inverse conformational selection model will extend also to lipid binding to targets other than proteins, such as drugs, RNA, and viruses.

High-resolution structure of a naturally red-shifted LOV domain

LOV domains are widespread photosensory modules that have also found applications in fluorescence microscopy, optogenetics, and light-driven generation of reactive oxygen species. Many of these applications require stable proteins with altered spectra. Here, we report a flavin-based fluorescent protein CisFbFP derived from Chloroflexus islandicus LOV domain-containing protein. We show that CisFbFP is thermostable, and its absorption and fluorescence spectra are red-shifted for ∼6 nm, which has not been observed for other cysteine-substituted natural LOV domains. We also provide a crystallographic structure of CisFbFP at the resolution of 1.2 Å that reveals alterations in the active site due to replacement of conservative asparagine with a serine. Finally, we discuss the possible effects of presence of cis-proline in the Aβ-Bβ loop on the protein's structure and stability. The findings provide the basis for engineering and color tuning of LOV-based tools for molecular biology.

Nitrate- and Nitrite-Sensing Histidine Kinases: Function, Structure, and Natural Diversity

Under anaerobic conditions, bacteria may utilize nitrates and nitrites as electron acceptors. Sensitivity to nitrous compounds is achieved via several mechanisms, some of which rely on sensor histidine kinases (HKs). The best studied nitrate- and nitrite-sensing HKs (NSHKs) are NarQ and NarX from Escherichia coli. Here, we review the function of NSHKs, analyze their natural diversity, and describe the available structural information. In particular, we show that around 6000 different NSHK sequences forming several distinct clusters may now be found in genomic databases, comprising mostly the genes from Beta- and Gammaproteobacteria as well as from Bacteroidetes and Chloroflexi, including those from anaerobic ammonia oxidation (annamox) communities. We show that the architecture of NSHKs is mostly conserved, although proteins from Bacteroidetes lack the HAMP and GAF-like domains yet sometimes have PAS. We reconcile the variation of NSHK sequences with atomistic models and pinpoint the structural elements important for signal transduction from the sensor domain to the catalytic module over the transmembrane and cytoplasmic regions spanning more than 200 Å.

Role of hydrogen bond alternation and charge transfer states in photoactivation of the Orange Carotenoid Protein

Here, we propose a possible photoactivation mechanism of a 35-kDa blue light-triggered photoreceptor, the Orange Carotenoid Protein (OCP), suggesting that the reaction involves the transient formation of a protonated ketocarotenoid (oxocarbenium cation) state. Taking advantage of engineering an OCP variant carrying the Y201W mutation, which shows superior spectroscopic and structural properties, it is shown that the presence of Trp201 augments the impact of one critical H-bond between the ketocarotenoid and the protein. This confers an unprecedented homogeneity of the dark-adapted OCP state and substantially increases the yield of the excited photoproduct S*, which is important for the productive photocycle to proceed. A 1.37 Å crystal structure of OCP Y201W combined with femtosecond time-resolved absorption spectroscopy, kinetic analysis, and deconvolution of the spectral intermediates, as well as extensive quantum chemical calculations incorporating the effect of the local electric field, highlighted the role of charge-transfer states during OCP photoconversion.

Yaroshevich et al. present a chemical reaction mechanism of a 35-kDa blue light-triggered photoreceptor, the Orange Carotenoid Protein (OCP). They find that photoactivation critically involves the transient formation of a protonated ketocarotenoid (oxocarbenium cation) state. This study suggests the role of charge-transfer states during OCP photoconversion.

The molecular basis of spectral tuning in blue- and red-shifted flavin-binding fluorescent proteins

Photoactive biological systems modify the optical properties of their chromophores, known as spectral tuning. Determining the molecular origin of spectral tuning is instrumental for understanding the function and developing applications of these biomolecules. Spectral tuning in flavin-binding fluorescent proteins (FbFPs), an emerging class of fluorescent reporters, is limited by their dependency on protein-bound flavins, whose structure and hence electronic properties cannot be altered by mutation. A blue-shifted variant of the plant-derived improved light, oxygen, voltage FbFP has been created by introducing a lysine within the flavin-binding pocket, but the molecular basis of this shift remains unconfirmed. We here structurally characterize the blue-shifted improved light, oxygen, voltage variant and construct a new blue-shifted CagFbFP protein by introducing an analogous mutation. X-ray structures of both proteins reveal displacement of the lysine away from the chromophore and opening up of the structure as instrumental for the blue shift. Site saturation mutagenesis and high-throughput screening yielded a red-shifted variant, and structural analysis revealed that the lysine side chain of the blue-shifted variant is stabilized close to the flavin by a secondary mutation, accounting for the red shift. Thus, a single additional mutation in a blue-shifted variant is sufficient to generate a red-shifted FbFP. Using spectroscopy, X-ray crystallography, and quantum mechanics molecular mechanics calculations, we provide a firm structural and functional understanding of spectral tuning in FbFPs. We also show that the identified blue- and red-shifted variants allow for two-color microscopy based on spectral separation. In summary, the generated blue- and red-shifted variants represent promising new tools for application in life sciences.

Insights into the mechanisms of light-oxygen-voltage domain color tuning from a set of high-resolution X-ray structures

Light-oxygen-voltage (LOV) domains are widespread photosensory modules that can be used in fluorescence microscopy, optogenetics and controlled production of reactive oxygen species. All of the currently known LOV domains have absorption maxima in the range of ~440 to ~450 nm, and it is not clear whether they can be shifted significantly using mutations. Here, we have generated a panel of LOV domain variants by mutating the key chromophore-proximal glutamine aminoacid of a thermostable flavin based fluorescent protein CagFbFP (Gln148) to asparagine, aspartate, glutamate, histidine, lysine and arginine. Absorption spectra of all of the mutants are blue-shifted, with the maximal shift of 8 nm observed for the Q148H variant. While CagFbFP and its Q148N/D/E variants are not sensitive to pH, Q148H/K/R reveal a moderate red shift induced byacidic pH. To gain further insight, we determined high resolution crystal structures of all of the mutants studied at the resolutions from 1.07 Å for Q148D to 1.63 Å for Q148R. Whereas in some of the variants, the aminoacid 148 remains in the vicinity of the flavin, in Q148K, Q148R and partially Q148D, the C-terminus of the protein unlatches and the side chain of the residue 148 is reoriented away from the chromophore. Our results explain the absence of color shifts from replacing Gln148 with charged aminoacids and pave the way for rational design of color-shifted flavin based fluorescent proteins.

Rational Design of a Split Flavin-Based Fluorescent Reporter

Protein-fragment complementation assays are used ubiquitously for probing protein-protein interactions. Most commonly, the reporter protein is split in two parts, which are then fused to the proteins of interest and can reassemble and provide a readout if the proteins of interest interact with each other. The currently known split fluorescent proteins either can be used only in aerobic conditions and assemble irreversibly, or require addition of exogenous chromophores, which complicates the design of experiments. In recent years, light-oxygen-voltage (LOV) domains of several photoreceptor proteins have been developed into flavin-based fluorescent proteins (FbFPs) that, under some circumstances, can outperform commonly used fluorescent proteins such as GFP. Here, we show that CagFbFP, a small thermostable FbFP based on a LOV domain-containing protein from Chloroflexus aggregans, can serve as a split fluorescent reporter. We use the available genetic and structural information to identify three loops between the conserved secondary structure elements, Aβ-Bβ, Eα-Fα, and Hβ-Iβ, that tolerate insertion of flexible poly-Gly/Ser segments and eventually splitting. We demonstrate that the designed split pairs, when fused to interacting proteins, are fluorescent in vivo in E. coli and human cells and have low background fluorescence. Our results enable probing protein-protein interactions in anaerobic conditions without using exogenous fluorophores and provide a basis for further development of LOV and PAS (Per-Arnt-Sim) domain-based fluorescent reporters and optogenetic tools.

Principal component analysis highlights the influence of temperature, curvature and cholesterol on conformational dynamics of lipids

Membrane lipids are inherently highly dynamic molecules. Currently, it is difficult to probe the structures of individual lipids experimentally at the timescales corresponding to atomic motions, and consequently molecular dynamics simulations are used widely. In our previous work, we have introduced the principal component analysis (PCA) as a convenient framework for comprehensive quantitative description of lipid motions. Here, we present a newly developed open source script, PCAlipids, which automates the analysis and allows us to refine the approach and test its limitations. We use PCAlipids to determine the influence of temperature, cholesterol and curvature on individual lipids, and show that the most prominent lipid tail scissoring motion is strongly affected by these factors and allows tracking of phase transition. Addition of cholesterol affects the conformations and selectively changes the dynamics of lipid molecules, impacting the large-amplitude motions. Introduction of curvature biases the conformational ensembles towards more extended structures. We hope that the developed approach will be useful for understanding the molecular basis of different processes occurring in lipid membrane systems and will stimulate development of complementary experimental techniques probing the conformations of individual lipid molecules.

Molecular mechanism of light-driven sodium pumping

The light-driven sodium-pumping rhodopsin KR2 from Krokinobacter eikastus is the only non-proton cation active transporter with demonstrated potential for optogenetics. However, the existing structural data on KR2 correspond exclusively to its ground state, and show no sodium inside the protein, which hampers the understanding of sodium-pumping mechanism. Here we present crystal structure of the O-intermediate of the physiologically relevant pentameric form of KR2 at the resolution of 2.1 Å, revealing a sodium ion near the retinal Schiff base, coordinated by N112 and D116 of the characteristic NDQ triad. We also obtained crystal structures of D116N and H30A variants, conducted metadynamics simulations and measured pumping activities of putative pathway mutants to demonstrate that sodium release likely proceeds alongside Q78 towards the structural sodium ion bound between KR2 protomers. Our findings highlight the importance of pentameric assembly for sodium pump function, and may be used for rational engineering of enhanced optogenetic tools.

Sensor Histidine Kinase NarQ Activates via Helical Rotation, Diagonal Scissoring, and Eventually Piston-Like Shifts

Membrane-embedded sensor histidine kinases (HKs) and chemoreceptors are used ubiquitously by bacteria and archaea to percept the environment, and are often crucial for their survival and pathogenicity. The proteins can transmit the signal from the sensor domain to the catalytic kinase domain reliably over the span of several hundreds of angstroms, and regulate the activity of the cognate response regulator proteins, with which they form two-component signaling systems (TCSs). Several mechanisms of transmembrane signal transduction in TCS receptors have been proposed, dubbed (swinging) piston, helical rotation, and diagonal scissoring. Yet, despite decades of studies, there is no consensus on whether these mechanisms are common for all TCS receptors. Here, we extend our previous work on Escherichia coli nitrate/nitrite sensor kinase NarQ. We determined a crystallographic structure of the sensor-TM-HAMP fragment of the R50S mutant, which, unexpectedly, was found in a ligand-bound-like conformation, despite an inability to bind nitrate. Subsequently, we reanalyzed the structures of the ligand-free and ligand-bound NarQ and NarX sensor domains, and conducted extensive molecular dynamics simulations of ligand-free and ligand-bound wild type and mutated NarQ. Based on the data, we show that binding of nitrate to NarQ causes, first and foremost, helical rotation and diagonal scissoring of the α-helices at the core of the sensor domain. These conformational changes are accompanied by a subtle piston-like motion, which is amplified by a switch in the secondary structure of the linker between the sensor and TM domains. We conclude that helical rotation, diagonal scissoring, and piston are simply different degrees of freedom in coiled-coil proteins and are not mutually exclusive in NarQ, and likely in other nitrate sensors and TCS proteins as well.

High-resolution structural insights into the heliorhodopsin family

Rhodopsins are the most abundant light-harvesting proteins. A new family of rhodopsins, heliorhodopsins (HeRs), has recently been discovered. Unlike in the known rhodopsins, in HeRs the N termini face the cytoplasm. The function of HeRs remains unknown. We present the structures of the bacterial HeR-48C12 in two states at the resolution of 1.5 Å, which highlight its remarkable difference from all known rhodopsins. The interior of HeR's extracellular part is completely hydrophobic, while the cytoplasmic part comprises a cavity (Schiff base cavity [SBC]) surrounded by charged amino acids and containing a cluster of water molecules, presumably being a primary proton acceptor from the Schiff base. At acidic pH, a planar triangular molecule (acetate) is present in the SBC. Structure-based bioinformatic analysis identified 10 subfamilies of HeRs, suggesting their diverse biological functions. The structures and available data suggest an enzymatic activity of HeR-48C12 subfamily and their possible involvement in fundamental redox biological processes.

Assembly of Spinach Chloroplast ATP Synthase Rotor Ring Protein-Lipid Complex

Rotor ATPases are large multisubunit membrane protein complexes found in all kingdoms of life. The membrane parts of these ATPases include a ring-like assembly, so-called c-ring, consisting of several subunits c, plugged by a patch of phospholipids. In this report, we use a nature-inspired approach to model the assembly of the spinach (Spinacia oleracea) c₁₄ ring protein-lipid complex, where partially assembled oligomers are pulled toward each other using a biasing potential. The resulting assemblies contain 23 to 26 encapsulated plug lipids, general position of which corresponds well to experimental maps. However, best fit to experimental data is achieved with 15 to 17 lipids inside the c-ring. In all of the simulations, the lipids from one leaflet (loop side of the c subunit) are ordered and static, whereas the lipids from the other leaflet are disordered and dynamic. Spontaneous permeation of water molecules toward Glu61 at the active site is also observed. The presented assembly approach is expected to be generalizable to other protein complexes with encapsulated lipid patches.

Structure and mechanisms of sodium-pumping KR2 rhodopsin

Rhodopsins are the most universal biological light-energy transducers and abundant phototrophic mechanisms that evolved on Earth and have a remarkable diversity and potential for biotechnological applications. Recently, the first sodium-pumping rhodopsin KR2 from Krokinobacter eikastus was discovered and characterized. However, the existing structures of KR2 are contradictory, and the mechanism of Na⁺ pumping is not yet understood. Here, we present a structure of the cationic (non H⁺) light-driven pump at physiological pH in its pentameric form. We also present 13 atomic structures and functional data on the KR2 and its mutants, including potassium pumps, which show that oligomerization of the microbial rhodopsin is obligatory for its biological function. The studies reveal the structure of KR2 at nonphysiological low pH where it acts as a proton pump. The structure provides new insights into the mechanisms of microbial rhodopsins and opens the way to a rational design of novel cation pumps for optogenetics.

A thermostable flavin-based fluorescent protein from Chloroflexus aggregans: a framework for ultra-high resolution structural studies

A new thermostable fluorescent protein is shown to be a promising model for ultra-high resolution structural studies of LOV domains and for application as a fluorescent reporter.

Efficient non-cytotoxic fluorescent staining of halophiles

Research on halophilic microorganisms is important due to their relation to fundamental questions of survival of living organisms in a hostile environment. Here we introduce a novel method to stain halophiles with MitoTracker fluorescent dyes in their growth medium. The method is based on membrane-potential sensitive dyes, which were originally used to label mitochondria in eukaryotic cells. We demonstrate that these fluorescent dyes provide high staining efficiency and are beneficial for multi-staining purposes due to the spectral range covered (from orange to deep red). In contrast with other fluorescent dyes used so far, MitoTracker does not affect growth rate, and remains in cells after several washing steps and several generations in cell culture. The suggested dyes were tested on three archaeal (Hbt. salinarum, Haloferax sp., Halorubrum sp.) and two bacterial (Salicola sp., Halomonas sp.) strains of halophilic microorganisms. The new staining approach provides new insights into biology of Hbt. salinarum. We demonstrated the interconversion of rod-shaped cells of Hbt. salinarium to spheroplasts and submicron-sized spheres, as well as the cytoplasmic integrity of giant rod Hbt. salinarum species. By expanding the variety of tools available for halophile detection, MitoTracker dyes overcome long-standing limitations in fluorescence microscopy studies of halophiles.

Microbial Rhodopsins

Microbial rhodopsins (MRs) are a large family of photoactive membrane proteins, found in microorganisms belonging to all kingdoms of life, with new members being constantly discovered. Among the MRs are light-driven proton, cation and anion pumps, light-gated cation and anion channels, and various photoreceptors. Due to their abundance and amenability to studies, MRs served as model systems for a great variety of biophysical techniques, and recently found a great application as optogenetic tools. While the basic aspects of microbial rhodopsins functioning have been known for some time, there is still a plenty of unanswered questions. This chapter presents and summarizes the available knowledge, focusing on the functional and structural studies.

Transmembrane Signal Transduction in Two-Component Systems: Piston, Scissoring, or Helical Rotation?

Allosteric and transmembrane (TM) signaling are among the major questions of structural biology. Here, we review and discuss signal transduction in four-helical TM bundles, focusing on histidine kinases and chemoreceptors found in two-component systems. Previously, piston, scissors, and helical rotation have been proposed as the mechanisms of TM signaling. We discuss theoretically possible conformational changes and examine the available experimental data, including the recent crystallographic structures of nitrate/nitrite sensor histidine kinase NarQ and phototaxis system NpSRII:NpHtrII. We show that TM helices can flex at multiple points and argue that the various conformational changes are not mutually exclusive, and often are observed concomitantly, throughout the TM domain or in its part. The piston and scissoring motions are the most prominent motions in the structures, but more research is needed for definitive conclusions.

Effects of Coarse Graining and Saturation of Hydrocarbon Chains on Structure and Dynamics of Simulated Lipid Molecules

Molecular dynamics simulations are used extensively to study the processes on biological membranes. The simulations can be conducted at different levels of resolution: all atom (AA), where all atomistic details are provided; united atom (UA), where hydrogen atoms are treated inseparably of corresponding heavy atoms; and coarse grained (CG), where atoms are grouped into larger particles. Here, we study the behavior of model bilayers consisting of saturated and unsaturated lipids DOPC, SOPC, OSPC and DSPC in simulations performed using all atom CHARMM36 and coarse grained Martini force fields. Using principal components analysis, we show that the structural and dynamical properties of the lipids are similar, both in AA and CG simulations, although the unsaturated molecules are more dynamic and favor more extended conformations. We find that CG simulations capture 75 to 100% of the major collective motions, overestimate short range ordering, result in more flexible molecules and 5–7 fold faster sampling. We expect that the results reported here will be useful for comprehensive quantitative comparisons of simulations conducted at different resolution levels and for further development and improvement of CG force fields.

Inward H+ pump xenorhodopsin: Mechanism and alternative optogenetic approach

Generation of an electrochemical proton gradient is the first step of cell bioenergetics. In prokaryotes, the gradient is created by outward membrane protein proton pumps. Inward plasma membrane native proton pumps are yet unknown. We describe comprehensive functional studies of the representatives of the yet noncharacterized xenorhodopsins from Nanohaloarchaea family of microbial rhodopsins. They are inward proton pumps as we demonstrate in model membrane systems, Escherichia coli cells, human embryonic kidney cells, neuroblastoma cells, and rat hippocampal neuronal cells. We also solved the structure of a xenorhodopsin from the nanohalosarchaeon Nanosalina (NsXeR) and suggest a mechanism of inward proton pumping. We demonstrate that the NsXeR is a powerful pump, which is able to elicit action potentials in rat hippocampal neuronal cells up to their maximal intrinsic firing frequency. Hence, inwardly directed proton pumps are suitable for light-induced remote control of neurons, and they are an alternative to the well-known cation-selective channelrhodopsins.

Mechanism of transmembrane signaling by sensor histidine kinases

One of the major and essential classes of transmembrane (TM) receptors, present in all domains of life, is sensor histidine kinases, parts of two-component signaling systems (TCSs). The structural mechanisms of TM signaling by these sensors are poorly understood. We present crystal structures of the periplasmic sensor domain, the TM domain, and the cytoplasmic HAMP domain of the Escherichia coli nitrate/nitrite sensor histidine kinase NarQ in the ligand-bound and mutated ligand-free states. The structures reveal that the ligand binding induces rearrangements and pistonlike shifts of TM helices. The HAMP domain protomers undergo leverlike motions and convert these pistonlike motions into helical rotations. Our findings provide the structural framework for complete understanding of TM TCS signaling and for development of antimicrobial treatments targeting TCSs.

Fast iodide-SAD phasing for high-throughput membrane protein structure determination

We describe a fast, easy, and potentially universal method for the de novo solution of the crystal structures of membrane proteins via iodide-single-wavelength anomalous diffraction (I-SAD). The potential universality of the method is based on a common feature of membrane proteins-the availability at the hydrophobic-hydrophilic interface of positively charged amino acid residues with which iodide strongly interacts. We demonstrate the solution using I-SAD of four crystal structures representing different classes of membrane proteins, including a human G protein-coupled receptor (GPCR), and we show that I-SAD can be applied using data collection strategies based on either standard or serial x-ray crystallography techniques.

New Insights on Signal Propagation by Sensory Rhodopsin II/Transducer Complex

The complex of two membrane proteins, sensory rhodopsin II (NpSRII) with its cognate transducer (NpHtrII), mediates negative phototaxis in halobacteria N. pharaonis. Upon light activation NpSRII triggers a signal transduction chain homologous to the two-component system in eubacterial chemotaxis. Here we report on crystal structures of the ground and active M-state of the complex in the space group I2₁2₁2₁. We demonstrate that the relative orientation of symmetrical parts of the dimer is parallel (“U”-shaped) contrary to the gusset-like (“V”-shaped) form of the previously reported structures of the NpSRII/NpHtrII complex in the space group P2₁2₁2, although the structures of the monomers taken individually are nearly the same. Computer modeling of the HAMP domain in the obtained “V”- and “U”-shaped structures revealed that only the “U”-shaped conformation allows for tight interactions of the receptor with the HAMP domain. This is in line with existing data and supports biological relevance of the “U” shape in the ground state. We suggest that the “V”-shaped structure may correspond to the active state of the complex and transition from the “U” to the “V”-shape of the receptor-transducer complex can be involved in signal transduction from the receptor to the signaling domain of NpHtrII.

Principal Component Analysis of Lipid Molecule Conformational Changes in Molecular Dynamics Simulations

Molecular dynamics simulations of lipid bilayers are ubiquitous nowadays. Usually, either global properties of the bilayer or some particular characteristics of each lipid molecule are evaluated in such simulations, but the structural properties of the molecules as a whole are rarely studied. Here, we show how a comprehensive quantitative description of conformational space and dynamics of a single lipid molecule can be achieved via the principal component analysis (PCA). We illustrate the approach by analyzing and comparing simulations of DOPC bilayers obtained using eight different force fields: all-atom generalized AMBER, CHARMM27, CHARMM36, Lipid14, and Slipids and united-atom Berger, GROMOS43A1-S3, and GROMOS54A7. Similarly to proteins, most of the structural variance of a lipid molecule can be described by only a few principal components. These major components are similar in different simulations, although there are notable distinctions between the older and newer force fields and between the all-atom and united-atom force fields. The DOPC molecules in the simulations generally equilibrate on the time scales of tens to hundreds of nanoseconds. The equilibration is the slowest in the GAFF simulation and the fastest in the Slipids simulation. Somewhat unexpectedly, the equilibration in the united-atom force fields is generally slower than in the all-atom force fields. Overall, there is a clear separation between the more variable previous generation force fields and significantly more similar new generation force fields (CHARMM36, Lipid14, Slipids). We expect that the presented approaches will be useful for quantitative analysis of conformations and dynamics of individual lipid molecules in other simulations of lipid bilayers.

An Approach to Heterologous Expression of Membrane Proteins. The Case of Bacteriorhodopsin

Heterologous overexpression of functional membrane proteins is a major bottleneck of structural biology. Bacteriorhodopsin from Halobium salinarum (bR) is a striking example of the difficulties in membrane protein overexpression. We suggest a general approach with a finite number of steps which allows one to localize the underlying problem of poor expression of a membrane protein using bR as an example. Our approach is based on constructing chimeric proteins comprising parts of a protein of interest and complementary parts of a homologous protein demonstrating advantageous expression. This complementary protein approach allowed us to increase bR expression by two orders of magnitude through the introduction of two silent mutations into bR coding DNA. For the first time the high quality crystals of bR expressed in E. Coli were obtained using the produced protein. The crystals obtained with in meso nanovolume crystallization diffracted to 1.67 Å.

Crystal structure of a light-driven sodium pump

Recently, the first known light-driven sodium pumps, from the microbial rhodopsin family, were discovered. We have solved the structure of one of them, Krokinobacter eikastus rhodopsin 2 (KR2), in the monomeric blue state and in two pentameric red states, at resolutions of 1.45 Å and 2.2 and 2.8 Å, respectively. The structures reveal the ion-translocation pathway and show that the sodium ion is bound outside the protein at the oligomerization interface, that the ion-release cavity is capped by a unique N-terminal α-helix and that the ion-uptake cavity is unexpectedly large and open to the surface. Obstruction of the cavity with the mutation G263F imparts KR2 with the ability to pump potassium. These results pave the way for the understanding and rational design of cation pumps with new specific properties valuable for optogenetics.

Structural and functional investigation of flavin binding center of the NqrC subunit of sodium-translocating NADH:quinone oxidoreductase from Vibrio harveyi

Na⁺-translocating NADH:quinone oxidoreductase (NQR) is a redox-driven sodium pump operating in the respiratory chain of various bacteria, including pathogenic species. The enzyme has a unique set of redox active prosthetic groups, which includes two covalently bound flavin mononucleotide (FMN) residues attached to threonine residues in subunits NqrB and NqrC. The reason of FMN covalent bonding in the subunits has not been established yet. In the current work, binding of free FMN to the apo-form of NqrC from Vibrio harveyi was studied showing very low affinity of NqrC to FMN in the absence of its covalent bonding. To study structural aspects of flavin binding in NqrC, its holo-form was crystallized and its 3D structure was solved at 1.56 Å resolution. It was found that the isoalloxazine moiety of the FMN residue is buried in a hydrophobic cavity and that its pyrimidine ring is squeezed between hydrophobic amino acid residues while its benzene ring is extended from the protein surroundings. This structure of the flavin-binding pocket appears to provide flexibility of the benzene ring, which can help the FMN residue to take the bended conformation and thus to stabilize the one-electron reduced form of the prosthetic group. These properties may also lead to relatively weak noncovalent binding of the flavin. This fact along with periplasmic location of the FMN-binding domains in the vast majority of NqrC-like proteins may explain the necessity of the covalent bonding of this prosthetic group to prevent its loss to the external medium.

Crystal structure of Escherichia coli-expressed Haloarcula marismortui bacteriorhodopsin I in the trimeric form

Bacteriorhodopsins are a large family of seven-helical transmembrane proteins that function as light-driven proton pumps. Here, we present the crystal structure of a new member of the family, Haloarcula marismortui bacteriorhodopsin I (HmBRI) D94N mutant, at the resolution of 2.5 Å. While the HmBRI retinal-binding pocket and proton donor site are similar to those of other archaeal proton pumps, its proton release region is extended and contains additional water molecules. The protein's fold is reinforced by three novel inter-helical hydrogen bonds, two of which result from double substitutions relative to Halobacterium salinarum bacteriorhodopsin and other similar proteins. Despite the expression in Escherichia coli and consequent absence of native lipids, the protein assembles as a trimer in crystals. The unique extended loop between the helices D and E of HmBRI makes contacts with the adjacent protomer and appears to stabilize the interface. Many lipidic hydrophobic tail groups are discernible in the membrane region, and their positions are similar to those of archaeal isoprenoid lipids in the crystals of other proton pumps, isolated from native or native-like sources. All these features might explain the HmBRI properties and establish the protein as a novel model for the microbial rhodopsin proton pumping studies.

Low-dose X-ray radiation induces structural alterations in proteins

X-ray-radiation-induced alterations to protein structures are still a severe problem in macromolecular crystallography. One way to avoid the influence of radiation damage is to reduce the X-ray dose absorbed by the crystal during data collection. However, here it is demonstrated using the example of the membrane protein bacteriorhodopsin (bR) that even a low dose of less than 0.06 MGy may induce structural alterations in proteins. This dose is about 500 times smaller than the experimental dose limit which should ideally not be exceeded per data set (i.e. 30 MGy) and 20 times smaller than previously detected specific radiation damage at the bR active site. To date, it is the lowest dose at which radiation modification of a protein structure has been described. Complementary use was made of high-resolution X-ray crystallography and online microspectrophotometry to quantitatively study low-dose X-ray-induced changes. It is shown that structural changes of the protein correlate with the spectroscopically observed formation of the so-called bR orange species. Evidence is provided for structural modifications taking place at the protein active site that should be taken into account in crystallographic studies which aim to elucidate the molecular mechanisms of bR function.

Fourier modal method for relief gratings with oblique boundary conditions

Oblique boundary conditions are introduced in the Fourier modal method at each slice of the staircase decomposition of an arbitrary profile of a dielectric corrugation grating. The precision and convergence improvement are demonstrated by comparison with reference methods.