Insight into the molecular switch mechanism of human Rab5a from molecular dynamics simulations

An in-silico study of FIKK9.5 protein of Plasmodium falciparum for identification of therapeutics

The FIKK protein family, encompassing 21 serine-threonine protein kinases, is a distinctive cluster exclusive to the Apicomplexa phylum. Predominantly located in Plasmodium falciparum which is a malarial parasite, with a solitary gene identified in a distinct apicomplexan species, this family derives its nomenclature from - phenylalanine, isoleucine, lysine, lysine (FIKK), a conserved amino acid motif. Integral to the parasite's life cycle and consequential to malaria pathogenesis, the absence of orthologous proteins in eukaryotic organisms designates it as a promising antimalarial drug target. Among the FIKKs, FIKK9.5 plays a pivotal role in the parasite's development within red blood cells (RBCs). This investigation acquired the three-dimensional structure of FIKK9.5 and its ligands through extensive database searches and literature review. Computational screening of natural phytochemicals derived from plants traditionally used in antimalarial remedies was conducted by employing the Glide docking suite. AutoDock Vina was utilized to discern the inhibitor exhibiting optimal binding affinity. Subsequently, Molecular Dynamics (MD) simulations employing GROMACS validated Rufigallol as the most potent inhibitory compound against FIKK9.5. The robustness of the protein-ligand complex was scrutinized through a 200 nanosecond molecular dynamics (MD) trajectory. Trajectory analysis and determination of binding free energies were accomplished using MM-GBSA and MM-PBSA approaches. The ligand-binding exhibited sustained stability throughout the simulation, manifesting an approximate binding free energy of -25.5986 kcal/mol. This comprehensive computational study lays the groundwork for potential experimental validation in the laboratory, paving the way for the development of novel therapeutics targeting FIKK9.5 in the pursuit of innovative antimalarial.

Unorthodox regulation of the MglA Ras-like GTPase controlling polarity in Myxococcus xanthus

Motile cells have developed a large array of molecular machineries to actively change their direction of movement in response to spatial cues from their environment. In this process, small GTPases act as molecular switches and work in tandem with regulators and sensors of their guanine nucleotide status (GAP, GEF, GDI and effectors) to dynamically polarize the cell and regulate its motility. In this review, we focus on Myxococcus xanthus as a model organism to elucidate the function of an atypical small Ras GTPase system in the control of directed cell motility. M. xanthus cells direct their motility by reversing their direction of movement through a mechanism involving the redirection of the motility apparatus to the opposite cell pole. The reversal frequency of moving M. xanthus cells is controlled by modular and interconnected protein networks linking the chemosensory-like frizzy (Frz) pathway - that transmits environmental signals - to the downstream Ras-like Mgl polarity control system - that comprises the Ras-like MglA GTPase protein and its regulators. Here, we discuss how variations in the GTPase interactome landscape underlie single-cell decisions and consequently, multicellular patterns.

Unravelling Binding of Human Serum Albumin with Galantamine: Spectroscopic, Calorimetric, and Computational Approaches

Human serum albumin (HSA), an abundant plasma protein, binds to various ligands, acting as a transporter for numerous endogenous and exogenous substances. Galantamine (GAL), an alkaloid, treats cognitive decline in mild to moderate Alzheimer's disease and other memory impairments. A vital step in pharmacological profiling involves the interaction of plasma protein with the drugs, and this serves as an essential platform for pharmaceutical industry advancements. This study is carried out to understand the binding mechanism of GAL with HSA using computational and experimental approaches. Molecular docking revealed that GAL preferentially occupies Sudlow's site I, i.e., binds to subdomain IIIA. The results unveiled that GAL binding does not induce any conformational change in HSA and hence does not compromise the functionality of HSA. Molecular dynamics simulation (250 ns) deciphered the stability of the HSA-GAL complex. We performed the fluorescence binding and isothermal titration calorimetry (ITC) to analyze the actual binding of GAL with HSA. The results suggested that GAL binds to HSA with a significant binding affinity. ITC measurements also delineated thermodynamic parameters associated with the binding of GAL to HSA. Altogether, the present study deciphers the binding mechanism of GAL with HSA.

A Distinct Motif in a Prokaryotic Small Ras-Like GTPase Highlights Unifying Features of Walker B Motifs in P-Loop NTPases

A hallmark of the catalytically essential Walker B motif of P-loop NTPases is the presence of an acidic residue (aspartate/glutamate) for efficient Mg²⁺ coordination. Although the Walker B motif has been identified in well-studied examples of P-loop NTPases, its identity is ambiguous in many families, for example, in the prokaryotic small Ras-like GTPase family of MglA. MglA, belonging to TRAFAC class of P-loop NTPases, possesses a threonine at the position equivalent to Walker B aspartate in eukaryotic Ras-like GTPases. To resolve the identity of the Walker B residue in MglA, we carried out a comprehensive analysis of Mg²⁺ coordination on P-loop NTPase structures. Atoms in the octahedral coordination of Mg²⁺ and their interactions comprise a network including water molecules, Walker A, Walker B and switch motifs of P-loop NTPases. Based on the conserved geometry of Mg²⁺ coordination, we confirm that a conserved aspartate functions as the Walker B residue of MglA, and validate it through mutagenesis and biochemical characterization. Location of the newly identified aspartate is spatially equivalent to the Walker B residue of the ASCE division of P-loop NTPases. Furthermore, similar to the allosteric regulation of the Walker B aspartate conformation in MglA, we identify protein families in which large conformational changes involving Walker B motif potentially function as allosteric regulators. The study unravels conserved features of Mg²⁺ coordination among divergent families of P-loop NTPases, especially between ancient Ras-like GTPases and ASCE family of ATPases. The conserved geometric features provide a foundation for design of nucleotide-hydrolyzing enzymes.

Identification and structure-activity relationship (SAR) studies of carvacrol derivatives as potential anti-malarial against Plasmodium falciparum falcipain-2 protease

In an effort to develop a potent anti-malarial agent against Plasmodium falciparum, a structure-guided virtual screening using an in-house library comprising 652 compounds was performed. By docking studies, we identified two compounds (JMI-105 and JMI-346) which formed significant non-covalent interactions and fit well in the binding pocket of PfFP-2. We affirmed this observation by MD simulation studies. As evident by the biochemical analysis, such as enzyme inhibition assay, Surface Plasmon Resonance (SPR), live-cell imaging and hemozoin inhibition, JMI-105 and JMI-346 at 25 µM concentration showed an inhibitory effect on purified PfFP-2. JMI-105 and JMI-346 inhibited the growth of CQ^S (3D7; IC₅₀ = 8.8 and 13 µM) and CQ^R (RKL-9; IC₅₀ = 14.3 and 33 µM) strains of P. falciparum. Treatment with compounds resulted in defect in parasite growth and development. No significant hemolysis or cytotoxicity towards human cells was observed suggesting that these molecules are non-toxic. We pursued, structural optimization on JMI-105 and in the process, SAR oriented derivatives (5a-5l) were synthesized and evaluated for growth inhibition potential. JMI-105 significantly decreased parasitemia and prolonged host survival in a murine model with P. berghei ANKA infection. The compounds (JMI-105 and JMI-346) against PfFP-2 have the potential to be used as an anti-malarial agent.

An Insightful 10-year Recollection Since the Emergence of the 5-steps Rule

OBJECTIVE

One of the most challenging and also the most difficult problems is how to formulate a biological sequence with a vector but considerably keep its sequence order information.

METHODS

To address such a problem, the approach of Pseudo Amino Acid Components or PseAAC has been developed.

RESULTS AND CONCLUSION

It has become increasingly clear via the 10-year recollection that the aforementioned proposal has been indeed very powerful.

Differential levels of Alpha-1-inhibitor III, Immunoglobulin heavy chain variable region, and Hypertrophied skeletal muscle protein GTF3 in rat mammary tumorigenesis

Early detection of breast cancer can be best facilitated by the development of precancerous markers. Serum proteins being the sensitive signatures, can be the ideal choice. We previously demonstrated the reduced levels of two serum proteins at a very early stage of tumorigenesis in a breast cancer model, developed in Wistar rats by 7,12-dimethylbenz[a]anthracene (DMBA) administration. Here we report the dysregulation of three more proteins in the serum collected at another early stage (15 weeks) of tumorigenesis in the same model. The proteins were identified (as Alpha-1-inhibitor III (Mug1), Immunoglobulin heavy chain variable region (IGHV), and Hypertrophied skeletal muscle protein GTF3) by MALDI-TOF MS after the screening and fingerprinting of serum samples by one-dimensional (1D) and two-dimensional (2D) electrophoresis respectively. Relative expression analysis of corresponding genes was also carried out, and the results were found as supporting the proteomic findings. In addition, the candidate proteins of the study and their corresponding ribonucleic acids (RNAs) were subjected to homology modelling and docking (using softwares like MODELLER, 3dRNA, Autodock4.0, and GROMACS etc), which revealed the binding sites for carcinogen (DMBA) and its nature of interaction with proteins and RNAs. Moreover, the network analysis by GeneMANIA unraveled the protein/gene functional network in which Mug1, IGHV, and GTF3 are involved. Based on the significant protein and gene expression alterations in early tumorigenesis, these proteins may prove very effective in search for biomarkers for the early detection of mammary cancer. Further, these proteins can also be tried as targets for chemotherapy.

Progress in the Mechanism and Clinical Application of Cilostazol

Cilostazol is a unique platelet inhibitor that has been used clinically for more than 20 years. As a phosphodiesterase type III inhibitor, cilostazol is capable of reversible inhibition of platelet aggregation and vasodilation, has antiproliferative effects, and is widely used in the treatment of peripheral arterial disease, cerebrovascular disease, percutaneous coronary intervention, etc. This article briefly reviews the pharmacological mechanisms and clinical application of cilostazol.

Evaluation of ellagic acid as an inhibitor of sphingosine kinase 1: A targeted approach towards anticancer therapy

Sphingosine kinase 1 (SphK1) is one of the central enzymes of sphingolipid metabolism whose high expression level is presumed to be correlated with cancer and other inflammatory diseases. Using a virtual screening approach and in vitro studies, we have identified the ellagic acid (EA), a dietary polyphenol, as a potent inhibitor of SphK1. Molecular docking study has suggested a strong binding affinity of EA to the SphK1. Fluorescence binding and isothermal titration calorimetry (ITC) measurements has also indicated an appreciable binding affinity. Kinase inhibition assay revealed an excellent inhibitory action of EA towards SphK1 (IC₅₀ = 0.74 ± 0.06 μM). Cell viability studies point towards the antiproliferative effects of EA on lung cancer cell line (A549) without affecting human embryonic kidney cells (HEK293). Binding and inhibition mechanism of EA was unveiled by docking analysis of SphK1-EA complex. EA binds to the SphK1 and forms several interactions with catalytically important residues of ATP-binding pocket. Structural stability and dynamics analysis of SphK1-EA complex during 100 ns molecular dynamic simulation studies suggested that EA forms a stable complex with SphK1 without inducing any significant conformational shift. Taken together, our study suggests that EA can be utilized as a chemical prototype to develop potent therapeutics targeting SphK1-associated pathologies.

Molecular interaction studies on ellagic acid for its anticancer potential targeting pyruvate dehydrogenase kinase 3

Pyruvate dehydrogenase kinase 3 (PDK3) plays a central role in the cancer metabolic switch through the reversible phosphorylation of pyruvate dehydrogenase complex thereby blocking the entry of pyruvate for its catabolism into the TCA cycle, and thus it is considered as an important drug target for various types of cancers. We have successfully expressed full length human PDK3 and investigated its interaction mechanism with dietary polyphenols in the search for potential inhibitors. Molecular docking analysis revealed that the selected compounds preferentially bind to the ATP-binding pocket of PDK3 and interact with functionally important residues. In silico observations were further complemented by experimental measurements of the fluorescence quenching of PDK3 and confirmed with the isothermal titration calorimetry measurements. Ellagic acid (EA) significantly binds and inhibits the kinase activity of PDK3. In vitro cytotoxicity and the anti-proliferative properties of EA were evaluated by MTT assay. Conformational dynamics of the EA-PDK3 complex during molecular dynamics simulation revealed that a stable complex was maintained by a significant number of hydrogen bonds throughout the 100 ns trajectories. In conclusion, EA may be considered as a promising molecule for PDK3 inhibition and could be exploited as a lead molecule against PDK3 associated diseases.

Recognition and stabilization of geranylgeranylated human Rab5 by the GDP Dissociation Inhibitor (GDI)

The small GTPase Rab5 is the key regulator of early endosomal fusion. It is post-translationally modified by covalent attachment of two geranylgeranyl (GG) chains to adjacent cysteine residues of the C-terminal hypervariable region (HVR). The GDP dissociation inhibitor (GDI) recognizes membrane-associated Rab5(GDP) and serves to release it into the cytoplasm where it is kept in a soluble state. A detailed new structural and dynamic model for human Rab5(GDP) recognition and binding with human GDI at the early endosome membrane and in its dissociated state is presented. In the cytoplasm, the GDI protein accommodates the GG chains in a transient hydrophobic binding pocket. In solution, two different binding modes of the isoprenoid chains inserted into the hydrophobic pocket of the Rab5(GDP):GDI complex can be identified. This equilibrium between the two states helps to stabilize the protein-protein complex in solution. Interprotein contacts between the Rab5 switch regions and characteristic patches of GDI residues from the Rab binding platform (RBP) and the C-terminus coordinating region (CCR) reveal insight on the formation of such a stable complex. GDI binding to membrane-anchored Rab5(GDP) is initially mediated by the solvent accessible switch regions of the Rab-specific RBP. Formation of the membrane-associated Rab5(GDP):GDI complex induces a GDI reorientation to establish additional interactions with the Rab5 HVR. These results allow to devise a detailed structural model for the process of extraction of GG-Rab5(GDP) by GDI from the membrane and the dissociation from targeting factors and effector proteins prior to GDI binding.

Hydrogen bond analysis of the EGFR-ErbB3 heterodimer related to non-small cell lung cancer and drug resistance

Lung cancer is the predominant cause of cancer deaths on a worldwide scale. A mutation in the epidermal growth factor receptor (EGFR) can cause non-small cell lung cancer (NSCLC). The L858R one-point mutation in exon 21 in EGFR is the most prevalent in NSCLC. For over 60% of EGFR-muted NSCLC, another mutation T790M can cause drug resistance. In this paper, we consider EGFR and ErbB3 heterodimers involving three structures of EGFR, wild-type, with L858R mutation, and with L858R and T790M mutations. We perform molecular dynamics (MD) simulations to analyze hydrogen bonds in all three instances. The hydrogen bonds contribute to the conformational stability of the protein and molecular recognition. Several other parameters are also investigated in the present study, which reveals significant changes in the dimer at different levels of mutation. The knowledge and results obtained from this study lead to useful insight into the mechanism of NSCLC drug resistance.

Exploration of synthetic multifunctional amides as new therapeutic agents for Alzheimer's disease through enzyme inhibition, chemoinformatic properties, molecular docking and dynamic simulation insights

A new series of multifunctional amides has been synthesized having moderate enzyme inhibitory potentials and mild cytotoxicity. 2-Furyl(1-piperazinyl)methanone (1) was coupled with 3,5-dichloro-2-hydroxybenzenesulfonyl chloride (2) to form {4-[(3,5-dichloro-2-hydroxyphenyl)sulfonyl]-1-piperazinyl}(2-furyl)methanone (3). Different elecrophiles were synthesized by the reaction of various un/substituted anilines (4a-o) with 2-bromoacetylbromide (5), 2‑bromo‑N-(un/substituted-phenyl)acetamides (6a-o). Further, equimolar ratios of 3 and 6a-o were allowed to react in the presence of K₂CO₃ in acetonitrile to form desired multifunctional amides (7a-o). The structural confirmation of all the synthesized compounds was carried out by their EI-MS, IR, ¹H NMR and ¹³C NMR spectral data. Enzyme inhibition activity was performed against acetyl and butyrylcholinestrase enzymes, whereby 7e showed very good activity having IC₅₀ value of 5.54 ± 0.03 and 9.15 ± 0.01 μM, respectively, relative to eserine, a reference standard. Hemolytic activity of the molecules was checked to asertain their cytotoxicity towards red blood cell membrance and it was observed that most of the compounds were not toxic up to certain range. Moreover, chemoinformatic protepties and docking simulation results also showed the significance of 7e as compared to other compounds. Based on in vitro and in silico analysis 7e could be used as a template for the development of new drugs against Alzheimer's disease.

Structural Changes Due to Antagonist Binding in Ligand Binding Pocket of Androgen Receptor Elucidated Through Molecular Dynamics Simulations

When a small molecule binds to the androgen receptor (AR), a conformational change can occur which impacts subsequent binding of co-regulator proteins and DNA. In order to accurately study this mechanism, the scientific community needs a crystal structure of the Wild type AR (WT-AR) ligand binding domain, bound with antagonist. To address this open need, we leveraged molecular docking and molecular dynamics (MD) simulations to construct a structure of the WT-AR ligand binding domain bound with antagonist bicalutamide. The structure of mutant AR (Mut-AR) bound with this same antagonist informed this study. After molecular docking analysis pinpointed the suitable binding orientation of a ligand in AR, the model was further optimized through 1 μs of MD simulations. Using this approach, three molecular systems were studied: (1) WT-AR bound with agonist R1881, (2) WT-AR bound with antagonist bicalutamide, and (3) Mut-AR bound with bicalutamide. Our structures were very similar to the experimentally determined structures of both WT-AR with R1881 and Mut-AR with bicalutamide, demonstrating the trustworthiness of this approach. In our model, when WT-AR is bound with bicalutamide, Val716/Lys720/Gln733, or Met734/Gln738/Glu897 move and thus disturb the positive and negative charge clumps of the AF2 site. This disruption of the AF2 site is key for understanding the impact of antagonist binding on subsequent co-regulator binding. In conclusion, the antagonist induced structural changes in WT-AR detailed in this study will enable further AR research and will facilitate AR targeting drug discovery.

Potential entry inhibitors of the envelope protein (E2) of Chikungunya virus: in silico structural modeling, docking and molecular dynamic studies

Chikungunya fever is an arboviral infection caused by the Chikungunya virus (CHIKV) and is transmitted by Aedes mosquito. The envelope protein (E2) of Chikungunya virus is involved in attachment of virion with the host cell. The present study was conceptualized to determine the structure of E2 protein of CHIKV and to identify the potential viral entry inhibitors. The secondary and tertiary structure of E2 protein was determined using bioinformatics tools. The mutational analysis of the E2 protein suggested that mutations may stabilize or de-stabilize the structure which may affect the structure-function relationship. In silico screening of various compounds from different databases identified two lead molecules i.e. phenothiazine and bafilomycin. Molecular docking and MD simulation studies of the E2 protein and compound complexes was carried out. This analysis revealed that bafilomycin has high docking score and thus high binding affinity with E2 protein suggesting stable protein-ligand interaction. Further, MD simulations suggested that both the compounds were stabilizing E2 protein. Thus, bafilomycin and phenothiazine may be considered as the lead compounds in terms of potential entry inhibitor for CHIKV. Further, these results should be confirmed by comprehensive cell culture, cytotoxic assays and animal experiments. Certain derivatives of phenothiazines can also be explored in future studies for entry inhibitors against CHIKV. The present investigation thus provides insight into protein structural dynamics of the envelope protein of CHIKV. In addition the study also provides information on the dynamics of interaction of E2 protein with entry inhibitors that will contribute towards structure based drug design.

3D pharmacophore-based virtual screening, docking and density functional theory approach towards the discovery of novel human epidermal growth factor receptor-2 (HER2) inhibitors

Human epidermal growth factor receptor 2 (HER2) is one of the four members of the epidermal growth factor receptor (EGFR) family and is expressed to facilitate cellular proliferation across various tissue types. Therapies targeting HER2, which is a transmembrane glycoprotein with tyrosine kinase activity, offer promising prospects especially in breast and gastric/gastroesophageal cancer patients. Persistence of both primary and acquired resistance to various routine drugs/antibodies is a disappointing outcome in the treatment of many HER2 positive cancer patients and is a challenge that requires formulation of new and improved strategies to overcome the same. Identification of novel HER2 inhibitors with improved therapeutics index was performed with a highly correlating (r=0.975) ligand-based pharmacophore model (Hypo1) in this study. Hypo1 was generated from a training set of 22 compounds with HER2 inhibitory activity and this well-validated hypothesis was subsequently used as a 3D query to screen compounds in a total of four databases of which two were natural product databases. Further, these compounds were analyzed for compliance with Veber's drug-likeness rule and optimum ADMET parameters. The selected compounds were then subjected to molecular docking and Density Functional Theory (DFT) analysis to discern their molecular interactions at the active site of HER2. The findings thus presented would be an important starting point towards the development of novel HER2 inhibitors using well-validated computational techniques.

iPPBS-Opt: A Sequence-Based Ensemble Classifier for Identifying Protein-Protein Binding Sites by Optimizing Imbalanced Training Datasets

Knowledge of protein-protein interactions and their binding sites is indispensable for in-depth understanding of the networks in living cells. With the avalanche of protein sequences generated in the postgenomic age, it is critical to develop computational methods for identifying in a timely fashion the protein-protein binding sites (PPBSs) based on the sequence information alone because the information obtained by this way can be used for both biomedical research and drug development. To address such a challenge, we have proposed a new predictor, called iPPBS-Opt, in which we have used: (1) the K-Nearest Neighbors Cleaning (KNNC) and Inserting Hypothetical Training Samples (IHTS) treatments to optimize the training dataset; (2) the ensemble voting approach to select the most relevant features; and (3) the stationary wavelet transform to formulate the statistical samples. Cross-validation tests by targeting the experiment-confirmed results have demonstrated that the new predictor is very promising, implying that the aforementioned practices are indeed very effective. Particularly, the approach of using the wavelets to express protein/peptide sequences might be the key in grasping the problem's essence, fully consistent with the findings that many important biological functions of proteins can be elucidated with their low-frequency internal motions. To maximize the convenience of most experimental scientists, we have provided a step-by-step guide on how to use the predictor's web server (http://www.jci-bioinfo.cn/iPPBS-Opt) to get the desired results without the need to go through the complicated mathematical equations involved.

Identification of Glucose-Binding Pockets in Human Serum Albumin Using Support Vector Machine and Molecular Dynamics Simulations

Human Serum Albumin (HSA) has been suggested to be an alternate biomarker to the existing Hemoglobin-A1c (HbA1c) marker for glycemic monitoring. Development and usage of HSA as an alternate biomarker requires the identification of glycation sites, or equivalently, glucose-binding pockets. In this work, we combine molecular dynamics simulations of HSA and the state-of-art machine learning method Support Vector Machine (SVM) to predict glucose-binding pockets in HSA. SVM uses the three dimensional arrangement of atoms and their chemical properties to predict glucose-binding ability of a pocket. Feature selection reveals that the arrangement of atoms and their chemical properties within the first 4Å from the centroid of the pocket play an important role in the binding of glucose. With a 10-fold cross validation accuracy of 84 percent, our SVM model reveals seven new potential glucose-binding sites in HSA of which two are exposed only during the dynamics of HSA. The predictions are further corroborated using docking studies. These findings can complement studies directed towards the development of HSA as an alternate biomarker for glycemic monitoring.

Identification of potential CCR5 inhibitors through pharmacophore-based virtual screening, molecular dynamics simulation and binding free energy analysis

Employing the combined strategy to identify novel CCR5 inhibitors and provide a basis for rational drug design.

Modelling the molecular mechanism of protein-protein interactions and their inhibition: CypD-p53 case study

Cyclophilin D (CypD) is an important regulatory protein involved in mitochondrial membrane permeability transition and cell death. Further, the mitochondrial CypD-p53 axis is an important contributor to necroptosis, a form of programmed necrosis, involved in various cardiovascular and neurological disorders. The CypD ligand, Cyclosporin A (CsA), was identified as an inhibitor of this interaction. In this study, using computational methods, we have attempted to model the CypD-p53 interaction in order to delineate their mode of binding and also to disclose the molecular mechanism, by means of which CsA interferes with this interaction. It was observed that p53 binds at the CsA-binding site of CypD. The knowledge obtained from this modelling was employed to identify novel CypD inhibitors through structure-based methods. Further, the identified compounds were tested by a similar strategy, adopted during the modelling process. This strategy could be applied to study the mechanism of protein-protein interaction (PPI) inhibition and to identify novel PPI inhibitors.

Identification of protein-protein binding sites by incorporating the physicochemical properties and stationary wavelet transforms into pseudo amino acid composition

With the explosive growth of protein sequences entering into protein data banks in the post-genomic era, it is highly demanded to develop automated methods for rapidly and effectively identifying the protein-protein binding sites (PPBSs) based on the sequence information alone. To address this problem, we proposed a predictor called iPPBS-PseAAC, in which each amino acid residue site of the proteins concerned was treated as a 15-tuple peptide segment generated by sliding a window along the protein chains with its center aligned with the target residue. The working peptide segment is further formulated by a general form of pseudo amino acid composition via the following procedures: (1) it is converted into a numerical series via the physicochemical properties of amino acids; (2) the numerical series is subsequently converted into a 20-D feature vector by means of the stationary wavelet transform technique. Formed by many individual "Random Forest" classifiers, the operation engine to run prediction is a two-layer ensemble classifier, with the 1st-layer voting out the best training data-set from many bootstrap systems and the 2nd-layer voting out the most relevant one from seven physicochemical properties. Cross-validation tests indicate that the new predictor is very promising, meaning that many important key features, which are deeply hidden in complicated protein sequences, can be extracted via the wavelets transform approach, quite consistent with the facts that many important biological functions of proteins can be elucidated with their low-frequency internal motions. The web server of iPPBS-PseAAC is accessible at http://www.jci-bioinfo.cn/iPPBS-PseAAC , by which users can easily acquire their desired results without the need to follow the complicated mathematical equations involved.

Homology modeling, functional annotation and comparative genomics of outer membrane protein H of Pasteurella multocida

Pasteurella multocida is an important pathogen of animals and humans. Outer Membrane Protein (Omp) H is a major conserved protein in the envelope of P. multocida and has been commonly targeted as a protective antigen. However, not much is known about its structure and function due to the difficulties that are typically associated with obtaining sufficient amounts of purified prokaryotic transmembrane proteins. The present work is aimed at studying the OmpH using an in silico approach and consolidate the findings in light of existing experimental evidences. Our study describes the first 3D model of the P. multocida OmpH obtained through a combination of several in silico modeling approaches. From our results, OmpH of P. multocida could be classified as a homotrimeric, 16 stranded, β-barrel porin involved in the non-specific transport of small, hydrophilic molecules, serving essential osmoregulatory function. Moreover, very small homologous sequences could be identified in the host proteome, strengthening the probability of a successful OmpH-based vaccine against the pathogen with remote chances of cross-reaction to host proteins.

Exploring the structure and conformational landscape of human leptin. A molecular dynamics approach

Leptin is a hormone that regulates energy homeostasis, inflammation, hematopoiesis and immune response, among other functions (Houseknecht et al., 1998; Zhang et al., 1995; Paz-Filho et al., 2010). To obtain its crystallographic structure, it was necessary to substitute a tryptophan for a glutamic acid at position 100, thus creating a mutant leptin that has been reported to have biological activity comparable to the activity of the wild type but that crystallizes more readily. Here, we report a comparative study of the conformational space of WT and W100E leptin using molecular dynamics simulations performed at 300, 400, and 500 K. We detected differences between the interactions of the two proteins with local and distal effects, resulting in changes in the conformation, accessible surface area, compactness, electrostatic potential and dynamic behavior. Additionally, the series of unfolding events that occur when leptin is subjected to high temperature differs for the two constructs. We observed that both proteins are mostly unstructured after 20 ns of MD simulation at 500 K. However, WT leptin maintains a significant amount of secondary structure in helix α2, while the most stable region of W100E leptin is helix α3. Furthermore, we found that the region between residues 25 and 42 might adopt interconverting secondary structures ranging from α-helices and random coils to β-strand structures. Thus, this region can be considered an intrinsically disordered region. This atomistic description supports our understanding of leptin signaling and consequently might facilitate the use of leptin in treatments for the pathophysiologies in which it is implicated.

Computational assay of Zanamivir binding affinity with original and mutant influenza neuraminidase 9 using molecular docking

Based upon molecular docking, this study aimed to find notable in silico neuraminidase 9 (NA9) point mutations of the avian influenza A H7N9 virus that possess a Zanamivir resistant property and to determine the lead compound capable of inhibiting these NA9 mutations. Seven amino acids (key residues) at the binding site of neuraminidase 9 responsible for Zanamivir-NA9 direct interactions were identified and 72 commonly occurring mutant NA9 versions were created using the Sybyl-X 2.0 software. The docking scores obtained after Zanamivir was bound to all mutant molecules of NA9 revealed 3 notable mutations R292W, R118P, and R292K that could greatly reduce the binding affinity of the medicine. These 3 mutant NA9 versions were then bound to each of 154 different molecules chosen from 5 groups of compounds to determine which molecule(s) might be capable of inhibiting mutant neuraminidase 9, leading to the discovery of the lead compound of potent mutant NA9 inhibitors. This compound, together with other mutations occurring to NA9 identified in the study, would be used as data for further research regarding neuraminidase inhibitors and synthesizing new viable medications used in the fight against the virus.

Insight into a molecular interaction force supporting peptide backbones and its implication to protein loops and folding

Although not being classified as the most fundamental protein structural elements like α-helices and β-strands, the loop segment may play considerable roles for protein stability, flexibility, and dynamic activity. Meanwhile, the protein loop is also quite elusive; i.e. its interactions with the other parts of protein as well as its own shape-maintaining forces have still remained as a puzzle or at least not quite clear yet. Here, we report a molecular force, the so-called polar hydrogen-π interaction (Hp-π), which may play an important role in supporting the backbones of protein loops. By conducting the potential energy surface scanning calculations on the quasi π-plane of peptide bond unit, we have observed the following intriguing phenomena: (1) when the polar hydrogen atom of a peptide unit is perpendicularly pointing to the π-plane of other peptide bond units, a remarkable Hp-π interaction occurs; (2) the interaction is distance and orientation dependent, acting in a broad space, and belonging to the 'point-to-plane' one. The molecular force reported here may provide useful interaction concepts and insights into better understanding the loop's unique stability and flexibility feature, as well as the driving force of the protein global folding.

Thermostable chitinase II from Thermomyces lanuginosus SSBP: Cloning, structure prediction and molecular dynamics simulations

Thermomyces lanuginosus is a thermophilic fungus that produces large number of industrially-significant enzymes owing to their inherent stability at high temperatures and wide range of pH optima, including thermostable chitinases that have not been fully characterized. Here, we report cloning, characterization and structure prediction of a gene encoding thermostable chitinase II. Sequence analysis revealed that chitinase II gene encodes a 343 amino acid protein of molecular weight 36.65kDa. Our study reports that chitinase II exhibits a well-defined TIM-barrel topology with an eight-stranded α/β domain. Structural analysis and molecular docking studies suggested that Glu176 is essential for enzyme activity. Folding studies of chitinase II using molecular dynamics simulations clearly demonstrated that the stability of the protein was evenly distributed at 350K.

A comparative study of structural and conformational properties of casein kinase-1 isoforms: insights from molecular dynamics and principal component analysis

Wnt signaling pathway regulates several developmental processes in human; however recently this pathway has been associated with development of different types of cancers. Casein kinase-1 (CK1) constitutes a family of serine-threonine protein kinase; various members of this family participate in Wnt signal transduction pathway and serve as molecular switch to this pathway. Among the known six isoforms of CK1, in human, at least three isoforms (viz. alpha, delta and epsilon) have been reported as oncogenic. The development of common therapeutics against these kinases is an arduous task; unless we have the detailed information of their tertiary structures and conformational properties. In the present work, the dynamical and conformational properties for each of three isoforms of CK1 are explored through molecular dynamics (MD) simulations. The conformational space distribution of backbone atoms is evaluated using principal component analysis of MD data, which are further validated on the basis of potential energy surface. Based on these analytics, it is suggested that conformational subspace shifts upon binding to ligands and guides the kinase action of CK1 isoforms. Further, this paper as a first effort to concurrently study all the three isoforms of CK1 provides structural basis for development of common anticancer therapeutics against three isoforms of CK1.

QSAR prediction of HIV-1 protease inhibitory activities using docking derived molecular descriptors

In this study, application of a new hybrid docking-quantitative structure activity relationship (QSAR) methodology to model and predict the HIV-1 protease inhibitory activities of a series of newly synthesized chemicals is reported. This hybrid docking-QSAR approach can provide valuable information about the most important chemical and structural features of the ligands that affect their inhibitory activities. Docking studies were used to find the actual conformations of chemicals in active site of HIV-1 protease. Then the molecular descriptors were calculated from these conformations. Multiple linear regression (MLR) and least square support vector machine (LS-SVM) were used as QSAR models, respectively. The obtained results reveal that statistical parameters of the LS-SVM model are better than the MLR model, which indicate that there are some non-linear relations between selected molecular descriptors and anti-HIV activities of interested chemicals. The correlation coefficient (R), root mean square error (RMSE) and average absolute error (AAE) for LS-SVM are: R=0.988, RMSE=0.207 and AAE=0.145 for the training set, and R=0.965, RMSE=0.403 and AAE=0.338 for the test set. Leave one out cross validation test was used for assessment of the predictive power and validity of models which led to cross-validation correlation coefficient QUOTE of 0.864 and 0.850 and standardized predicted relative error sum of squares (SPRESS) of 0.553 and 0.581 for LS-SVM and MLR models, respectively.

Bioassay directed isolation and biological evaluation of compounds isolated from Rubus fairholmianus Gard

The in vitro and in silico analysis of Rubus fairholmianus acetone extract for antioxidant, antiproliferative, and anti-inflammatory activity led to the isolation of six compounds. Amongst all the six isolated compounds tested, 1-(2-hydroxyphenyl)-4-methylpentan-1-one (compound 1) and 2-[(3-methylbutoxy) carbonyl] benzoic acid (compound 2) were found to be more active in inhibiting BRCA and COX target proteins, which also showed the better results for DPPH and ABTS radical scavenging assays. The promising results of this investigation emphasize the importance of using R. fairholmianus in the treatment of radical generated disorders mainly cancer and other inflammatory diseases.

Sequence-specific flexibility organization of splicing flanking sequence and prediction of splice sites in the human genome

More and more reported results of nucleosome positioning and histone modifications showed that DNA structure play a well-established role in splicing. In this study, a set of DNA geometric flexibility parameters originated from molecular dynamics (MD) simulations were introduced to discuss the structure organization around splice sites at the DNA level. The obtained profiles of specific flexibility/stiffness around splice sites indicated that the DNA physical-geometry deformation could be used as an alternative way to describe the splicing junction region. In combination with structural flexibility as discriminatory parameter, we developed a hybrid computational model for predicting potential splicing sites. And the better prediction performance was achieved when the benchmark dataset evaluated. Our results showed that the mechanical deformability character of a splice junction is closely correlated with both the splice site strength and structural information in its flanking sequences.

In SilicoInvestigation of Cytochrome P450 2C9 in relation to Aging Using Traditional Chinese Medicine

Cytochrome P450 2C9 (CYP2C9) metabolizes dehydroepiandrosterone-sulfate (DHEA-S), but in elderly people the amount of DHEA-S remaining after CYP2C9 metabolization may be insufficient for optimal health. A prediction model, molecular docking, and molecular dynamics were used to screen the Traditional Chinese Medicine (TCM) database to determine molecular compounds that may inhibit CYP2C9. The candidate compounds apocynoside(I), 4-methoxymagndialdehyde, and prunasin have higher Dock Scores, and prediction bioactivity than warfarin (the control drug). The interaction between 4-methoxymagndialdehyde and CYP2C9 is more intense than with other TCM compounds, but the simulation is longer. In these compounds, apocynoside(I) and prunasin have a greater number of pathways for their flexible structure, but these structures create weak interactions. These candidate compounds, which are known to have antioxidation and hypolipidemic functions that have an indirect effect on the aging process, can be extracted from traditional Chinese medicines. Thus, these candidate compounds may become CYP2C9 inhibitors and play an important role in providing optimal health in the elderly.

Molecular modeling study for inhibition mechanism of human chymase and its application in inhibitor design

Human chymase catalyzes the hydrolysis of peptide bonds. Three chymase inhibitors with very similar chemical structures but highly different inhibitory profiles towards the hydrolase function of chymase were selected with the aim of elucidating the origin of disparities in their biological activities. As a substrate (angiotensin-I) bound crystal structure is not available, molecular docking was performed to dock the substrate into the active site. Molecular dynamics simulations of chymase complexes with inhibitors and substrate were performed to calculate the binding orientation of inhibitors and substrate as well as to characterize conformational changes in the active site. The results elucidate details of the 3D chymase structure as well as the importance of K40 in hydrolase function. Binding mode analysis showed that substitution of a heavier Cl atom at the phenyl ring of most active inhibitor produced a great deal of variation in its orientation causing the phosphinate group to interact strongly with residue K40. Dynamics simulations revealed the conformational variation in region of V36-F41 upon substrate and inhibitor binding induced a shift in the location of K40 thus changing its interactions with them. Chymase complexes with the most active compound and substrate were used for development of a hybrid pharmacophore model which was applied in databases screening. Finally, hits which bound well at the active site, exhibited key interactions and favorable electronic properties were identified as possible inhibitors for chymase. This study not only elucidates inhibitory mechanism of chymase inhibitors but also provides key structural insights which will aid in the rational design of novel potent inhibitors of the enzyme. In general, the strategy applied in the current study could be a promising computational approach and may be generally applicable to drug design for other enzymes.

Structural origins for the loss of catalytic activities of bifunctional human LTA4H revealed through molecular dynamics simulations

Human leukotriene A4 hydrolase (hLTA4H), which is the final and rate-limiting enzyme of arachidonic acid pathway, converts the unstable epoxide LTA4 to a proinflammatory lipid mediator LTB4 through its hydrolase function. The LTA4H is a bi-functional enzyme that also exhibits aminopeptidase activity with a preference over arginyl tripeptides. Various mutations including E271Q, R563A, and K565A have completely or partially abolished both the functions of this enzyme. The crystal structures with these mutations have not shown any structural changes to address the loss of functions. Molecular dynamics simulations of LTA4 and tripeptide complex structures with functional mutations were performed to investigate the structural and conformation changes that scripts the observed differences in catalytic functions. The observed protein-ligand hydrogen bonds and distances between the important catalytic components have correlated well with the experimental results. This study also confirms based on the structural observation that E271 is very important for both the functions as it holds the catalytic metal ion at its location for the catalysis and it also acts as N-terminal recognition residue during peptide binding. The comparison of binding modes of substrates revealed the structural changes explaining the importance of R563 and K565 residues and the required alignment of substrate at the active site. The results of this study provide valuable information to be utilized in designing potent hLTA4H inhibitors as anti-inflammatory agents.

Design novel dual agonists for treating type-2 diabetes by targeting peroxisome proliferator-activated receptors with core hopping approach

Owing to their unique functions in regulating glucose, lipid and cholesterol metabolism, PPARs (peroxisome proliferator-activated receptors) have drawn special attention for developing drugs to treat type-2 diabetes. By combining the lipid benefit of PPAR-alpha agonists (such as fibrates) with the glycemic advantages of the PPAR-gamma agonists (such as thiazolidinediones), the dual PPAR agonists approach can both improve the metabolic effects and minimize the side effects caused by either agent alone, and hence has become a promising strategy for designing effective drugs against type-2 diabetes. In this study, by means of the powerful “core hopping” and “glide docking” techniques, a novel class of PPAR dual agonists was discovered based on the compound GW409544, a well-known dual agonist for both PPAR-alpha and PPAR-gamma modified from the farglitazar structure. It was observed by molecular dynamics simulations that these novel agonists not only possessed the same function as GW409544 did in activating PPAR-alpha and PPAR-gamma, but also had more favorable conformation for binding to the two receptors. It was further validated by the outcomes of their ADME (absorption, distribution, metabolism, and excretion) predictions that the new agonists hold high potential to become drug candidates. Or at the very least, the findings reported here may stimulate new strategy or provide useful insights for discovering more effective dual agonists for treating type-2 diabetes. Since the “core hopping” technique allows for rapidly screening novel cores to help overcome unwanted properties by generating new lead compounds with improved core properties, it has not escaped our notice that the current strategy along with the corresponding computational procedures can also be utilized to find novel and more effective drugs for treating other illnesses.

Insights into the mutation-induced HHH syndrome from modeling human mitochondrial ornithine transporter-1

Human mitochondrial ornithine transporter-1 is reported in coupling with the hyperornithinemia-hyperammonemia-homocitrullinuria (HHH) syndrome, which is a rare autosomal recessive disorder. For in-depth understanding of the molecular mechanism of the disease, it is crucially important to acquire the 3D structure of human mitochondrial ornithine transporter-1. Since no such structure is available in the current protein structure database, we have developed it via computational approaches based on the recent NMR structure of human mitochondrial uncoupling protein (Berardi MJ, Chou JJ, et al. Nature 2011, 476:109–113). Subsequently, we docked the ligand L-ornithine into the computational structure to search for the favorable binding mode. It was observed that the binding interaction for the most favorable binding mode is featured by six remarkable hydrogen bonds between the receptor and ligand, and that the most favorable binding mode shared the same ligand-binding site with most of the homologous mitochondrial carriers from different organisms, implying that the ligand-binding sites are quite conservative in the mitochondrial carriers family although their sequences similarity is very low with 20% or so. Moreover, according to our structural analysis, the relationship between the disease-causing mutations of human mitochondrial ornithine transporter-1 and the HHH syndrome can be classified into the following three categories: (i) the mutation occurs in the pseudo-repeat regions so as to change the region of the protein closer to the mitochondrial matrix; (ii) the mutation is directly affecting the substrate binding pocket so as to reduce the substrate binding affinity; (iii) the mutation is located in the structural region closer to the intermembrane space that can significantly break the salt bridge networks of the protein. These findings may provide useful insights for in-depth understanding of the molecular mechanism of the HHH syndrome and developing effective drugs against the disease.

Novel inhibitor design for hemagglutinin against H1N1 influenza virus by core hopping method

The worldwide spread of H1N1 avian influenza and the increasing reports about its resistance to the current drugs have made a high priority for developing new anti-influenza drugs. Owing to its unique function in assisting viruses to bind the cellular surface, a key step for them to subsequently penetrate into the infected cell, hemagglutinin (HA) has become one of the main targets for drug design against influenza virus. To develop potent HA inhibitors, the ZINC fragment database was searched for finding the optimal compound with the core hopping technique. As a result, the Neo6 compound was obtained. It has been shown through the subsequent molecular docking studies and molecular dynamic simulations that Neo6 not only assumes more favorable conformation at the binding pocket of HA but also has stronger binding interaction with its receptor. Accordingly, Neo6 may become a promising candidate for developing new and more powerful drugs for treating influenza. Or at the very least, the findings reported here may provide useful insights to stimulate new strategy in this area.