Computational systems approach for drug target discovery

Importance of the field: The shift in focus from ligand based design approaches to target based discovery over the last two to three decades has been a major milestone in drug discovery research. Currently, it is witnessing another major paradigm shift by leaning towards the holistic systems based approaches rather the reductionist single molecule based methods. The effect of this new trend is likely to be felt strongly in terms of new strategies for therapeutic intervention, new targets individually and in combinations, and design of specific and safer drugs. Computational modeling and simulation form important constituents of new-age biology because they are essential to comprehend the large-scale data generated by high-throughput experiments and to generate hypotheses, which are typically iterated with experimental validation. Areas covered in this review: This review focuses on the repertoire of systems-level computational approaches currently available for target identification. The review starts with a discussion on levels of abstraction of biological systems and describes different modeling methodologies that are available for this purpose. The review then focuses on how such modeling and simulations can be applied for drug target discovery. Finally, it discusses methods for studying other important issues such as understanding targetability, identifying target combinations and predicting drug resistance, and considering them during the target identification stage itself. What the reader will gain: The reader will get an account of the various approaches for target discovery and the need for systems approaches, followed by an overview of the different modeling and simulation approaches that have been developed. An idea of the promise and limitations of the various approaches and perspectives for future development will also be obtained. Take home message: Systems thinking has now come of age enabling a `bird's eye view' of the biological systems under study, at the same time allowing us to `zoom in', where necessary, for a detailed description of individual components. A number of different methods available for computational modeling and simulation of biological systems can be used effectively for drug target discovery.

Cytotoxicity of half sandwich ruthenium(II) complexes with strong hydrogen bond acceptor ligands and their mechanism of action

Neutral and cationic organometallic ruthenium(II) piano stool complexes of the type [(eta(6)-cymene)R-uCl(X)(Y)] (complexes R1-R8) has been synthesized and characterized. In cationic complexes, X, Y is either a eta(2) phosphorus ligand such as 1,1-bis(diphenylphosphino)methane (DPPM) and 1,2-bis(diphenylphosphino)ethane (DPPE) or partially oxidized ligands such as 1,2-bis(diphenylphosphino)methane monooxide (DPPMO) and 1,2-bis(diphenylphosphino)ethane monooxide (DPPEO) which are strong hydrogen bond acceptors. In neutral complexes. X is chloride and Y is a monodentate phosphorous donor. Complexes with DPPM and DPPMO ligands ([(eta(6)-cymene)Ru(eta(2)-DPPM)Cl]PF6 (R2), [(eta(6)-cymene)Ru(eta(2)-DPPMO)Cl]PF6 (R3), [(eta(6)-cymene)Ru(eta(1)-DPPM)Cl-2] (R5) and [(eta(6)-cymene)Ru(eta(1)-DPPMO)Cl-2] (R6) show good cytotoxicity. Growth inhibition study of several human cancer cell lines by these complexes has been carried out. Mechanistic studies for R5 and R6 show that inhibition of cancer cell growth involves both cell cycle arrest and apoptosis induction. Using an apoptosis PCR array, we identified the sets of antiapoptotic genes that were down regulated and pro-apoptotic genes that were up regulated. These complexes were also found to be potent metastasis inhibitors as they prevented cell invasion through matrigel. The complexes were shown to bind DNA in a non intercalative fashion and cause unwinding of plasmid DNA in cell-free medium by competitive ethidium bromide binding, viscosity measurements, thermal denaturation and gel mobility shift assays.

Mycobacterial Stress Regulation: The Dps ``Twin Sister'' Defense Mechanism and Structure-Function Relationship

In this work, we have tried to emphasize the connection between mycobacterial growth and regulation of gene expression. Utilization of multiple carbon sources and diauxic growth helps bacteria to regulate gene expression at an optimum level so that the inhospitable conditions encountered during nutrient depletion can be circumvented. These aspects will be discussed with respect to mycobacterial growth in subsequent sections. Identification and characterization of genes induced under such conditions is helpful to understand the physiology of the bacterium. Although it is necessary to compare the total expression profile of proteins as they transit from vegetative growth to stationary phase, at times a lot of insights can be deciphered from the expression pattern of one or two proteins. We have compared the protein expression and sigma factor selectivity of two such proteins in M. smegmatis to understand the differential regulation of genes playing diverse function in the same species. Some newer insights on the structure and function of one of the Dps proteins are also explained.

X-ray crystallographic and NMR studies of pantothenate synthetase provide insights into the mechanism of homotropic inhibition by pantoate

The structural basis for the homotropic inhibition of pantothenate synthetase by the substrate pantoate was investigated by X-ray crystallography and high-resolution NMR spectroscopic methods. The tertiary structure of the dimeric N-terminal domain of Escherichia coli pantothenate synthetase, determined by X-ray crystallography to a resolution of 1.7 Å, showed a second molecule of pantoate bound in the ATP-binding pocket. Pantoate binding to the ATP-binding site induced large changes in structure, mainly for backbone and side chain atoms of residues in the ATP binding HXGH(34–37) motif. Sequence-specific NMR resonance assignments and solution secondary structure of the dimeric N-terminal domain, obtained using samples enriched in 2H, 13C, and 15N, indicated that the secondary structural elements were conserved in solution. Nitrogen-15 edited two-dimensional solution NMR chemical shift mapping experiments revealed that pantoate, at 10 mm, bound at these two independent sites. The solution NMR studies unambiguously demonstrated that ATP stoichiometrically displaced pantoate from the ATP-binding site. All NMR and X-ray studies were conducted at substrate concentrations used for enzymatic characterization of pantothenate synthetase from different sources [Jonczyk R & Genschel U (2006) J Biol Chem 281, 37435–37446]. As pantoate binding to its canonical site is structurally conserved, these results demonstrate that the observed homotropic effects of pantoate on pantothenate biosynthesis are caused by competitive binding of this substrate to the ATP-binding site. The results presented here have implications for the design and development of potential antibacterial and herbicidal agents.

Mechanism of failure in a free-standing Pt-aluminide bond coat during tensile testing at room temperature

The room temperature (RT) tensile behaviour of a free-standing high activity Pt-aluminide bond coat has been evaluated by microtensile testing technique. The coating had a typical three-layer microstructure. The stress-strain plot for the free-standing coating was linear, indicating the coating to be brittle at RT. Different fracture features were observed across the coating layers, namely quasi-cleavage in the outer layer and inner interdiffusion zone, and cleavage in the intermediate layer. By employing interrupted tensile test and observing the cross-sectional microstructure of the tested specimens, it was determined that failure of the microtensile samples occurred by the initiation of a single crack in the intermediate layer of the coating and its subsequent inside-out propagation. Such a mechanism of failure has been explained in terms of the fracture features observed across the sample thickness. This mechanism of failure is consistent with fracture toughness values of the individual coating layers. (C) 2009 Elsevier B.V. All rights reserved.

An adaptive drug delivery design using neural networks for effective treatment of infectious diseases: A simulation study

An adaptive drug delivery design is presented in this paper using neural networks for effective treatment of infectious diseases. The generic mathematical model used describes the coupled evolution of concentration of pathogens, plasma cells, antibodies and a numerical value that indicates the relative characteristic of a damaged organ due to the disease under the influence of external drugs. From a system theoretic point of view, the external drugs can be interpreted as control inputs, which can be designed based on control theoretic concepts. In this study, assuming a set of nominal parameters in the mathematical model, first a nonlinear controller (drug administration) is designed based on the principle of dynamic inversion. This nominal drug administration plan was found to be effective in curing "nominal model patients" (patients whose immunological dynamics conform to the mathematical model used for the control design exactly. However, it was found to be ineffective in curing "realistic model patients" (patients whose immunological dynamics may have off-nominal parameter values and possibly unwanted inputs) in general. Hence, to make the drug delivery dosage design more effective for realistic model patients, a model-following adaptive control design is carried out next by taking the help of neural networks, that are trained online. Simulation studies indicate that the adaptive controller proposed in this paper holds promise in killing the invading pathogens and healing the damaged organ even in the presence of parameter uncertainties and continued pathogen attack. Note that the computational requirements for computing the control are very minimal and all associated computations (including the training of neural networks) can be carried out online. However it assumes that the required diagnosis process can be carried out at a sufficient faster rate so that all the states are available for control computation.

A Bayesian incentive compatible mechanism for decentralized supply chain formation

In this paper we consider a decentralized supply chain formation problem for linear multi-echelon supply chains when the managers of the individual echelons are autonomous, rational, and intelligent. At each echelon, there is a choice of service providers and the specific problem we solve is that of determining a cost-optimal mix of service providers so as to achieve a desired level of end-to-end delivery performance. The problem can be broken up into two sub-problems following a mechanism design approach: (1) Design of an incentive compatible mechanism to elicit the true cost functions from the echelon managers; (2) Formulation and solution of an appropriate optimization problem using the true cost information. In this paper we propose a novel Bayesian incentive compatible mechanism for eliciting the true cost functions. This improves upon existing solutions in the literature which are all based on the classical Vickrey-Clarke-Groves mechanisms, requiring significant incentives to be paid to the echelon managers for achieving dominant strategy incentive compatibility. The proposed solution, which we call SCF-BIC (Supply Chain Formation with Bayesian Incentive Compatibility), significantly reduces the cost of supply chain formation. We illustrate the efficacy of the proposed methodology using the example of a three echelon manufacturing supply chain.

Prospecting for alternate sources of shikimic acid, a precursor of Tamiflu, a bird-flu drug

Shikimic acid, more commonly known by its anionic form, shikimate, is an important intermediate compound of the ‘shikimate pathway’ in plants and microorganisms1. It is the principal precursor for the synthesis of aromatic amino acids, phenylalanine, tryptophan and tyrosine and other compounds such as alkaloids, phenolics and phenyl propanoids2. It is used extensively as a chiral building block for the synthesis of a number of compounds in both pharmaceutical and cosmetic industries3. In the recent past, the focus on shikimic acid has increased since it is the key precursor for the synthesis of Tamiflu, the only drug against avian flu caused by the H5N1 virus4,5. Shikimic acid is converted to a diethyl ketal intermediate, which is then reduced in two steps to an epoxide that is finally transformed to Tamiflu6.

Structural Insights into the acyl-intermediates of plasmodium falciparum fatty acid synthesis pathway; the mechanism of expansion of acyl carrier protein core

Acyl carrier protein (ACP) plays a central role in fatty acid biosynthesis. However, the molecular machinery that mediates its function is not yet fully understood. Therefore, structural studies were carried out on the acyl-ACP intermediates of Plasmodium falciparum using NMR as a spectroscopic probe. Chemical shift perturbation studies put forth a new picture of the interaction of ACP molecule with the acyl chain, namely, the hydrophobic core can protect up to 12 carbon units, and additional carbons protrude out from the top of the hydrophobic cavity. The latter hypothesis stems from chemical shift changes observed in C-alpha and C-beta of Ser-37 in tetradecanoyl-ACP. C-13, N-15-Double-filtered nuclear Overhauser effect (NOE) spectroscopy experiments further substantiate the concept; in octanoyl (C-8)- and dodecanoyl (C-12)-ACP, a long range NOE is observed within the phosphopantetheine arm, suggesting an arch-like conformation. This NOE is nearly invisible in tetradecanoyl (C-14)-ACP, indicating a change in conformation of the prosthetic group. Furthermore, the present study provides insights into the molecular mechanism of ACP expansion, as revealed from a unique side chain-to-backbone hydrogen bond between two fairly conserved residues, Ile-55 HN and Glu-48 O. The backbone amide of Ile-55 HN reports a pK(a) value for the carboxylate, similar to 1.9 pH units higher than model compound value, suggesting strong electrostatic repulsion between helix II and helix III. Charge-charge repulsion between the helices in combination with thrust from inside due to acyl chain would energetically favor the separation of the two helices. Helix III has fewer structural restraints and, hence, undergoes major conformational change without altering the overall-fold of P. falciparum ACP.

Design of mechanism-based inhibitors of transthyretin amyloidosis:Studies with biphenyl ethers and new structural templates

Transthyretin (TTR), a tetrameric thyroxine (T4) carrier protein, is associated with a variety of amyloid diseases. In this study, we explore the potential of biphenyl ethers (BPE), which are shown to interact with a high affinity to its T4 binding site thereby preventing its aggregation and fibrillogenesis. They prevent fibrillogenesis by stabilizing the tetrameric ground state of transthyretin. Additionally, we identify two new structural templates (2-(5-mercapto-[1,3,4]oxadiazol-2-yl)-phenol and 2,3,6-trichloro-N-(4H-[1,2,4]triazol-3-yl) represented as compounds 11 and 12, respectively, throughout the manuscript) exhibiting the ability to arrest TTR amyloidosis. The dissociation constants for the binding of BPEs and compound 11 and 12 to TTR correlate with their efficacies of inhibiting amyloidosis. They also have the ability to inhibit the elongation of intermediate fibrils as well as show nearly complete (> 90%) disruption of the preformed fibrils. The present study thus establishes biphenyl ethers and compounds 11 and 12 as very potent inhibitors of TTR fibrillization and inducible cytotoxicity.

DNA looping and translocation provide an optimal cleavage mechanism for the type III restriction enzymes

EcoP15I is a type III restriction enzyme that requires two recognition sites in a defined orientation separated by up to 3.5 kbp to efficiently cleave DNA. The mechanism through which site- bound EcoP15I enzymes communicate between the two sites is unclear. Here, we use atomic force microscopy to study EcoP15I-DNA pre-cleavage complexes. From the number and size distribution of loops formed, we conclude that the loops observed do not result from translocation, but are instead formed by a contact between site- bound EcoP15I and a nonspecific region of DNA. This conclusion is confirmed by a theoretical polymer model. It is further shown that translocation must play some role, because when translocation is blocked by a Lac repressor protein, DNA cleavage is similarly blocked. On the basis of these results, we present a model for restriction by type III restriction enzymes and highlight the similarities between this and other classes of restriction enzymes.

Ultrasound-Triggered Controlled Drug Delivery and Biosensing Using Silica Nanotubes

Silica nanotubes (SNTs) have been demonstrated here as a versatile host for controlled drug delivery and biosensing. The sol-gel template synthesized SNTs have a slow rate of drug release. Application of an external stimulus in the form of ultrasound to or chemical functionalization of synthesized SNT results in higher yield of drug release as well as yield of drug release varying linearly with time. In case of controlled drug delivery triggered by ultrasound, drug yield as function of time is found to be heavily dependent on the ultrasound impulse protocol. Impulses of shorter duration (similar to 0.5 min) and shorter time intervals between successive impulses resulted in higher drug yields. Confinement of hemoglobin (Hb) inside nanometer sized channels of SNT does not have any detrimental effect on the native protein structure and function. Observance of significant enhancement in direct electron transfer of Hb makes the SNTs also promising for application in biosensors.

Three-dimensional structure of heat shock protein 90 from Plasmodium falciparum: molecular modelling approach to rational drug design against malaria

We have recently implicated heat shock protein 90 from Plasmodium falciparum (PfHsp90) as a potential drug target against malaria. Using inhibitors specific to the nucleotide binding domain of Hsp90, we have shown potent growth inhibitory effects on development of malarial parasite in human erythrocytes. To gain better understanding of the vital role played by PfHsp90 in parasite growth, we have modeled its three dimensional structure using recently described full length structure of yeast Hsp90. Sequence similarity found between PfHsp90 and yeast Hsp90 allowed us to model the core structure with high confidence. The superimposition of the predicted structure with that of the template yeast Hsp90 structure reveals an RMSD of 3.31 angstrom. The N-terminal and middle domains showed the least RMSD (1.76 angstrom) while the more divergent C-terminus showed a greater RMSD (2.84 angstrom) with respect to the template. The structure shows overall conservation of domains involved in nucleotide binding, ATPase activity, co-chaperone binding as well as inter-subunit interactions. Important co-chaperones known to modulate Hsp90 function in other eukaryotes are conserved in malarial parasite as well. An acidic stretch of amino acids found in the linker region, which is uniquely extended in PfHsp90 could not be modeled in this structure suggesting a flexible conformation. Our results provide a basis to compare the overall structure and functional pathways dependent on PfHsp90 in malarial parasite. Further analysis of differences found between human and parasite Hsp90 may make it possible to design inhibitors targeted specifically against malaria.

Triosephosphate isomerase from Plasmodium falciparum:. the crystal structure provides insights into antimalarial drug design

Background: Malaria caused by the parasite Plasmodium falciparum is a major public health concern. The parasite lacks a functional tricarboxylic acid cycle, making glycolysis its sole energy source. Although parasite enzymes have been considered as potential antimalarial drug targets, little is known about their structural biology. Here we report the crystal structure of triosephosphate isomerase (TIM) from P. falciparum at 2.2 Angstrom resolution. Results: The crystal structure of P. falciparum TIM (PfTIM), expressed in Escherichia coli, was determined by the molecular replacement method using the structure of trypanosomal TIM as the starting model. Comparison of the PfTIM structure with other TIM structures, particularly human TIM, revealed several differences, In most TIMs the residue at position 183 is a glutamate but in PtTIM it is a leucine, This leucine residue is completely exposed and together with the surrounding positively charged patch, may be responsible for binding TIM to the erythrocyte membrane. Another interesting feature is the occurrence of a cysteine residue at the dimer interface of PfTIM (Cys13), in contrast to human TIM where this residue is a methionine. Finally, residue 96 of human TIM (Ser96), which occurs near the active site, has been replaced by phenylalanine in PfTIM.