Paul W. Fenimore

Slaving: Solvent fluctuations dominate protein dynamics and functions

Protein motions are essential for function. Comparing protein processes with the dielectric fluctuations of the surrounding solvent shows that they fall into two classes: nonslaved and slaved. Nonslaved processes are independent of the solvent motions; their rates are determined by the protein conformation and vibrational dynamics. Slaved processes are tightly coupled to the solvent; their rates have approximately the same temperature dependence as the rate of the solvent fluctuations, but they are smaller. Because the temperature dependence is determined by the activation enthalpy, we propose that the solvent is responsible for the activation enthalpy, whereas the protein and the hydration shell control the activation entropy through the energy landscape. Bond formation is the prototype of nonslaved processes; opening and closing of channels are quintessential slaved motions. The prevalence of slaved motions highlights the importance of the environment in cells and membranes for the function of proteins.

A unified model of protein dynamics

Protein functions require conformational motions. We show here that the dominant conformational motions are slaved by the hydration shell and the bulk solvent. The protein contributes the structure necessary for function. We formulate a model that is based on experiments, insights from the physics of glass-forming liquids, and the concepts of a hierarchically organized energy landscape. To explore the effect of external fluctuations on protein dynamics, we measure the fluctuations in the bulk solvent and the hydration shell with broadband dielectric spectroscopy and compare them with internal fluctuations measured with the Mössbauer effect and neutron scattering. The result is clear. Large-scale protein motions are slaved to the fluctuations in the bulk solvent. They are controlled by the solvent viscosity, and are absent in a solid environment. Internal protein motions are slaved to the beta fluctuations of the hydration shell, are controlled by hydration, and are absent in a dehydrated protein. The model quantitatively predicts the rapid increase of the mean-square displacement above ≈200 K, shows that the external beta fluctuations determine the temperature- and time-dependence of the passage of carbon monoxide through myoglobin, and explains the nonexponential time dependence of the protein relaxation after photodissociation.

Protein folding is slaved to solvent motions

Proteins, the workhorses of living systems, are constructed from chains of amino acids, which are synthesized in the cell based on the instructions of the genetic code and then folded into working proteins. The time for folding varies from microseconds to hours. What controls the folding rate is hotly debated. We postulate here that folding has the same temperature dependence as the α-fluctuations in the bulk solvent but is much slower. We call this behavior slaving. Slaving has been observed in folded proteins: Large-scale protein motions follow the solvent fluctuations with rate coefficient kα but can be slower by a large factor. Slowing occurs because large-scale motions proceed in many small steps, each determined by kα. If conformational motions of folded proteins are slaved, so a fortiori must be the motions during folding. The unfolded protein makes a Brownian walk in the conformational space to the folded structure, with each step controlled by kα. Because the number of conformational substates in the unfolded protein is extremely large, the folding rate coefficient, kf, is much smaller than kα. The slaving model implies that the activation enthalpy of folding is dominated by the solvent, whereas the number of steps nf = kα/kf is controlled by the number of accessible substates in the unfolded protein and the solvent. Proteins, however, undergo not only α- but also β-fluctuations. These additional fluctuations are local protein motions that are essentially independent of the bulk solvent fluctuations and may be relevant at late stages of folding.

SARS outbreaks in Ontario, Hong Kong and Singapore: the role of diagnosis and isolation as a control mechanism

Abstract In this article we use global and regional data from the SARS epidemic in conjunction with a model of susceptible, exposed, infective, diagnosed, and recovered classes of people (“SEIJR”) to extract average properties and rate constants for those populations. The model is fitted to data from the Ontario (Toronto) in Canada, Hong Kong in China and Singapore outbreaks and predictions are made based on various assumptions and observations, including the current effect of isolating individuals diagnosed with SARS. The epidemic dynamics for Hong Kong and Singapore appear to be different from the dynamics in Toronto, Ontario. Toronto shows a very rapid increase in the number of cases between March 31st and April 6th, followed by a significant slowing in the number of new cases. We explain this as the result of an increase in the diagnostic rate and in the effectiveness of patient isolation after March 26th. Our best estimates are consistent with SARS eventually being contained in Toronto, although the time of containment is sensitive to the parameters in our model. It is shown that despite the empirically modeled heterogeneity in transmission, SARS’ average reproductive number is 1.2, a value quite similar to that computed for some strains of influenza (J. Math. Biol. 27 (1989) 233). Although it would not be surprising to see levels of SARS infection higher than 10% in some regions of the world (if unchecked), lack of data and the observed heterogeneity and sensitivity of parameters prevent us from predicting the long-term impact of SARS. The possibility that 10 or more percent of the world population at risk could eventually be infected with the virus in conjunction with a mortality rate of 3–7% or more, and indications of significant improvement in Toronto support the stringent measures that have been taken to isolate diagnosed cases.

Myoglobin: the hydrogen atom of biology and a paradigm of complexity.

In a tour de force in this issue of PNAS, Bourgeois et al. (1) have used 2.2-ns x-ray pulses to observe the motion of carbon monoxide (CO) through myoglobin (Mb) and the relaxation of the protein from 3.2 ns to 3 ms after photodissociation. This work follows the pioneering experiments of Moffat and collaborators (2). It demonstrates how far advances in x-ray sources and computers have moved the field of protein structure determination since the path-breaking work of Kendrew et al . and Perutz (3, 4). The recent breakthrough shows how careful studies of proteins, in particular of Mb, impact many different fields. Mb is a monomeric protein that gives muscle its red color. Thirty years ago the textbook function of Mb, storage of dioxygen at the heme iron, was considered to be simple, fully understood, and consequently boring. Mb was essentially written off as a topic of serious research. Since then, the situation has changed: Mb is no longer fully understood. It plays roles other than O2 storage, serves as a prototype for complex systems, and yields insight into the chemistry and physics of soft matter and of chemical reactions. Mb consists of 153 amino acids that fold into a structure that is ≈3 nm in diameter, as depicted in Fig. 1 a . Fig. 1 b gives a schematic cross section through Mb that shows the active center: a heme group with a central iron atom. Surrounding the heme group are five cavities, the heme cavity and four cavities denoted by Xe1 to Xe4 (5). The amino acids lining the xenon cavities are much more conserved than other amino acids in mammalian Mb; they are thus likely to be important for function. A major part of Mb is taken up by amino acids that do not appear to …

Model Parameters and Outbreak Control for SARS

Tool for estimating basic reproductive number for the SARS outbreak suggests need for multiple methods of control.

Hydration, slaving and protein function

Protein dynamics is crucial for protein function. Proteins in living systems are not isolated, but operate in networks and in a carefully regulated environment. Understanding the external control of protein dynamics is consequently important. Hydration and solvent viscosity are among the salient properties of the environment. Dehydrated proteins and proteins in a rigid environment do not function properly. It is consequently important to understand the effect of hydration and solvent viscosity in detail. We discuss experiments that separate the two effects. These experiments have predominantly been performed with wild-type horse and sperm whale myoglobin, using the binding of carbon monoxide over a broad range of temperatures as a tool. The experiments demonstrate that data taken only in the physiological temperature range are not sufficient to understand the effect of hydration and solvent on protein relaxation and function. While the actual data come from myoglobin, it is expected that the results apply to most or all globular proteins.

Protein dynamics and function: Insights from the energy landscape and solvent slaving

Protein motions are complex and a good way to describe them is in terms of a very high‐dimensional conformation space. We give here a simple explanation of the conformation space and the energy landscape, the conformational motions and protein reactions, based on an analogy to a traffic problem. The analogy provides insight into the slaving of protein processes to bulk solvent fluctuations, in both the native and unfolded states.

Filovirus RefSeq Entries: Evaluation and Selection of Filovirus Type Variants, Type Sequences, and Names

Sequence determination of complete or coding-complete genomes of viruses is becoming common practice for supporting the work of epidemiologists, ecologists, virologists, and taxonomists. Sequencing duration and costs are rapidly decreasing, sequencing hardware is under modification for use by non-experts, and software is constantly being improved to simplify sequence data management and analysis. Thus, analysis of virus disease outbreaks on the molecular level is now feasible, including characterization of the evolution of individual virus populations in single patients over time. The increasing accumulation of sequencing data creates a management problem for the curators of commonly used sequence databases and an entry retrieval problem for end users. Therefore, utilizing the data to their fullest potential will require setting nomenclature and annotation standards for virus isolates and associated genomic sequences. The National Center for Biotechnology Information’s (NCBI’s) RefSeq is a non-redundant, curated database for reference (or type) nucleotide sequence records that supplies source data to numerous other databases. Building on recently proposed templates for filovirus variant naming [ ()////-], we report consensus decisions from a majority of past and currently active filovirus experts on the eight filovirus type variants and isolates to be represented in RefSeq, their final designations, and their associated sequences.

Designing and Testing Broadly-Protective Filoviral Vaccines Optimized for Cytotoxic T-Lymphocyte Epitope Coverage

We report the rational design and in vivo testing of mosaic proteins for a polyvalent pan-filoviral vaccine using a computational strategy designed for the Human Immunodeficiency Virus type 1 (HIV-1) but also appropriate for Hepatitis C virus (HCV) and potentially other diverse viruses. Mosaics are sets of artificial recombinant proteins that are based on natural proteins. The recombinants are computationally selected using a genetic algorithm to optimize the coverage of potential cytotoxic T lymphocyte (CTL) epitopes. Because evolutionary history differs markedly between HIV-1 and filoviruses, we devised an adapted computational technique that is effective for sparsely sampled taxa; our first significant result is that the mosaic technique is effective in creating high-quality mosaic filovirus proteins. The resulting coverage of potential epitopes across filovirus species is superior to coverage by any natural variants, including current vaccine strains with demonstrated cross-reactivity. The mosaic cocktails are also robust: mosaics substantially outperformed natural strains when computationally tested against poorly sampled species and more variable genes. Furthermore, in a computational comparison of cross-reactive potential a design constructed prior to the Bundibugyo outbreak performed nearly as well against all species as an updated design that included Bundibugyo. These points suggest that the mosaic designs would be more resilient than natural-variant vaccines against future Ebola outbreaks dominated by novel viral variants. We demonstrate in vivo immunogenicity and protection against a heterologous challenge in a mouse model. This design work delineates the likely requirements and limitations on broadly-protective filoviral CTL vaccines.