I am working in the field of theoretical and computational (bio)chemsitry, with particular emphasis on 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 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 physically motivated backbone and tertiary potentials with bioinformatics-based local interactions that are parameterized using experimental information (e.g. from 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 (currentely at Rice University) and is used by a wider scientific community today.
Offical 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 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 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, as one of my ongoing postdoctoral projects in Voth group, we are developing a method for construction of dynamic CG models that accurately reproduce both equilibrium and dynamic properties. This is 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.
My latest work in Voth group is devoted to mesoscopic modeling of membrane-protein systems using discrete quasi-particle approach. It has been shown before that such approaches allow for 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 recently formulated makes it possible to target principally new problems and phenomena on 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 will be available to anyone interested in mesoscale simulations of membranes-protein systems after the publication of the first paper that introduces the method and that is currently in works. 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 fundamental understanding of molecular mechanisms behind vital biological processes, such as cell signaling and bacterial resistance. My group will develop novel and efficient computational methods that will enable large-scale and long-time simulations of complex biological systems containing tens or hundreds of proteins and biological membranes. The predictions made by those methods will promote the development of new drugs and therapeutic treatments.
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" (working paper)
Z. Jarin, A. Davtyan, F.-C. Tsai, J. Grime, P. Bassereau, G. A. Voth, "Organization of I-BAR Proteins on Tubular and Vesicular Membranes" (working paper)
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