Source author record

Eugene Shakhnovich

Eugene Shakhnovich appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

7works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

7 published item(s)

preprint2023arXiv

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 tradeoffs between barriers and absorption wavelengths, and open-source our data and software to accelerate research in photopharmacology.

preprint2020arXiv

Liquid-liquid microphase separation leads to formation of membraneless organelles

Proteins and nucleic acids can spontaneously self-assemble into membraneless droplet-like compartments, both in vitro and in vivo. A key component of these droplets are multi-valent proteins that possess several adhesive domains with specific interaction partners (whose number determines total valency of the protein) separated by disordered regions. Here, using multi-scale simulations we show that such proteins self-organize into micro-phase separated droplets of various sizes as opposed to the Flory-like macro-phase separated equilibrium state of homopolymers or equilibrium physical gels. We show that the micro-phase separated state is a dynamic outcome of the interplay between two competing processes: a diffusion-limited encounter between proteins, and the dynamics within small clusters that results in exhaustion of available valencies whereby all specifically interacting domains find their interacting partners within smaller clusters, leading to arrested phase separation. We first model these multi-valent chains as bead-spring polymers with multiple adhesive domains separated by semi-flexible linkers and use Langevin Dynamics (LD) to assess how key timescales depend on the molecular properties of associating polymers. Using the time-scales from LD simulations, we develop a coarse-grained kinetic model to study this phenomenon at longer times. Consistent with LD simulations, the macro-phase separated state was only observed at high concentrations and large interaction valencies. Further, in the regime where cluster sizes approach macro-phase separation, the condensed phase becomes dynamically solid-like, suggesting that it might no longer be biologically functional. Therefore, the micro-phase separated state could be a hallmark of functional droplets formed by proteins with the sticker-spacer architecture.

preprint2016arXiv

Benchmarking inverse statistical approaches for protein structure and design with exactly solvable models

Inverse statistical approaches to determine protein structure and function from Multiple Sequence Alignments (MSA) are emerging as powerful tools in computational biology. However the underlying assumptions of the relationship between the inferred effective Potts Hamiltonian and real protein structure and energetics remain untested so far. Here we use lattice protein model (LP) to benchmark those inverse statistical approaches. We build MSA of highly stable sequences in target LP structures, and infer the effective pairwise Potts Hamiltonians from those MSA. We find that inferred Potts Hamiltonians reproduce many important aspects of 'true' LP structures and energetics. Careful analysis reveals that effective pairwise couplings in inferred Potts Hamiltonians depend not only on the energetics of the native structure but also on competing folds; in particular, the coupling values reflect both positive design (stabilization of native conformation) and negative design (destabilization of competing folds). In addition to providing detailed structural information, the inferred Potts models used as protein Hamiltonian for design of new sequences are able to generate with high probability completely new sequences with the desired folds, which is not possible using independent-site models. Those are remarkable results as the effective LP Hamiltonians used to generate MSA are not simple pairwise models due to the competition between the folds. Our findings elucidate the reasons for the success of inverse approaches to the modelling of proteins from sequence data, and their limitations.

preprint2009arXiv

Interplay between pleiotropy and secondary selection determines rise and fall of mutators in stress response

Dramatic rise of mutators has been found to accompany adaptation of bacteria in response to many kinds of stress. Two views on the evolutionary origin of this phenomenon emerged: the pleiotropic hypothesis positing that it is a byproduct of environmental stress or other specific stress response mechanisms and the second order selection which states that mutators hitchhike to fixation with unrelated beneficial alleles. Conventional population genetics models could not fully resolve this controversy because they are based on certain assumptions about fitness landscape. Here we address this problem using a microscopic multiscale model, which couples physically realistic molecular descriptions of proteins and their interactions with population genetics of carrier organisms without assuming any a priori fitness landscape. We found that both pleiotropy and second order selection play a crucial role at different stages of adaptation: the supply of mutators is provided through destabilization of error correction complexes or fluctuations of production levels of prototypic mismatch repair proteins (pleiotropic effects), while rise and fixation of mutators occur when there is a sufficient supply of beneficial mutations in replication-controlling genes. This general mechanism assures a robust and reliable adaptation of organisms to unforeseen challenges. This study highlights physical principles underlying physical biological mechanisms of stress response and adaptation.

preprint2007arXiv

A first-principles model of early evolution: Emergence of gene families, species and preferred protein folds

In this work we develop a microscopic physical model of early evolution, where phenotype,organism life expectancy, is directly related to genotype, the stability of its proteins in their native conformations which can be determined exactly in the model. Simulating the model on a computer, we consistently observe the Big Bang scenario whereby exponential population growth ensues as soon as favorable sequence-structure combinations (precursors of stable proteins) are discovered. Upon that, random diversity of the structural space abruptly collapses into a small set of preferred proteins. We observe that protein folds remain stable and abundant in the population at time scales much greater than mutation or organism lifetime, and the distribution of the lifetimes of dominant folds in a population approximately follows a power law. The separation of evolutionary time scales between discovery of new folds and generation of new sequences gives rise to emergence of protein families and superfamilies whose sizes are power-law distributed, closely matching the same distributions for real proteins. On the population level we observe emergence of species, subpopulations which carry similar genomes. Further we present a simple theory that relates stability of evolving proteins to the sizes of emerging genomes. Together, these results provide a microscopic first principles picture of how first gene families developed in the course of early evolution

preprint2007arXiv

Universality and diversity of folding mechanics for three-helix bundle proteins

In this study we evaluate, at full atomic detail, the folding processes of two small helical proteins, the B domain of protein A and the Villin headpiece. Folding kinetics are studied by performing a large number of ab initio Monte Carlo folding simulations using a single transferable all-atom potential. Using these trajectories, we examine the relaxation behavior, secondary structure formation, and transition-state ensembles (TSEs) of the two proteins and compare our results with experimental data and previous computational studies. To obtain a detailed structural information on the folding dynamics viewed as an ensemble process, we perform a clustering analysis procedure based on graph theory. Moreover, rigorous pfold analysis is used to obtain representative samples of the TSEs and a good quantitative agreement between experimental and simulated Fi-values is obtained for protein A. Fi-values for Villin are also obtained and left as predictions to be tested by future experiments. Our analysis shows that two-helix hairpin is a common partially stable structural motif that gets formed prior to entering the TSE in the studied proteins. These results together with our earlier study of Engrailed Homeodomain and recent experimental studies provide a comprehensive, atomic-level picture of folding mechanics of three-helix bundle proteins.

preprint2000arXiv

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, even if the quality of the solvent with respect to the neutral monomers is good, a long chain will collapse into a globule. For weakly charged chains, the interior of this globule is semi-dilute. This paper considers mainly theta-solvents, and we calculate the dynamic charge density correlation function g(k,t) in the interior of the globules, 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 R be the blob size, and let T be the characteristic time scale at the blob level. Define the dimensionless wave vector q = R k, and the dimensionless time s = t/T. We find that for q<1, corresponding to length scales larger than the blob size, the charge density fluctuations relax according to g(q,s) = q^2(1-s^(1/2)) at short times s < 1, and according to g(q,s) = q^2 s^(-1/2) at intermediate times 1 < s < q^(-4). We expect these results to be valid for wave vectors q > 0.1, where entanglements are unimportant.