Quadratic Optimization-Based Nonlinear Control for Protein Conformation Prediction

Predicting the final folded structure of protein molecules and simulating their folding pathways is of crucial importance for designing viral drugs and studying diseases such as Alzheimer’s at the molecular level. To this end, this letter investigates the problem of protein conformation prediction under the constraint of avoiding high-entropy-loss routes during folding. Using the well-established kinetostatic compliance (KCM)-based nonlinear dynamics of a protein molecule, this letter formulates the protein conformation prediction as a pointwise optimal control synthesis problem cast as a quadratic program (QP). It is shown that the KCM torques in the protein folding literature can be utilized for defining a reference vector field for the QP-based control generation problem. The resulting kinetostatic control torque inputs will be close to the KCM-based reference vector field and guaranteed to be constrained by a predetermined bound; hence, high-entropy-loss routes during folding are avoided while the energy of the molecule is decreased.

[1]  Constantinos Mavroidis,et al.  Analysis and Design of Protein Based Nanodevices: Challenges and Opportunities in Mechanical Design , 2005 .

[2]  Mark W. Spong,et al.  The control of robot manipulators with bounded input , 1986 .

[3]  Horea T. Ilies,et al.  Residue Level Three-dimensional Workspace Maps for Conformational Trajectory Planning of Proteins , 2009, Int. J. Robotics Res..

[4]  H. Berendsen,et al.  Essential dynamics of proteins , 1993, Proteins.

[5]  Mohammed AlQuraishi A watershed moment for protein structure prediction , 2020, Nature.

[6]  Anna Tramontano,et al.  Critical assessment of methods of protein structure prediction (CASP) — round x , 2014, Proteins.

[7]  Stephen P. Boyd,et al.  Fast Evaluation of Quadratic Control-Lyapunov Policy , 2011, IEEE Transactions on Control Systems Technology.

[8]  Dianne P. O'Leary,et al.  HOPE: A Homotopy Optimization Method for Protein Structure Prediction , 2005, J. Comput. Biol..

[9]  Toshiyuki Ohtsuka,et al.  Solutions to the Hamilton-Jacobi Equation With Algebraic Gradients , 2011, IEEE Transactions on Automatic Control.

[10]  Aaron D. Ames,et al.  Sufficient conditions for the Lipschitz continuity of QP-based multi-objective control of humanoid robots , 2013, 52nd IEEE Conference on Decision and Control.

[11]  R. Barnard An optimal-aim control strategy for nonlinear regulation systems , 1975 .

[12]  Erik I. Verriest,et al.  Graceful Transitions between Periodic Walking Gaits of Fully Actuated Bipedal Robots , 2020, 2020 American Control Conference (ACC).

[14]  M. Ikeguchi,et al.  Computational Methods for Configurational Entropy Using Internal and Cartesian Coordinates. , 2016, Journal of chemical theory and computation.

[15]  Kazem Kazerounian,et al.  Hydrogen Bonds and Kinematic Mobility of Protein Molecules , 2010 .

[16]  H. Rabitz,et al.  Interdiction of Protein Folding for Therapeutic Drug Development in SARS CoV-2 , 2020, The journal of physical chemistry. B.

[17]  Khalid Latif,et al.  Nano-Kinematics for Analysis Of Protein Molecules , 2005 .

[18]  Morad Behandish,et al.  Protofold II: Enhanced Model and Implementation for Kinetostatic Protein Folding , 2015, ArXiv.

[19]  Frank L. Lewis,et al.  Fixed-Final-Time-Constrained Optimal Control of Nonlinear Systems Using Neural Network HJB Approach , 2007, IEEE Transactions on Neural Networks.

[20]  Michael P Eastwood,et al.  A common, avoidable source of error in molecular dynamics integrators. , 2007, The Journal of chemical physics.

[21]  Alexei V. Finkelstein,et al.  Protein Physics: A Course of Lectures , 2002 .

[22]  Paulo Tabuada,et al.  Control Barrier Function Based Quadratic Programs for Safety Critical Systems , 2016, IEEE Transactions on Automatic Control.

[23]  Burak Erman,et al.  Prediction of Optimal Folding Routes of Proteins That Satisfy the Principle of Lowest Entropy Loss: Dynamic Contact Maps and Optimal Control , 2010, PloS one.

[24]  Combining Optimal Control Theory and Molecular Dynamics for Protein Folding , 2012, PloS one.

[25]  Kazem Kazerounian,et al.  On the rotational operators in protein structure simulations. , 2003, Protein engineering.

[26]  A. Iserles A First Course in the Numerical Analysis of Differential Equations: Gaussian elimination for sparse linear equations , 2008 .

[27]  Pouya Tavousi,et al.  On the Systematic Design and Analysis of Artificial Molecular Machines , 2016 .

[28]  Chang-Hee Won,et al.  Optimal control using an algebraic method for control-affine non-linear systems , 2007, Int. J. Control.

[29]  James S. Thorp,et al.  A model-referenced controller for stabilizing large transient swings in power systems , 1976 .

[30]  A. Liwo,et al.  Physics-based protein-structure prediction using a hierarchical protocol based on the UNRES force field: assessment in two blind tests. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[31]  M. Ediger,et al.  Brownian dynamics simulations of local motions in polyisoprene , 1991 .

[32]  Yanran Ding,et al.  qpSWIFT: A Real-Time Sparse Quadratic Program Solver for Robotic Applications , 2019, IEEE Robotics and Automation Letters.

[33]  K. Kazerounian,et al.  Protofold: A Successive Kinetostatic Compliance Method for Protein Conformation Prediction , 2005 .

[34]  Mark W. Spong,et al.  Integral line-of-sight path following control of magnetic helical microswimmers subject to step-out frequencies , 2021, Autom..

[35]  Roland Siegwart,et al.  LQR-Assisted Whole-Body Control of a Wheeled Bipedal Robot With Kinematic Loops , 2020, IEEE Robotics and Automation Letters.

[36]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[37]  Anna Tramontano,et al.  Critical assessment of methods of protein structure prediction—Round VII , 2007, Proteins.

[38]  O. Mangasarian,et al.  The Fritz John Necessary Optimality Conditions in the Presence of Equality and Inequality Constraints , 1967 .

[39]  K. Dill,et al.  Folding rates and low-entropy-loss routes of two-state proteins. , 2003, Journal of molecular biology.

[40]  A Caflisch,et al.  Computer simulations of protein folding by targeted molecular dynamics , 2000, Proteins.

[41]  Qiang Shao,et al.  Effects of turn stability and side-chain hydrophobicity on the folding of β-structures. , 2010, Journal of molecular biology.