Molecular Dynamic Simulation Reveals Damaging Impact of RAC1 F28L Mutation in the Switch I Region

Ras-related C3 botulinum toxin substrate 1 (RAC1) is a plasma membrane-associated small GTPase which cycles between the active GTP-bound and inactive GDP-bound states. There is wide range of evidences indicating its active participation in inducing cancer-associated phenotypes. RAC1 F28L mutation (RAC^F28L) is a fast recycling mutation which has been implicated in several cancer associated cases. In this work we have performed molecular docking and molecular dynamics simulation (~0.3 μs) to investigate the conformational changes occurring in the mutant protein. The RMSD, RMSF and NHbonds results strongly suggested that the loss of native conformation in the Switch I region in RAC1 mutant protein could be the reason behind its oncogenic transformation. The overall results suggested that the mutant protein attained compact conformation as compared to the native. The major impact of mutation was observed in the Switch I region which might be the crucial reason behind the loss of interaction between the guanine ring and F28 residue.     

Relationship between a point mutation S97C in CK1δ protein and its affect on ATP-binding affinity

CK1δ (Casein kinase I isoform delta) is a member of CK1 kinase family protein that mediates neurite outgrowth and the function as brain-specific microtubule-associated protein. ATP binding kinase domain of CK1δ is essential for regulating several key cell cycle signal transduction pathways. Mutation in CK1δ protein is reported to cause cancers and affects normal brain development. S97C mutation in kinase domain of CK1δ protein has been involved to induce breast cancer and ductal carcinoma. We performed molecular docking studies to examine the effect of this mutation on its ATP-binding affinity. Further, we conducted molecular dynamics simulations to understand the structural consequences of S97C mutation over the kinase domain of CK1δ protein. Docking results indicated the loss of ATP-binding affinity of mutant structure, which were further rationalized by molecular dynamics simulations, where a notable loss in 3-D conformation of CK1δ kinase domain was observed in mutant as compared to native. Our results explained the underlying molecular mechanism behind the observed cancer associated phenotype caused by S97C mutation in CK1δ protein.     

PhyA gene product of Aspergillus ficuum and Peniophora lycii produces dissimilar phytases

PhyA gene products of Aspergillus ficuum (AF) and Peniophora lycii (PL) as expressed in industrial strains of Aspergillus niger and Aspergillus oryzae, respectively, were purified to homogeneity and then characterized for both physical and biochemical properties. The PL phytase is 26 amino acid residues shorter than the AF phytase. Dynamic light scattering studies indicate that the active AF phytase is a monomer while the PL phytase is a dimer. While both of the phytases retained four identical glycosylatable Asn residues, unique glycosylation sites, six for PL and seven for AF phytase, were observed. Global alignment of both the phytases has shown 38% sequence homology between the two proteins. At 58 degrees C and pH 5.0, the PL phytase gave a specific activity of 22,000 nKat/mg as opposed to about 3000 nKat/mg for AF phytase. However, the AF phytase is more thermostable than its counterpart PL phytase at 65 degrees C. Also, AF phytase is more stable at pH 7.5 than the PL phytase. The two phytases differed in K(m) for phytate, K(i) for myo-inositol hexasulfate (MIHS), and pH optima profile. Despite similarities in the active site sequences, the two phytases show remarkable differences in turnover number, pH optima profile, stability at higher temperature, and alkaline pH. These biochemical differences indicate that phytases from ascomycete and basidiomycete fungi may have evolved to degrade phytate in different environments.     

Site-directed mutagenesis of disulfide bridges in Aspergillus niger NRRL 3135 phytase (PhyA), their expression in Pichia pastoris and catalytic characterization

Earlier studies have established the importance of five disulfide bridges (DBs) in Aspergillus niger phytase. In this study, the relative importance of each of the individual disulfide bridge is determined by its removal by site-directed mutagenesis of specific cysteines in the cloned A. niger phyA gene. Individually, these mutant phytases were expressed in a Pichia expression system and their product purified and characterized. The removal of disulfide bridge 2 yielded a mutant phytase with a complete loss of catalytic activity. The other disulfide mutants displayed a broad array of altered catalytic properties including a lower optimum temperature from 58°C to 53°C for bridge number 1, 37°C for bridge number 3 and 4, and 42°C for bridge number 5. The pH versus activity profile was also modified in the DB mutants. The pH profile of the wild-type phytase was modified by the DB mutations. In bridge number 1, 3, and 4, the second peak at pH 2.5 was abolished, and in bridge number 5, the peak at pH 5.0 was abolished completely leaving only the pH 2.5. While the K _m was not affected drastically, the turnover number was lowered significantly in bridge number 3, 4, and 5.     

Nature of cation-π interactions and their role in structural stability of immunoglobulin proteins

Cation-π interactions are known to be important contributors to protein stability and ligand-protein interactions. In this study, we have analyzed the influence of cation-π interactions in single chain immunoglobulin proteins. We observed 87 cation-π interactions in a data set of 33 proteins. These interactions are mainly formed by long-range contacts, and there is preference of Arg over Lys in these interactions. Arg-Tyr interactions are predominant among the various pairs analyzed. Despite the scarcity of interactions involving Trp, the average energy for Trp-cation interactions is quite high. This information suggests that the cation-π interactions involving Trp might be of high relevance to the proteins. Secondary structure analysis reveals that cation-π interactions are formed preferably between residues in which at least one is in β-strand. Proteins having β-strand regions have the highest number of cation-π interaction-forming residues.     

Exploring the role of cation-pi interactions in glycoproteins lipid-binding proteins and RNA-binding proteins

We have analyzed and compared the influence of cation-pi interactions in glycoproteins (GPs), lipid-binding proteins (LBPs) and RNA-binding proteins (RBPs) in this study. We observed that all the proteins included in the study had profound cation-pi interactions. There is an average of one energetically significant cation-pi interaction for every 71 residues in GPs, for every 58 residues in LBPs and for every 64 residues in RBPs. Long-range contacts are predominant in all the three types of proteins studied. The pair-wise cation-pi interaction energy between the positively charged and aromatic residues shows that Arg-Trp pair energy was the strongest among all six possible pairs in all the three types of proteins studied. There were considerable differences in the preference of cation-pi interacting residues to different secondary structure elements and ASA and these might contribute to differences in biochemical functions of GPs, LBPs and RBPs. It was interesting to note that all the five residues involved in cation-pi interactions were found to have stabilization centers in GPs, LBPs and RBPs. Majority of the cation-pi interacting residues investigated in the present study had a conservation score of 6, the cutoff value used to identify the stabilizing residues. A small percentage of cation-pi interacting residues were also present as stabilizing residues. The cation-pi interaction-forming residues play an important role in the structural stability of in GPs, LBPs and RBPs. The results obtained in this study will be helpful in further understanding the stability, specificity and differences in the biochemical functions of GPs, LBPs and RBPs.     

Investigations on unconventional hydrogen bonds in RNA binding proteins: the role of CH...O=C interactions

We have investigated the roles played by C-H...O=C interactions in RNA binding proteins. There was an average of 78 CH...O=C interactions per protein and also there was an average of one significant CH...O=C interactions for every 6 residues in the 59 RNA binding proteins studied. Main chain-Main chain (MM) CH...O=C interactions are the predominant type of interactions in RNA binding proteins. The donor atom contribution to CH...O=C interactions was mainly from aliphatic residues. The acceptor atom contribution for MM CH...O=C interactions was mainly from Val, Phe, Leu, Ile, Arg and Ala. The secondary structure preference analysis of CH...O=C interacting residues showed that, Arg, Gln, Glu and Tyr preferred to be in helix, while Ala, Asp, Cys, Gly, Ile, Leu, Lys, Met, Phe, Trp and Val preferred to be in strand conformation. Most of the CH...O=C interacting polar amino acid residues were solvent exposed while, majority of the CH...O=C interacting non polar residues were excluded from the solvent. Long and medium-range CH...O=C interactions are the predominant type of interactions in RNA binding proteins. More than 50% of CH...O=C interacting residues had a higher conservation score. Significant percentage of CH...O=C interacting residues had one or more stabilization centers. Sixty-six percent of the theoretically predicted stabilizing residues were also involved in CH...O=C interactions and hence these residues may also contribute additional stability to RNA binding proteins.     

Identification and in silico analysis of functional SNPs of the BRCA1 gene

Single-nucleotide polymorphisms (SNPs) play a major role in the understanding of the genetic basis of many complex human diseases. Also, the genetics of human phenotype variation could be understood by knowing the functions of these SNPs. It is still a major challenge to identify the functional SNPs in a disease-related gene. In this work, we have analyzed the genetic variation that can alter the expression and the function of the BRCA1 gene using computational methods. Of the total 477 SNPs, 65 were found to be nonsynonymous (ns) SNPs. Among the 14 SNPs in the untranslated region, 4 were found in the 5' and 10 were found in the 3' untranslated region (UTR). It was found that 16.9% of the nsSNPs were damaging, by both the SIFT and the PolyPhen servers. The UTR Resource tool suggested that 2 of 4 SNPs in the 5' UTR and 3 of 10 SNPs in the 3' UTR might change the protein expression levels. We identified major mutations from proline to serine at positions 1776 and 1812 of the native protein of the BRCA1 gene. From a comparison of the stabilizing residues of the native and mutant proteins, we propose that an nsSNP (rs1800751) could be an important candidate for the breast cancer caused by the BRCA1 gene.     

Computation of non-covalent interactions in TNF proteins and interleukins

The roles played by the non-covalent interactions have been investigated for a set of six TNF proteins and nine Interleukins. The stabilizing residues have been identified by a consensus approach using the concepts of available surface area, medium and long-range interactions and conservation of amino acid residues. The cation-pi interactions have been computed based on a geometric approach such as distance and energy criteria. We identified an average of 1 energetically significant cation-pi interactions in every 94 residues in TNF proteins and 1 in every 62 residues in Interleukins. In TNF proteins, the cationic groups Lys preferred to be in helix while Arg preferred to be in strand regions while in Interleukins the Arg residues preferred to be in helix and Lys preferred to be in strand regions. From the available surface area calculations, we found that, almost all the cation and pi residues in TNF proteins and Interleukins were either in buried or partially buried regions and none of them in the exposed regions. Medium and long-range interactions were predominant in both TNF proteins and Interleukins. It was observed that the percentage of stabilizing centers were more in TNF proteins as compared to the Interleukins, while the percentage of conserved residues were more in Interleukins than in TNF proteins. In the stabilizing residues Lys was observed to be a stabilizing residue in both TNF proteins and Interleukins. Among the aromatic group, Phe was seen to be a stabilizing residue in both TNF and Interleukins. We suggest that this study on the computation of cation-pi interactions in TNF proteins and Interleukins would be very helpful in further understanding the structure, stability and functional similarity of these proteins.     

Synthesis,anticoagulant and PIVKA-II induced by new 4-hydroxycoumarin derivatives

The action of the coumarin-type drugs and related compounds is reviewed to their VKOR antagonistic effects. In our study, twenty 3-pyridinyl, pyrimidinyl and pyrazolyl-4-hydroxycoumarin derivatives were synthesized. A comparative in vivo (CT, PT determination) and in vitro (measurement of PIVKA-II levels) anticoagulant study with respect to warfarin showed that the synthesized compounds have different anticoagulant activities, the most prospective compounds were the 3-pyrazolyl-4-hydroxycoumarin derivatives.