Computationally guided identification of Akt-3, a serine/threonine kinase inhibitors: Insights from homology modelling,structure-based screening,molecular dynamics and quantum mechanical calculations

Abstract

Chemical entities targeting kinase signalling pathways serve as a potential strategy to combat malignancies. Protein Kinase B or Akt is a validated target for various malignancies and Akt3 remains the least explored isoform among all its isoforms. Initially, homology modelling technique was used for generating protein structure and further validation was performed using molecular dynamics simulation and Ramachandran plot. The validated protein structure was then subjected for active site analysis which led to identification of active site residues based on metrics provided by site score. The important residues in binding site were identified as Thr81, Asp271 and Asp289 for binding energetics and inhibition. Subsequently, virtual screening methodologies were used for identification of novel hits for inhibition of Protein Kinase B or Akt3. This led to the identification of two hits, i.e. thiophene derivative and thieno-pyridine derivative which were selected on the basis of their binding affinity and drug likeliness. These identified hits were subjected for molecular dynamics simulations, quantum mechanical and synthetic accessibility studies. The role of crucial residues in binding site stood validated as suggested by molecular dynamics simulations studies.

Communicated by Ramaswamy H. Sarma     

R516Q mutation in Melanoma differentiation-associated protein 5 (MDA5) and its pathogenic role towards rare Singleton-Merten syndrome; a signature associated molecular dynamics study

Singleton-Merten syndrome, a critical and rare multifactorial disorder that is closely linked to R516Q mutation in MDA5 protein associated with an enhanced interferon response in the affected individual. In the present study, we provide conclusive key evidence on R516Q mutation and their connectivity towards sequence-structural basis dysfunction of MDA5 protein. Among the various mutations, we found R516Q is the most pathogenic mutation based on mutational signature Q-A-[RE]-G-R-[GA]-R-A-[ED]-[DE]-S-[ST]-Y-[TSAV]-L-V designed from our work. Further, we derived a distant ortholog for this mutational signature from which we identified 343 intra-residue interactions that fall communally in the position required to maintain the structural and functional integration of protein architecture. This identification served us to understand the critical role of hot spots in residual aggregation that holds a native form of folding conformation in the functional region. In addition, the long-range molecular dynamics simulation demarcated the residual dependencies of conformational transition in distinct regions (_L29360–370_α18, _α19380–410_L31, _α21430–480_L33-α22-L35 and _α24510–520_L38) occurring upon R516Q mutation. Together, our results emphasise that the dislocation of functional hot spots Pro229, Arg414, Val498, Met510, Ala513, Gly515 and Arg516 in MDA5 protein which is important for interior structural packing and fold arrangements. In a nutshell, our findings are perfectly conceded with other experimental reports and will have potential implications in immune therapeutical advancement for rare singleton-merten syndrome.     

Anisotropy and anharmonicity of atomic fluctuations in proteins: analysis of a molecular dynamics simulation

Positional probability density functions (pdf) for the atomic fluctuations are determined from a molecular dynamics simulation for hen egg-white lysozyme. Most atoms are found to have motions that are highly anisotropic but only slightly anharmonic. The largest deviations from harmonic motion are in the direction of the largest rms fluctuations in the local principal axis frame. Backbone atoms tend to be more nearly harmonic than sidechain atoms. The atoms with the largest anharmonicities tend to have pdfs with multiple peaks, each of which is close to harmonic. Several model pdfs are evaluated on the basis of how well they fit probability densities from the dynamics simulations when parameterized in terms of the moments of the distribution. Gram-Charlier and Edgeworth perturbation expansions, which have been successful in describing the motions of small molecules in crystals, are shown to be inadequate for the distributions found in the dynamics of proteins. Multipeaked distribution functions are found to be more appropriate.     

Conformational dynamics of a short antigenic peptide in its free and antibody bound forms gives insight into the role of β‐turns in peptide immunogenicity

Earlier immunological experiments with a synthetic 36‐residue peptide (75‐110) from Influenza hemagglutinin have been shown to elicit anti‐peptide antibodies (Ab) which could cross‐react with the parent protein. In this article, we have studied the conformational features of a short antigenic (Ag) peptide (⁹⁸YPYDVPDYASLRS¹¹⁰) from Influenza hemagglutinin in its free and antibody (Ab) bound forms with molecular dynamics simulations using GROMACS package and OPLS‐AA/L all‐atom force field at two different temperatures (293 K and 310 K). Multiple simulations for the free Ag peptide show sampling of ordered conformations and suggest different conformational preferences of the peptide at the two temperatures. The free Ag samples a conformation crucial for Ab binding (β‐turn formed by “DYAS” sequence) with greater preference at 310 K while, it samples a native‐like conformation with relatively greater propensity at 293 K. The sequence “DYAS” samples β‐turn conformation with greater propensity at 310 K as part of the hemagglutinin protein also. The bound Ag too samples the β‐turn involving “DYAS” sequence and in addition it also samples a β‐turn formed by the sequence “YPYD” at its N‐terminus, which seems to be induced upon binding to the Ab. Further, the bound Ag displays conformational flexibility at both 293 K and 310 K, particularly at terminal residues. The implications of these results for peptide immunogenicity and Ag–Ab recognition are discussed. Proteins 2015; 83:1352–1367. © 2015 Wiley Periodicals, Inc.     

Identification of a conserved 8 aa insert in the PIP5K protein in the Saccharomycetaceae family of fungi and the molecular dynamics simulations and structural analysis to investigate its potential functional role

Homologs of the phosphatidylinositol‐4‐phosphate‐5‐kinase (PIP5K), which controls a multitude of essential cellular functions, contain a 8 aa insert in a conserved region that is specific for the Saccharomycetaceae family of fungi. Using structures of human PIP4K proteins as templates, structural models were generated of the Saccharomyces cerevisiae and human PIP5K proteins. In the modeled S. cerevisiae PIP5K, the 8 aa insert forms a surface exposed loop, present on the same face of the protein as the activation loop of the kinase domain. Electrostatic potential analysis indicates that the residues from 8 aa conserved loop form a highly positively charged surface patch, which through electrostatic interaction with the anionic portions of phospholipid head groups, is expected to play a role in the membrane interaction of the yeast PIP5K. To unravel this prediction, molecular dynamics (MD) simulations were carried out to examine the binding interaction of PIP5K, either containing or lacking the conserved signature insert, with two different membrane lipid bilayers. The results from MD studies provide insights concerning the mechanistic of interaction of PIP5K with lipid bilayer, and support the contention that the identified 8 aa conserved insert in fungal PIP5K plays an important role in the binding of this protein with membrane surface. Proteins 2017; 85:1454–1467. © 2017 Wiley Periodicals, Inc.     

Molecular and structural analysis of C4-specific PEPC isoform from Pennisetum glaucum plays a role in stress adaptation

As an activator of adenylate cyclase, the neuropeptide Pituitary Adenylate Cyclase Activating Peptide (PACAP) impacts levels of cyclic AMP, a key second messenger available in brain cells. PACAP is involved in certain adult behaviors. To elucidate PACAP interactions, a compendium of microarrays representing mRNA expression in the adult mouse whole brain was pooled from the Phenogen database for analysis. A regulatory network was computed based on mutual information between gene pairs using gene expression data across the compendium. Clusters among genes directly linked to PACAP, and probable interactions between corresponding proteins were computed. Database “experts” affirmed some of the inferred relationships. The findings suggest ADCY7 is probably the adenylate cyclase isoform most relevant to PACAP's action. They also support intervening roles for kinases including GSK3B, PI 3-kinase, SGK3 and AMPK. Other high-confidence interactions are hypothesized for future testing. This new information has implications for certain behavioral and other disorders.