Computer modelling studies of ribonuclease A-pyrimidine nucleotide complexes

Different modes of binding of pyrimidine monophosphates 2'-UMP, 3'-UMP, 2'-CMP and 3'-CMP to ribonuclease (RNase) A are studied by energy minimization in torsion angle and subsequently in Cartesian coordinate space. The results are analysed in the light of primary binding sites. The hydrogen bonding pattern brings out roles for amino acids such as Asn44 and Ser123 apart from the well known active site residues viz., His12,Lys41,Thr45 and His119. Amino acid segments 43-45 and 119-121 seem to be guiding the ligand binding by forming a pocket. Many of the active site charged residues display considerable movement upon nucleotide binding.

Computer modeling studies on the binding of 2',5'-linked dinucleoside phosphates to ribonuclease T1-influence of subsite interactions on the substrate specificity

The modes of binding of Gp(2',5')A, Gp(2',5')C, Gp(2',5')G and Gp(2',5')U to RNase T1 have been determined by computer modelling studies. All these dinucleoside phosphates assume extended conformations in the active site leading to better interactions with the enzyme. The 5'-terminal guanine of all these ligands is placed in the primary base binding site of the enzyme in an orientation similar to that of 2'-GMP in the RNase T1-2'-GMP complex. The 2'-terminal purines are placed close to the hydrophobic pocket formed by the residues Gly71, Ser72, Pro73 and Gly74 which occur in a loop region. However, the orientation of the 2'-terminal pyrimidines is different from that of 2'-terminal purines. This perhaps explains the higher binding affinity of the 2',5'-linked guanine dinucleoside phosphates with 2'-terminal purines than those with 2'-terminal pyrimidines. A comparison of the binding of the guanine dinucleoside phosphates with 2',5'- and 3',5'-linkages suggests significant differences in the ribose pucker and hydrogen bonding interactions between the catalytic residues and the bound nucleoside phosphate implying that 2',5'-linked dinucleoside phosphates may not be the ideal ligands to probe the role of the catalytic amino acid residues. A change in the amino acid sequence in the surface loop region formed by the residues Gly71 to Gly74 drastically affects the conformation of the base binding subsite, and this may account for the inactivity of the enzyme with altered sequence i.e., with Pro, Gly and Ser at positions 71 to 73 respectively. These results thus suggest that in addition to recognition and catalytic sites, interactions at the loop regions which constitute the subsite for base binding are also crucial in determining the substrate specificity.

Effect of peptide-based captopril analogues on angiotensin converting enzyme activity and peroxynitrite-mediated tyrosine nitration

Angiotensin converting enzyme (ACE) regulates the blood pressure by converting angiotensin I to angiotensin II and bradykinin to bradykinin 1-7. These two reactions elevate the blood pressure as angiotensin II and bradykinin are vasoconstrictory and vasodilatory hormones, respectively. Therefore, inhibition of ACE is an important strategy for the treatment of hypertension. The natural substrates of ACE, i.e., angiotensin II and bradykinin, contain a Pro-Phe motif near the site of hydrolysis. Therefore, there may be a Pro-Phe binding pocket at the active site of ACE, which may facilitate the substrate binding. In view of this, we have synthesized a series of thiol-and selenol-containing dipeptides and captopril analogues and studied their ACE inhibition activities. This study reveals that both the selenol or thiol moiety and proline residues are essential for ACE inhibition. Although the introduction of a Phe residue to captopril and its selenium analogue considerably reduces the inhibitory effect, there appears to be a Phe binding pocket at the active site of ACE.

DEPTH: a web server to compute depth and predict small-molecule binding cavities in proteins

Depth measures the extent of atom/residue burial within a protein. It correlates with properties such as protein stability, hydrogen exchange rate, protein-protein interaction hot spots, post-translational modification sites and sequence variability. Our server, DEPTH, accurately computes depth and solvent-accessible surface area (SASA) values. We show that depth can be used to predict small molecule ligand binding cavities in proteins. Often, some of the residues lining a ligand binding cavity are both deep and solvent exposed. Using the depth-SASA pair values for a residue, its likelihood to form part of a small molecule binding cavity is estimated. The parameters of the method were calibrated over a training set of 900 high-resolution X-ray crystal structures of single-domain proteins bound to small molecules (molecular weight < 1.5 KDa). The prediction accuracy of DEPTH is comparable to that of other geometry-based prediction methods including LIGSITE, SURFNET and Pocket-Finder (all with Matthew's correlation coefficient of similar to 0.4) over a testing set of 225 single and multi-chain protein structures. Users have the option of tuning several parameters to detect cavities of different sizes, for example, geometrically flat binding sites. The input to the server is a protein 3D structure in PDB format. The users have the option of tuning the values of four parameters associated with the computation of residue depth and the prediction of binding cavities. The computed depths, SASA and binding cavity predictions are displayed in 2D plots and mapped onto 3D representations of the protein structure using Jmol. Links are provided to download the outputs. Our server is useful for all structural analysis based on residue depth and SASA, such as guiding site-directed mutagenesis experiments and small molecule docking exercises, in the context of protein functional annotation and drug discovery.

Analysis of laminates with ply drops

Thickness tapered laminates obtained by terminating a certain number of plies contain resin-rich areas called ‘resin pockets’ near ply drops, where high stress concentrations exist. Study of the effects of ply drops and resin pockets on the tensile behaviour of tapered laminates considering certain important parameters like taper angle, the number of plies dropped, and the fiber orientation is reported here. Estimation of the tensile strength of tapered laminates necessitates accurate determination of the state of stress near the ply-drop region, which is, in general, three-dimensional (3-D) in nature. Recognising the fact that full 3-D finite-element analysis becomes computationally exorbitant, special layered 3-D finite-element analysis is carried out. Laminates with ply drops along only one direction are analysed to elicit the nature of the local bending effects occurring near the ply drops. Complete 3-D Tsai–Wu criterion considering all the six stress components is used to obtain a quick and comparative assessment of the tensile strength of these laminates. High stress concentration zones are identified and the effects of number of plies dropped at a station and resin pocket geometry are illustrated. The mechanism of load transfer near ply drops and the local bending that occurs are described. Susceptibility of ply drop zones to the onset and subsequent growth of delaminations is also brought out.

Identification and characterization of cysteine proteinases of Trypanosoma evansi

Trypanosoma evansi is a causative agent of `surra', a common haemoprotozoan disease of livestock in India causing high morbidity and mortality in disease endemic areas. The proteinases released by live and dead trypanosomes entail immunosuppression in the infected host, which immensely contribute in disease pathogenesis. Cysteine proteinases are identified in the infectious cycle of trypanosomes such as cruzain from Trypanosoma cruzi, rhodesain or brucipain from Trypanosoma brucei rhodesiense and congopain from Trypanosoma congelense. These enzymes localised in lysosome-like organelles, flagellar pocket and on cell surface, which play a critical role in the life cycle of protozoan parasites, viz. in host invasion, nutrition and alteration of the host immune response. The paper describes the identification of cysteine proteinases of T. evansi lysate, activity profile at different pH optima and inhibition pattern using a specific inhibitor, besides the polypeptide profile of an antigen. Eight proteinases of T. evansi were identified in the molecular weight (MW) ranges of 28-170 kDa using gelatin substrate-polyacrylamide gel electrophoresis (GS-PAGE), and of these proteinases, six were cysteine proteinases, as they were inhibited by L-3-carboxy-2,3-transepoxypropionyl-lecuylamido (4-guanidino)-butane (E-64), a specific inhibitor. These proteolytic enzymes were most reactive in acidic pH between 3.0 and 5.5 in the presence of dithiothreitol and completely inactive at alkaline pH 10.0. Similarly, the GS-PAGE profile of the serum samples of rats infected with T. evansi revealed strong proteolytic activity only at the 28-kDa zone at pH 5.5, while no proteolytic activity was observed in serum samples of uninfected rats. Further, the other zones of clearance, which were evident in T. evansi antigen zymogram, could not be observed in the serum samples of rats infected with T. evansi. The polypeptide pattern of the whole cell lysate antigen revealed 12-15 polypeptide bands ranging from 28 to 81 kDa along with five predominant polypeptides bands (MW of 81, 66, 62, 55 and 45 kDa), which were immunoreactive with hyperimmune serum (HIS) and serum of experimentally infected rabbits with T. evansi infection. The immunoblot recognised antibodies in experimentally infected rabbits and against HIS as well, corresponding to the zone of clearances at lower MW ranges (28-41 kDa), which may be attributed to the potential of these proteinases in the diagnosis of T. evansi infection. Since these thiol-dependent enzymes are most active in acidic pH and considering their inhibition characteristics, these data suggest that they resemble to the mammalian lysosomal cathepsin B and L.

Functional regulation of PVBV Nuclear Inclusion protein-a protease activity upon interaction with Viral Protein genome-linked and phosphorylation

Regulation of NIa-Pro is crucial for polyprotein processing and hence, for successful infection of potyviruses. We have examined two novel mechanisms that could regulate NIa-Pro activity. Firstly, the influence of VPg domain on the proteolytic activity of NIa-Pro was investigated. It was shown that the turnover number of the protease increases when these two domains interact (as: two-fold; trans: seven-fold) with each other. Secondly, the protease activity of NIa-Pro could also be modulated by phosphorylation at Ser129. A mutation of this residue either to aspartate (phosphorylation-mimic) or alanine (phosphorylation-deficient) drastically reduces the protease activity. Based on these observations and molecular modeling studies, we propose that interaction with VPg as well as phosphorylation of Ser129 could relay a signal through Trp143 present at the protein surface to the active site pocket by subtle conformational changes, thus modulating protease activity of NIa-Pro. (C) 2011 Elsevier Inc. All rights reserved.

Sustainability of apparent slip in micro-channel flows

In the present work, we experimentally study the flow of water over textured hydrophobic surfaces in a micro-channel. Shear stress measurements are done along with direct visualization of trapped air pockets on the hydrophobic surface. The trapped air pockets on such surfaces are known to be responsible for apparent slip at these surfaces and hence in significant drag reduction. In typical circumstances, the apparent slip reduces over time as seen, for example, from our shear stress measurements. This implies that the drag reduction will not be sustained. We have performed extensive visualizations of the trapped air pockets while varying flow parameters like the flow rate and the pressure. We present here direct visualizations that show that under some conditions, the air pockets can grow with time. The variation of the air pocket size with time is found to change qualitatively and quantitatively as the flow rate is varied. These measured changes in the air pocket size with time have a direct bearing on the sustainability of apparent slip in micro-channel flows.

Crystal structures and binding studies of atovaquone and its derivatives with cytochrome bc(1): a molecular basis for drug design

Crystal structure of trans-atovaquone (antimalarial drug), its polymorph and its stereoisomer (cis) along with five other derivatives with different functional groups have been analyzed. Based on the conformational features of these compounds and the characteristics of the nature of intermolecular interactions, valuable insights into the atomistic details of protein-inhibitor interactions have been derived by docking studies. Atovaquone and its derivatives pack in the crystal lattice using intermolecular O-H center dot center dot center dot O hydrogen bond dimer motifs supported by surrogate weak interactions including C-H center dot center dot center dot O and C-H center dot center dot center dot Cl hydrogen bonds. The docking results of these molecules with cytochrome bc(1) show preferences to form N-H center dot center dot center dot O, O-H center dot center dot center dot O and O-H center dot center dot center dot Cl hydrogen bonds. The involvement of halogen atoms in the binding pocket appears to be significant and is contrary to the theoretically predicted mechanism of protein-ligand docking reported earlier based on mimicking experimental binding results of stigmatellin with cytochrome bc(1). The significance of subtle energy factors controlled by weak intermolecular interactions appears to play a major role in drug binding.

A multi-level multi-scale approach to study essential genes in Mycobacterium tuberculosis

Background: The set of indispensable genes that are required by an organism to grow and sustain life are termed as essential genes. There is a strong interest in identification of the set of essential genes, particularly in pathogens, not only for a better understanding of the pathogen biology, but also for identifying drug targets and the minimal gene set for the organism. Essentiality is inherently a systems property and requires consideration of the system as a whole for their identification. The available experimental approaches capture some aspects but each method comes with its own limitations. Moreover, they do not explain the basis for essentiality in most cases. A powerful prediction method to recognize this gene pool including rationalization of the known essential genes in a given organism would be very useful. Here we describe a multi-level multi-scale approach to identify the essential gene pool in a deadly pathogen, Mycobacterium tuberculosis. Results: The multi-level workflow analyses the bacterial cell by studying (a) genome-wide gene expression profiles to identify the set of genes which show consistent and significant levels of expression in multiple samples of the same condition, (b) indispensability for growth by using gene expression integrated flux balance analysis of a genome-scale metabolic model, (c) importance for maintaining the integrity and flow in a protein-protein interaction network and (d) evolutionary conservation in a set of genomes of the same ecological niche. In the gene pool identified, the functional basis for essentiality has been addressed by studying residue level conservation and the sub-structure at the ligand binding pockets, from which essential amino acid residues in that pocket have also been identified. 283 genes were identified as essential genes with high-confidence. An agreement of about 73.5% is observed with that obtained from the experimental transposon mutagenesis technique. A large proportion of the identified genes belong to the class of intermediary metabolism and respiration. Conclusions: The multi-scale, multi-level approach described can be generally applied to other pathogens as well. The essential gene pool identified form a basis for designing experiments to probe their finer functional roles and also serve as a ready shortlist for identifying drug targets.

Conserved water-mediated recognition and dynamics of NAD(+) (carboxamide group) to hIMPDH enzyme: water mimic approach toward the design of isoform-selective inhibitor

Inosine monophosphate dehydrogenase (IMPDH) enzyme involves in GMP biosynthesis pathway. Type I hIMPDH is expressed at lower levels in all cells, whereas type II is especially observed in acute myelogenous leukemia, chronic myelogenous leukemia cancer cells, and 10 ns simulation of the IMP-NAD(+) complex structures (PDB ID. 1B3O and 1JCN) have revealed the presence of a few conserved hydrophilic centers near carboxamide group of NAD(+). Three conserved water molecules (W1, W, and W1 `) in di-nucleotide binding pocket of enzyme have played a significant role in the recognition of carboxamide group (of NAD(+)) to D274 and H93 residues. Based on H-bonding interaction of conserved hydrophilic (water molecular) centers within IMP-NAD(+)-enzyme complexes and their recognition to NAD(+), some covalent modification at carboxamide group of di-nucleotide (NAD(+)) has been made by substituting the -CONH(2)group by -CONHNH2 (carboxyl hydrazide group) using water mimic inhibitor design protocol. The modeled structure of modified ligand may, though, be useful for the development of antileukemic agent or it could be act as better inhibitor for hIMPDH-II.

PLIC: protein-ligand interaction clusters

Most of the biological processes are governed through specific protein-ligand interactions. Discerning different components that contribute toward a favorable protein-ligand interaction could contribute significantly toward better understanding protein function, rationalizing drug design and obtaining design principles for protein engineering. The Protein Data Bank (PDB) currently hosts the structure of similar to 68 000 protein-ligand complexes. Although several databases exist that classify proteins according to sequence and structure, a mere handful of them annotate and classify protein-ligand interactions and provide information on different attributes of molecular recognition. In this study, an exhaustive comparison of all the biologically relevant ligand-binding sites (84 846 sites) has been conducted using PocketMatch: a rapid, parallel, in-house algorithm. PocketMatch quantifies the similarity between binding sites based on structural descriptors and residue attributes. A similarity network was constructed using binding sites whose PocketMatch scores exceeded a high similarity threshold (0.80). The binding site similarity network was clustered into discrete sets of similar sites using the Markov clustering (MCL) algorithm. Furthermore, various computational tools have been used to study different attributes of interactions within the individual clusters. The attributes can be roughly divided into (i) binding site characteristics including pocket shape, nature of residues and interaction profiles with different kinds of atomic probes, (ii) atomic contacts consisting of various types of polar, hydrophobic and aromatic contacts along with binding site water molecules that could play crucial roles in protein-ligand interactions and (iii) binding energetics involved in interactions derived from scoring functions developed for docking. For each ligand-binding site in each protein in the PDB, site similarity information, clusters they belong to and description of site attributes are provided as a relational database-protein-ligand interaction clusters (PLIC).

Characterizing the pocketome of Mycobacterium tuberculosis and application in rationalizing polypharmacological target selection

Polypharmacology is beginning to emerge as an important concept in the field of drug discovery. However, there are no established approaches to either select appropriate target sets or design polypharmacological drugs. Here, we propose a structural-proteomics approach that utilizes the structural information of the binding sites at a genome-scale obtained through in-house algorithms to characterize the pocketome, yielding a list of ligands that can participate in various biochemical events in the mycobacterial cell. The pocket-type space is seen to be much larger than the sequence or fold-space, suggesting that variations at the site-level contribute significantly to functional repertoire of the organism. All-pair comparisons of binding sites within Mycobacterium tuberculosis (Mtb), pocket-similarity network construction and clustering result in identification of binding-site sets, each containing a group of similar binding sites, theoretically having a potential to interact with a common set of compounds. A polypharmacology index is formulated to rank targets by incorporating a measure of druggability and similarity to other pockets within the proteome. This study presents a rational approach to identify targets with polypharmacological potential along with possible drugs for repurposing, while simultaneously, obtaining clues on lead compounds for use in new drug-discovery pipelines.

Hydrodynamic instability and shear layer effect on the response of an acoustically excited laminar premixed flame

Combustion instabilities can cause serious problems which limit the operating envelope of low-emission lean premixed combustion systems. Predicting the onset of combustion instability requires a description of the unsteady heat release driving the instability, i.e., the heat release response transfer function of the system. This study focuses on the analysis of fully coupled two-way interactions between a disturbance field and a laminar premixed flame that incorporates gas expansion effects by solving the conservation equations of a compressible fluid. Results of the minimum and maximum flame front deflections are presented to underline the impact of the hydrodynamic instability on the flame and the shear layer effect on the initial flame front wrinkling which is increased at decreasing gas expansion. These phenomena influence the magnitude of the burning area and burning area rate response of the flame at lower frequency excitation more drastically than reduced-order model (ROM) predictions even for low temperature ratios. It is shown that the general trend of the flame response magnitudes can be well captured at higher frequency excitation, where stretch effects are dominant. The phase response is influenced by the DL mechanism, which cannot be captured by the ROM, and by the resulting discrepancy in the flame pocket formation and annihilation process at the flame tip. (C) 2014 The Combustion Institute. Published by Elsevier Inc. All rights reserved,

On the impact of a concave nosed axisymmetric body on a free surface

We report on an experimental study of the vertical impact of a concave nosed axisymmetric body on a free surface. Previous studies have shown that bodies with a convex nose, like a sphere, produce a well defined splash with a relatively large cavity behind the model. In contrast, we find that with a concave nose, there is hardly a splash and the cavity extent is greatly reduced. This may be explained by the fact that in the concave nosed case, the initial impact is between a confined air pocket and the free surface unlike in the convex nosed case. From measurements of the unsteady pressure in the concave nose portion, we show that in this case, the maximum pressures are significantly lower than the classically expected ``water hammer'' pressures and also lower than those generally measured on other geometries. Thus, the presence of an air pocket in the case of a concave nosed body adds an interesting dimension to the classical problem of impact of solid bodies on to a free surface. (C) 2015 AIP Publishing LLC.