Reducing the flexibility of retinal restores a wild-type-like photocycle in bacteriorhodopsin mutants defective in protein–retinal coupling

The thermal re-isomerization of retinal from the 13-cis to the all-trans state is a key step in the final stages of the photocycle of the light-driven proton pump, bacteriorhodopsin. This step is greatly slowed upon replacement of Leu-93, a residue in van der Waals contact with retinal. The most likely role of this key interaction is that it restricts the flexibility of retinal. To test this hypothesis, we have exchanged native retinal in Leu-93 mutants with bridged retinal analogs that render retinal less flexible by restricting free rotation around either the C10—C11 (9,11-bridged retinal) or C12—C13 (11,13-bridged retinal) single bonds. The effect of the analogs on the photocycle was then determined spectroscopically by taking advantage of the previous finding that the decay of the O intermediate in the Leu-93 mutants provides a convenient marker for retinal re-isomerization. Time-resolved spectroscopic studies showed that both retinal analogs resulted in a dramatic acceleration of the photocycling time by increasing the rate of decay of the O intermediate. In particular, exchange of native retinal in the Leu-93 → Ala mutant with the 9,11-bridged retinal resulted in an acceleration of the decay of the O intermediate to a rate similar to that seen in wild-type bacteriorhodopsin. We conclude that the protein-induced restriction of conformational flexibility in retinal is a key structural requirement for efficient protein–retinal coupling in the bacteriorhodopsin photocycle.

Cell locomotion and focal adhesions are regulated by substrate flexibility

Responses of cells to mechanical properties of the adhesion substrate were examined by culturing normal rat kidney epithelial and 3T3 fibroblastic cells on a collagen-coated polyacrylamide substrate that allows the flexibility to be varied while maintaining a constant chemical environment. Compared with cells on rigid substrates, those on flexible substrates showed reduced spreading and increased rates of motility or lamellipodial activity. Microinjection of fluorescent vinculin indicated that focal adhesions on flexible substrates were irregularly shaped and highly dynamic whereas those on firm substrates had a normal morphology and were much more stable. Cells on flexible substrates also contained a reduced amount of phosphotyrosine at adhesion sites. Treatment of these cells with phenylarsine oxide, a tyrosine phosphatase inhibitor, induced the formation of normal, stable focal adhesions similar to those on firm substrates. Conversely, treatment of cells on firm substrates with myosin inhibitors 2,3-butanedione monoxime or KT5926 caused the reduction of both vinculin and phosphotyrosine at adhesion sites. These results demonstrate the ability of cells to survey the mechanical properties of their surrounding environment and suggest the possible involvement of both protein tyrosine phosphorylation and myosin-generated cortical forces in this process. Such response to physical parameters likely represents an important mechanism of cellular interaction with the surrounding environment within a complex organism.

Prion protein: Evolution caught en route

The prion protein displays a unique structural ambiguity in that it can adopt multiple stable conformations under physiological conditions. In our view, this puzzling feature resulted from a sudden environmental change in evolution when the prion, previously an integral membrane protein, got expelled into the extracellular space. Analysis of known vertebrate prions unveils a primordial transmembrane protein encrypted in their sequence, underlying this relocalization hypothesis. Apparently, the time elapsed since this event was insufficient to create a “minimally frustrated” sequence in the new milieu, probably due to the functional constraints set by the importance of the very flexibility that was created in the relocalization. This scenario may explain why, in a structural sense, the prion protein is still en route toward becoming a foldable globular protein.

Identification of interdomain sequences promoting the intronless evolution of a bacterial protein family.

In the evolution of eukaryotic genes, introns are believed to have played a major role in increasing the probability of favorable duplication events, chance recombinations, and exon shuffling resulting in functional hybrid proteins. As a rule, prokaryotic genes lack introns, and the examples of prokaryotic introns described do not seem to have contributed to gene evolution by exon shuffling. Still, certain protein families in modern bacteria evolve rapidly by recombination of genes, duplication of functional domains, and as shown for protein PAB of the anaerobic bacterial species Peptostreptococcus magnus, by the shuffling of an albumin-binding protein module from group C and G streptococci. Characterization of a protein PAB-related gene in a P. magnus strain with less albumin-binding activity revealed that the shuffled module was missing. Based on this fact and observations made when comparing gene sequences of this family of bacterial surface proteins interacting with albumin and/or immunoglobulin, a model is presented that can explain how this rapid intronless evolution takes place. A new kind of genetic element is introduced: the recer sequence promoting interdomain, in frame recombination and acting as a structure-less flexibility-promoting spacer in the corresponding protein. The data presented also suggest that antibiotics could represent the selective pressure behind the shuffling of protein modules in P. magnus, a member of the indigenous bacterial flora.

Internal molecular motions of bacteriorhodopsin: hydration-induced flexibility studied by quasielastic incoherent neutron scattering using oriented purple membranes.

Quasielastic incoherent neutron scattering from hydrogen atoms, which are distributed nearly homogeneously in biological molecules, allows the investigation of diffusive motions occurring on the pico- to nanosecond time scale. A quasielastic incoherent neutron scattering study was performed on the integral membrane protein bacteriorhodopsin (BR), which is a light-driven proton pump in Halobacterium salinarium. BR is embedded in lipids, forming patches in the cell membrane of the organism, which are the so called purple membranes (PMs). Measurements were carried out at room temperature on oriented PM-stacks hydrated at two different levels (low hydration, h = 0.03 g of D2O per g of PM; high hydration, h = 0.28 g of D2O per g of PM) using time-of-flight spectrometers. From the measured spectra, different diffusive components were identified and analyzed with respect to the influence of hydration. This study supports the idea that a decrease in hydration results in an appreciable decrease in internal molecular flexibility of the protein structure. Because it is known from studies on the function of BR that the pump activity is reduced if the hydration level of the protein is insufficient, we conclude that the observed diffusive motions are essential for the function of this protein. A detailed analysis and classification of the different kinds of diffusive motions, predominantly occurring in PMs under physiological conditions, is presented.

Repacking protein cores with backbone freedom: structure prediction for coiled coils.

Progress in homology modeling and protein design has generated considerable interest in methods for predicting side-chain packing in the hydrophobic cores of proteins. Present techniques are not practically useful, however, because they are unable to model protein main-chain flexibility. Parameterization of backbone motions may represent a general and efficient method to incorporate backbone relaxation into such fixed main-chain models. To test this notion, we introduce a method for treating explicitly the backbone motions of alpha-helical bundles based on an algebraic parameterization proposed by Francis Crick in 1953 [Crick, F. H. C. (1953) Acta Crystallogr. 6, 685-689]. Given only the core amino acid sequence, a simple calculation can rapidly reproduce the crystallographic main-chain and core side-chain structures of three coiled coils (one dimer, one trimer, and one tetramer) to within 0.6-A root-mean-square deviations. The speed of the predictive method [approximately 3 min per rotamer choice on a Silicon Graphics (Mountain View, CA) 4D/35 computer] permits it to be used as a design tool.

Mechanistic Importance of Redox Potentials and Conformational Flexibility in Electron-Transferring Flavoproteins

The mitochondrial matrix flavoproteins electron transfer flavoprotein (ETF) and electron transfer flavoprotein-ubiquinone oxidoreductase (ETF-QO) are responsible for linking fatty acid β-oxidation with the main mitochondrial respiratory chain. Electrons derived from flavoprotein dehydrogenases are transferred sequentially through ETF and ETF-QO to ubiquinone and then into the respiratory chain via complex III. In this study, the effects of changes in ETF-QO redox potentials on its activity and the conformational flexibility of ETF were investigated. ETF-QO contains one [4Fe-4S]2+,1+ and one flavin adenine dinucleotide (FAD). In the porcine protein, threonine 367 is hydrogen bonded to N1 and O2 of the flavin ring of the FAD. The analogous site in Rhodobacter sphaeroides ETF-QO is asparagine 338. Mutations N338T and N338A were introduced into the R. sphaeroides protein by site-directed mutagenesis to determine the impact of hydrogen bonding at this site on redox potentials and activity. FAD redox potentials were measured by potentiometric titration probed by electron paramagnetic resonance (EPR) spectroscopy. The N338T and N338A mutations lowered the midpoint potentials, which resulted in a decrease in the quinone reductase activity and negligible impact on disproportionation of ETF1e- catalyzed by ETF-QO. These observations indicate that the FAD is involved in electron transfer to ubiquinone, but not in electron transfer from ETF to ETF-QO. Therefore it is proposed that the iron-sulfur cluster is the immediate acceptor from ETF. It has been proposed that the αII domain of ETF is mobile, allowing promiscuous interactions with structurally different partners. Double electron-electron resonance (DEER) was used to measure the distance between spin labels at various sites and an enzymatically reduced FAD cofactor in Paracoccus denitrificans ETF. Two or three interspin distance distributions were observed for spin-labels in the αI (A43C) and βIII (A111C) domains, but only one is observed for a label in the βII (A210C) domain. This suggests that the αII domain adopts several stable conformations which may correspond to a closed/inactive conformation and an open/active conformation. An additional mutation, E162A, was introduced to increase the mobility of the αII domain. The E162A mutation doubled the activity compared to wild-type and caused the distance distributions to become wider. The DEER method has the potential to characterize conformational changes in ETF that occur when it interacts with various redox partners.

Prediction of protein B-factor profiles

The polypeptide backbones and side chains of proteins are constantly moving due to thermal motion and the kinetic energy of the atoms. The B-factors of protein crystal structures reflect the fluctuation of atoms about their average positions and provide important information about protein dynamics. Computational approaches to predict thermal motion are useful for analyzing the dynamic properties of proteins with unknown structures. In this article, we utilize a novel support vector regression (SVR) approach to predict the B-factor distribution (B-factor profile) of a protein from its sequence. We explore schemes for encoding sequences and various settings for the parameters used in SVR. Based on a large dataset of high-resolution proteins, our method predicts the B-factor distribution with a Pearson correlation coefficient (CC) of 0.53. In addition, our method predicts the B-factor profile with a CC of at least 0.56 for more than half of the proteins. Our method also performs well for classifying residues (rigid vs. flexible). For almost all predicted B-factor thresholds, prediction accuracies (percent of correctly predicted residues) are greater than 70%. These results exceed the best results of other sequence-based prediction methods. (C) 2005 Wiley-Liss, Inc.

Prediction of protein continuum secondary structure with probabilistic models based on NMR solved structures

Background: The structure of proteins may change as a result of the inherent flexibility of some protein regions. We develop and explore probabilistic machine learning methods for predicting a continuum secondary structure, i.e. assigning probabilities to the conformational states of a residue. We train our methods using data derived from high-quality NMR models. Results: Several probabilistic models not only successfully estimate the continuum secondary structure, but also provide a categorical output on par with models directly trained on categorical data. Importantly, models trained on the continuum secondary structure are also better than their categorical counterparts at identifying the conformational state for structurally ambivalent residues. Conclusion: Cascaded probabilistic neural networks trained on the continuum secondary structure exhibit better accuracy in structurally ambivalent regions of proteins, while sustaining an overall classification accuracy on par with standard, categorical prediction methods.

Structure and dynamics of ASC2, a pyrin domain-only protein that regulates inflammatory signaling

Pyrin domain (PYD)-containing proteins are key components of pathways that regulate inflammation, apoptosis, and cytokine processing. Their importance is further evidenced by the consequences of mutations in these proteins that give rise to autoimmune and hyperinflammatory syndromes. PYDs, like other members of the death domain ( DD) superfamily, are postulated to mediate homotypic interactions that assemble and regulate the activity of signaling complexes. However, PYDs are presently the least well characterized of all four DD subfamilies. Here we report the three-dimensional structure and dynamic properties of ASC2, a PYD-only protein that functions as a modulator of multidomain PYD-containing proteins involved in NF-KB and caspase-1 activation. ASC2 adopts a six-helix bundle structure with a prominent loop, comprising 13 amino acid residues, between helices two and three. This loop represents a divergent feature of PYDs from other domains with the DD fold. Detailed analysis of backbone N-15 NMR relaxation data using both the Lipari-Szabo model-free and reduced spectral density function formalisms revealed no evidence of contiguous stretches of polypeptide chain with dramatically increased internal motion, except at the extreme N and C termini. Some mobility in the fast, picosecond to nanosecond timescale, was seen in helix 3 and the preceding alpha 2-alpha 3 loop, in stark contrast to the complete disorder seen in the corresponding region of the NALP1 PYD. Our results suggest that extensive conformational flexibility in helix 3 and the alpha 2-alpha 3 loop is not a general feature of pyrin domains. Further, a transition from complete disorder to order of the alpha 2-alpha 3 loop upon binding, as suggested for NALP1, is unlikely to be a common attribute of pyrin domain interactions.

The MK2/3 cascade regulates AMPAR trafficking and cognitive flexibility

The interplay between long-term potentiation and long-term depression (LTD) is thought to be involved in learning and memory formation. One form of LTD expressed in the hippocampus is initiated by the activation of the group 1 metabotropic glutamate receptors (mGluRs). Importantly, mGluRs have been shown to be critical for acquisition of new memories and for reversal learning, processes that are thought to be crucial for cognitive flexibility. Here we provide evidence that MAPK-activated protein kinases 2 and 3 (MK2/3) regulate neuronal spine morphology, synaptic transmission and plasticity. Furthermore, mGluR-LTD is impaired in the hippocampus of MK2/3 double knockout (DKO) mice, an observation that is mirrored by deficits in endocytosis of GluA1 subunits. Consistent with compromised mGluR-LTD, MK2/3 DKO mice have distinctive deficits in hippocampal-dependent spatial reversal learning. These novel findings demonstrate that the MK2/3 cascade plays a strategic role in controlling synaptic plasticity and cognition. © 2014 Macmillan Publishers Limited. All rights reserved.

Structural flexibility and oxygen diffusion pathways in monomeric fluorescent proteins

Fluorescent proteins are valuable tools as biochemical markers for studying cellular processes. Red fluorescent proteins (RFPs) are highly desirable for in vivo applications because they absorb and emit light in the red region of the spectrum where cellular autofluorescence is low. The naturally occurring fluorescent proteins with emission peaks in this region of the spectrum occur in dimeric or tetrameric forms. The development of mutant monomeric variants of RFPs has resulted in several novel FPs known as mFruits. Though oxygen is required for maturation of the chromophore, it is known that photobleaching of FPs is oxygen sensitive, and oxygen-free conditions result in improved photostabilities. Therefore, understanding oxygen diffusion pathways in FPs is important for both photostabilites and maturation of the chromophores. We used molecular dynamics calculations to investigate the protein barrel fluctuations in mCherry, which is one of the most useful monomeric mFruit variants, and its GFP homolog citrine. We employed implicit ligand sampling and locally enhanced sampling to determine oxygen pathways from the bulk solvent into the mCherry chromophore in the interior of the protein. The pathway contains several oxygen hosting pockets, which were identified by the amino acid residues that form the pocket. We calculated the free-energy of an oxygen molecule at points along the path. We also investigated an RFP variant known to be significantly less photostable than mCherry and find much easier oxygen access in this variant. We showed that oxygen pathways can be blocked or altered, and barrel fluctuations can be reduced by strategic amino acid substitutions. The results provide a better understanding of the mechanism of molecular oxygen access into the fully folded mCherry protein barrel and provide insight into the photobleaching process in these proteins.

Protein and Drug Design Algorithms Using Improved Biophysical Modeling

This thesis focuses on the development of algorithms that will allow protein design calculations to incorporate more realistic modeling assumptions. Protein design algorithms search large sequence spaces for protein sequences that are biologically and medically useful. Better modeling could improve the chance of success in designs and expand the range of problems to which these algorithms are applied. I have developed algorithms to improve modeling of backbone flexibility (DEEPer) and of more extensive continuous flexibility in general (EPIC and LUTE). I’ve also developed algorithms to perform multistate designs, which account for effects like specificity, with provable guarantees of accuracy (COMETS), and to accommodate a wider range of energy functions in design (EPIC and LUTE).

Detection of African Swine Fever Virus Antibodies in Serum and Oral Fluid Specimens Using a Recombinant Protein 30 (p30) Dual Matrix Indirect ELISA

In the absence of effective vaccine(s), control of African swine fever caused by African swine fever virus (ASFV) must be based on early, efficient, cost-effective detection and strict control and elimination strategies. For this purpose, we developed an indirect ELISA capable of detecting ASFV antibodies in either serum or oral fluid specimens. The recombinant protein used in the ELISA was selected by comparing the early serum antibody response of ASFV-infected pigs (NHV-p68 isolate) to three major recombinant polypeptides (p30, p54, p72) using a multiplex fluorescent microbead-based immunoassay (FMIA). Non-hazardous (non-infectious) antibody-positive serum for use as plate positive controls and for the calculation of sample-to-positive (S:P) ratios was produced by inoculating pigs with a replicon particle (RP) vaccine expressing the ASFV p30 gene. The optimized ELISA detected anti-p30 antibodies in serum and/or oral fluid samples from pigs inoculated with ASFV under experimental conditions beginning 8 to 12 days post inoculation. Tests on serum (n = 200) and oral fluid (n = 200) field samples from an ASFV-free population demonstrated that the assay was highly diagnostically specific. The convenience and diagnostic utility of oral fluid sampling combined with the flexibility to test either serum or oral fluid on the same platform suggests that this assay will be highly useful under the conditions for which OIE recommends ASFV antibody surveillance, i.e., in ASFV-endemic areas and for the detection of infections with ASFV isolates of low virulence.