Constrained Mutational Sampling of Amino Acids in HIV-1 Protease Evolution.

The evolution of HIV-1 protein sequences should be governed by a combination of factors including nucleotide mutational probabilities, the genetic code, and fitness. The impact of these factors on protein sequence evolution is interdependent, making it challenging to infer the individual contribution of each factor from phylogenetic analyses alone. We investigated the protein sequence evolution of HIV-1 by determining an experimental fitness landscape of all individual amino acid changes in protease. We compared our experimental results to the frequency of protease variants in a publicly available data set of 32,163 sequenced isolates from drug-naïve individuals. The most common amino acids in sequenced isolates supported robust experimental fitness, indicating that the experimental fitness landscape captured key features of selection acting on protease during viral infections of hosts. Amino acid changes requiring multiple mutations from the likely ancestor were slightly less likely to support robust experimental fitness than single mutations, consistent with the genetic code favoring chemically conservative amino acid changes. Amino acids that were common in sequenced isolates were predominantly accessible by single mutations from the likely protease ancestor. Multiple mutations commonly observed in isolates were accessible by mutational walks with highly fit single mutation intermediates. Our results indicate that the prevalence of multiple-base mutations in HIV-1 protease is strongly influenced by mutational sampling.

CRISPR-Cas9-mediated saturated mutagenesis screen predicts clinical drug resistance with improved accuracy.

Developing tools to accurately predict the clinical prevalence of drug-resistant mutations is a key step toward generating more effective therapeutics. Here we describe a high-throughput CRISPR-Cas9-based saturated mutagenesis approach to generate comprehensive libraries of point mutations at a defined genomic location and systematically study their effect on cell growth. As proof of concept, we mutagenized a selected region within the leukemic oncogene BCR-ABL1 Using bulk competitions with a deep-sequencing readout, we analyzed hundreds of mutations under multiple drug conditions and found that the effects of mutations on growth in the presence or absence of drug were critical for predicting clinically relevant resistant mutations, many of which were cancer adaptive in the absence of drug pressure. Using this approach, we identified all clinically isolated BCR-ABL1 mutations and achieved a prediction score that correlated highly with their clinical prevalence. The strategy described here can be broadly applied to a variety of oncogenes to predict patient mutations and evaluate resistance susceptibility in the development of new therapeutics.

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.

Quantifying and understanding the fitness effects of protein mutations: Laboratory versus nature.

The last decade has seen a growing number of experiments aimed at systematically mapping the effects of mutations in different proteins, and of attempting to correlate their biophysical and biochemical effects with organismal fitness. While insightful, systematic laboratory measurements of fitness effects present challenges and difficulties. Here, we discuss the limitations associated with such measurements, and in particular the challenge of correlating the effects of mutations at the single protein level ("protein fitness") with their effects on organismal fitness. A variety of experimental setups are used, with some measuring the direct effects on protein function and others monitoring the growth rate of a model organism carrying the protein mutants. The manners by which fitness effects are calculated and presented also vary, and the conclusions, including the derived distributions of fitness effects of mutations, vary accordingly. The comparison of the effects of mutations in the laboratory to the natural protein diversity, namely to amino acid changes that have fixed in the course of millions of years of evolution, is also debatable. The results of laboratory experiments may, therefore, be less relevant to understanding long-term inter-species variations yet insightful with regard to short-term polymorphism, for example, in the study of the effects of human SNPs.

An atomic-resolution view of neofunctionalization in the evolution of apicomplexan lactate dehydrogenases.

Malate and lactate dehydrogenases (MDH and LDH) are homologous, core metabolic enzymes that share a fold and catalytic mechanism yet possess strict specificity for their substrates. In the Apicomplexa, convergent evolution of an unusual LDH from MDH produced a difference in specificity exceeding 12 orders of magnitude. The mechanisms responsible for this extraordinary functional shift are currently unknown. Using ancestral protein resurrection, we find that specificity evolved in apicomplexan LDHs by classic neofunctionalization characterized by long-range epistasis, a promiscuous intermediate, and few gain-of-function mutations of large effect. In canonical MDHs and LDHs, a single residue in the active-site loop governs substrate specificity: Arg102 in MDHs and Gln102 in LDHs. During the evolution of the apicomplexan LDH, however, specificity switched via an insertion that shifted the position and identity of this 'specificity residue' to Trp107f. Residues far from the active site also determine specificity, as shown by the crystal structures of three ancestral proteins bracketing the key duplication event. This work provides an unprecedented atomic-resolution view of evolutionary trajectories creating a nascent enzymatic function.

Viewing protein fitness landscapes through a next-gen lens.

High-throughput sequencing has enabled many powerful approaches in biological research. Here, we review sequencing approaches to measure frequency changes within engineered mutational libraries subject to selection. These analyses can provide direct estimates of biochemical and fitness effects for all individual mutations across entire genes (and likely compact genomes in the near future) in genetically tractable systems such as microbes, viruses, and mammalian cells. The effects of mutations on experimental fitness can be assessed using sequencing to monitor time-dependent changes in mutant frequency during bulk competitions. The impact of mutations on biochemical functions can be determined using reporters or other means of separating variants based on individual activities (e.g., binding affinity for a partner molecule can be interrogated using surface display of libraries of mutant proteins and isolation of bound and unbound populations). The comprehensive investigation of mutant effects on both biochemical function and experimental fitness provide promising new avenues to investigate the connections between biochemistry, cell physiology, and evolution. We summarize recent findings from systematic mutational analyses; describe how they relate to a field rich in both theory and experimentation; and highlight how they may contribute to ongoing and future research into protein structure-function relationships, systems-level descriptions of cell physiology, and population-genetic inferences on the relative contributions of selection and drift.

Quantitative measurement of epidermal growth factor receptor-mitogen-activated protein kinase signal transduction using a nine-plex, peptide-based immunoassay.

Aberrant epidermal growth factor receptor (EGFR, ErbB1) signaling is implicated in cell transformation, motility, and invasion in a variety of cell types, and EGFR is the target of several anticancer drugs. However, the kinetics of EGFR signaling and the individual contributions of site-specific phosphorylation events remain largely unknown. A peptide-based, multiplex immunoassay approach was developed to simultaneously measure both total and phosphorylated protein in a single sample. The approach involves the proteolytic digestion of proteins prior to the isolation and quantitation of site-specific phosphorylation events within an individual protein. Quantitation of phosphorylated and total proteins, in picomolar to nanomolar concentrations, were interpolated from standard curves generated with synthetic peptides that correspond to the peptide targets used in the immunoassays. In this study, a bead-based, nine-plex immunoassay measuring total and phosphorylated protein was constructed to measure temporal, site-specific phosphorylation of key members of the EGFR pathway (ErbB1 receptor, MEK1, MEK2, ERK1, and ERK2) in A431 cells stimulated with epidermal growth factor. The effect of MEK inhibition on this pathway was determined using a known MEK kinase inhibitor, SL327. The results reported herein are the first quantitative measurements of site-specific phosphorylation events and total proteins in a single sample, at the same time representing a new paradigm for standardized protein and phosphorylation analysis using multiplexed, peptide-based, sandwich immunoassays.