<p>Protein conformational transitions are fundamental to the functions of many proteins, and computational methods are valuable for elucidating the transitions that are not readily accessible by experimental techniques. Here we developed accelerated sampling methods to calculate optimized all-atom protein conformational transition paths. Adaptively biased path optimization (ABPO) is a computational simulation method to optimize the conformational transition path between two states. We first examined the transition paths of three systems with relatively simple transitions. The ways to define reduced variables were explored and transition paths were built at convergence of the optimizations. We constructed the all-atom conformational transition path between the active and the inactive states of the Src kinase domain. The C helix rotation was identified as the main free energy barrier in the all‑atom system, and the intermediate conformations and key interactions along the transition path were analyzed. This is the first demonstration of the robustness of a computational method for calculating protein conformational transitions without restraints to a specified path. We also evaluated protein‑peptide interactions using both molecular dynamics simulations and peptide docking. Long unbiased simulations were used to evaluate Src‑SSP interactions and complex stability in both implicit and explicit solvent. Molecular docking was used to build possible protein‑peptide interaction models, using both Src regulatory domain SH2 and the kinase domain. Possible Src‑SSP complexes were built as the first Src‑substrate complex structure models.</p>

No description available.

abstract: In a typical living cell, millions to billions of proteins—nanomachines that fluctuate and cycle among many conformational states—convert available free energy into mechanochemical work. A fundamental goal of biophysics is to ascertain how 3D protein structures encode specific functions, such as catalyzing chemical reactions or transporting nutrients into a cell. Protein dynamics span femtosecond timescales (i.e., covalent bond oscillations) to large conformational transition timescales in, and beyond, the millisecond regime (e.g., glucose transport across a phospholipid bilayer). Actual transition events are fast but rare, occurring orders of magnitude faster than typical metastable equilibrium waiting times. Equilibrium molecular dynamics (EqMD) can capture atomistic detail and solute-solvent interactions, but even microseconds of sampling attainable nowadays still falls orders of magnitude short of transition timescales, especially for large systems, rendering observations of such "rare events" difficult or effectively impossible. Advanced path-sampling methods exploit reduced physical models or biasing to produce plausible transitions while balancing accuracy and efficiency, but quantifying their accuracy relative to other numerical and experimental data has been challenging. Indeed, new horizons in elucidating protein function necessitate that present methodologies be revised to more seamlessly and quantitatively integrate a spectrum of methods, both numerical and experimental. In this dissertation, experimental and computational methods are put into perspective using the enzyme adenylate kinase (AdK) as an illustrative example. We introduce Path Similarity Analysis (PSA)—an integrative computational framework developed to quantify transition path similarity. PSA not only reliably distinguished AdK transitions by the originating method, but also traced pathway differences between two methods back to charge-charge interactions (neglected by the stereochemical model, but not the all-atom force field) in several conserved salt bridges. Cryo-electron microscopy maps of the transporter Bor1p are directly incorporated into EqMD simulations using MD flexible fitting to produce viable structural models and infer a plausible transport mechanism. Conforming to the theme of integration, a short compendium of an exploratory project—developing a hybrid atomistic-continuum method—is presented, including initial results and a novel fluctuating hydrodynamics model and corresponding numerical code. / Dissertation/Thesis / Doctoral Dissertation Physics 2017