Mechanical models of proteins

textIn general, this dissertation is concerned with modeling of mechanical behavior of protein molecules. In particular, we focus on coarse-grained models, which bridge the gap in time and length scale between the atomistic simulation and biological processes. The dissertation presents three independent studies involving such models. The first study is concerned with a rigorous coarse-graining method for dynamics of linear systems. In this method, as usual, the conformational space of the original atomistic system is divided into master and slave degrees of freedom. Under the assumption that the characteristic timescales of the masters are slower than those of the slaves, the method results in Langevin-type equations of motion governed by an effective potential of mean force. In addition, coarse-graining introduces hydrodynamic-like coupling among the masters as well as non-trivial inertial effects. Application of our method to the long-timescale part of the relaxation spectra of proteins shows that such dynamic coupling is essential for reproducing their relaxation rates and modes. The second study is concerned with calibration of elastic network models based on the so-called B-factors, obtained from x-ray crystallographic measurements. We show that a proper calibration procedure must account for rigid-body motion and constraints imposed by the crystalline environment on the protein. These fundamental aspects of protein dynamics in crystals are often ignored in currently used elastic network models, leading to potentially erroneous network parameters. We develop an elastic network model that properly takes rigid-body motion and crystalline constraints into account. This model reveals that B-factors are dominated by rigid-body motion rather than deformation, and therefore B-factors are poorly suited for identifying elastic properties of protein molecules. Furthermore, it turns out that B-factors for a benchmark set of three hundred and thirty protein molecules can be well approximated by assuming that the protein molecules are rigid. The third study is concerned with the polymer mediated interaction between two planar surfaces. In particular, we consider the case where a thin polymer layer bridges two parallel plates. We consider two models of monodisperse and polydisperse for the polymer layer and obtain an analytical expression for the force-distance relationship of the two plates.Engineering Mechanic

The Role of Mutations in Protein Structural Dynamics and Function: A Multi-scale Computational Approach

abstract: Proteins are a fundamental unit in biology. Although proteins have been extensively studied, there is still much to investigate. The mechanism by which proteins fold into their native state, how evolution shapes structural dynamics, and the dynamic mechanisms of many diseases are not well understood. In this thesis, protein folding is explored using a multi-scale modeling method including (i) geometric constraint based simulations that efficiently search for native like topologies and (ii) reservoir replica exchange molecular dynamics, which identify the low free energy structures and refines these structures toward the native conformation. A test set of eight proteins and three ancestral steroid receptor proteins are folded to 2.7Å all-atom RMSD from their experimental crystal structures. Protein evolution and disease associated mutations (DAMs) are most commonly studied by in silico multiple sequence alignment methods. Here, however, the structural dynamics are incorporated to give insight into the evolution of three ancestral proteins and the mechanism of several diseases in human ferritin protein. The differences in conformational dynamics of these evolutionary related, functionally diverged ancestral steroid receptor proteins are investigated by obtaining the most collective motion through essential dynamics. Strikingly, this analysis shows that evolutionary diverged proteins of the same family do not share the same dynamic subspace. Rather, those sharing the same function are simultaneously clustered together and distant from those functionally diverged homologs. This dynamics analysis also identifies 77% of mutations (functional and permissive) necessary to evolve new function. In silico methods for prediction of DAMs rely on differences in evolution rate due to purifying selection and therefore the accuracy of DAM prediction decreases at fast and slow evolvable sites. Here, we investigate structural dynamics through computing the contribution of each residue to the biologically relevant fluctuations and from this define a metric: the dynamic stability index (DSI). Using DSI we study the mechanism for three diseases observed in the human ferritin protein. The T30I and R40G DAMs show a loss of dynamic stability at the C-terminus helix and nearby regulatory loop, agreeing with experimental results implicating the same regulatory loop as a cause in cataracts syndrome.Dissertation/ThesisPh.D. Physics 201

Without magic bullets: the biological basis for public health interventions against protein folding disorders

Protein folding disorders of aging like Alzheimer's and Parkinson's diseases currently present intractable medical challenges. 'Small molecule' interventions - drug treatments - often have, at best, palliative impact, failing to alter disease course. The design of individual or population level interventions will likely require a deeper understanding of protein folding and its regulation than currently provided by contemporary 'physics' or culture-bound medical magic bullet models. Here, a topological rate distortion analysis is applied to the problem of protein folding and regulation that is similar in spirit to Tlusty's (2010a) elegant exploration of the genetic code. The formalism produces large-scale, quasi-equilibrium 'resilience' states representing normal and pathological protein folding regulation under a cellular-level cognitive paradigm similar to that proposed by Atlan and Cohen (1998) for the immune system. Generalization to long times produces diffusion models of protein folding disorders in which epigenetic or life history factors determine the rate of onset of regulatory failure, in essence, a premature aging driven by familiar synergisms between disjunctions of resource allocation and need in the context of socially or physiologically toxic exposures and chronic powerlessness at individual and group scales. Application of an HPA axis model is made to recent observed differences in Alzheimer's onset rates in White and African American subpopulations as a function of an index of distress-proneness

Single-domain protein folding: a multi-faceted problem

We review theoretical approaches, experiments and numerical simulations that have been recently proposed to investigate the folding problem in single-domain proteins. From a theoretical point of view, we emphasize the energy landscape approach. As far as experiments are concerned, we focus on the recent development of single-molecule techniques. In particular, we compare the results obtained with two main techniques: single protein force measurements with optical tweezers and single-molecule fluorescence in studies on the same protein (RNase H). This allows us to point out some controversial issues such as the nature of the denatured and intermediate states and possible folding pathways. After reviewing the various numerical simulation techniques, we show that on-lattice protein-like models can help to understand many controversial issues.Comment: 26 pages, AIP Conference Proceeding

Biopolymer Unfolding as a Process of Biased Activated Barrier Crossing

The biopolymer unfolding process can be conceptualized as the thermally activated crossing of a barrier in some energy landscape; the real-world unfolding dynamics are in correspondence with the rate of transit over the barrier. In this dissertation, we report on four projects that exploit this analogy to advance our understanding of so-called pulling experiments, in which a biopolymer chain is encouraged to unfold by direct application of forces to its two ends. Regardless of the specific experimental realization, standard pulling techniques give rise to a biased barrier crossing problem. The applied force has the effect of tilting the energy landscape. The height and shape of the barrier are altered by the tilt, and the rate of transit over the barrier is highly sensitive to such changes

Contextualized genome-scale model unveils high-order metabolic effects of the specific growth rate and oxygenation level in recombinant Pichia pastoris

Pichia pastoris is recognized as a biotechnological workhorse for recombinant protein expression. The metabolic performance of this microorganism depends on genetic makeup and culture conditions, amongst which the specific growth rate and oxygenation level are critical. Despite their importance, only their individual effects have been assessed so far, and thus their combined effects and metabolic consequences still remain to be elucidated. In this work, we present a comprehensive framework for revealing high-order (i.e., individual and combined) metabolic effects of the above parameters in glucose-limited continuous cultures of P. pastoris, using thaumatin production as a case study. Specifically, we employed a rational experimental design to calculate statistically significant metabolic effects from multiple chemostat data, which were later contextualized using a refined and highly predictive genome-scale metabolic model of this yeast under the simulated conditions. Our results revealed a negative effect of the oxygenation on the specific product formation rate (thaumatin), and a positive effect on the biomass yield. Notably, we identified a novel positive combined effect of both the specific growth rate and oxygenation level on the specific product formation rate. Finally, model predictions indicated an opposite relationship between the oxygenation level and the growth-associated maintenance energy (GAME) requirement, suggesting a linear GAME decrease of 0.56 mmol ATP/g per each 1% increase in oxygenation level, which translated into a 44% higher metabolic cost under low oxygenation compared to high oxygenation. Overall, this work provides a systematic framework for mapping high-order metabolic effects of different culture parameters on the performance of a microbial cell factory. Particularly in this case, it provided valuable insights about optimal operational conditions for protein production in P. pastoris

Optimization of Enzymatic Production of Oligopeptides from Apricot Almonds Meal with Neutrase and N120P

Neutrase 0.8L and N120P proteases were used for oligopeptide production from apricot almonds meal, and response surface design was carried out to optimize the effect of hydrolysis conditions on hydrolysis degree (DH) and oligopeptide yield rate. Four independent variables were used to optimize the hydrolysis process: hydrolysis temperature (X1), enzyme-to substrate ratio (E/S) (X2), substrate concentration (X3) and reaction time (X4). Statistical analysis indicated that the four variables, quadratic terms of X1, X3, and X4, and the interaction terms with X1 had a significant (p < 0.05) effect on DH. The yield rate was also significantly affected by the four variables and quadratic terms of X1, X2 and X4. Two mathematical models with high determination coefficient were obtained and could be employed to optimize protein hydrolysis. The optimal hydrolysis conditions were determined as follows: hydrolysis temperature 52.5 °C; enzyme-to-substrate ratio (E/S) 7200 U/g; substrate concentration 2%; reaction time 173 min. The initial pH 6.5 and Neutrase-to-N120P dosage ratio 2:1 were fixed in this study according to the preliminary research. Under these conditions, the experimental DH and yield rate were 34.10 ± 5.25% and 72.42 ± 2.27%, respectively

PH-dependent aggregation in intrinsically disordered proteins is determined by charge and lipophilicity

Protein aggregation is associated with an increasing number of human disorders and premature aging. Moreover, it is a central concern in the manufacturing of recombinant proteins for biotechnological and therapeutic applications. Nevertheless, the unique architecture of protein aggregates is also exploited by nature for functional purposes, from bacteria to humans. The relevance of this process in health and disease has boosted the interest in understanding and controlling aggregation, with the concomitant development of a myriad of algorithms aimed to predict aggregation propensities. However, most of these programs are blind to the protein environment and, in particular, to the influence of the pH. Here, we developed an empirical equation to model the pH-dependent aggregation of intrinsically disordered proteins (IDPs) based on the assumption that both the global protein charge and lipophilicity depend on the solution pH. Upon its parametrization with a model IDP, this simple phenomenological approach showed unprecedented accuracy in predicting the dependence of the aggregation of both pathogenic and functional amyloidogenic IDPs on the pH. The algorithm might be useful for diverse applications, from large-scale analysis of IDPs aggregation properties to the design of novel reversible nanofibrillar materials