Homology modeling and metabolism prediction of human carboxylesterase-2 using docking analyses by GriDock: a parallelized tool based on AutoDock 4.0.

Metabolic problems lead to numerous failures during clinical trials, and much effort is now devoted to developing in silico models predicting metabolic stability and metabolites. Such models are well known for cytochromes P450 and some transferases, whereas less has been done to predict the activity of human hydrolases. The present study was undertaken to develop a computational approach able to predict the hydrolysis of novel esters by human carboxylesterase hCES2. The study involved first a homology modeling of the hCES2 protein based on the model of hCES1 since the two proteins share a high degree of homology (congruent with 73%). A set of 40 known substrates of hCES2 was taken from the literature; the ligands were docked in both their neutral and ionized forms using GriDock, a parallel tool based on the AutoDock4.0 engine which can perform efficient and easy virtual screening analyses of large molecular databases exploiting multi-core architectures. Useful statistical models (e.g., r (2) = 0.91 for substrates in their unprotonated state) were calculated by correlating experimental pK(m) values with distance between the carbon atom of the substrate's ester group and the hydroxy function of Ser228. Additional parameters in the equations accounted for hydrophobic and electrostatic interactions between substrates and contributing residues. The negatively charged residues in the hCES2 cavity explained the preference of the enzyme for neutral substrates and, more generally, suggested that ligands which interact too strongly by ionic bonds (e.g., ACE inhibitors) cannot be good CES2 substrates because they are trapped in the cavity in unproductive modes and behave as inhibitors. The effects of protonation on substrate recognition and the contrasting behavior of substrates and products were finally investigated by MD simulations of some CES2 complexes.

Isolation, homology modeling and renal effects of a C-type natriuretic peptide from the venom of the Brazilian yellow scorpion (Tityus serrulatus)

Mammalian natriuretic peptides (NPs) have been extensively investigated for use as therapeutic agents in the treatment of cardiovascular diseases. Here, we describe the isolation, sequencing and tridimensional homology modeling of the first C-type natriuretic peptide isolated from scorpion venom. In addition, its effects on the renal function of rats and on the mRNA expression of natriuretic peptide receptors in the kidneys are delineated. Fractionation of Tityusserrulatus venom using chromatographic techniques yielded a peptide with a molecular mass of 2190.64Da, which exhibited the pattern of disulfide bridges that is characteristic of a C-type NP (TsNP, T. serrulatus Natriuretic Peptide). In the isolated perfused rat kidney assay, treatment with two concentrations of TsNP (0.03 and 0.1μg/mL) increased the perfusion pressure, glomerular filtration rate and urinary flow. After 60min of treatment at both concentrations, the percentages of sodium, potassium and chloride transport were decreased, and the urinary cGMP concentration was elevated. Natriuretic peptide receptor-A (NPR-A) mRNA expression was down regulated in the kidneys treated with both concentrations of TsNP, whereas NPR-B, NPR-C and CG-C mRNAs were up regulated at the 0.1μg/mL concentration. In conclusion, this work describes the isolation and modeling of the first natriuretic peptide isolated from scorpion venom. In addition, examinations of the renal actions of TsNP indicate that its effects may be related to the activation of NPR-B, NPR-C and GC-C. © 2013 Elsevier Ltd.

Isolation and characterization of a cone snail protease with homology to CRISP proteins of the pathogenesis-related protein superfamily

The pathogenesis-related (PR) protein superfamily is widely distributed in the animal, plant, and fungal kingdoms and is implicated in human brain tumor growth and plant pathogenesis. The precise biological activity of PR proteins, however, has remained elusive. Here we report the characterization, cloning and structural homology modeling of Tex31 from the venom duct of Conus textile. Tex31 was isolated to >95% purity by activity-guided fractionation using a para-nitroanilide substrate based on the putative cleavage site residues found in the propeptide precursor of conotoxin TxVIA. Tex31 requires four residues including a leucine N-terminal of the cleavage site for efficient substrate processing. The sequence of Tex31 was determined using two degenerate PCR primers designed from N-terminal and tryptic digest Edman sequences. A BLAST search revealed that Tex31 was a member of the PR protein superfamily and most closely related to the CRISP family of mammalian proteins that have a cysteine-rich C-terminal tail. A homology model constructed from two PR proteins revealed that the likely catalytic residues in Tex31 fall within a structurally conserved domain found in PR proteins. Thus, it is possible that other PR proteins may also be substrate-specific proteases.

Structural characterization of respiratory syncytial virus fusion inhibitor escape mutants: homology model of the F protein and a syncytium formation assay

Respiratory syncytial virus (RSV) is a ubiquitous human pathogen and the leading cause of lower respiratory tract infections in infants. Infection of cells and subsequent formation of syncytia occur through membrane fusion mediated by the RSV fusion protein (RSV-F). A novel in vitro assay of recombinant RSV-F function has been devised and used to characterize a number of escape mutants for three known inhibitors of RSV-F that have been isolated. Homology modeling of the RSV-F structure has been carried out on the basis of a chimera derived from the crystal structures of the RSV-F core and a fragment from the orthologous fusion protein from Newcastle disease virus (NDV). The structure correlates well with the appearance of RSV-F in electron micrographs, and the residues identified as contributing to specific binding sites for several monoclonal antibodies are arranged in appropriate solvent-accessible clusters. The positions of the characterized resistance mutants in the model structure identify two promising regions for the design of fusion inhibitors. (C) 2003 Elsevier Science (USA). All rights reserved.

Functional analysis and structural modeling of human APOBEC3G reveal the role of evolutionarily conserved elements in the inhibition of human immunodeficiency virus type 1 infection and Alu transposition.

Retroelements are important evolutionary forces but can be deleterious if left uncontrolled. Members of the human APOBEC3 family of cytidine deaminases can inhibit a wide range of endogenous, as well as exogenous, retroelements. These enzymes are structurally organized in one or two domains comprising a zinc-coordinating motif. APOBEC3G contains two such domains, only the C terminal of which is endowed with editing activity, while its N-terminal counterpart binds RNA, promotes homo-oligomerization, and is necessary for packaging into human immunodeficiency virus type 1 (HIV-1) virions. Here, we performed a large-scale mutagenesis-based analysis of the APOBEC3G N terminus, testing mutants for (i) inhibition of vif-defective HIV-1 infection and Alu retrotransposition, (ii) RNA binding, and (iii) oligomerization. Furthermore, in the absence of structural information on this domain, we used homology modeling to examine the positions of functionally important residues and of residues found to be under positive selection by phylogenetic analyses of primate APOBEC3G genes. Our results reveal the importance of a predicted RNA binding dimerization interface both for packaging into HIV-1 virions and inhibition of both HIV-1 infection and Alu transposition. We further found that the HIV-1-blocking activity of APOBEC3G N-terminal mutants defective for packaging can be almost entirely rescued if their virion incorporation is forced by fusion with Vpr, indicating that the corresponding region of APOBEC3G plays little role in other aspects of its action against this pathogen. Interestingly, residues forming the APOBEC3G dimer interface are highly conserved, contrasting with the rapid evolution of two neighboring surface-exposed amino acid patches, one targeted by the Vif protein of primate lentiviruses and the other of yet-undefined function.

Molecular Modeling Approaches for Determining Gene Function: application to a Putative Poly-A Binding Protein from Leishmania amazonensis (LaPABP)

The great expansion in the number of genome sequencing projects has revealed the importance of computational methods to speed up the characterization of unknown genes. These studies have been improved by the use of three dimensional information from the predicted proteins generated by molecular modeling techniques. In this work, we disclose the structure-function relationship of a gene product from Leishmania amazonensis by applying molecular modeling and bioinformatics techniques. The analyzed sequence encodes a 159 aminoacids polypeptide (estimated 18 kDa) and was denoted LaPABP for its high homology with poly-A binding proteins from trypanosomatids. The domain structure, clustering analysis and a three dimensional model of LaPABP, basically obtained by homology modeling on the structure of the human poly-A binding protein, are described. Based on the analysis of the electrostatic potential mapped on the model's surface and conservation of intramolecular contacts responsible for folding stabilization we hypothesize that this protein may have less avidity to RNA than it's L. major counterpart but still account for a significant functional activity in the parasite. The model obtained will help in the design of mutagenesis experiments aimed to elucidate the mechanism of gene expression in trypanosomatids and serve as a starting point for its exploration as a potential source of targets for a rational chemotherapy.

Molecule modeling of human pentameric alpha(7) neuronal nicotinic acetylcholine receptor and its interaction with its agonist and competitive antagonist

The nicotinic Acetylcholine Receptor (nAChR) is the major class of neurotransmitter receptors that is involved in many neurodegenerative conditions such as schizophrenia, Alzheimer's and Parkinson's diseases. The N-terminal region or Ligand Binding Domain (LBD) of nAChR is located at pre- and post-synaptic nervous system, which mediates synaptic transmission. nAChR acts as the drug target for agonist and competitive antagonist molecules that modulate signal transmission at the nerve terminals. Based on Acetylcholine Binding Protein (AChBP) from Lymnea stagnalis as the structural template, the homology modeling approach was carried out to build three dimensional model of the N-terminal region of human alpha(7)nAChR. This theoretical model is an assembly of five alpha(7) subunits with 5 fold axis symmetry, constituting a channel, with the binding picket present at the interface region of the subunits. alpha-netlrotoxin is a potent nAChR competitive antagonist that readily blocks the channel resulting in paralysis. The molecular interaction of alpha-Bungarotoxin, a long chain alpha-neurotoxin from (Bungarus multicinctus) and human alpha(7)nAChR seas studied. Agonists such as acetylcholine, nicotine, which are used in it diverse array of biological activities, such as enhancements of cognitive performances, were also docked with the theoretical model of human alpha(7)nAChR. These docked complexes were analyzed further for identifying the crucial residues involved in interaction. These results provide the details of interaction of agonists and competitive antagonists with three dimensional model of the N-terminal region of human alpha(7)nAChR and thereby point to the design of novel lead compounds.

Advanced Protein Modeling Method: Benchmarking and Applications in Computer-Aided Drug Discovery

Development of homology modeling methods will remain an area of active research. These methods aim to develop and model increasingly accurate three-dimensional structures of yet uncrystallized therapeutically relevant proteins e.g. Class A G-Protein Coupled Receptors. Incorporating protein flexibility is one way to achieve this goal. Here, I will discuss the enhancement and validation of the ligand-steered modeling, originally developed by Dr. Claudio Cavasotto, via cross modeling of the newly crystallized GPCR structures. This method uses known ligands and known experimental information to optimize relevant protein binding sites by incorporating protein flexibility. The ligand-steered models were able to model, reasonably reproduce binding sites and the co-crystallized native ligand poses of the β2 adrenergic and Adenosine 2A receptors using a single template structure. They also performed better than the choice of template, and crude models in a small scale high-throughput docking experiments and compound selectivity studies. Next, the application of this method to develop high-quality homology models of Cannabinoid Receptor 2, an emerging non-psychotic pain management target, is discussed. These models were validated by their ability to rationalize structure activity relationship data of two, inverse agonist and agonist, series of compounds. The method was also applied to improve the virtual screening performance of the β2 adrenergic crystal structure by optimizing the binding site using β2 specific compounds. These results show the feasibility of optimizing only the pharmacologically relevant protein binding sites and applicability to structure-based drug design projects.

Enzymatic characterization and homology model of a catalytically active recombinant West Nile virus NS3 protease

West Nile Virus (WNV) is a mosquito-borne flavivirus with a rapidly expanding global distribution. Infection causes severe neurological disease and fatalities in both human and animal hosts. The West Nile viral protease (NS2B-NS3) is essential for post-translational processing in host-infected cells of a viral polypeptide precursor into structural and functional viral proteins, and its inhibition could represent a potential treatment for viral infections. This article describes the design, expression, and enzymatic characterization of a catalytically active recombinant WNV protease, consisting of a 40-residue component of cofactor NS2B tethered via a noncleavable nonapeptide (G(4)SG(4)) to the N-terminal 184 residues of NS3. A chromogenic assay using synthetic para-nitroanilide (pNA) hexapeptide substrates was used to identify optimal enzyme-processing conditions (pH 9.5, I < 0.1 M, 30% glycerol, 1 mM CHAPS), preferred substrate cleavage sites, and the first competitive inhibitor (Ac-FASGKR- H, IC50 &SIM; 1 μM). A putative three-dimensional structure of WNV protease, created through homology modeling based on the crystal structures of Dengue-2 and Hepatitis C NS3 viral proteases, provides some valuable insights for structure-based design of potent and selective inhibitors of WNV protease.

Structural aspects of the p.P222Q homozygous mutation of HSD3B2 gene in a patient with congenital adrenal hyperplasia

Type II 3&#946;-hydroxysteroid dehydrogenase/&#916;5-&#916;4-isomerase (3&#946;-HSD2), encoded by the HSD3B2 gene, is a key enzyme involved in the biosynthesis of all the classes of steroid hormones. Deleterious mutations in the HSD3B2 gene cause the classical deficiency of 3&#946;-HSD2, which is a rare autosomal recessive disease that leads to congenital adrenal hyperplasia (CAH). CAH is the most frequent cause of ambiguous genitalia and adrenal insufficiency in newborn infants with variable degrees of salt losing. Here we report the molecular and structural analysis of the HSD3B2 gene in a 46,XY child, who was born from consanguineous parents, and presented with ambiguous genitalia and salt losing. The patient carries a homozygous nucleotide c.665C>A change in exon 4 that putatively substitutes the proline at codon 222 for glutamine. Molecular homology modeling of normal and mutant 3&#946;-HSD2 enzymes emphasizes codon 222 as an important residue for the folding pattern of the enzyme and validates a suitable model for analysis of new mutations.

Molecular determinants of improved cathepsin B inhibition by new cystatins obtained by DNA shuffling

Background: Cystatins are inhibitors of cysteine proteases. The majority are only weak inhibitors of human cathepsin B, which has been associated with cancer, Alzheimer's disease and arthritis. Results: Starting from the sequences of oryzacystatin-1 and canecystatin-1, a shuffling library was designed and a hybrid clone obtained, which presented higher inhibitory activity towards cathepsin B. This clone presented two unanticipated point mutations as well as an N-terminal deletion. Reversing each point mutation independently or both simultaneously abolishes the inhibitory activity towards cathepsin B. Homology modeling together with experimental studies of the reverse mutants revealed the likely molecular determinants of the improved inhibitory activity to be related to decreased protein stability. Conclusion: A combination of experimental approaches including gene shuffling, enzyme assays and reverse mutation allied to molecular modeling has shed light upon the unexpected inhibitory properties of certain cystatin mutants against Cathepsin B. We conclude that mutations disrupting the hydrophobic core of phytocystatins increase the flexibility of the N-terminus, leading to an increase in inhibitory activity. Such mutations need not affect the inhibitory site directly but may be observed distant from it and manifest their effects via an uncoupling of its three components as a result of increased protein flexibility.

Computer-aided Drug Design of Novel PLA(2) Inhibitor Candidates for Treatment of Snakebite

Phospholipases A(2) (PLA(2)) are enzymes commonly found in snake venoms from Viperidae and Elaphidae families, which are major components thereof. Many plants are used in traditional medicine its active agents against various effects induced by snakebite. This article presents the PLA(2) BthTX-I structure prediction based on homology modeling. In addition, we have performed virtual screening in a large database yielding a set of potential bioactive inhibitors. A flexible docking program was used to investigate the interactions between the receptor and the new ligands. We have performed molecular interaction fields (MIFs) calculations with the phospholipase model. Results confirm the important role of Lys49 for binding ligands and suggest three additional residues as well. We have proposed a theoretically nontoxic, drug-like, and potential novel BthTX-I inhibitor. These calculations have been used to guide the design of novel phospholipase inhibitors as potential lead compounds that may be optimized for future treatment of snakebite victims as well as other human diseases in which PLA(2) enzymes are involved.

The AF-1 Domain of the Orphan Nuclear Receptor NOR-1 Mediates Trans-activation, Coactivator Recruitment, and Activation by the Purine Anti-metabolite 6-Mercaptopurine

NOR-1/NR4A3 is an orphan member of the nuclear hormone receptor superfamily. NOR-1 and its close relatives Nurr1 and Nur77 are members of the NR4A subgroup of nuclear receptors. Members of the NR4A subgroup are induced through multiple signal transduction pathways. They have been implicated in cell proliferation, differentiation, T-cell apoptosis, chondrosarcomas, neurological disorders, inflammation, and atherogenesis. However, the mechanism of transcriptional activation, coactivator recruitment, and agonist-mediated activation remain obscure. Hence, we examined the molecular basis of NOR-1-mediated activation. We observed that NOR-1 trans-activates gene expression in a cell- and target-specific manner; moreover, it operates in an activation function (AF)-1-dependent manner. The N-terminal AF-1 domain delimited to between amino acids 1 and 112, preferentially recruits the steroid receptor coactivator (SRC). Furthermore, SRC-2 modulates the activity of the AF-1 domain but not the C-terminal ligand binding domain (LBD). Homology modeling indicated that the NOR-1 LBD was substantially different from that of hRORbeta, a closely related AF-2-dependent receptor. In particular, the hydrophobic cleft characteristic of nuclear receptors was replaced with a very hydrophilic surface with a distinct topology. This observation may account for the inability of this nuclear receptor LBD to efficiently mediate cofactor recruitment and transcriptional activation. In contrast, the N-terminal AF-1 is necessary for cofactor recruitment and can independently conscript coactivators. Finally, we demonstrate that the purine anti-metabolite 6-mercaptopurine, a widely used antineoplastic and anti-inflammatory drug, activates NOR-1 in an AF-1-dependent manner. Additional 6-mercaptopurine analogs all efficiently activated NOR-1, suggesting that the signaling pathways that modulate proliferation via inhibition of de novo purine and/or nucleic acid biosynthesis are involved in the regulation NR4A activity. We hypothesize that the NR4A subgroup mediates the genotoxic stress response and suggest that this subgroup may function as sensors that respond to genotoxicity.

The fourth extracellular loop of the alpha subunit of Na,K-ATPase. Functional evidence for close proximity with the second extracellular loop.

Na,K-ATPase is the main active transport system that maintains the large gradients of Na(+) and K(+) across the plasma membrane of animal cells. The crystal structure of a K(+)-occluding conformation of this protein has been recently published, but the movements of its different domains allowing for the cation pumping mechanism are not yet known. The structure of many more conformations is known for the related calcium ATPase SERCA, but the reliability of homology modeling is poor for several domains with low sequence identity, in particular the extracellular loops. To better define the structure of the large fourth extracellular loop between the seventh and eighth transmembrane segments of the alpha subunit, we have studied the formation of a disulfide bond between pairs of cysteine residues introduced by site-directed mutagenesis in the second and the fourth extracellular loop. We found a specific pair of cysteine positions (Y308C and D884C) for which extracellular treatment with an oxidizing agent inhibited the Na,K pump function, which could be rapidly restored by a reducing agent. The formation of the disulfide bond occurred preferentially under the E2-P conformation of Na,K-ATPase, in the absence of extracellular cations. Using recently published crystal structure and a distance constraint reproducing the existence of disulfide bond, we performed an extensive conformational space search using simulated annealing and showed that the Tyr(308) and Asp(884) residues can be in close proximity, and simultaneously, the SYGQ motif of the fourth extracellular loop, known to interact with the extracellular domain of the beta subunit, can be exposed to the exterior of the protein and can easily interact with the beta subunit.