Structure determination of high-energy states in a dynamic protein ensemble.

Macromolecular function frequently requires that proteins change conformation into high-energy states1-4. However, methods for solving the structures of these functionally essential, lowly populated states are lacking. Here we develop a method for high-resolution structure determination of minorly populated states by coupling NMR spectroscopy-derived pseudocontact shifts5 (PCSs) with Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion6 (PCS-CPMG). Our approach additionally defines the corresponding kinetics and thermodynamics of high-energy excursions, thereby characterizing the entire free-energy landscape. Using a large set of simulated data for adenylate kinase (Adk), calmodulin and Src kinase, we find that high-energy PCSs accurately determine high-energy structures (with a root mean squared deviation of less than 3.5 angström). Applying our methodology to Adk during catalysis, we find that the high-energy excursion involves surprisingly small openings of the AMP and ATP lids. This previously unresolved high-energy structure solves a longstanding controversy about conformational interconversions that are rate-limiting for catalysis. Primed for either substrate binding or product release, the high-energy structure of Adk suggests a two-step mechanism combining conformational selection to this state, followed by an induced-fit step into a fully closed state for catalysis of the phosphoryl-transfer reaction. Unlike other methods for resolving high-energy states, such as cryo-electron microscopy and X-ray crystallography, our solution PCS-CPMG approach excels in cases involving domain rearrangements of smaller systems (less than 60 kDa) and populations as low as 0.5%, and enables the simultaneous determination of protein structure, kinetics and thermodynamics while proteins perform their function.

Extant Sequence Reconstruction: The Accuracy of Ancestral Sequence Reconstructions Evaluated by Extant Sequence Cross-Validation.

Ancestral sequence reconstruction (ASR) is a phylogenetic method widely used to analyze the properties of ancient biomolecules and to elucidate mechanisms of molecular evolution. Despite its increasingly widespread application, the accuracy of ASR is currently unknown, as it is generally impossible to compare resurrected proteins to the true ancestors. Which evolutionary models are best for ASR? How accurate are the resulting inferences? Here we answer these questions using a cross-validation method to reconstruct each extant sequence in an alignment with ASR methodology, a method we term "extant sequence reconstruction" (ESR). We thus can evaluate the accuracy of ASR methodology by comparing ESR reconstructions to the corresponding known true sequences. We find that a common measure of the quality of a reconstructed sequence, the average probability, is indeed a good estimate of the fraction of correct amino acids when the evolutionary model is accurate or overparameterized. However, the average probability is a poor measure for comparing reconstructions from different models, because, surprisingly, a more accurate phylogenetic model often results in reconstructions with lower probability. While better (more predictive) models may produce reconstructions with lower sequence identity to the true sequences, better models nevertheless produce reconstructions that are more biophysically similar to true ancestors. In addition, we find that a large fraction of sequences sampled from the reconstruction distribution may have fewer errors than the single most probable (SMP) sequence reconstruction, despite the fact that the SMP has the lowest expected error of all possible sequences. Our results emphasize the importance of model selection for ASR and the usefulness of sampling sequence reconstructions for analyzing ancestral protein properties. ESR is a powerful method for validating the evolutionary models used for ASR and can be applied in practice to any phylogenetic analysis of real biological sequences. Most significantly, ESR uses ASR methodology to provide a general method by which the biophysical properties of resurrected proteins can be compared to the properties of the true protein.

A gustatory receptor paralogue controls rapid warmth avoidance in Drosophila.

Behavioural responses to temperature are critical for survival, and animals from insects to humans show strong preferences for specific temperatures. Preferred temperature selection promotes avoidance of adverse thermal environments in the short term and maintenance of optimal body temperatures over the long term, but its molecular and cellular basis is largely unknown. Recent studies have generated conflicting views of thermal preference in Drosophila, attributing importance to either internal or peripheral warmth sensors. Here we reconcile these views by showing that thermal preference is not a singular response, but involves multiple systems relevant in different contexts. We found previously that the transient receptor potential channel TRPA1 acts internally to control the slowly developing preference response of flies exposed to a shallow thermal gradient. We now find that the rapid response of flies exposed to a steep warmth gradient does not require TRPA1; rather, the gustatory receptor GR28B(D) drives this behaviour through peripheral thermosensors. Gustatory receptors are a large gene family, widely studied in insect gustation and olfaction, and are implicated in host-seeking by insect disease vectors, but have not previously been implicated in thermosensation. At the molecular level, GR28B(D) misexpression confers thermosensitivity upon diverse cell types, suggesting that it is a warmth sensor. These data reveal a new type of thermosensory molecule and uncover a functional distinction between peripheral and internal warmth sensors in this tiny ectotherm reminiscent of thermoregulatory systems in larger, endothermic animals. The use of multiple, distinct molecules to respond to a given temperature, as observed here, may facilitate independent tuning of an animal's distinct thermosensory responses.

Functional and Structural Resilience of the Active Site Loop in the Evolution of Plasmodium Lactate Dehydrogenase.

The malarial pathogen Plasmodium falciparum ( Pf) is a member of the Apicomplexa, which independently evolved a highly specific lactate dehydrogenase (LDH) from an ancestral malate dehydrogenase (MDH) via a five-residue insertion in a key active site loop. PfLDH is widely considered an attractive drug target because of its unique active site. The conservation of the apicomplexan loop suggests that a precise insertion sequence was required for the evolution of LDH specificity. Aside from a single critical tryptophan, W107f, the functional and structural roles of residues in the loop are currently unknown. Here we show that the loop is remarkably robust to mutation, as activity is resilient to radical perturbations of both loop identity and length. Thus, alternative insertions could have evolved LDH specificity as long as they contained a tryptophan in the proper location. PfLDH likely has great potential to develop resistance to drugs designed to target its distinctive active site loop.

Expression, purification, and spectral tuning of RhoGC, a retinylidene/guanylyl cyclase fusion protein and optogenetics tool from the aquatic fungus Blastocladiella emersonii.

RhoGC is a rhodopsin (Rho)-guanylyl cyclase (GC) gene fusion molecule that is central to zoospore phototaxis in the aquatic fungus Blastocladiella emersonii It has generated considerable excitement because of its demonstrated potential as a tool for optogenetic manipulation of cell-signaling pathways involving cyclic nucleotides. However, a reliable method for expressing and purifying RhoGC is currently lacking. We present here an expression and purification system for isolation of the full-length RhoGC protein expressed in HEK293 cells in detergent solution. The protein exhibits robust light-dependent guanylyl cyclase activity, whereas a truncated form lacking the 17- to 20-kDa N-terminal domain is completely inactive under identical conditions. Moreover, we designed several RhoGC mutants to increase the utility of the protein for optogenetic studies. The first class we generated has altered absorption spectra designed for selective activation by different wavelengths of light. Two mutants were created with blue-shifted (E254D, λmax = 390 nm; D380N, λmax = 506 nm) and one with red-shifted (D380E, λmax = 533 nm) absorption maxima relative to the wild-type protein (λmax = 527 nm). We also engineered a double mutant, E497K/C566D, that changes the enzyme to a specific, light-stimulated adenylyl cyclase that catalyzes the formation of cAMP from ATP. We anticipate that this expression/purification system and these RhoGC mutants will facilitate mechanistic and structural exploration of this important enzyme.

Purification and Characterization of RhoPDE, a Retinylidene/Phosphodiesterase Fusion Protein and Potential Optogenetic Tool from the Choanoflagellate Salpingoeca rosetta.

RhoPDE is a type I rhodopsin/phosphodiesterase gene fusion product from the choanoflagellate Salpingoeca rosetta. The gene was discovered around the time that a similar type I rhodopsin/guanylyl cyclase fusion protein, RhoGC, was shown to control phototaxis of an aquatic fungus through a cGMP signaling pathway. RhoPDE has potential as an optogenetic tool catalyzing the hydrolysis of cyclic nucleotides. Here we provide an expression and purification system for RhoPDE, as well as a crystal structure of the C-terminal phosphodiesterase catalytic domain. We show that RhoPDE contains an even number of transmembrane segments, with N- and C-termini both located on the cytoplasmic surface of the cell membrane. The purified protein exhibits an absorption maximum at 490 nm in the dark state, which shifts to 380 nm upon exposure to light. The protein acts as a cGMP-selective phosphodiesterase. However, the activity does not appear to be modulated by light. The protein is also active with cAMP as a substrate, but with a roughly 5-7-fold lower kcat. A truncation consisting solely of the phosphodiesterase domain is also active with a kcat for cGMP roughly 6-9-fold lower than that of the full-length protein. The isolated PDE domain was crystallized, and the X-ray structure showed the protein to be a dimer similar to human PDE9. We anticipate that the purification system introduced here will enable further structural and biochemical experiments to improve our understanding of the function and mechanism of this unique fusion protein.

Relocating the Active-Site Lysine in Rhodopsin: 2. Evolutionary Intermediates.

The visual pigment rhodopsin is a G protein-coupled receptor that covalently binds its retinal chromophore via a Schiff base linkage to an active-site Lys residue in the seventh transmembrane helix. Although this residue is strictly conserved among all type II retinylidene proteins, we found previously that the active-site Lys in bovine rhodopsin (Lys296) can be moved to three other locations (G90K, T94K, S186K) while retaining the ability to form a pigment with retinal and to activate transducin in a light-dependent manner [ Devine et al. ( 2013 ) Proc. Natl. Acad. Sci. USA 110 , 13351 - 13355 ]. Because the active-site Lys is not functionally constrained to be in helix seven, it is possible that it could relocate within the protein, most likely via an evolutionary intermediate with two active-site Lys. Therefore, in this study we characterized potential evolutionary intermediates with two Lys in the active site. Four mutant rhodopsins were prepared in which the original Lys296 was left untouched and a second Lys residue was substituted for G90K, T94K, S186K, or F293K. All four constructs covalently bind 11-cis-retinal, form a pigment, and activate transducin in a light-dependent manner. These results demonstrate that rhodopsin can tolerate a second Lys in the retinal binding pocket and suggest that an evolutionary intermediate with two Lys could allow migration of the Schiff base Lys to a position other than the observed, highly conserved location in the seventh TM helix. From sequence-based searches, we identified two groups of natural opsins, insect UV cones and neuropsins, that contain Lys residues at two positions in their active sites and also have intriguing spectral similarities to the mutant rhodopsins studied here.

A formal test of the theory of universal common ancestry.

Universal common ancestry (UCA) is a central pillar of modern evolutionary theory. As first suggested by Darwin, the theory of UCA posits that all extant terrestrial organisms share a common genetic heritage, each being the genealogical descendant of a single species from the distant past. The classic evidence for UCA, although massive, is largely restricted to 'local' common ancestry-for example, of specific phyla rather than the entirety of life-and has yet to fully integrate the recent advances from modern phylogenetics and probability theory. Although UCA is widely assumed, it has rarely been subjected to formal quantitative testing, and this has led to critical commentary emphasizing the intrinsic technical difficulties in empirically evaluating a theory of such broad scope. Furthermore, several researchers have proposed that early life was characterized by rampant horizontal gene transfer, leading some to question the monophyly of life. Here I provide the first, to my knowledge, formal, fundamental test of UCA, without assuming that sequence similarity implies genetic kinship. I test UCA by applying model selection theory to molecular phylogenies, focusing on a set of ubiquitously conserved proteins that are proposed to be orthologous. Among a wide range of biological models involving the independent ancestry of major taxonomic groups, the model selection tests are found to overwhelmingly support UCA irrespective of the presence of horizontal gene transfer and symbiotic fusion events. These results provide powerful statistical evidence corroborating the monophyly of all known life.

Analysis of Drosophila TRPA1 reveals an ancient origin for human chemical nociception.

Chemical nociception, the detection of tissue-damaging chemicals, is important for animal survival and causes human pain and inflammation, but its evolutionary origins are largely unknown. Reactive electrophiles are a class of noxious compounds humans find pungent and irritating, such as allyl isothiocyanate (in wasabi) and acrolein (in cigarette smoke). Diverse animals, from insects to humans, find reactive electrophiles aversive, but whether this reflects conservation of an ancient sensory modality has been unclear. Here we identify the molecular basis of reactive electrophile detection in flies. We demonstrate that Drosophila TRPA1 (Transient receptor potential A1), the Drosophila melanogaster orthologue of the human irritant sensor, acts in gustatory chemosensors to inhibit reactive electrophile ingestion. We show that fly and mosquito TRPA1 orthologues are molecular sensors of electrophiles, using a mechanism conserved with vertebrate TRPA1s. Phylogenetic analyses indicate that invertebrate and vertebrate TRPA1s share a common ancestor that possessed critical characteristics required for electrophile detection. These findings support emergence of TRPA1-based electrophile detection in a common bilaterian ancestor, with widespread conservation throughout vertebrate and invertebrate evolution. Such conservation contrasts with the evolutionary divergence of canonical olfactory and gustatory receptors and may relate to electrophile toxicity. We propose that human pain perception relies on an ancient chemical sensor conserved across approximately 500 million years of animal evolution.

Relocating the active-site lysine in rhodopsin and implications for evolution of retinylidene proteins.

Type I and type II rhodopsins share several structural features including a G protein-coupled receptor fold and a highly conserved active-site Lys residue in the seventh transmembrane segment of the protein. However, the two families lack significant sequence similarity that would indicate common ancestry. Consequently, the rhodopsin fold and conserved Lys are widely thought to have arisen from functional constraints during convergent evolution. To test for the existence of such a constraint, we asked whether it were possible to relocate the highly conserved Lys296 in the visual pigment bovine rhodopsin. We show here that the Lys can be moved to three other locations in the protein while maintaining the ability to form a pigment with 11-cis-retinal and activate the G protein transducin in a light-dependent manner. These results contradict the convergent hypothesis and support the homology of type I and type II rhodopsins by divergent evolution from a common ancestral protein.

An empirical test of convergent evolution in rhodopsins.

Rhodopsins are photochemically reactive membrane proteins that covalently bind retinal chromophores. Type I rhodopsins are found in both prokaryotes and eukaryotic microbes, whereas type II rhodopsins function as photoactivated G-protein coupled receptors (GPCRs) in animal vision. Both rhodopsin families share the seven transmembrane α-helix GPCR fold and a Schiff base linkage from a conserved lysine to retinal in helix G. Nevertheless, rhodopsins are widely cited as a striking example of evolutionary convergence, largely because the two families lack detectable sequence similarity and differ in many structural and mechanistic details. Convergence entails that the shared rhodopsin fold is so especially suited to photosensitive function that proteins from separate origins were selected for this architecture twice. Here we show, however, that the rhodopsin fold is not required for photosensitive activity. We engineered functional bacteriorhodopsin variants with novel folds, including radical noncircular permutations of the α-helices, circular permutations of an eight-helix construct, and retinal linkages relocated to other helices. These results contradict a key prediction of convergence and thereby provide an experimental attack on one of the most intractable problems in molecular evolution: how to establish structural homology for proteins devoid of discernible sequence similarity.

Structural basis of ligand specificity and channel activation in an insect gustatory receptor.

Gustatory receptors (GRs) are critical for insect chemosensation and are potential targets for controlling pests and disease vectors, making their structural investigation a vital step toward such applications. We present structures of Bombyx mori Gr9 (BmGr9), a fructose-gated cation channel, in agonist-free and fructose-bound states. BmGr9 forms a tetramer similar to distantly related insect odorant receptors (ORs). Upon fructose binding, BmGr9's channel gate opens through helix S7b movements. In contrast to ORs, BmGr9's ligand-binding pocket, shaped by a kinked helix S4 and a shorter extracellular S3-S4 loop, is larger and solvent accessible in both agonist-free and fructose-bound states. Also, unlike ORs, fructose binding by BmGr9 involves helix S5 and a pocket lined with aromatic and polar residues. Structure-based sequence alignments reveal distinct patterns of ligand-binding pocket residue conservation in GR subfamilies associated with different ligand classes. These data provide insight into the molecular basis of GR ligand specificity and function.

Crystal Structure of Caryolan-1-ol Synthase, a Sesquiterpene Synthase Catalyzing an Initial Anti-Markovnikov Cyclization Reaction.

In a continuing effort to understand reaction mechanisms of terpene synthases catalyzing initial anti-Markovnikov cyclization reactions, we solved the X-ray crystal structure of (+)-caryolan-1-ol synthase (CS) from Streptomyces griseus , with and without an inactive analog of the FPP substrate, 2-fluorofarnesyl diphosphate (2FFPP), bound in the active site of the enzyme. The CS-2FFPP complex was solved to 2.65 Å resolution and showed the ligand in a linear, elongated orientation, incapable of undergoing the initial cyclization event to form a bond between carbons C1 and C11. Intriguingly, the apo CS structure (2.2 Å) also had electron density in the active site, in this case density that was well fit with a curled-up tetraethylene glycol molecule presumably recruited from the crystallization medium. The density was also well fit by a molecule of farnesene suggesting that the structure may mimic an intermediate along the reaction coordinate. The curled-up conformation of tetraethylene glycol was accompanied by dramatic rotamer shifts among active-site residues. Most notably, W56 was observed to undergo a 90° rotation between the 2FFPP complex and apo-enzyme structures, suggesting that it contributes to steric interactions that help curl the tetraethylene glycol molecule in the active site, and by extension perhaps also a derivative of the FPP substrate in the normal course of the cyclization reaction. In support of this proposal, the CS W56L variant lost the ability to cyclize the FPP substrate and produced only the linear terpene products farnesol and α- and ß-farnesene.

Likelihood-based classification of cryo-EM images using FREALIGN.

We describe an implementation of maximum likelihood classification for single particle electron cryo-microscopy that is based on the FREALIGN software. Particle alignment parameters are determined by maximizing a joint likelihood that can include hierarchical priors, while classification is performed by expectation maximization of a marginal likelihood. We test the FREALIGN implementation using a simulated dataset containing computer-generated projection images of three different 70S ribosome structures, as well as a publicly available dataset of 70S ribosomes. The results show that the mixed strategy of the new FREALIGN algorithm yields performance on par with other maximum likelihood implementations, while remaining computationally efficient.

Optimal simultaneous superpositioning of multiple structures with missing data.

MOTIVATION: Superpositioning is an essential technique in structural biology that facilitates the comparison and analysis of conformational differences among topologically similar structures. Performing a superposition requires a one-to-one correspondence, or alignment, of the point sets in the different structures. However, in practice, some points are usually 'missing' from several structures, for example, when the alignment contains gaps. Current superposition methods deal with missing data simply by superpositioning a subset of points that are shared among all the structures. This practice is inefficient, as it ignores important data, and it fails to satisfy the common least-squares criterion. In the extreme, disregarding missing positions prohibits the calculation of a superposition altogether. RESULTS: Here, we present a general solution for determining an optimal superposition when some of the data are missing. We use the expectation-maximization algorithm, a classic statistical technique for dealing with incomplete data, to find both maximum-likelihood solutions and the optimal least-squares solution as a special case. AVAILABILITY AND IMPLEMENTATION: The methods presented here are implemented in THESEUS 2.0, a program for superpositioning macromolecular structures. ANSI C source code and selected compiled binaries for various computing platforms are freely available under the GNU open source license from http://www.theseus3d.org. CONTACT: dtheobald@brandeis.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Structure of an insect gustatory receptor.

Gustatory Receptors (GRs) are critical for insect chemosensation and are potential targets for controlling pests and disease vectors. However, GR structures have not been experimentally determined. We present structures of Bombyx mori Gr9 (BmGr9), a fructose-gated cation channel, in agonist-free and fructose-bound states. BmGr9 forms a tetramer similar to distantly related insect Olfactory Receptors (ORs). Upon fructose binding, BmGr9's ion channel gate opens through helix S7b movements. In contrast to ORs, BmGR9's ligand-binding pocket, shaped by a kinked helix S4 and a shorter extracellular S3-S4 loop, is larger and solvent accessible in both agonist-free and fructose-bound states. Also unlike ORs, fructose binding by BmGr9 involves helix S5 and a binding pocket lined with aromatic and polar residues. Structure-based sequence alignments reveal distinct patterns of ligand-binding pocket residue conservation in GR subfamilies associated with distinct ligand classes. These data provide insight into the molecular basis of GR ligand specificity and function.

Bayesian machine learning analysis of single-molecule fluorescence colocalization images.

Multi-wavelength single-molecule fluorescence colocalization (CoSMoS) methods allow elucidation of complex biochemical reaction mechanisms. However, analysis of CoSMoS data is intrinsically challenging because of low image signal-to-noise ratios, non-specific surface binding of the fluorescent molecules, and analysis methods that require subjective inputs to achieve accurate results. Here, we use Bayesian probabilistic programming to implement Tapqir, an unsupervised machine learning method that incorporates a holistic, physics-based causal model of CoSMoS data. This method accounts for uncertainties in image analysis due to photon and camera noise, optical non-uniformities, non-specific binding, and spot detection. Rather than merely producing a binary 'spot/no spot' classification of unspecified reliability, Tapqir objectively assigns spot classification probabilities that allow accurate downstream analysis of molecular dynamics, thermodynamics, and kinetics. We both quantitatively validate Tapqir performance against simulated CoSMoS image data with known properties and also demonstrate that it implements fully objective, automated analysis of experiment-derived data sets with a wide range of signal, noise, and non-specific binding characteristics.

A thermodynamic model of protein structure evolution explains empirical amino acid substitution matrices.

Proteins evolve under a myriad of biophysical selection pressures that collectively control the patterns of amino acid substitutions. These evolutionary pressures are sufficiently consistent over time and across protein families to produce substitution patterns, summarized in global amino acid substitution matrices such as BLOSUM, JTT, WAG, and LG, which can be used to successfully detect homologs, infer phylogenies, and reconstruct ancestral sequences. Although the factors that govern the variation of amino acid substitution rates have received much attention, the influence of thermodynamic stability constraints remains unresolved. Here we develop a simple model to calculate amino acid substitution matrices from evolutionary dynamics controlled by a fitness function that reports on the thermodynamic effects of amino acid mutations in protein structures. This hybrid biophysical and evolutionary model accounts for nucleotide transition/transversion rate bias, multi-nucleotide codon changes, the number of codons per amino acid, and thermodynamic protein stability. We find that our theoretical model accurately recapitulates the complex yet universal pattern observed in common global amino acid substitution matrices used in phylogenetics. These results suggest that selection for thermodynamically stable proteins, coupled with nucleotide mutation bias filtered by the structure of the genetic code, is the primary driver behind the global amino acid substitution patterns observed in proteins throughout the tree of life.

Fast determination of the optimal rotational matrix for macromolecular superpositions.

Finding the rotational matrix that minimizes the sum of squared deviations between two vectors is an important problem in bioinformatics and crystallography. Traditional algorithms involve the inversion or decomposition of a 3 x 3 or 4 x 4 matrix, which can be computationally expensive and numerically unstable in certain cases. Here, we present a simple and robust algorithm to rapidly determine the optimal rotation using a Newton-Raphson quaternion-based method and an adjoint matrix. Our method is at least an order of magnitude more efficient than conventional inversion/decomposition methods, and it should be particularly useful for high-throughput analyses of molecular conformations.

Accurate structural correlations from maximum likelihood superpositions.

The cores of globular proteins are densely packed, resulting in complicated networks of structural interactions. These interactions in turn give rise to dynamic structural correlations over a wide range of time scales. Accurate analysis of these complex correlations is crucial for understanding biomolecular mechanisms and for relating structure to function. Here we report a highly accurate technique for inferring the major modes of structural correlation in macromolecules using likelihood-based statistical analysis of sets of structures. This method is generally applicable to any ensemble of related molecules, including families of nuclear magnetic resonance (NMR) models, different crystal forms of a protein, and structural alignments of homologous proteins, as well as molecular dynamics trajectories. Dominant modes of structural correlation are determined using principal components analysis (PCA) of the maximum likelihood estimate of the correlation matrix. The correlations we identify are inherently independent of the statistical uncertainty and dynamic heterogeneity associated with the structural coordinates. We additionally present an easily interpretable method ("PCA plots") for displaying these positional correlations by color-coding them onto a macromolecular structure. Maximum likelihood PCA of structural superpositions, and the structural PCA plots that illustrate the results, will facilitate the accurate determination of dynamic structural correlations analyzed in diverse fields of structural biology.