Evolution favors protein mutational robustness in sufficiently large populations

BACKGROUND: An important question is whether evolution favors properties such as mutational robustness or evolvability that do not directly benefit any individual, but can influence the course of future evolution. Functionally similar proteins can differ substantially in their robustness to mutations and capacity to evolve new functions, but it has remained unclear whether any of these differences might be due to evolutionary selection for these properties. RESULTS: Here we use laboratory experiments to demonstrate that evolution favors protein mutational robustness if the evolving population is sufficiently large. We neutrally evolve cytochrome P450 proteins under identical selection pressures and mutation rates in populations of different sizes, and show that proteins from the larger and thus more polymorphic population tend towards higher mutational robustness. Proteins from the larger population also evolve greater stability, a biophysical property that is known to enhance both mutational robustness and evolvability. The excess mutational robustness and stability is well described by existing mathematical theories, and can be quantitatively related to the way that the proteins occupy their neutral network. CONCLUSIONS: Our work is the first experimental demonstration of the general tendency of evolution to favor mutational robustness and protein stability in highly polymorphic populations. We suggest that this phenomenon may contribute to the mutational robustness and evolvability of viruses and bacteria that exist in large populations

How Protein Stability and New Functions Trade Off

Numerous studies have noted that the evolution of new enzymatic specificities is accompanied by loss of the protein's thermodynamic stability (ΔΔG), thus suggesting a tradeoff between the acquisition of new enzymatic functions and stability. However, since most mutations are destabilizing (ΔΔG>0), one should ask how destabilizing mutations that confer new or altered enzymatic functions relative to all other mutations are. We applied ΔΔG computations by FoldX to analyze the effects of 548 mutations that arose from the directed evolution of 22 different enzymes. The stability effects, location, and type of function-altering mutations were compared to ΔΔG changes arising from all possible point mutations in the same enzymes. We found that mutations that modulate enzymatic functions are mostly destabilizing (average ΔΔG = +0.9 kcal/mol), and are almost as destabilizing as the “average” mutation in these enzymes (+1.3 kcal/mol). Although their stability effects are not as dramatic as in key catalytic residues, mutations that modify the substrate binding pockets, and thus mediate new enzymatic specificities, place a larger stability burden than surface mutations that underline neutral, non-adaptive evolutionary changes. How are the destabilizing effects of functional mutations balanced to enable adaptation? Our analysis also indicated that many mutations that appear in directed evolution variants with no obvious role in the new function exert stabilizing effects that may compensate for the destabilizing effects of the crucial function-altering mutations. Thus, the evolution of new enzymatic activities, both in nature and in the laboratory, is dependent on the compensatory, stabilizing effect of apparently “silent” mutations in regions of the protein that are irrelevant to its function

Synonymous Genes Explore Different Evolutionary Landscapes

The evolutionary potential of a gene is constrained not only by the amino acid sequence of its product, but by its DNA sequence as well. The topology of the genetic code is such that half of the amino acids exhibit synonymous codons that can reach different subsets of amino acids from each other through single mutation. Thus, synonymous DNA sequences should access different regions of the protein sequence space through a limited number of mutations, and this may deeply influence the evolution of natural proteins. Here, we demonstrate that this feature can be of value for manipulating protein evolvability. We designed an algorithm that, starting from an input gene, constructs a synonymous sequence that systematically includes the codons with the most different evolutionary perspectives; i.e., codons that maximize accessibility to amino acids previously unreachable from the template by point mutation. A synonymous version of a bacterial antibiotic resistance gene was computed and synthesized. When concurrently submitted to identical directed evolution protocols, both the wild type and the recoded sequence led to the isolation of specific, advantageous phenotypic variants. Simulations based on a mutation isolated only from the synthetic gene libraries were conducted to assess the impact of sub-functional selective constraints, such as codon usage, on natural adaptation. Our data demonstrate that rational design of synonymous synthetic genes stands as an affordable improvement to any directed evolution protocol. We show that using two synonymous DNA sequences improves the overall yield of the procedure by increasing the diversity of mutants generated. These results provide conclusive evidence that synonymous coding sequences do experience different areas of the corresponding protein adaptive landscape, and that a sequence's codon usage effectively constrains the evolution of the encoded protein

Structure and Age Jointly Influence Rates of Protein Evolution

What factors determine a protein's rate of evolution are actively debated. Especially unclear is the relative role of intrinsic factors of present-day proteins versus historical factors such as protein age. Here we study the interplay of structural properties and evolutionary age, as determinants of protein evolutionary rate. We use a large set of one-to-one orthologs between human and mouse proteins, with mapped PDB structures. We report that previously observed structural correlations also hold within each age group – including relationships between solvent accessibility, designabililty, and evolutionary rates. However, age also plays a crucial role: age modulates the relationship between solvent accessibility and rate. Additionally, younger proteins, despite being less designable, tend to evolve faster than older proteins. We show that previously reported relationships between age and rate cannot be explained by structural biases among age groups. Finally, we introduce a knowledge-based potential function to study the stability of proteins through large-scale computation. We find that older proteins are more stable for their native structure, and more robust to mutations, than younger ones. Our results underscore that several determinants, both intrinsic and historical, can interact to determine rates of protein evolution

High-Precision, In Vitro Validation of the Sequestration Mechanism for Generating Ultrasensitive Dose-Response Curves in Regulatory Networks

Our ability to recreate complex biochemical mechanisms in designed, artificial systems provides a stringent test of our understanding of these mechanisms and opens the door to their exploitation in artificial biotechnologies. Motivated by this philosophy, here we have recapitulated in vitro the “target sequestration” mechanism used by nature to improve the sensitivity (the steepness of the input/output curve) of many regulatory cascades. Specifically, we have employed molecular beacons, a commonly employed optical DNA sensor, to recreate the sequestration mechanism and performed an exhaustive, quantitative study of its key determinants (e.g., the relative concentrations and affinities of probe and depletant). We show that, using sequestration, we can narrow the pseudo-linear range of a traditional molecular beacon from 81-fold (i.e., the transition from 10% to 90% target occupancy spans an 81-fold change in target concentration) to just 1.5-fold. This narrowing of the dynamic range improves the sensitivity of molecular beacons to that equivalent of an oligomeric, allosteric receptor with a Hill coefficient greater than 9. Following this we have adapted the sequestration mechanism to steepen the binding-site occupancy curve of a common transcription factor by an order of magnitude over the sensitivity observed in the absence of sequestration. Given the success with which the sequestration mechanism has been employed by nature, we believe that this strategy could dramatically improve the performance of synthetic biological systems and artificial biosensors

Search for a Technicolor omega_T Particle in Events with a Photon and a b-quark Jet at CDF

If the Technicolor omega_T particle exists, a likely decay mode is omega_T -> gamma pi_T, followed by pi_T -> bb-bar, yielding the signature gamma bb-bar. We have searched 85 pb^-1 of data collected by the CDF experiment at the Fermilab Tevatron for events with a photon and two jets, where one of the jets must contain a secondary vertex implying the presence of a b quark. We find no excess of events above standard model expectations. We express the result of an exclusion region in the M_omega_T - M_pi_T mass plane.Comment: 14 pages, 2 figures. Available from the CDF server (PS with figs): http://www-cdf.fnal.gov/physics/pub98/cdf4674_omega_t_prl_4.ps FERMILAB-PUB-98/321-

Measurement of the B0 anti-B0 oscillation frequency using l- D*+ pairs and lepton flavor tags

The oscillation frequency Delta-md of B0 anti-B0 mixing is measured using the partially reconstructed semileptonic decay anti-B0 -> l- nubar D*+ X. The data sample was collected with the CDF detector at the Fermilab Tevatron collider during 1992 - 1995 by triggering on the existence of two lepton candidates in an event, and corresponds to about 110 pb-1 of pbar p collisions at sqrt(s) = 1.8 TeV. We estimate the proper decay time of the anti-B0 meson from the measured decay length and reconstructed momentum of the l- D*+ system. The charge of the lepton in the final state identifies the flavor of the anti-B0 meson at its decay. The second lepton in the event is used to infer the flavor of the anti-B0 meson at production. We measure the oscillation frequency to be Delta-md = 0.516 +/- 0.099 +0.029 -0.035 ps-1, where the first uncertainty is statistical and the second is systematic.Comment: 30 pages, 7 figures. Submitted to Physical Review

Search for New Particles Decaying to top-antitop in proton-antiproton collisions at squareroot(s)=1.8 TeV

We use 106 \ipb of data collected with the Collider Detector at Fermilab to search for narrow-width, vector particles decaying to a top and an anti-top quark. Model independent upper limits on the cross section for narrow, vector resonances decaying to \ttbar are presented. At the 95% confidence level, we exclude the existence of a leptophobic \zpr boson in a model of topcolor-assisted technicolor with mass M_{\zpr} $<$ 480 \gev for natural width $\Gamma$ = 0.012 M_{\zpr}, and M_{\zpr} $<$ 780 \gev for $\Gamma$ = 0.04 M_{\zpr}.Comment: The CDF Collaboration, submitted to PRL 25-Feb-200

A Study of B0 -> J/psi K(*)0 pi+ pi- Decays with the Collider Detector at Fermilab

We report a study of the decays B0 -> J/psi K(*)0 pi+ pi-, which involve the creation of a u u-bar or d d-bar quark pair in addition to a b-bar -> c-bar(c s-bar) decay. The data sample consists of 110 1/pb of p p-bar collisions at sqrt{s} = 1.8 TeV collected by the CDF detector at the Fermilab Tevatron collider during 1992-1995. We measure the branching ratios to be BR(B0 -> J/psi K*0 pi+ pi-) = (8.0 +- 2.2 +- 1.5) * 10^{-4} and BR(B0 -> J/psi K0 pi+ pi-) = (1.1 +- 0.4 +- 0.2) * 10^{-3}. Contributions to these decays are seen from psi(2S) K(*)0, J/psi K0 rho0, J/psi K*+ pi-, and J/psi K1(1270)