I am working in the field of theoretical and computational (bio)chemsitry, with particular emphasis on the development and application of multiscale computational techniques. Such techniques allow for multi-resolution investigation of molecular systems using all-atom, coarse-grained (CG), and mesoscale simulations. While all-atom simulation methods provide invaluable means for studying molecular systems, the available computational resources limit their direct applicability to many important biological processes due to the large length and time scales on which those processes occur. The CG methods provide a molecular-level description but are significantly more efficient due to the reduction in the number of degrees of freedom that represent the system. This is achieved by grouping some atoms into CG beads and formally eliminating others that implicitly influence the effective CG interactions. The mesoscale methods operate beyond molecular resolution, using continuum-level theories. Combined with systematic discretization techniques, they provide powerful tools for studying biological processes on submicrometer length scales. An example of such description can be the quasi-particle representation of a membrane, where each particle represents a patch containing hundreds or thousands of lipids. The clever combination of those methods and the link between them is a key for addressing many difficult biological and materials science problems.
Associative Memory, Water Mediated, Structure, and Energy Model for Molecular Dynamics (AWSEM-MD) is a CG model that is capable of predicting protein native state and protein-protein association using a fraction of computational time that would have been required using all-atom simulations. AWSEM combines the physically motivated backbone and tertiary potentials with bioinformatics-based local interactions that are parameterized using experimental information (e.g. from the PDB database). Thus, AWSEM-MD can be used for de novo structure prediction as well as for various dynamic and kinetic studies in the absence of a priory known structure of the protein of interest or its homologs. This method was developed during my doctoral studies at the University of Maryland under the supervision of Prof. Garegin A. Papoian in collaboration with Prof. Peter G. Wolynes's group (currently at Rice University) and is used by a wider scientific community today.
Official website: http://www.awsem-md.org
Direct source link: https://github.com/adavtyan/awsemmd
Protein folding in vivo
In vitro experiments and corresponding simulations in a dilute environment taught us a lot about how proteins fold and interact with each other. However, those conditions are drastically different from crowded and rapidly changing environment inside the cell. Part of my doctoral work was devoted to the investigation of protein folding under fluctuating environmental conditions. Using a combination of simple structure-based CG simulations and stochastic theories, we have shown that both periodic and colored random noise can accelerate protein folding if the characteristic time scale of the noise coincides with the folding time, consistent with the description of the stochastic resonance phenomenon. Our theoretical and computational prediction was confirmed experimentally by Prof. Martin Gruebele’s group at the University of Illinois, Urbana-Champaign. We continue working in this direction to investigate the implications of this finding on various biological processes.
Dynamic Force Matching
Bottom-up CGing approaches make it possible to extend the applicability of atomistic simulations to larger length and time scales, thus allowing the study of larger systems on time scales that otherwise would not be possible to achieve with the computational resources at hand. While such CG models are usually successful in reproducing equilibrium properties of the all-atom system (such as radial distribution function), they often fail to reproduce dynamic properties, such as self-diffusion rates and other transport coefficients. More specifically, faster diffusion is typically exhibited by the CG system due to the absence of the fluctuation forces from the movement of individual atoms that were removed or reduced into effective CG sites. In collaboration with Prof. Hans C. Andersen from Stanford University, we have developed a method for the construction of dynamic CG models that accurately reproduce both equilibrium and dynamic properties. This was achieved by reintroducing the fluctuation forces into the "non-dynamic" CG model in a manner that does not disrupt any of its equilibrium properties and does not affect the model's efficiency in a notable way.
One aspect of my work was devoted to mesoscopic modeling of membrane-protein systems using a discrete quasi-particle approach. It has been shown before that such approaches allow for the investigation of large-scale transformation of plasma membranes (or membrane remodeling) under the influence of curvature generating proteins. The Mesoscopic Membrane with Proteins Model (MesM-P) that we formulated makes it possible to target principally new problems and phenomena on the experimental length and time scales. MesM-P has been implemented in a highly scalable and computationally efficient manner as part of LAMMPS open-source package and is available to anyone interested in mesoscale simulations of membranes-protein systems. Similar quasi-continuum approaches can be used to study a wide range of biological problems, including polymer melts (e.g. of actin or DNA) and active nematic systems driven by motors.
Source code available at https://github.com/adavtyan/MesM-P
I am interested in the fundamental understanding of molecular mechanisms behind vital biological processes, such as innate and adaptive immunity, gene regulation, and metabolism. I am developing novel and efficient computational methods to enable large-scale and long-time simulations of complex biological systems containing tens or hundreds of proteins, nucleic acids, and biological membranes. The predictions made by those methods will promote the fundamental understanding as well as the development of new therapeutic strategies and drugs.
A. Davtyan, M. Y. Tsai, F. Bai, and P. G. Wolynes, “Rational Design of Drugs Against Oligomer-mediated Neurodegeneration Using Novel Ensemble of Oligomer Structures of Amyloid β”, (working paper)
Y. Xu, K. Knapp, K.N. Le, N. P. Schafer, M. S. Safari, A. Davtyan, P. G. Wolynes, P. G. Vekilov, “Frustrated peptide chains at the fibril tip control the kinetics of growth of amyloid-β fibrils”, PNAS, 2021, 118 (38), e2110995118.
N. Sobhani, D.R. Tardiel-Cyril, A. Davtyan, D. Generali, R. Roudi, Y. Li, “CTLA-4 in Regulatory T Cells for Cancer Immunotherapy”, Cancers, 2021, 13(6), 1440.
D. S. Yang, A. Saeedi, A. Davtyan, M. Fathi, M. B. Sherman, M. S. Safari, A. Klindziuk, M. C. Barton, N. Varadarajan, A. B. Kolomeisky, P. G. Vekilov, “Mesoscopic protein-rich clusters host the nucleation of mutant p53 amyloid fibrils”, PNAS, 2021, 118(10), e2015618118.
W. Lu, C. Bueno, N. P. Schafer, J. Moller, S. Jin, X. Chen, M. Chen, X. Gu, A. Davtyan, J. J. de Pablo, P. G. Wolynes, “OpenAWSEM with Open3SPN2: A fast, flexible, and accessible framework for large-scale coarse-grained biomolecular simulations”, PLoS Comput. Biol., 2021, 17(2), e1008308.
D. Krepel, A. Davtyan, N. P. Schafer, P. G. Wolynes, J. N. Onuchic, “Braiding topology and the energy landscape of chromosome organization proteins”, PNAS, 2020, 117(3), 1468-1477
S. Jin, V. G. Contessoto, M. Chen, N.P. Schafer, W. Lu, X. Chen, C. Bueno, A. Hajitaheri, B. J. Sirovetz, A. Davtyan, G. A. Papoian, M. Y. Tsai, P. G. Wolynes, “AWSEM-Suite: a protein structure prediction server based on template-guided, coevolutionary-enhanced optimized folding landscapes”, Nucleic Acids Res., 2020.
Z. Jarin, F.-C. Tsai, A. Davtyan, A. J. Pak, P. Bassereau, G. A. Voth, “Unusual Organization of I-BAR Proteins on Tubular and Vesicular Membranes”, Biophys. J., 2019, 117(3), 553-562.
A. Davtyan, A. B. Kolomeisky “Theoretical Insights into Mechanisms of Stochastic Gating in Channel-Facilitated Molecular Transport”, J. Chem. Phys., 2019, 150(12), 124111
H. H. Katkar, A. Davtyan, A. E. P. Durumeric, G. M. Hocky, A. Schramm, M. Enrique, G. A. Voth, “Insights into the Cooperative Nature of ATP Hydrolysis in Actin Filaments”, Biophys. J., 2018, 115(8), 1589-1602
A. Davtyan, M. Simunovic, G. A. Voth, “The Mesoscopic Membrane with Proteins (MesM-P) Model”, J. Chem. Phys., 2017, 147, 044101
A. Davtyan, G. A. Voth, H. C. Andersen, “Dynamic Force Matching: Construction of Dynamical Coarse-grained Models for Simultaneous Capturing of Translational and Rotational Diffusion Rates”, J. Chem. Phys., 2016, 145, 224107
A. Davtyan, M. Simunovic, G. A. Voth, “Multiscale Simulations of Protein Facilitated Membrane Remodeling”, J. Struct. Biol., 2016, 196(1), 57-63
A. Davtyan, M. Platkov, M. Gruebele, G. A. Papoian, “Stochastic Resonance in Protein Folding Dynamics”, ChemPhysChem, 2016, 17(9), 1305-1313
K. Dave, A. Davtyan, M. Platkov, G. A. Papoian, M. Gruebele, “Environmental Fluctuations and Stochastic Resonance in Protein Folding”, ChemPhysChem, 2016, 17(9), 1341-1348
A. Davtyan, J. F. Dama, G. A. Voth, H. C. Andersen, “Dynamic Force Matching: A Method for Constructing Dynamical Coarse-Grained Models with Realistic Time Dependence”, J. Chem. Phys., 2015, 142, 154104
A. Davtyan, J. F. Dama, A. V. Sinitskiy, G. A. Voth, “The Theory of Ultra-Coarse-Graining. 2. Numerical Implementation”, J. Chem. Theory Comput., 2014, 10(12), 5265-5275
W. Zheng, N. Schafer, A. Davtyan, G. A. Papoian, P. G. Wolynes , “Predictive Energy Landscapes for Protein-Protein Association”, PNAS, 2012, 109(47), 19244-19249
A. Davtyan, W. Zheng, N. Schafer, C. Clementi, P. G. Wolynes and G. A Papoian, “AWSEM-MD: Protein Structure Prediction Using Coarse-grained Physical Potential and Bioinformatically Based Local Structure Biasing”, J. Phys. Chem. B., 2012, 116(29), 8494-8503