Enzymatic metabolons dramatically enhance metabolic fluxes of low-efficiency biochemical reactions

In this work, we investigate how spatial proximity of enzymes belonging to the same pathway (metabolon) affects metabolic flux. Using off-lattice Langevin dynamics simulations in tandem with a stochastic reaction-diffusion protocol and a semi-analytical reaction-diffusion model, we systematically explored how strength of protein-protein interactions, catalytic efficiency, and protein-ligand interactions affect metabolic flux through the metabolon. Formation of a metabolon leads to a greater speedup for longer pathways and especially for reaction-limited enzymes, whereas, for fully optimized diffusion-limited enzymes, the effect is negligible. Notably, specific cluster architectures are not a prerequisite for enhancing reaction flux. Simulations uncover the crucial role of optimal nonspecific protein-ligand interactions in enhancing catalytic efficiency of a metabolon. Our theory implies, and bioinformatics analysis confirms, that longer catalytic pathways are enriched in less optimal enzymes, whereas most diffusion-limited enzymes populate shorter pathways. Our findings point toward a plausible evolutionary strategy where enzymes compensate for less-than-optimal efficiency by increasing their local concentration in the clustered state.

Mapping the Space of Photoswitchable Ligands and Photodruggable Proteins with Computational Modeling

Light-activated drugs are a promising way to localize biological activity and minimize side effects. However, their development is complicated by the numerous photophysical and biological properties that must be simultaneously optimized. To accelerate the design of photoactive drugs, we describe a procedure that combines ligand-protein docking with chemical property prediction based on machine learning (ML). We apply this procedure to 58 proteins and 9000 photo-drug candidates based on azobenzene cis-trans isomerism. We find that most proteins display a preference for trans isomers over cis and that the binding affinities of nominally active/inactive pairs are in fact highly correlated. These findings have significant value for photopharmacology research, and reinforce the need for virtual screening to identify compounds with rare desirable properties. Further, we combine our procedure with quantum chemical validation to identify promising candidates for the photoactive inhibition of PARP1, an enzyme that is over-expressed in cancer cells. The top compounds are predicted to have long-lived active forms, differential bioactivity, and absorption in the near-infrared therapeutic window.

Phosphorylation sites are evolutionary checkpoints against liquid-solid transition in protein condensates

Assemblies of multivalent RNA-binding protein fused in sarcoma (FUS) can exist in the functional liquid-like state as well as less dynamic and potentially toxic amyloid- and hydrogel-like states. How could then cells form liquid-like condensates while avoiding their transformation to amyloids? Here, we show how posttranslational phosphorylation can provide a "handle" that prevents liquid-solid transition of intracellular condensates containing FUS. Using residue-specific coarse-grained simulations, for 85 different mammalian FUS sequences, we show how the number of phosphorylation sites and their spatial arrangement affect intracluster dynamics preventing conversion to amyloids. All atom simulations further confirm that phosphorylation can effectively reduce the β-sheet propensity in amyloid-prone fragments of FUS. A detailed evolutionary analysis shows that mammalian FUS PLDs are enriched in amyloid-prone stretches compared to control neutrally evolved sequences, suggesting that mammalian FUS proteins evolved to self-assemble. However, in stark contrast to proteins that do not phase-separate for their function, mammalian sequences have phosphosites in close proximity to these amyloid-prone regions. These results suggest that evolution uses amyloid-prone sequences in prion-like domains to enhance phase separation of condensate proteins while enriching phosphorylation sites in close proximity to safeguard against liquid-solid transitions.

Thermal Half-Lives of Azobenzene Derivatives: Virtual Screening Based on Intersystem Crossing Using a Machine Learning Potential

Molecular photoswitches are the foundation of light-activated drugs. A key photoswitch is azobenzene, which exhibits trans-cis isomerism in response to light. The thermal half-life of the cis isomer is of crucial importance, since it controls the duration of the light-induced biological effect. Here we introduce a computational tool for predicting the thermal half-lives of azobenzene derivatives. Our automated approach uses a fast and accurate machine learning potential trained on quantum chemistry data. Building on well-established earlier evidence, we argue that thermal isomerization proceeds through rotation mediated by intersystem crossing, and incorporate this mechanism into our automated workflow. We use our approach to predict the thermal half-lives of 19,000 azobenzene derivatives. We explore trends and trade-offs between barriers and absorption wavelengths, and open-source our data and software to accelerate research in photopharmacology.

Excited state non-adiabatic dynamics of large photoswitchable molecules using a chemically transferable machine learning potential

Light-induced chemical processes are ubiquitous in nature and have widespread technological applications. For example, photoisomerization can allow a drug with a photo-switchable scaffold such as azobenzene to be activated with light. In principle, photoswitches with desired photophysical properties like high isomerization quantum yields can be identified through virtual screening with reactive simulations. In practice, these simulations are rarely used for screening, since they require hundreds of trajectories and expensive quantum chemical methods to account for non-adiabatic excited state effects. Here we introduce a diabatic artificial neural network (DANN), based on diabatic states, to accelerate such simulations for azobenzene derivatives. The network is six orders of magnitude faster than the quantum chemistry method used for training. DANN is transferable to azobenzene molecules outside the training set, predicting quantum yields for unseen species that are correlated with experiment. We use the model to virtually screen 3100 hypothetical molecules, and identify novel species with high predicted quantum yields. The model predictions are confirmed using high-accuracy non-adiabatic dynamics. Our results pave the way for fast and accurate virtual screening of photoactive compounds.

The physics of liquid-to-solid transitions in multi-domain protein condensates

Many RNA-binding proteins (RBPs) that assemble into membraneless organelles have a common architecture including disordered prion-like domain (PLD) and folded RNA-binding domain (RBD). An enrichment of PLD within the condensed phase gives rise to formation, on longer time scales, of amyloid-like fibrils (aging). In this study, we employ coarse-grained Langevin dynamics simulations to explore the physical basis for the structural diversity in condensed phases of multi-domain RBPs. We discovered a highly cooperative first-order transition between disordered structures and an ordered phase whereby chains of PLD organize in fibrils with high nematic orientational order. An interplay between homodomain (PLD-PLD) and heterodomain (PLD-RBD) interactions results in variety of structures with distinct spatial architectures. Interestingly, the different structural phases also exhibit vastly different intracluster dynamics of proteins, with diffusion coefficients 5 times (disordered structures) to 50 times (ordered structures) lower than that of the dilute phase. Cooperativity of this liquid-solid transition makes fibril formation highly malleable to mutations or post-translational modifications. Our results provide a mechanistic understanding of how multi-domain RBPs could form assemblies with distinct structural and material properties.

Development of antibacterial compounds that constrain evolutionary pathways to resistance

Antibiotic resistance is a worldwide challenge. A potential approach to block resistance is to simultaneously inhibit WT and known escape variants of the target bacterial protein. Here, we applied an integrated computational and experimental approach to discover compounds that inhibit both WT and trimethoprim (TMP) resistant mutants of E. coli dihydrofolate reductase (DHFR). We identified a novel compound (CD15-3) that inhibits WT DHFR and its TMP resistant variants L28R, P21L and A26T with IC₅₀ 50–75 µM against WT and TMP-resistant strains. Resistance to CD15-3 was dramatically delayed compared to TMP in in vitro evolution. Whole genome sequencing of CD15-3-resistant strains showed no mutations in the target folA locus. Rather, gene duplication of several efflux pumps gave rise to weak (about twofold increase in IC₅₀) resistance against CD15-3. Altogether, our results demonstrate the promise of strategy to develop evolution drugs - compounds which constrain evolutionary escape routes in pathogens.

Effect of RNA on Morphology and Dynamics of Membraneless Organelles

Membraneless organelles (MLOs) are spatiotemporally regulated structures that concentrate multivalent proteins or RNA, often in response to stress. The proteins enriched within MLOs are often classified as high-valency "scaffolds" or low-valency "clients", with the former being associated with a phase-separation promoting role. In this study, we employ a minimal model for P-body components, with a defined protein-protein interaction network, to study their phase separation at biologically realistic low protein concentrations. Without RNA, multivalent proteins can assemble into solid-like clusters only in the regime of high concentration and stable interactions. RNA molecules promote cluster formation in an RNA-length-dependent manner, even in the regime of weak interactions and low protein volume fraction. Our simulations reveal that long RNA chains act as superscaffolds that stabilize large RNA-protein clusters by recruiting low-valency proteins within them while also ensuring functional "liquid-like" turnover of components. Our results suggest that RNA-mediated phase separation could be a plausible mechanism for spatiotemporally regulated phase separation in the cell.

Differential Enzyme Flexibility Probed Using Solid-State Nanopores

Enzymes and motor proteins are dynamic macromolecules that coexist in a number of conformations of similar energies. Protein function is usually accompanied by a change in structure and flexibility, often induced upon binding to ligands. However, while measuring protein flexibility changes between active and resting states is of therapeutic significance, it remains a challenge. Recently, our group has demonstrated that breadth of signal amplitudes in measured electrical signatures as an ensemble of individual protein molecules is driven through solid-state nanopores and correlates with protein conformational dynamics. Here, we extend our study to resolve subtle flexibility variation in dihydrofolate reductase mutants from unlabeled single molecules in solution. We first demonstrate using a canonical protein system, adenylate kinase, that both size and flexibility changes can be observed upon binding to a substrate that locks the protein in a closed conformation. Next, we investigate the influence of voltage bias and pore geometry on the measured electrical pulse statistics during protein transport. Finally, using the optimal experimental conditions, we systematically study a series of wild-type and mutant dihydrofolate reductase proteins, finding a good correlation between nanopore-measured protein conformational dynamics and equilibrium bulk fluorescence probe measurements. Our results unequivocally demonstrate that nanopore-based measurements reliably probe conformational diversity in native protein ensembles.

Rational Design of Novel Allosteric Dihydrofolate Reductase Inhibitors Showing Antibacterial Effects on Drug-Resistant Escherichia coli Escape Variants

In drug discovery, systematic variations of substituents on a common scaffold and bioisosteric replacements are often used to generate diversity and obtain molecules with better biological effects. However, this could saturate the small-molecule diversity pool resulting in drug resistance. On the other hand, conventional drug discovery relies on targeting known pockets on protein surfaces leading to drug resistance by mutations of critical pocket residues. Here, we present a two-pronged strategy of designing novel drugs that target unique pockets on a protein's surface to overcome the above problems. Dihydrofolate reductase, DHFR, is a critical enzyme involved in thymidine and purine nucleotide biosynthesis. Several classes of compounds that are structural analogues of the substrate dihydrofolate have been explored for their antifolate activity. Here, we describe 10 novel small-molecule inhibitors of Escherichia coli DHFR, EcDHFR, belonging to the stilbenoid, deoxybenzoin, and chalcone family of compounds discovered by a combination of pocket-based virtual ligand screening and systematic scaffold hopping. These inhibitors show a unique uncompetitive or noncompetitive inhibition mechanism, distinct from those reported for all known inhibitors of DHFR, indicative of binding to a unique pocket distinct from either substrate or cofactor-binding pockets. Furthermore, we demonstrate that rescue mutants of EcDHFR, with reduced affinity to all known classes of DHFR inhibitors, are inhibited at the same concentration as the wild-type. These compounds also exhibit antibacterial activity against E. coli harboring the drug-resistant variant of DHFR. This discovery is the first report on a novel class of inhibitors targeting a unique pocket on EcDHFR.

Exploring the Mutational Robustness of Nucleic Acids by Searching Genotype Neighborhoods in Sequence Space

To assess the mutational robustness of nucleic acids, many genome- and protein-level studies have been performed, where nucleic acids are treated as genetic information carriers and transferrers. However, the molecular mechanisms through which mutations alter the structural, dynamic, and functional properties of nucleic acids are poorly understood. Here we performed a SELEX in silico study to investigate the fitness distribution of the l-Arm-binding aptamer genotype neighborhoods. Two novel functional genotype neighborhoods were isolated and experimentally verified to have comparable fitness as the wild-type. The experimental aptamer fitness landscape suggests the mutational robustness is strongly influenced by the local base environment and ligand-binding mode, whereas bases distant from the binding pocket provide potential evolutionary pathways to approach the global fitness maximum. Our work provides an example of successful application of SELEX in silico to optimize an aptamer and demonstrates the strong sensitivity of mutational robustness to the site of genetic variation.

Searching the Sequence Space for Potent Aptamers Using SELEX in Silico

To isolate functional nucleic acids that bind to defined targets with high affinity and specificity, which are known as aptamers, the systematic evolution of ligands by exponential enrichment (SELEX) methodology has emerged as the preferred approach. Here, we propose a computational approach, SELEX in silico, that allows the sequence space to be more thoroughly explored regarding binding of a certain target. Our approach consists of two steps: (i) secondary structure-based sequence screening, which aims to collect the sequences that can form a desired RNA motif as an enhanced initial library, followed by (ii) sequence enrichment regarding target binding by molecular dynamics simulation-based virtual screening. Our SELEX in silico method provided a practical computational solution to three key problems in aptamer sequence searching: design of nucleic acid libraries, knowledge of sequence enrichment, and identification of potent aptamers. Six potent theophylline-binding aptamers, which were isolated by SELEX in silico from a sequence space containing 4(13) sequences, were experimentally verified to bind theophylline with high affinity: Kd ranging from 0.16 to 0.52 μM, compared with the dissociation constant of the original aptamer-theophylline, 0.32 μM. These results demonstrate the significant potential of SELEX in silico as a new method for aptamer discovery and optimization.

Systems-level response to point mutations in a core metabolic enzyme modulates genotype-phenotype relationship

Linking the molecular effects of mutations to fitness is central to understanding evolutionary dynamics. Here, we establish a quantitative relation between the global effect of mutations on the E. coli proteome and bacterial fitness. We created E. coli strains with specific destabilizing mutations in the chromosomal folA gene encoding dihydrofolate reductase (DHFR) and quantified the ensuing changes in the abundances of 2,000+ E. coli proteins in mutant strains using tandem mass tags with subsequent LC-MS/MS. mRNA abundances in the same E. coli strains were also quantified. The proteomic effects of mutations in DHFR are quantitatively linked to phenotype: the SDs of the distributions of logarithms of relative (to WT) protein abundances anticorrelate with bacterial growth rates. Proteomes hierarchically cluster first by media conditions, and within each condition, by the severity of the perturbation to DHFR function. These results highlight the importance of a systems-level layer in the genotype-phenotype relationship.

The interface of protein structure, protein biophysics, and molecular evolution

Abstract The interface of protein structural biology, protein biophysics, molecular evolution, and molecular population genetics forms the foundations for a mechanistic understanding of many aspects of protein biochemistry. Current efforts in interdisciplinary protein modeling are in their infancy and the state-of-the art of such models is described. Beyond the relationship between amino acid substitution and static protein structure, protein function, and corresponding organismal fitness, other considerations are also discussed. More complex mutational processes such as insertion and deletion and domain rearrangements and even circular permutations should be evaluated. The role of intrinsically disordered proteins is still controversial, but may be increasingly important to consider. Protein geometry and protein dynamics as a deviation from static considerations of protein structure are also important. Protein expression level is known to be a major determinant of evolutionary rate and several considerations including selection at the mRNA level and the role of interaction specificity are discussed. Lastly, the relationship between modeling and needed high-throughput experimental data as well as experimental examination of protein evolution using ancestral sequence resurrection and in vitro biochemistry are presented, towards an aim of ultimately generating better models for biological inference and prediction.

Parallel folding pathways in the SH3 domain protein

The transition-state ensemble (TSE) is the set of protein conformations with an equal probability to fold or unfold. Its characterization is crucial for an understanding of the folding process. We determined the TSE of the src-SH3 domain protein by using extensive molecular dynamics simulations of the Go model and computing the folding probability of a generated set of TSE candidate conformations. We found that the TSE possesses a well-defined hydrophobic core with variable enveloping structures resulting from the superposition of three parallel folding pathways. The most preferred pathway agrees with the experimentally determined TSE, while the two least preferred pathways differ significantly. The knowledge of the different pathways allows us to design the interactions between amino acids that guide the protein to fold through the least preferred pathway. This particular design is akin to a circular permutation of the protein. The finding motivates the hypothesis that the different experimentally observed TSEs in homologous proteins and circular permutants may represent potentially available pathways to the wild-type protein.

The mechanical unfolding of ubiquitin through all-atom Monte Carlo simulation with a Go-type potential

The mechanical unfolding of proteins under a stretching force has an important role in living systems and is a logical extension of the more general protein folding problem. Recent advances in experimental methodology have allowed the stretching of single molecules, thus rendering this process ripe for computational study. We use all-atom Monte Carlo simulation with a Gō-type potential to study the mechanical unfolding pathway of ubiquitin. A detailed, robust, well-defined pathway is found, confirming existing results in this vein though using a different model. Additionally, we identify the protein's fundamental stabilizing secondary structure interactions in the presence of a stretching force and show that this fundamental stabilizing role does not persist in the absence of mechanical stress. The apparent success of simulation methods in studying ubiquitin's mechanical unfolding pathway indicates their potential usefulness for future study of the stretching of other proteins and the relationship between protein structure and the response to mechanical deformation.

On the origin and highly likely completeness of single-domain protein structures

The size and origin of the protein fold universe is of fundamental and practical importance. Analyzing randomly generated, compact sticky homopolypeptide conformations constructed in generic simplified and all-atom protein models, all have similar folds in the library of solved structures, the Protein Data Bank, and conversely, all compact, single-domain protein structures in the Protein Data Bank have structural analogues in the compact model set. Thus, both sets are highly likely complete, with the protein fold universe arising from compact conformations of hydrogen-bonded, secondary structures. Because side chains are represented by their Cbeta atoms, these results also suggest that the observed protein folds are insensitive to the details of side-chain packing. Sequence specificity enters both in fine-tuning the structure and thermodynamically stabilizing a given fold with respect to the set of alternatives. Scanning the models against a three-dimensional active-site library, close geometric matches are frequently found. Thus, the presence of active-site-like geometries also seems to be a consequence of the packing of compact, secondary structural elements. These results have significant implications for the evolution of protein structure and function.

Protein structure and evolutionary history determine sequence space topology

Understanding the observed variability in the number of homologs of a gene is a very important unsolved problem that has broad implications for research into coevolution of structure and function, gene duplication, pseudogene formation, and possibly for emerging diseases. Here, we attempt to define and elucidate some possible causes behind the observed irregularity in sequence space. We present evidence that sequence variability and functional diversity of a gene or fold family is influenced by quantifiable characteristics of the protein structure. These characteristics reflect the structural potential for sequence plasticity, i.e., the ability to accept mutation without losing thermodynamic stability. We identify a structural feature of a protein domain-contact density-that serves as a determinant of entropy in sequence space, i.e., the ability of a protein to accept mutations without destroying the fold (also known as fold designability). We show that (log) of average gene family size exhibits statistical correlation (R(2) > 0.9.) with contact density of its three-dimensional structure. We present evidence that the size of individual gene families are influenced not only by the designability of the structure, but also by evolutionary history, e.g., the amount of time the gene family was in existence. We further show that our observed statistical correlation between gene family size and contact density of the structure is valid on many levels of evolutionary divergence, i.e., not only for closely related sequence, but also for less-related fold and superfamily levels of homology.

ELISA: structure-function inferences based on statistically significant and evolutionarily inspired observations

The problem of functional annotation based on homology modeling is primary to current bioinformatics research. Researchers have noted regularities in sequence, structure and even chromosome organization that allow valid functional cross-annotation. However, these methods provide a lot of false negatives due to limited specificity inherent in the system. We want to create an evolutionarily inspired organization of data that would approach the issue of structure-function correlation from a new, probabilistic perspective. Such organization has possible applications in phylogeny, modeling of functional evolution and structural determination. ELISA (Evolutionary Lineage Inferred from Structural Analysis, http://romi.bu.edu/elisa) is an online database that combines functional annotation with structure and sequence homology modeling to place proteins into sequence-structure-function "neighborhoods". The atomic unit of the database is a set of sequences and structural templates that those sequences encode. A graph that is built from the structural comparison of these templates is called PDUG (protein domain universe graph). We introduce a method of functional inference through a probabilistic calculation done on an arbitrary set of PDUG nodes. Further, all PDUG structures are mapped onto all fully sequenced proteomes allowing an easy interface for evolutionary analysis and research into comparative proteomics. ELISA is the first database with applicability to evolutionary structural genomics explicitly in mind.

AVAILABILITY

The database is available at http://romi.bu.edu/elisa.

Dynamic charge-density correlation function in weakly charged polyampholyte globules

We study solutions of statistically neutral polyampholyte chains containing a large fraction of neutral monomers. It is known that such solutions phase separate at very low concentrations, even if the quality of the solvent with respect to the neutral monomers is good. The precipitate is semidilute if the chains are weakly charged. This paper considers straight theta solvents and good solvents, and we calculate the dynamic charge density correlation function g(k,t) in the precipitate, using the quadratic approximation to the Martin-Siggia-Rose generating functional. It is convenient to express the results in terms of dimensionless space and time variables. Let xi be the blob size, and let tau be the characteristic time scale at the blob level. Define the dimensionless wave vector q=xik, and the dimensionless time s=t/tau. In the regime q<1, corresponding to length scales larger than the blob size, and 1<s<q(-4), corresponding to time scales in between the blob relaxation time and the relaxation time at scale q(-1), we find that the charge density fluctuations relax according to the power law g(q,s) approximately q(2)s(-1/2). This relaxation is qualitatively different from that of a neutral semidilute polymer solution. We expect our results to be valid for wave vectors q>0.1, where entanglements are unimportant.

Protein folding theory: from lattice to all-atom models

This review focuses on recent advances in understanding protein folding kinetics in the context of nucleation theory. We present basic concepts such as nucleation, folding nucleus, and transition state ensemble and then discuss recent advances and challenges in theoretical understanding of several key aspects of protein folding kinetics. We cover recent topology-based approaches as well as evolutionary studies and molecular dynamics approaches to determine protein folding nucleus and analyze other aspects of folding kinetics. Finally, we briefly discuss successful all-atom Monte-Carlo simulations of protein folding and conclude with a brief outlook for the future.

Evolutionary conservation of the folding nucleus

Here, we present statistical analysis of conservation profiles in families of homologous sequences for nine proteins whose folding nucleus was determined by protein engineering methods. We show that in all but one protein (AcP) folding nucleus residues are significantly more conserved than the rest of the protein. Two aspects of our study are especially important: (i) grouping of amino acid residues into classes according to their physical-chemical properties and (ii) proper normalization of amino acid probabilities that reflects the fact that evolutionary pressure to conserve some amino acid types may itself affect concentration of various amino acid types in protein families. Neglect of any of those two factors may make physical and biological "signals" from conservation profiles disappear.

Factors that affect the folding ability of proteins

The folding ability of a heteropolymer model for proteins subject to Monte Carlo dynamics on a simple cubic lattice is shown to be strongly correlated with the stability of the native state. We consider a number of estimates of the stability that can be determined without simulation, including the energy gap between the native state and the structurally dissimilar part of the spectrum (Z score) and, for sequences with fully compact native states, the gap in energy between the native and first excited fully compact states. These estimates are found to be more robust predictors of folding ability than a parameter sigma that requires simulation for its evaluation: sigma = 1 - Tf/Ttheta, where Tf is the temperature at which the fluctuation of an order parameter is at its maximum and Ttheta is the temperature at which the specific heat is at its maximum. We show that the interpretation of Ttheta as the collapse transition temperature is not correct in general and that the correlation between sigma and the folding ability arises from the fact that sigma is related to the energy gap (Z score).

A strategy for detecting the conservation of folding-nucleus residues in protein superfamilies

BACKGROUND

Nucleation-growth theory predicts that fast-folding peptide sequences fold to their native structure via structures in a transition-state ensemble that share a small number of native contacts (the folding nucleus). Experimental and theoretical studies of proteins suggest that residues participating in folding nuclei are conserved among homologs. We attempted to determine if this is true in proteins with highly diverged sequences but identical folds (superfamilies).

RESULTS

We describe a strategy based on comparisons of residue conservation in natural superfamily sequences with simulated sequences (generated with a Monte-Carlo sequence design strategy) for the same proteins. The basic assumptions of the strategy were that natural sequences will conserve residues needed for folding and stability plus function, the simulated sequences contain no functional conservation, and nucleus residues make native contacts with each other. Based on these assumptions, we identified seven potential nucleus residues in ubiquitin superfamily members. Non-nucleus conserved residues were also identified; these are proposed to be involved in stabilizing native interactions. We found that all superfamily members conserved the same potential nucleus residue positions, except those for which the structural topology is significantly different.

CONCLUSIONS

Our results suggest that the conservation of the nucleus of a specific fold can be predicted by comparing designed simulated sequences with natural highly diverged sequences that fold to the same structure. We suggest that such a strategy could be used to help plan protein folding and design experiments, to identify new superfamily members, and to subdivide superfamilies further into classes having a similar folding mechanism.

Conserved residues and the mechanism of protein folding

Experimental and simulation studies show that small monomeric proteins fold in one kinetic step, which entails overcoming the free-energy barrier between the unfolded and the native protein through a transition state. Two models of transition state formation have been proposed: a 'nonspecific' one in which it depends on the formation of a sufficient number of native-like contacts regardless of what amino acids are involved, and a 'specific' one, in which it depends on formation of a specific subset of the native structure (a folding nucleus). The latter requires that some amino acids form most of their contacts in the transition state, whereas others only do so on reaching the native conformation. If so, mutations affecting the stability of the transition state nucleus should have a greater effect on the folding kinetics than mutations elsewhere, and the residues involved should be evolutionarily conserved. Lattice-model simulations and experiments suggest that such mutations exist. Here we present a method for determining the folding nucleus of a protein with known structure with two-state folding kinetics. This method is based on the alignment of many sequences designed to fold into the native conformation of a protein to identify the positions where amino acids are most conserved in designed sequences. The method is applied to chymotrypsin inhibitor 2 (CI2), a protein whose transition state has been previously studied by protein engineering. The involvement of residues in folding nucleus of CI2 is clearly correlated with their conservation in design, and the residues forming the nucleus are highly conserved in 23 natural sequences homologous to CI2.