Spatio-temporal coordination among functional residues in protein

The effect of phosphorylation on the conformational dynamics and allostery of the association of death-associated protein kinase with calmodulin

Protein phosphorylation plays an important role in the signal transduction and is capable of regulation of cell activity. The death-associated protein kinase 1 (DAPK1), as a Ser/Thr kinase, interacts with calmodulin (CaM) to regulate apoptotic and autophagic signaling. Autophosphorylation of DAPK1 at Ser308 located at the autoregulatory domain (ARD) blocks CaM binding and inhibits kinase catalytic activity. However, the mechanism underlying the influence of Ser308 phosphorylation (pS308) on the DAPK1 activity remains unclear. Here, we performed multiple, microsecond length molecular dynamics (MD) simulations, the molecular mechanics generalized Born/surface area (MM-GBSA) binding free energy calculations, principal component analysis, and dynamic cross-correlation analysis to unravel the conformational dynamics and allostery of the DAPK1 - CaM interaction triggered by the pS308 at the ARD. MD simulations showed that pS308 affected the conformational stability of the DAPK1 - CaM complex. Further energetic and structural exploration revealed that pS308 weakened the association of the phosphorylated DAPK1 to CaM, which lowered the susceptibility of DAPK1 to be activated by CaM. This result can provide mechanistic insights into the molecular underpinning through which the DAPK1 kinase activity is modulated by the auto-phosphorylation.Communicated by Ramaswamy H. Sarma.

Uncovering allostery and regulation in SORCIN through molecular dynamics simulations

Soluble resistance-related calcium-binding protein or Sorcin is an allosteric, calcium-binding Penta-EF hand (PEF) family protein implicated in multi-drug resistant cancers. Sorcin is known to bind chemotherapeutic molecules such as Doxorubicin. This study uses in-silico molecular dynamics simulations to explore the dynamics and allosteric behavior of Sorcin in the context of Ca2+ uptake and Doxorubicin binding. The results show that Ca2+ binding induces large, but reversible conformational changes in the Sorcin structure which manifest as rigid body reorientations that preserve the local secondary structure. A reciprocal allosteric handshake centered around the EF5 hand is found to be key in Sorcin dimer formation and stabilization. Binding of Doxorubicin results in rearrangement of allosteric communities which disrupts long-range allosteric information transfer from the N-terminal domain to the middle lobe. However, this binding does not result in secondary structure destabilization. Sorcin does not appear to have a distinct Ca2+ activated mode of Doxorubicin binding.Communicated by Ramaswamy H. Sarma.

Kinetic coherence underlies the dynamics of disordered proteins

The dynamics of two proteins of similar size, the globular lysozyme and the intrinsically disordered Huntingtin interacting protein, has been simulated in three states resembling a globule, a pre-molten globule, and a molten globule. A coherence time τ has been defined, measuring the delay in the display of a stochastic behaviour after a perturbation of the system. This time has been computed for two sets of collective variables: the projection of the phase point onto the positions and momenta subspaces (τ_r and τ_p), and the principal components (PCs) of positions q and momenta π produced by a covariance analysis in these subspaces (τ_q and τ_π). In all states τ_p ≈ 3.5τ_r, and τ_π ≈ 3.5τ_q. The coherence times of individual PCs, τ^(l)_q and τ^(l)_π, have also been computed, and τ^(l)_π > τ^(l)_q in all states. The prevalence of τ_p over τ_r, or of τ_π over τ_q, drives the dynamics of the protein over a time range of ≈1–2 ps; moreover, a hidden synchronism appears to raise the momenta subspace's coherence above that of its individual PCs. In the transition of lysozyme to the molten globule the τ^(l)_q decrease but, unexpectedly, the τ^(l)_π increase; after this transition τ_p ≈ 5τ_r and τ_π ≈ 5τ_q. A gain of kinetic coherence accompanies thus the loss of structural coherence caused by the denaturation of the protein in the transition from globule to molten globule. The increase of the τ^(l)_π does not take place in the analogous transition of the Huntingtin protein. These results are compared with those of a similar analysis performed on three pseudo-proteins designed by scrambling the primary sequence of the Huntingtin interacting protein, and on two oligopeptides. The hidden synchronism appears to be a generic property of these polypeptides. The τ^(l)_π spectrum is similar in denaturated and in intrinsically disordered biomolecules; but the gain of kinetic coherence as a result of denaturation seems to be a specific property of the biologically functional lysozyme.

In the phase space of a globular or intrinsically disordered protein, the momenta's dynamics is less chaotic than the coordinates' dynamics. When a protein is denaturated, a gain in kinetic coherence accompanies the loss of structural coherence.

Microscopic insight to specificity of metal ion cofactor in DNA cleavage by restriction endonuclease EcoRV

Restriction endonucleases protect bacterial cells against bacteriophage infection by cleaving the incoming foreign DNA into fragments. In presence of Mg²⁺ ions, EcoRV is able to cleave the DNA but not in presence of Ca²⁺ , although the protein binds to DNA in presence of both metal ions. We make an attempt to understand this difference using conformational thermodynamics. We calculate the changes in conformational free energy and entropy of conformational degrees of freedom, like DNA base pair steps and dihedral angles of protein residues in Mg²⁺ (A)-EcoRV-DNA complex compared to Ca²⁺ (S)-EcoRV-DNA complex using all-atom molecular dynamics (MD) trajectories of the complexes. We find that despite conformational stability and order in both complexes, the individual degrees of freedom behave differently in the presence of two different metal ions. The base pairs in cleavage region are highly disordered in Ca²⁺ (S)-EcoRV-DNA compared to Mg²⁺ (A)-EcoRV-DNA. One of the acidic residues ASP90, coordinating to the metal ion in the vicinity of the cleavage site, is conformationally destabilized and disordered, while basic residue LYS92 gets conformational stability and order in Ca²⁺ (S) bound complex than in Mg²⁺ (A) bound complex. The enhanced fluctuations hinder placement of the metal ion in the vicinity of the scissile phosphate of DNA. Similar loss of conformational stability and order in the cleavage region is observed by the replacement of the metal ion. Considering the placement of the metal ion near scissile phosphate as requirement for cleavage action, our results suggest that the changes in conformational stability and order of the base pair steps and the protein residues lead to cofactor sensitivity of the enzyme. Our method based on fluctuations of microscopic conformational variables can be applied to understand enzyme activities in other protein-DNA systems.

In-silico studies on conformational stability of flagellin-receptor complexes

Flagellin is a protein, responsible for virulent activities of bacteria. The host cell surface receptor protein TLR5 is known to interact with flagellin in order to activate immune response. However, the underlying microscopic details of this immune response are still elusive. In this study, we report on conformational stability of flagellin of two different organisms known as fliC and flaD in bilayer with reference to water. We find that both the flagellin is conformationally more stable in bilayer than in water. We also observe that fliC-TLR5 and flaD-TLR5 complexes are conformationally stable when the extracellular domain of the protein binds to conserved D1 domain of both fliC and flaD, although the binding interface between fliC-TLR5 and flaD-TLR5 is not identical. Our studies suggest that this might lead to differences in coreceptor bindings involved in immune response and thus have potential application in pharmaceutical developments. AbbreviationsA2Aadenosine receptorDPPCdipalmitoyl phosphatidylcholineecdextracellular domainecl2extracellular loop 2eLRRextracellular Leucine rich repeat domainflaDflagellin of Vibrio choleraefliCflagellin of Salmonella typhimuriumHPVhyper-variableMDmolecular dynamicsRMSDroot means squared deviationTIRtoll-interleukin receptorTLR5toll like receptor 5VPAC1vasoactive intestinal peptide receptorCommunicated by Ramaswamy H. Sarma.

Allostery in Orai1 binding to calmodulin revealed from conformational thermodynamics

Here, we study microscopic mechanism of complex formation between Ca²⁺-bound calmodulin (holoCaM) and Orai1 that regulates Ca²⁺-dependent inactivation process in eukaryotic cells. We compute conformational thermodynamic changes in holoCaM with respect to complex of Orai1 bound to C-terminal domain of holoCaM using histograms of dihedral angles of the proteins over trajectories from molecular dynamics simulations. Our analysis shows that the N-terminal domain residues L4, T5, Q41, N42, T44 and E67 of holoCaM get destabilized and disordered due to Orai1 binding to C-terminal domain of calmodulin affect the N-terminal domain residues. Among these residues, polar T44, having maximum destabilization and disorder via backbone fluctuations, shows the largest change in solvent exposure. This suggests that N-terminal domain is allosterically regulated via T44 by the binding of Orai1 to the C-terminal domain.