Lewyn Li

Topological determinants of protein folding

The folding of many small proteins is kinetically a two-state process that represents overcoming the major free-energy barrier. A kinetic characteristic of a conformation, its probability to descend to the native state domain in the amount of time that represents a small fraction of total folding time, has been introduced to determine to which side of the free-energy barrier a conformation belongs. However, which features make a protein conformation on the folding pathway become committed to rapidly descending to the native state has been a mystery. Using two small, well characterized proteins, CI2 and C-Src SH3, we show how topological properties of protein conformations determine their kinetic ability to fold. We use a macroscopic measure of the protein contact network topology, the average graph connectivity, by constructing graphs that are based on the geometry of protein conformations. We find that the average connectivity is higher for conformations with a high folding probability than for those with a high probability to unfold. Other macroscopic measures of protein structural and energetic properties such as radius of gyration, rms distance, solvent-accessible surface area, contact order, and potential energy fail to serve as predictors of the probability of a given conformation to fold.

Kinetics, thermodynamics and evolution of non-native interactions in a protein folding nucleus.

A lattice model with side chains was used to investigate protein folding with computer simulations. In this model, we rigorously demonstrate the existence of a specific folding nucleus. This nucleus contains specific interactions not present in the native state that, when weakened, slow folding but do not change protein stability. Such a decoupling of folding kinetics from thermodynamics has been observed experimentally for real proteins. From our results, we conclude that specific non-native interactions in the transition state would give rise to φ-values that are negative or larger than unity. Furthermore, we demonstrate that residue Ile 34 in src SH3, which has been shown to be kinetically, but not thermodynamically, important, is universally conserved in proteins with the SH3 fold. This is a clear example of evolution optimizing the folding rate of a protein independent of its stability and function.

Constructing, verifying, and dissecting the folding transition state of chymotrypsin inhibitor 2 with all-atom simulations

Experimentally, protein engineering and φ-value analysis is the method of choice to characterize the structure in folding transition state ensemble (TSE) of any protein. Combining experimental φ values and computer simulations has led to a deeper understanding of how proteins fold. In this report, we construct the TSE of chymotrypsin inhibitor 2 from published φ values. Importantly, we verify, by means of multiple independent simulations, that the conformations in the TSE have a probability of ≈0.5 to reach the native state rapidly, so the TSE consists of true transition states. This finding validates the use of transition state theory underlying all φ-value analyses. Also, we present a method to dissect and study the TSE by generating conformations that have a disrupted α-helix (α-disrupted states) or disordered β-strands 3 and 4 (β-disrupted states). Surprisingly, the α-disrupted states have a stronger tendency to fold than the β-disrupted states, despite the higher φ values for the α-helix in the TSE. We give a plausible explanation for this result and discuss its implications on protein folding and design. Our study shows that, by using both experiments and computer simulations, we can gain many insights into protein folding.

Amino acids determining enzyme-substrate specificity in prokaryotic and eukaryotic protein kinases

The binding between a PK and its target is highly specific, despite the fact that many different PKs exhibit significant sequence and structure homology. There must be, then, specificity-determining residues (SDRs) that enable different PKs to recognize their unique substrate. Here we use and further develop a computational procedure to discover putative SDRs (PSDRs) in protein families, whereby a family of homologous proteins is split into orthologous proteins, which are assumed to have the same specificity, and paralogous proteins, which have different specificities. We reason that PSDRs must be similar among orthologs, whereas they must necessarily be different among paralogs. Our statistical procedure and evolutionary model identifies such residues by discriminating a functional signal from a phylogenetic one. As case studies we investigate the prokaryotic two-component system and the eukaryotic AGC (i.e., cAMP-dependent PK, cGMP-dependent PK, and PKC) PKs. Without using experimental data, we predict PSDRs in prokaryotic and eukaryotic PKs, and suggest precise mutations that may convert the specificity of one PK to another. We compare our predictions with current experimental results and obtain considerable agreement with them. Our analysis unifies much of existing data on PK specificity. Finally, we find PSDRs that are outside the active site. Based on our results, as well as structural and biochemical characterizations of eukaryotic PKs, we propose the testable hypothesis of “specificity via differential activation” as a way for the cell to control kinase specificity.

Somatic Cell Fusions Reveal Extensive Heterogeneity in Basal-like Breast Cancer

Basal-like and luminal breast tumors have distinct clinical behavior and molecular profiles, yet the underlying mechanisms are poorly defined. To interrogate processes that determine these distinct phenotypes and their inheritance pattern, we generated somatic cell fusions and performed integrated genetic and epigenetic (DNA methylation and chromatin) profiling. We found that the basal-like trait is generally dominant and is largely defined by epigenetic repression of luminal transcription factors. Definition of super-enhancers highlighted a core program common in luminal cells but a high degree of heterogeneity in basal-like breast cancers that correlates with clinical outcome. We also found that protein extracts of basal-like cells are sufficient to induce a luminal-to-basal phenotypic switch, implying a trigger of basal-like autoregulatory circuits. We determined that KDM6A might be required for luminal-basal fusions, and we identified EN1, TBX18, and TCF4 as candidate transcriptional regulators of the luminal-to-basal switch. Our findings highlight the remarkable epigenetic plasticity of breast cancer cells.

Structural insights into the effects of 2'-5' linkages on the RNA duplex.

Significance The nonenzymatic replication of RNA is thought to have been a critical step in the emergence of simple cellular life from prebiotic chemistry. However, the chemical copying of RNA templates generates product strands that contain 2′-5′ backbone linkages and normal 3′-5′ linkages. Our recent finding that RNAs with such mixed backbones can still fold into functional structures raised the question of how RNA accommodates the presence of 2′-5′ linkages. Here we use X-ray crystallography and molecular dynamics simulations to reveal how 3′-5′–linked RNA duplexes accommodate interspersed 2′-5′ linkages. The diminished thermal and chemical stability of such RNA duplexes reflects local structural changes, but compensatory changes result in a global RNA duplex structure with relatively minor alterations. The mixture of 2′-5′ and 3′-5′ linkages generated during the nonenzymatic replication of RNA has long been regarded as a central problem for the origin of the RNA world. However, we recently observed that both a ribozyme and an RNA aptamer retain considerable functionality in the presence of prebiotically plausible levels of linkage heterogeneity. To better understand the RNA structure and function in the presence of backbone linkage heterogeneity, we obtained high-resolution X-ray crystal structures of a native 10-mer RNA duplex (1.32 Å) and two variants: one containing one 2′-5′ linkage per strand (1.55 Å) and one containing three such linkages per strand (1.20 Å). We found that RNA duplexes adjust their local structures to accommodate the perturbation caused by 2′-5′ linkages, with the flanking nucleotides buffering the disruptive effects of the isomeric linkage and resulting in a minimally altered global structure. Although most 2′-linked sugars were in the expected 2′-endo conformation, some were partially or fully in the 3′-endo conformation, suggesting that the energy difference between these conformations was relatively small. Our structural and molecular dynamic studies also provide insight into the diminished thermal and chemical stability of the duplex state associated with the presence of 2′-5′ linkages. Our results contribute to the view that a low level of 2′-5′ substitution would not have been fatal in an early RNA world and may in contrast have been helpful for both the emergence of nonenzymatic RNA replication and the early evolution of functional RNAs.

Chromatin immunoprecipitation from fixed clinical tissues reveals tumor-specific enhancer profiles

Extensive cross-linking introduced during routine tissue fixation of clinical pathology specimens severely hampers chromatin immunoprecipitation followed by next-generation sequencing (ChIP-seq) analysis from archived tissue samples. This limits the ability to study the epigenomes of valuable, clinically annotated tissue resources. Here we describe fixed-tissue chromatin immunoprecipitation sequencing (FiT-seq), a method that enables reliable extraction of soluble chromatin from formalin-fixed paraffin-embedded (FFPE) tissue samples for accurate detection of histone marks. We demonstrate that FiT-seq data from FFPE specimens are concordant with ChIP-seq data from fresh-frozen samples of the same tumors. By using multiple histone marks, we generate chromatin-state maps and identify cis-regulatory elements in clinical samples from various tumor types that can readily allow us to distinguish between cancers by the tissue of origin. Tumor-specific enhancers and superenhancers that are elucidated by FiT-seq analysis correlate with known oncogenic drivers in different tissues and can assist in the understanding of how chromatin states affect gene regulation.

ChiLin: a comprehensive ChIP-seq and DNase-seq quality control and analysis pipeline

BackgroundTranscription factor binding, histone modification, and chromatin accessibility studies are important approaches to understanding the biology of gene regulation. ChIP-seq and DNase-seq have become the standard techniques for studying protein-DNA interactions and chromatin accessibility respectively, and comprehensive quality control (QC) and analysis tools are critical to extracting the most value from these assay types. Although many analysis and QC tools have been reported, few combine ChIP-seq and DNase-seq data analysis and quality control in a unified framework with a comprehensive and unbiased reference of data quality metrics.ResultsChiLin is a computational pipeline that automates the quality control and data analyses of ChIP-seq and DNase-seq data. It is developed using a flexible and modular software framework that can be easily extended and modified. ChiLin is ideal for batch processing of many datasets and is well suited for large collaborative projects involving ChIP-seq and DNase-seq from different designs. ChiLin generates comprehensive quality control reports that include comparisons with historical data derived from over 23,677 public ChIP-seq and DNase-seq samples (11,265 datasets) from eight literature-based classified categories. To the best of our knowledge, this atlas represents the most comprehensive ChIP-seq and DNase-seq related quality metric resource currently available. These historical metrics provide useful heuristic quality references for experiment across all commonly used assay types. Using representative datasets, we demonstrate the versatility of the pipeline by applying it to different assay types of ChIP-seq data. The pipeline software is available open source at https://github.com/cfce/chilin.ConclusionChiLin is a scalable and powerful tool to process large batches of ChIP-seq and DNase-seq datasets. The analysis output and quality metrics have been structured into user-friendly directories and reports. We have successfully compiled 23,677 profiles into a comprehensive quality atlas with fine classification for users.

Relating helix tilt in a bilayer to lipid disorder: a mean-field theory.

We present a mean-field theory relating the helix tilt angle in a bilayer to lipid disorder. The theory provides a method to compare the rotational barriers for different helices in lipid bilayers. The results suggest that the helix tilt angle is strongly affected by both the hydrophobicity of the helix and the average lipid disorder. This leads us to point out future experiments that could shed light on lipid-protein interactions.

Directed nucleation and growth by balancing local supersaturation and substrate/nucleus lattice mismatch

Significance The energy barrier for a classical heterogeneous crystal nucleation can be controlled by the energy contributions from the substrate/nucleus interface and local supersaturation. Exerting control over crystal growth thus requires modifying either one of these terms. We here introduce a strategy to modulate the contributions of both parameters simultaneously using substrates containing different crystal structures of calcium carbonate. Based on a theoretical analysis, we program both the positioning and growth direction of carbonate salts on preselected polymorphs. These findings may hold relevance for understanding, mimicking, and ultimately expanding upon nature’s mineralization strategies and for developing functional microscale materials. Controlling nucleation and growth is crucial in biological and artificial mineralization and self-assembly processes. The nucleation barrier is determined by the chemistry of the interfaces at which crystallization occurs and local supersaturation. Although chemically tailored substrates and lattice mismatches are routinely used to modify energy landscape at the substrate/nucleus interface and thereby steer heterogeneous nucleation, strategies to combine this with control over local supersaturations have remained virtually unexplored. Here we demonstrate simultaneous control over both parameters to direct the positioning and growth direction of mineralizing compounds on preselected polymorphic substrates. We exploit the polymorphic nature of calcium carbonate (CaCO3) to locally manipulate the carbonate concentration and lattice mismatch between the nucleus and substrate, such that barium carbonate (BaCO3) and strontium carbonate (SrCO3) nucleate only on specific CaCO3 polymorphs. Based on this approach we position different materials and shapes on predetermined CaCO3 polymorphs in sequential steps, and guide the growth direction using locally created supersaturations. These results shed light on nature’s remarkable mineralization capabilities and outline fabrication strategies for advanced materials, such as ceramics, photonic structures, and semiconductors.