Molecular dynamics simulations of aptamer-binding reveal generalized allostery in thrombin

Thrombin is an attractive target for antithrombotic therapy due to its central role in thrombosis and hemostasis as well as its role in inducing tumor growth, metastasis, and tumor invasion. The thrombin-binding DNA aptamer (TBA), is under investigation for anticoagulant drugs. Although aptamer binding experiments have been revealed various effects on thrombin’s enzymatic activities, the detailed picture of the thrombin’s allostery from TBA binding is still unclear. To investigate thrombin’s response to the aptamer-binding at the molecular level, we compare the mechanical properties and free energy landscapes of the free and aptamer-bound thrombin using microsecond-scale all-atom GPU-based molecular dynamics simulations. Our calculations on residue fluctuations and coupling illustrate the allosteric effects of aptamer-binding at the atomic level, highlighting the exosite II, 60s, γ and the sodium loops, and the alpha helix region in the light chains involved in the allosteric changes. This level of details clarifies the mechanisms of previous experimentally demonstrated phenomena, and provides a prediction of the reduced autolysis rate after aptamer-binding. The shifts in thrombin’s ensemble of conformations and free energy surfaces after aptamer-binding demonstrate that the presence of bound-aptamer restricts the conformational freedom of thrombin suggesting that conformational selection, i.e. generalized allostery, is the dominant mechanism of thrombin-aptamer binding. The profound perturbation on thrombin’s mechanical and thermodynamic properties due to the aptamer-binding, which was revealed comprehensively as a generalized allostery in this work, may be exploited in further drug discovery and development.     

Harmonicity and anharmonicity in protein dynamics: A normal mode analysis and principal component analysis

A comparison of a normal mode analysis and principal component analysis of a 200-ps molecular dynamics trajectory of bovine pancreatic trypsin inhibitor in vacuum has been made in order to further elucidate the harmonic and anharmonic aspects in the dynamics of proteins. An anharmonicity factor is defined which measures the degree of anharmonicity in the modes, be they principal modes or normal modes, and it is shown that the principal mode system naturally divides into anharmonic modes with peak frequencies below 80 cm^?1, and harmonic modes with frequencies above this value. In general the larger the mean-square fluctuation of a principal mode, the greater the degree of anharmonicity in its motion. The anharmonic modes represent only 12% of the total number of variables, but account for 98% of the total mean-square fluctuation. The transitional nature of the anharmonic motion is demonstrated. The results strongly suggest that in a large subspace, the free energy surface, as probed by the simulation, is approximated by a multi-dimensional parabola which is just a resealed version of the parabola corresponding to the harmonic approximation to the conformational energy surface at a single minimum. After 200 ps, the resealing factor, termed the “normal mode resealing factor,” has apparently converged to a value whereby the mean-square fluctuation within the subspace is about twice that predicted by the normal mode analysis. © 1995 Wiley-Liss, Inc.     

Engineered Myosin VI Motors Reveal Minimal Structural Determinants of Directionality and Processivity

Myosins have diverse mechanical properties reflecting a range of cellular roles. A major challenge is to understand the structural basis for generating novel functions from a common motor core. Myosin VI (M6) is specialized for processive motion toward the (−) end of actin filaments. We have used engineered M6 motors to test and refine the “redirected power stroke” model for (−) end directionality and to explore poorly understood structural requirements for processive stepping. Guided by crystal structures and molecular modeling, we fused artificial lever arms to the catalytic head of M6 at several positions, retaining varying amounts of native structure. We found that an 18-residue α-helical insert is sufficient to reverse the directionality of the motor, with no requirement for any calmodulin light chains. Further, we observed robust processive stepping of motors with artificial lever arms, demonstrating that processivity can arise without optimizing lever arm composition or mechanics.     

Dynamic properties of some β-chain mutant hemoglobins

The thermal behavior of the Soret band relative to the carbonmonoxy derivatives of some β-chain mutant hemoglobins is studied in the temperature range 300–10 K and compared to that of wild-type carbonmonoxy hemoglobin. The band profile at various temperatures is modeled as a Voigt function that accounts for homogeneous broadening and for the coupling with high- and low-frequency vibrational modes, while inhomogeneous broadening is taken into account with a gaussian distribution of purely electronic transition frequencies. The various contributions to the overall bandwidth are singled out With this analysis and their temperature dependence, in turn, gives information on structural and dynamic properties of the system studied. In the wildtype and mutant hemoglobins, the values of homogeneous bandwidth and of the coupling constants to high-frequency vibrational modes are not modified with respect to natural human hemoglobin, thus indicating that the local electronic and vibrational properties of the heme–CO complex are not altered by the recombinant procedures. On the contrary, differences in the protein dynamic behavior are observed. The most relevant are those relative to the “polar isosteric” βVal-67(Ell) →Thr substitution, localized in the heme pocket, which results in decreased coupling with low-frequency modes and increased anharmonic motions. Mutations involving residue βLys-144(HC1) at the C-terminal and residue βCys-112(G14) at the α₁β₁ interface have a smaller effect consisting in an increased coupling with low-frequency modes. Mutations at the β-N-terminal and at the α₁β₂ interface have no effect on the dynamic properties of the heme pocket. © 1995 Wiley-Liss, Inc.     

Probing light chain mutation effects on thrombin via molecular dynamics simulations and machine learning

Thrombin is a key component for chemotherapeutic and antithrombotic therapy development. As the physiologic and pathologic roles of the light chain still remain vague, here, we continue previous efforts to understand the impacts of the disease-associated single deletion of LYS9 in the light chain. By combining supervised and unsupervised machine learning methodologies and more traditional structural analyses on data from 10 μs molecular dynamics simulations, we show that the conformational ensemble of the ΔK9 mutant is significantly perturbed. Our analyses consistently indicate that LYS9 deletion destabilizes both the catalytic cleft and regulatory functional regions and result in some conformational changes that occur in tens to hundreds of nanosecond scaled motions. We also reveal that the two forms of thrombin each prefer a distinct binding mode of a Na⁺ ion. We expand our understanding of previous experimental observations and shed light on the mechanisms of the LYS9 deletion associated bleeding disorder by providing consistent but more quantitative and detailed structural analyses than early studies in literature. With a novel application of supervised learning, i.e. the decision tree learning on the hydrogen bonding features in the wild-type and ΔK9 mutant forms of thrombin, we predict that seven pairs of critical hydrogen bonding interactions are significant for establishing distinct behaviors of wild-type thrombin and its ΔK9 mutant form. Our calculations indicate the LYS9 in the light chain has both localized and long-range allosteric effects on thrombin, supporting the opinion that light chain has an important role as an allosteric effector.