Interacting mechanism of ID3 HLH domain towards E2A/E12 transcription factor - An Insight through molecular dynamics and docking approach

Dynamozones are the most obvious sign of the evolution of conformational dynamics in HIV-1 protease

Proteins are not static but are flexible molecules that can adopt many different conformations. The HIV-1 protease is an important target for the development of therapies to treat AIDS, due to its critical role in the viral life cycle. We investigated several dynamics studies on the HIV-1 protease families to illustrate the significance of examining the dynamic behaviors and molecular motions for an entire understanding of their dynamics-structure-function relationships. Using computer simulations and principal component analysis approaches, the dynamics data obtained revealed that: (i) The flap regions are the most obvious sign of the evolution of conformational dynamics in HIV-1 protease; (ii) There are dynamic structural regions in some proteins that contribute to the biological function and allostery of proteins via appropriate flexibility. These regions are a clear sign of the evolution of conformational dynamics of proteins, which we call dynamozones. The flap regions are one of the most important dynamozones members that are critical for HIV-1 protease function. Due to the existence of other members of dynamozones in different proteins, we propose to consider dynamozones as a footprint of the evolution of the conformational dynamics of proteins.

Phylogenetic, structural, functional characterisation and effect of exogenous spermidine on rice (Oryza sativa) HAK transporters under salt stress

The HAK (High-affinity K+ ) family members mediate K+ transport that confers normal plant growth and resistance against unfavourable environmental conditions. Rice (Oryza sativa L.) HAK transporters have been extensively investigated for phylogenetic analyses with other plants species with very few of them functionally characterised. But very little information is known about their evolutionary aspects, overall structural, functional characterisation, and global expression pattern of the complete HAK family members in response to salt stress. In this study, 27 rice transporters were phylogenetically clustered with different dicot and monocot family members. Subsequently, the exon-intron structural patterns, conserved motif analyses, evolutionary divergence based different substitution matrix, orthologous-paralogous relationships were studied elaborately. Structural characterisations included a comparative study of secondary and tertiary structure, post-translational modifications, correspondence analyses, normal mode analyses, K+ /Na+ binding affinities of each of the OsHAK gene members. Global expression profile under salt stress showed clade-specific expression pattern of the proteins. Additionally, five OsHAK genes were chosen for further expression analyses in root and shoot tissues of two rice varieties during short-term salinity in the presence and absence of exogenous spermidine. All the information can be used as first-hand data for dissecting the administrative role of rice HAK transporters under various abiotic stresses.

Crocin Inhibits the Fibrillation of Human α-synuclein and Disassembles Mature Fibrils: Experimental Findings and Mechanistic Insights from Molecular Dynamics Simulation

The aggregation of human alpha-synuclein (hαS) is pivotally implicated in the development of most types of synucleinopathies. Molecules that can inhibit or reverse the aggregation process of amyloidogenic proteins have potential therapeutic value. The anti-aggregating activity of multiple carotenoid compounds has been reported over the past decades against a growing list of amyloidogenic polypeptides. Here, we aimed to determine whether crocin, the main carotenoid glycoside component of saffron, would inhibit hαS aggregation or could disassemble its preformed fibrils. By employing a series of biochemical and biophysical techniques, crocin was exhibited to inhibit hαS fibrillation in a dose-dependent fashion by stabilizing very early aggregation intermediates in off-pathway non-toxic conformations with little β-sheet content. We also observed that crocin at high concentrations could efficiently destabilize mature fibrils and disassemble them into seeding-incompetent intermediates by altering their β-sheet conformation and reshaping their structure. Our atomistic molecular dynamics (MD) simulations demonstrated that crocin molecules bind to both the non amyloid-β component (NAC) region and C-terminal domain of hαS. These interactions could thereby stabilize the autoinhibitory conformation of the protein and prevent it from adopting aggregation-prone structures. MD simulations further suggested that ligand molecules prefer to reside longitudinally along the fibril axis onto the edges of the inter-protofilament interface where they establish hydrogen and hydrophobic bonds with steric zipper stabilizing residues. These interactions turned out to destabilize hαS fibrils by altering the interstrand twist angles, increasing the rigidity of the fibril core, and elevating its radius of gyration. Our findings suggest the potential pharmaceutical implication of crocin in synucleinopathies.

Concerted motion of structure and active site charge is required for tyrosine aminotransferase activity in Leishmania parasite

Leishmania donovani tyrosine aminotransferase (LdTAT) is an essential enzyme that catalyzes the first step of amino acid catabolism. To understand LdTAT activity at different pH, molecular dynamics simulations were performed and trajectory and T-pad analysis pad were conducted. Fluorescence spectroscopy of LdTAT at various pH was measured to understand structural stability. UV studies on PLP were performed to determine the binding of the enzyme to cofactor PLP at different pH. The MD simulations showed that the structure of LdTAT was stable and no structural denaturation was observed at pH 2, 7 and 12. LdTAT exhibited the highest activity at pH -8 and fluorescent spectroscopy also corroborated by exhibiting the highest intensity at pH -8. Moreover, no structural denaturation was observed during the pH gradient. UV studies concluded that the aldimine bond forms only around neutral pH and redshift was observed on enzyme binding. From our observation, we hypothesize that the activity of LdTAT is a close interplay between the structure and charges of K286 and PLP. This study may provide significant insight into understanding parasitic enzymes like LdTAT during the life-cycle of Leishmania parasite. Knowledge of such enzyme mechanisms can pave the way for the design and delivery of enzyme-specific inhibitors.

Alanine mutation of the catalytic sites of Pantothenate Synthetase causes distinct conformational changes in the ATP binding region

The enzyme Pantothenate synthetase (PS) represents a potential drug target in Mycobacterium tuberculosis. Its X-ray crystallographic structure has demonstrated the significance and importance of conserved active site residues including His44, His47, Asn69, Gln72, Lys160 and Gln164 in substrate binding and formation of pantoyl adenylate intermediate. In the current study, molecular mechanism of decreased affinity of the enzyme for ATP caused by alanine mutations was investigated using molecular dynamics (MD) simulations and free energy calculations. A total of seven systems including wild-type + ATP, H44A + ATP, H47A + ATP, N69A + ATP, Q72A + ATP, K160A + ATP and Q164A + ATP were subjected to 50 ns MD simulations. Docking score, MM-GBSA and interaction profile analysis showed weak interactions between ATP (substrate) and PS (enzyme) in H47A and H160A mutants as compared to wild-type, leading to reduced protein catalytic activity. However, principal component analysis (PCA) and free energy landscape (FEL) analysis revealed that ATP was strongly bound to the catalytic core of the wild-type, limiting its movement to form a stable complex as compared to mutants. The study will give insight about ATP binding to the PS at the atomic level and will facilitate in designing of non-reactive analogue of pantoyl adenylate which will act as a specific inhibitor for PS.

An in-silico glimpse into the pH dependent structural changes of T7 RNA polymerase: a protein with simplicity

The capability of performing an array of functions with its single subunit structure makes T7 RNA polymerase (T7RNAP) as one of the simplest yet attractive target for various investigations ranging from structure determinations to several biological tests. In this study, with the help of molecular dynamics (MD) calculations and molecular docking, we investigated the effect of varying pH conditions on conformational flexibility of T7RNAP. We also studied its effect on the interactions with a well established inhibitor (heparin), substrate GTP and T7 promoter of T7RNAP. The simulation studies were validated with the help of three dimensional reconstructions of the polymerase at different pH environments using transmission electron microscopy and single particle analysis. On comparing the simulated structures, it was observed that the structure of T7RNAP changes considerably and interactions with its binding partners also changes as the pH shifts from basic to acidic. Further, it was observed that the C-terminal end plays a vital role in the inefficiency of the polymerase at low pH. Thus, this in-silico study may provide a significant insight into the structural investigations on T7RNAP as well as in designing potent inhibitors against it in varying pH environments.

Mechanistic insights into the activity of Ptf1-p48 (pancreas transcription factor 1a): probing the interactions levels of Ptf1-p48 with E2A-E47 (transcription factor E2-alpha) and ID3 (inhibitor of DNA binding 3)

Ptf1-p48 (Pancreas specific transcription factor 1a) is transcription regulatory protein known for the activation of exocrine specific genes. Downregulation of its expression formulates early stages of pancreatic adenocarcinoma as deduced by its association with oncogenic bHLH (Basic Helix-Loop-Helix) protein ID3 (Inhibitor of DNA binding 3) protein whose overexpression induces cytoplasmic mislocalization of Ptf1-p48. The precise mechanism and/or functional role of Ptf1-p48in promoting pancreatic cancer is vague. The structural features of the Ptf1-p48 and its dimerization with E47 (Transcription factor E2-alpha) and ID3 mediated by their HLH (Helix-Loop-Helix) domain were perceived through MD (Molecular Dynamics) simulations of 50 ns. The interactions formed by the HLH domain in both Ptf1-E47 and Ptf1-ID3 complexes are favored by the synergistic movement of their domain helices. Accordingly, in the Ptf1-E47 complex α7 of Ptf1-p48 and α1 helix of E47 along with the loop residues of their HLH domain exhibit transitions marked by inward movement toward each other and forms polar and charged interactions. In the Ptf1-ID3 complex, α8 of Ptf1-p48 moves toward the α3 helix of ID3 and forms hydrogen bonds. The interface analysis also reveals better interface in the Ptf1-p48 complex than the Ptf1-ID3 evident by energetics and number of hydrogen bonds. The interactions in each of these complexes, supported by angular displacement and mode vector analyzes, comprehensibly describe the considerable structural changes induced upon dimer formation. It thereby gives an insight into the interfaces that could help in designing of potential inhibitors for ID3 to curb the cancer cell growth.

A computational assessment of pH-dependent differential interaction of T7 lysozyme with T7 RNA polymerase

Background

T7 lysozyme (T7L), also known as N-acetylmuramoyl-L-alanine amidase, is a T7 bacteriophage gene product. It involves two functions: It can cut amide bonds in the bacterial cell wall and interacts with T7 RNA polymerase (T7RNAP) as a part of transcription inhibition. In this study, with the help of molecular dynamics (MD) calculations and computational interaction studies, we investigated the effect of varying pH conditions on conformational flexibilities of T7L and their influence on T7RNAP -T7L interactions.

Results

From the MD studies of the T7L at three different pH strengths viz. 5, neutral and 7.9 it was observed that T7L structure at pH 5 exhibited less stable nature with more residue level fluctuations, decrease of secondary structural elements and less compactness as compared to its counterparts: neutral pH and pH 7.9. The T-pad analysis of the MD trajectories identified local fluctuations in few residues that influenced the conformational differences in three pH strengths. From the docking of the minimum energy representative structures of T7L at different pH strengths (obtained from the free energy landscape analysis) with T7RNAP structures at same pH strengths, we saw strong interaction patterns at pH 7.9 and pH 5. The MD analysis of these complexes also confirmed the observations of docking study. From the combined in silico studies, it was observed that there are conformational changes in N-terminal and near helix 1 of T7L at different pH strengths, which are involved in the T7RNAP interaction, thereby varying the interaction pattern.

Conclusion

Since T7L has been used for developing novel therapeutics and T7RNAP one of the most biologically useful protein in both in-vitro and in vivo experiments, this in silico study of pH dependent conformational differences in T7L and the differential interaction with T7RNAP at different pH can provide a significant insight into the structural investigations on T7L and T7RNAP in varying pH environments.

Electronic supplementary material

The online version of this article (doi:10.1186/s12900-017-0077-9) contains supplementary material, which is available to authorized users.