Structure and Dynamics Study of LeuT Using the Markov State Model and Perturbation Response Scanning Reveals Distinct Ion Induced Conformational States

Conformational transitions of caspase-6 in substrate-induced activation process explored by perturbation-response scanning combined with targeted molecular dynamics

Caspase-6 participates in a series of neurodegenerative pathways, and has aroused widespread attentions as a promising molecular target for the treatment of neurodegeneration. Caspase-6 is a homodimer with 6 central-stranded β-sheets and 5 α-helices in each monomer. Previous crystallographic studies suggested that the 60′s, 90′s and 130′s helices of caspase-6 undergo a distinctive conformational transition upon substrate binding. Although the caspase-6 structures in apo and active states have been determined, the conformational transition process between the two states remains poorly understood. In this work, perturbation-response scanning (PRS) combined with targeted molecular dynamics (TMD) simulations was employed to unravel the atomistic mechanism of the dynamic conformational transitions underlying the substrate-induced activation process of caspase-6. The results showed that the conformational transition of caspase-6 from apo to active states is mainly characterized by structural rearrangements of the substrate-binding site as well as the conformational changes of 60′s and 130′s extended helices. The H-bond interactions between L1, 130′s helix and 90′s helix are proved to be key determinant factors for substrate-induced conformational transition. These findings provide valuable insights into the activation mechanism of caspase-6 as well as the molecular design of caspase-6 inhibitors.

Synchronous and Asynchronous Response in Dynamically Perturbed Proteins

We present a dynamic perturbation-response model of proteins based on the Gaussian Network Model, where a residue is perturbed periodically, and the dynamic response of other residues is determined. The model shows that periodic perturbation causes a synchronous response in phase with the perturbation and an asynchronous response that is out of phase. The asynchronous component results from the viscous effects of the solvent and other dispersive factors in the system. The model is based on the solution of the Langevin equation in the presence of solvent, noise, and perturbation. We introduce several novel ideas: The concept of storage and loss compliance of the protein and their dependence on structure and frequency; the amount of work lost and the residues that contribute significantly to the lost work; new dynamic correlations that result from perturbation; causality, that is, the response of j when i is perturbed is not equal to the response of i when j is perturbed. As examples, we study two systems, namely, bovine rhodopsin and the class of nanobodies. The general results obtained are (i) synchronous and asynchronous correlations depend strongly on the frequency of perturbation, their magnitude decreases with increasing frequency, (ii) time-delayed mean-squared fluctuations of residues have only synchronous components. Asynchronicity is present only in cross correlations, that is, correlations between different residues, (iii) perturbation of loop residues leads to a large dissipation of work, (iv) correlations satisfy the hypothesis of pre-existing pathways according to which information transfer by perturbation rides on already existing equilibrium correlations in the system, (v) dynamic perturbation can introduce a selective response in the system, where the perturbation of each residue excites different sets of responding residues, and (vi) it is possible to identify nondissipative residues whose perturbation does not lead to dissipation in the protein. Despite its simplicity, the model explains several features of allosteric manipulation.

Computational approaches for deciphering the equilibrium and kinetic properties of iron transport proteins

With the advances in three-dimensional structure determination techniques, high quality structures of the iron transport proteins transferrin and the bacterial ferric binding protein (FbpA) have been deposited in the past decade. These are proteins of relatively large size, and developments in hardware and software have only recently made it possible to study their dynamics using standard computational resources. We review computational techniques towards understanding the equilibrium and kinetic properties of iron transport proteins under different environmental conditions. At the level of detail that requires quantum chemical treatments, the octahedral geometry around iron has been scrutinized and it has been established that the iron coordinating tyrosines are in an unusual deprotonated state. At the atomistic level, both the N-lobe and the full bilobal structure of transferrin have been studied under varying conditions of pH, ionic strength and binding of other metal ions by molecular dynamics (MD) simulations. These studies have allowed questions to be answered, among others, on the function of second shell residues in iron release, the role of synergistic anions in preparing the active site for iron binding, and the differences between the kinetics of the N- and the C-lobe. MD simulations on FbpA have led to the detailed observation of the binding kinetics of phosphate to the apo form, and to the conformational preferences of the holo form under conditions mimicking the environmental niches provided by the periplasmic space. To study the dynamics of these proteins with their receptors, one must resort to coarse-grained methodologies, since these systems are prohibitively large for atomistic simulations. A study of the complex of human transferrin (hTf) with its pathogenic receptor by such methods has revealed a potential mechanistic explanation for the defense mechanism that arises in evolutionary warfare. Meanwhile, the motions in the transferrin receptor bound hTf have been shown to disfavor apo hTf dissociation, explaining why the two proteins remain in complex during the recycling process from the endosome to the cell surface. Open problems and possible technological applications related to metal ion binding-release in iron transport proteins that may be handled by hybrid use of quantum mechanical, MD and coarse-grained approaches are discussed.