This paper investigates the problem of prediction of protein molecule folding pathways under entropy-loss constraints by formulating a control synthesis problem whose solutions are obtained by solving large-scale quadratic programming (QP) optimizations with nonlinear constraints. The utilized non-iterative and computationally efficient algorithm, which is based on solving generalized eigenvalue problems, prevents an unpredictable and potentially large number of iterations at each protein conformation for computing the folding control inputs. The synthesized control inputs remain close to the renowned kinetostatic compliance method (KCM) reference vector field while satisfying proper quadratic inequality constraints that limit the rate of molecule entropy-loss during folding.
more »
« less
Prediction of Protein Folding Pathways under Entropy-Loss Constraints using Quadratic Programming-Based Nonlinear Control
This paper investigates the problem of prediction of protein molecule folding pathways under entropy-loss constraints by formulating a control synthesis problem whose solutions are obtained by solving large-scale quadratic programming (QP) optimizations with nonlinear constraints. The utilized non-iterative and computationally efficient algorithm, which is based on solving generalized eigenvalue problems, prevents an unpredictable and potentially large number of iterations at each protein conformation for computing the folding control inputs. The synthesized control inputs remain close to the renowned kinetostatic compliance method (KCM) reference vector field while satisfying proper quadratic inequality constraints that limit the rate of molecule entropy-loss during folding.
more »
« less
- Award ID(s):
- 2153744
- PAR ID:
- 10427744
- Publisher / Repository:
- HAL open science
- Date Published:
- Journal Name:
- 2023 American Control Conference (ACC)
- Format(s):
- Medium: X
- Location:
- 2023 American Control Conference (ACC), San Diego, CA
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract This paper addresses the problem of algorithmic prediction of protein folding pathways, namely, the transient three-dimensional conformations of protein molecules during folding, under constrained rates of entropy change. We formulate the physics-based prediction of folding pathways as a control synthesis problem, where the control inputs guide the protein folding simulations. These folding control inputs are obtained from largescale trust-region subproblems (TRS) utilizing a computationally efficient algorithm with no need for outer iterations. The proposed control synthesis approach, which leverages the solutions obtained from a special generalized eigenvalue problem, avoids potentially cumbersome and unpredictable iterative computations at each protein conformation. Moreover, the TRS-based control inputs align the closed-loop dynamics closely with the kinetostatic compliance method (KCM) reference vector field while satisfying ellipsoidal constraints on the folding control inputs. Finally, we provide conditions for existence and uniqueness of the resulting closed-loop solutions, which are the protein folding pathways under constraints on the rate of entropy change. Numerical simulations utilizing the KCM approach on protein backbones confirm the effectiveness of the proposed framework.more » « less
-
Understanding the process of protein unfolding plays a crucial role in various applications such as design of folding-based protein engines. Using the well-established kinetostatic compliance (KCM)-based method for modeling of protein conformation dynamics and a recent nonlinear control theoretic approach to KCM-based protein folding, this letter formulates protein unfolding as a destabilizing control analysis/synthesis problem. In light of this formulation, it is shown that the Chetaev instability framework can be used to investigate the KCM-based unfolding dynamics. In particular, a Chetaev function for analysis of unfolding dynamics under the effect of optical tweezers and a class of control Chetaev functions for synthesizing control inputs that elongate protein strands from their folded conformations are presented. Based on the presented control Chetaev function, an unfolding input is derived from the Artstein-Sontag universal formula and the results are compared against optical tweezer-based unfolding.more » « less
-
This paper proposes a sign gradient descent (SGD) algorithm for predicting the three-dimensional folded protein molecule structures under the kinetostatic compliance method (KCM). In the KCM framework, which can be used to simulate the range of motion of peptide-based nanorobots/nanomachines, protein molecules are modeled as a large number of rigid nano-linkages that form a kinematic mechanism under motion constraints imposed by chemical bonds while folding under the kinetostatic effect of nonlinear interatomic force fields. In a departure from the conventional successive kinetostatic fold compliance framework, the proposed SGD-based iterative algorithm in this paper results in convergence to the local minima of the free energy of protein molecules corresponding to their final folded conformations in a faster and more robust manner. KCM-based folding dynamics simulations of the backbone chains of protein molecules demonstrate the effectiveness of the proposed algorithm.more » « less
-
We consider a multi-agent linear quadratic optimal control problem. Due to communication constraints, the agents are required to quantize their local state measurements before communicating them to the rest of the team, thus resulting in a decentralized information structure. The optimal controllers are to be synthesized under this decentralized and quantized information structure. The agents are given a set of quantizers with varying quantization resolutions—higher resolution incurs higher communication cost and vice versa. The team must optimally select the quantizer to prioritize agents with ‘highquality’ information for optimizing the control performance under communication constraints. We show that there exist a sepatation between the optimal solution to the control problem and the choice of the optimal quantizer. We show that the optimal controllers are linear and the optimal selection of the quantizers can be determined by solving a linear program.more » « less
An official website of the United States government

