On the Validation of Protein Force Fields Based on Structural Criteria

Molecular dynamics simulations depend critically on the quality of the force field used to describe the interatomic interactions and the extent to which it has been validated for use in a specific application. Using a curated test set of 52 high-resolution structures, 39 derived from X-ray diffraction and 13 solved using NMR, we consider the extent to which different parameter sets of the GROMOS protein force field can be distinguished based on comparing a range of structural criteria, including the number of backbone hydrogen bonds, the number of native hydrogen bonds, polar and nonpolar solvent-accessible surface area, radius of gyration, the prevalence of secondary structure elements, J-coupling constants, nuclear Overhauser effect (NOE) intensities, positional root-mean-square deviations (RMSD), and the distribution of backbone ϕ and ψ dihedral angles. It is shown that while statistically significant differences between the average values of individual metrics could be detected, these were in general small. Furthermore, improvements in agreement in one metric were often offset by loss of agreement in another. The work establishes a framework and test set against which protein force fields can be validated. It also highlights the danger of inferring the relative quality of a given force field based on a small range of structural properties or small number of proteins.

Permeability of Dermatological Solutes through the Short Periodicity Phase of Human Stratum Corneum Lipid Bilayers

Determining the permeability of drug-like solutes through the densely packed and heterogeneous stratum corneum lipid layer presents a significant challenge. In this study, we employed umbrella sampling with a periodic weighing function applied to the center of mass of the lipid bilayers. Precise umbrella sampling was conducted with an interframe distance of 0.5 Å, spanning from the bilayer center to the water phase, and each frame was simulated for at least 20 ns. Autocorrelation functions, potential of mean force (PMF), and diffusivity profiles were analyzed for six solutes (testosterone, benzene, caffeine, ethanol, mannitol, and histidine). The results revealed that autocorrelations were dependent solely on the medium, whether water or lipid phase. Diffusivity and PMF profiles along the reaction coordinate were influenced by the hydrophilicity of the solute rather than its size. For hydrophobic solutes, the PMF curves exhibited a minimum at the bilayer center, while for hydrophilic solutes, the PMFs peaked at the bilayer center and lipid tails (where the lipid tails are not interacting with the cholesterol). Diffusivity curves were low at the bilayer center and water phase, with peaks observed at the headgroup or the boundary between fatty acid and cholesterol (1 nm from the bilayer center). The quantitative findings presented in this work hold significance for pharmacists and drug designers.

Development of a model for granule-bound starch synthase I activity using free-energy calculations

Starch is a branched polymer of glucose with two components, both of which have (1 → 4)-α linear links and (1 → 6)-α branch points: amylopectin, of high molecular weight with many short branches, and amylose, of lower molecular weight and only a few long-chain branches. Granule-bound starch synthase I (GBSSI) is one of the main enzymes controlling amylose synthesis and chain-length distribution. As production of different GBSSI mutants is time-consuming and laborious, molecular dynamics (MD) simulations are used here to predict the binding of different GBSSI mutants to a representative amylose fragment. The simulations were atomistic, with explicit solvent and docking, a method successfully used to understand the binding of wild-type GBSSI to amylose fragments. The binding of GBSSI to G5 (a pentasaccharide amylose fragment) is combined with free-energy calculations employing a thermodynamic integration method to predict the effects of mutations on enzyme activity. Ten GBSSI mutants with different enzyme activities were analyzed to find the structural and energy changes among different single amino-acid mutants and their possible relationship to starch characteristics. Comparing the structural changes and the relative binding free energy of G5 to the wild type GBSSI and GBSSI mutants, it was found that mutants with negative binding energy (lower than -2.0 kcal/mol) are more likely to have higher enzyme activity and amylose content compared to the wild type. This theoretical paper used simulations and robust free energy calculations to interpret in planta data with potential predictions as to what mutants might be generated to give desired properties. This study can be used to help develop grains with improved functional properties.

Probing Disaccharide Binding to Triplatin as Models for Tumor Cell Heparan Sulfate (GAG) Interactions

In this study, we have used [1H, 15N] NMR spectroscopy to investigate the interactions of the trinuclear platinum anticancer drug triplatin (1) (1,0,1/t,t,t or BBR3464) with site-specific sulfated and carboxylated disaccharides. Specifically, the disaccharides GlcNS(6S)-GlcA (I) and GlcNS(6S)-IdoA(2S) (II) are useful models of longer-chain glycosaminoglycans (GAGs) such as heparan sulfate (HS). For both the reactions of 15N-1 with I and II, equilibrium conditions were achieved more slowly (65 h) compared to the reaction with the monosaccharide GlcNS(6S) (9 h). The data suggest both carboxylate and sulfate binding of disaccharide I to the Pt with the sulfato species accounting for <1% of the total species at equilibrium. The rate constant for sulfate displacement of the aqua ligand (kL2) is 4 times higher than the analogous rate constant for carboxylate displacement (kL1). There are marked differences in the equilibrium concentrations of the chlorido, aqua, and carboxy-bound species for reactions with the two disaccharides, notably a significantly higher concentration of carboxylate-bound species for II, where sulfate-bound species were barely detectable. The trend mirrors that reported for the corresponding dinuclear platinum complex 1,1/t,t, where the rate constant for sulfate displacement of the aqua ligand was 3 times higher than that for acetate. Also similar to what we observed for the reactions of 1,1/t,t with the simple anions, aquation of the sulfato group is rapid, and the rate constant k-L2 is 3 orders of magnitude higher than that for displacement of the carboxylate (k-L1). Molecular dynamics calculations suggest that extra hydrogen-bonding interactions with the more sulfated disaccharide II may prevent or diminish sulfate binding of the triplatin moiety. The overall results suggest that Pt-O donor interactions should be considered in any full description of platinum complex cellular chemistry.

Understanding the Binding of Starch Fragments to Granule-Bound Starch Synthase

Granule-bound starch synthase (GBSS) plays a major role, that of chain elongation, in the biosynthesis of amylose, a starch component with mostly (1 → 4)-α connected long chains of glucose with a few (1 → 6)-α branch points. Chain-length distributions (CLDs) of amylose affect functional properties, which can be controlled by changing appropriate residues on granule-bound starch synthase (GBSS). Knowing the binding of GBSS and amylose at a molecular level can help better determine the key amino acids on GBSS that affect CLDs of amylose for subsequent use in molecular engineering. Atomistic molecular dynamics simulations with explicit solvent and docking approaches were used in this study to build a model of the binding between rice GBSS and amylose. Amylose fragments containing 3-12 linearly linked glucose units were built to represent the starch fragments. The stability of the complexes, interactions between GBSS and sugars, and difference in structure/conformation of bound and free starch fragments were analyzed. The study found that starch/amylose fragments with 5 or 6 glucose units were suitable for modeling starch binding to GBSS. The removal of an interdomain disulfide on GBSS was found to affect both GBSS and starch stability. Key residues that could affect the binding ability were also indicated. This model can help rationalize the design of mutants and suggest ways to make single-point mutations, which could be used to develop plants producing starches with improved functional properties.

Platinum complexes act as shielding agents against virus infection

We determine that the substitution-inert polynuclear platinum complex (PPC) TriplatinNC is an antiviral agent and protects cells from enterovirus 71 and human metapneumovirus infection. This protection occurs through the formation of adducts with cell-surface glycosaminoglycans. Our detailed mechanistic investigation demonstrates that TriplatinNC blocks viral entry by shielding cells from virus attack, opening new directions for metalloshielding antiviral drug development.

MyD88 TIR domain higher-order assembly interactions revealed by microcrystal electron diffraction and serial femtosecond crystallography

MyD88 and MAL are Toll-like receptor (TLR) adaptors that signal to induce pro-inflammatory cytokine production. We previously observed that the TIR domain of MAL (MALTIR) forms filaments in vitro and induces formation of crystalline higher-order assemblies of the MyD88 TIR domain (MyD88TIR). These crystals are too small for conventional X-ray crystallography, but are ideally suited to structure determination by microcrystal electron diffraction (MicroED) and serial femtosecond crystallography (SFX). Here, we present MicroED and SFX structures of the MyD88TIR assembly, which reveal a two-stranded higher-order assembly arrangement of TIR domains analogous to that seen previously for MALTIR. We demonstrate via mutagenesis that the MyD88TIR assembly interfaces are critical for TLR4 signaling in vivo, and we show that MAL promotes unidirectional assembly of MyD88TIR. Collectively, our studies provide structural and mechanistic insight into TLR signal transduction and allow a direct comparison of the MicroED and SFX techniques.

SARM1 is a metabolic sensor activated by an increased NMN/NAD+ ratio to trigger axon degeneration

Axon degeneration is a central pathological feature of many neurodegenerative diseases. Sterile alpha and Toll/interleukin-1 receptor motif-containing 1 (SARM1) is a nicotinamide adenine dinucleotide (NAD+)-cleaving enzyme whose activation triggers axon destruction. Loss of the biosynthetic enzyme NMNAT2, which converts nicotinamide mononucleotide (NMN) to NAD+, activates SARM1 via an unknown mechanism. Using structural, biochemical, biophysical, and cellular assays, we demonstrate that SARM1 is activated by an increase in the ratio of NMN to NAD+ and show that both metabolites compete for binding to the auto-inhibitory N-terminal armadillo repeat (ARM) domain of SARM1. We report structures of the SARM1 ARM domain bound to NMN and of the homo-octameric SARM1 complex in the absence of ligands. We show that NMN influences the structure of SARM1 and demonstrate via mutagenesis that NMN binding is required for injury-induced SARM1 activation and axon destruction. Hence, SARM1 is a metabolic sensor responding to an increased NMN/NAD+ ratio by cleaving residual NAD+, thereby inducing feedforward metabolic catastrophe and axonal demise.

Bisubstrate Ether-Linked Uridine-Peptide Conjugates as O-GlcNAc Transferase Inhibitors

The O-linked β-N-acetylglucosamine (O-GlcNAc) transferase (OGT) is a master regulator of installing O-GlcNAc onto serine or threonine residues on a multitude of target proteins. Numerous nuclear and cytosolic proteins of varying functional classes, including translational factors, transcription factors, signaling proteins, and kinases are OGT substrates. Aberrant O-GlcNAcylation of proteins is implicated in signaling in metabolic diseases such as diabetes and cancer. Selective and potent OGT inhibitors are valuable tools to study the role of OGT in modulating a wide range of effects on cellular functions. We report linear bisubstrate ether-linked uridine-peptide conjugates as OGT inhibitors with micromolar affinity. In vitro evaluation of the compounds revealed the importance of donor substrate, linker and acceptor substrate in the rational design of bisubstrate analogue inhibitors. Molecular dynamics simulations shed light on the binding of this novel class of inhibitors and rationalized the effect of amino acid truncation of acceptor peptide on OGT inhibition.

Simplified Novel Muraymycin Analogues; using a Serine Template Strategy for Linking Key Pharmacophores

The present status of antibiotic research requires the urgent invention of novel agents that act on multidrug-resistant bacteria. The World Health Organization has classified antibiotic-resistant bacteria into critical, high and medium priority according to the urgency of need for new antibiotics. Naturally occurring uridine-derived "nucleoside antibiotics" have shown promising activity against numerous priority resistant organisms by inhibiting the transmembrane protein MraY (translocase I), which is yet to be explored in a clinical context. The catalytic activity of MraY is an essential process for bacterial cell viability and growth including that of priority organisms. Muraymycins are one subclass of naturally occurring MraY inhibitors. Despite having potent antibiotic properties, the structural complexity of muraymycins advocates for simplified analogues as potential lead structures. Herein, we report a systematic structure-activity relationship (SAR) study of serine template-linked, simplified muraymycin-type analogues. This preliminary SAR lead study of serine template analogues successfully revealed that the complex structure of naturally occurring muraymycins could be easily simplified to afford bioactive scaffolds against resistant priority organisms. This study will pave the way for the development of novel antibacterial lead compounds based on a simplified serine template.

Automated Topology Builder Version 3.0: Prediction of Solvation Free Enthalpies in Water and Hexane

The ability of atomic interaction parameters generated using the Automated Topology Builder and Repository version 3.0 (ATB3.0) to predict experimental hydration free enthalpies (Δ G^water) and solvation free enthalpies in the apolar solvent hexane (Δ G^hexane) is presented. For a validation set of 685 molecules the average unsigned error (AUE) between Δ G^water values calculated using the ATB3.0 and experiment is 3.8 kJ·mol^-1. The slope of the line of best fit is 1.00, the intercept -1.0 kJ·mol^-1, and the R² 0.90. For the more restricted set of 239 molecules used to validate OPLS3 ( J. Chem. Theory Comput. 2016 , 12 , 281 - 296 , DOI: 10.1021/acs.jctc.5b00864 ) the AUE using the ATB3.0 is just 2.7 kJ·mol^-1 and the R² 0.93. A roadmap for further improvement of the ATB parameters is presented together with a discussion of the challenges of validating force fields against the available experimental data.

Predicting the Prevalence of Alternative Warfarin Tautomers in Solution

Warfarin, a widely used oral anticoagulant, is prescribed as a racemic mixture. Each enantiomer of neutral Warfarin can exist in 20 possible tautomeric states leading to complex pharmacokinetics and uncertainty as to the relevant species under different conditions. Here, the ability of alternative computational approaches to predict the preferred tautomeric form(s) of neutral Warfarin in different solvents is examined. It is shown that varying the method used to estimate the heat of formation in vacuum (direct or via homodesmic reactions), whether entropic corrections were included, and the method used to estimate the free enthalpy of solvation (i.e., PCM, COSMO, or SMD implicit models or explicit solvent) lead to large differences in the predicted rank and relative populations of the tautomers. In this case, only a combination of the enthalpy of formation using homodesmic reactions and explicit solvent to estimate the free enthalpy of solvation yielded results compatible with the available experimental data. The work also suggests that a small but significant subset of the possible Warfarin tautomers are likely to be physiologically relevant.

Molecular dynamics simulations reveal structural insights into inhibitor binding modes and functionality in human Group IIA phospholipase A2

Human Group IIA phospholipase A₂ (hGIIA) promotes inflammation in immune-mediated pathologies by regulating the arachidonic acid pathway through both catalysis-dependent and -independent mechanisms. The hGIIA crystal structure, both alone and inhibitor-bound, together with structures of closely related snake-venom-derived secreted phospholipase enzymes has been well described. However, differentiation of biological and nonbiological contacts and the relevance of structures determined from snake venom enzymes to human enzymes are not clear. We employed molecular dynamics (MD) and docking approaches to understand the binding of inhibitors that selectively or nonselectively block the catalysis-independent mechanism of hGIIA. Our results indicate that hGIIA behaves as a monomer in the solution environment rather than a dimer arrangement that is in the asymmetric unit of some crystal structures. The binding mode of a nonselective inhibitor, KH064, was validated by a combination of the experimental electron density and MD simulations. The binding mode of the selective pentapeptide inhibitor FLSYK to hGIIA was stipulated to be different to that of the snake venom phospholipases A₂ of Daboia russelli pulchella (svPLA₂ ). Our data suggest that the application of MD approaches to crystal structure data is beneficial in evaluating the robustness of conclusions drawn based on crystal structure data alone. Proteins 2017; 85:827-842. © 2016 Wiley Periodicals, Inc.

The B(OH)-NH Analog Is a Surrogate for the Amide Bond (CO-NH) in Peptides: An ab Initio Study

The conformational preferences of N-methyl-methylboronamide (NMB), a B(OH)-NH analog of the amide CO-NH in natural peptides, have been investigated at the Hartree-Fock; Becke's three-parameter exchange functional and the gradient-corrected functional of Lee, Yang, and Parr; and second-order Møller-Plesset levels of theory with the 6-31+G* basis set. The minima, saddle points, and rotation barriers on the potential energy surface of NMB have been located and the energy barriers estimated. Besides the global minimum, there are three local minima within 2.0 kcal mol(-)(1) of the global minimum characterized by specific ω and τ torsion values. The energy barriers for rotation about the "ω angle" are 16.4-18.8 kcal mol(-)(1) and are a consequence of the double-bond character of the B-N bond as revealed by natural bond orbitals calculations. The "ω angle" and the ω rotation barrier are nearly the same as those seen in natural peptides. The τ rotation barriers (B-O bond) are relatively low because of the single-bond character of the B-O bond. Ala-BON, the Ala-dipeptide derived from NMB, has been constructed as a model peptide to study the conformational preferences about the φ and ψ torsion angles. The study reveals a strong preference for α-helix, type-II β-turn, 2.27 ribbon, and antiparallel β-sheet conformations, and mirror images of both type-II β-turn and 2.27 ribbon motifs whose φ and ψ values fall in the "disfavored regions" of the Ramachandran map. Thus, the replacement of the carbonyl group by B-OH retains the geometry and barrier around the "ω angle" and induces a strong preference for regular secondary structure motifs and also structures with positive φ values. This makes the B(OH)-NH analog an important surrogate for the peptide bond, with the additional advantage of stability to proteolytic enzymes.

The ω, φ, and ψ Space of N-Hydroxy-N-methylacetamide and N-Acetyl-N '-hydroxy-N '-methylamide of Alanine and Their Boron Isosteres

The conformational space of N-hydroxy-N-methylacetamide [CH3-CO-N(OH)CH3, NMAOH] and its boron isostere [CH3-CO-B(OH)CH3, BMAOH] has been studied by quantum chemical methods. The potential energy surface of NMAOH and BMAOH has been built at the HF, B3LYP, and MP2 levels of theory with the 6-31+G* basis set. The minima and transition states for rotations about various torsional angles have been located, and the energy barriers have been estimated. The global minimum energy structure of both peptides exhibits an intramolecular hydrogen bond between the carbonyl oxygen and the hydroxyl group, imparting a conformational rigidity to the peptides. The omega rotation barrier is lower in the boron isostere than in NMAOH. The difference in the rotation barrier has been attributed to second-order orbital interactions, like negative hyperconjugation, as revealed by NBO calculations. In contrast, the rotation barrier around the torsion angle tau (torsion governing rotation about the N-O and B-O bonds) is relatively higher in the boron analogue. This difference is due to the double bond character in the B-O bond as opposed to the N-O bond which has the character of a single bond. As an extension, N-acetyl-N'-hydroxy-N'-methylamide of alanine (Ala-NOH) and its boron isostere (Ala-BOH) have been adopted as model peptides to study the conformational preferences about the φ and ψ torsion angles. The study reveals a strong preference for a Type I beta turn as well as inclinations for a left-handed alpha helix, for positive phi torsions, and for extended psi conformations for Ala-NOH; Ala-BOH, on the other hand, shows a leaning toward positive phi and extended psi, with no preference for any regular secondary structure motifs. The replacement of nitrogen by boron changes the electronic and conformational properties of the peptide, extending greater flexibility around the omega angle, a strong preference for positive phi values, and a shift in the site of nucleophilic attack from the carbonyl group to boron.

The characterization of modified starch branching enzymes: toward the control of starch chain-length distributions

Starch is a complex branched glucose polymer whose branch molecular weight distribution (the chain-length distribution, CLD) influences nutritionally important properties such as digestion rate. Chain-stopping in starch biosynthesis is by starch branching enzyme (SBE). Site-directed mutagenesis was used to modify SBEIIa from Zea mays (mSBEIIa) to produce mutants, each differing in a single conserved amino-acid residue. Products at different times from in vitro branching were debranched and the time evolution of the CLD measured by size-exclusion chromatography. The results confirm that Tyr352, Glu513, and Ser349 are important for mSBEIIa activity while Arg456 is important for determining the position at which the linear glucan is cut. The mutant mSBEIIa enzymes have different activities and suggest the length of the transferred chain can be varied by mutation. The work shows analysis of the molecular weight distribution can yield information regarding the enzyme branching sites useful for development of plants yielding starch with improved functionality.

Missing fragments: detecting cooperative binding in fragment-based drug design

The aim of fragment-based drug design (FBDD) is to identify molecular fragments that bind to alternate subsites within a given binding pocket leading to cooperative binding when linked. In this study, the binding of fragments to human phenylethanolamine N-methyltransferase is used to illustrate how (a) current protocols may fail to detect fragments that bind cooperatively, (b) theoretical approaches can be used to validate potential hits, and (c) apparent false positives obtained when screening against cocktails of fragments may in fact indicate promising leads.

Challenges in the determination of the binding modes of non-standard ligands in X-ray crystal complexes

Despite its central role in structure based drug design the determination of the binding mode (position, orientation and conformation in addition to protonation and tautomeric states) of small heteromolecular ligands in protein:ligand complexes based on medium resolution X-ray diffraction data is highly challenging. In this perspective we demonstrate how a combination of molecular dynamics simulations and free energy (FE) calculations can be used to correct and identify thermodynamically stable binding modes of ligands in X-ray crystal complexes. The consequences of inappropriate ligand structure, force field and the absence of electrostatics during X-ray refinement are highlighted. The implications of such uncertainties and errors for the validation of virtual screening and fragment-based drug design based on high throughput X-ray crystallography are discussed with possible solutions and guidelines.

CoMFA study of distamycin analogs binding to the minor-groove of DNA: a unified model for broad-spectrum activity

A 3D-QSAR analysis has been carried out by comparative molecular field analysis (CoMFA) on a series of distamycin analogs that bind to the DNA of drug-resistant bacterial strains MRSA, PRSP and VSEF. The structures of the molecules were derived from the X-ray structure of distamycin bound to DNA and were aligned using the Database alignment method in Sybyl. Statistically significant CoMFA models for each activity were generated. The CoMFA contours throw light on the structure activity relationship (SAR) and help to identify novel features that can be incorporated into the distamycin framework to improve the activity. Common contours have been gleaned from the three models to construct a unified model that explains the steric and electrostatic requirements for antimicrobial activity against the three resistant strains.

Design of inhibitors of the MurF enzyme of Streptococcus pneumoniae using docking, 3D-QSAR, and de Novo design

The biosynthetic pathway for formation of the bacterial cell wall (peptidoglycan) presents an attractive target for intervention. This is exploited by many of the clinically useful antibiotics, which inhibit enzymes involved in the later stages of peptidoglycan synthesis. MurF is one of the four amide bond-forming enzymes (d-alanyl-d-alanine ligating enzyme) that catalyzes the ATP-dependent formation of UDP-MurNAc-tripeptide. In the present study, several MurF inhibitors were docked into the active site of MurF to explore their binding modes and also to gain an insight into the crucial ligand-receptor interactions at the molecular level. The final selection of the "bioactive" conformation of every ligand was influenced by consensus scoring in which various independent scoring functions such as GoldScore, ChemScore, HINT score and X-CScore were employed. Subsequently, 3D-QSAR studies using comparative molecular field analysis (CoMFA) and the new approach comparative residue interaction analysis (CoRIA) have been carried out on the enzyme-inhibitor complexes obtained by docking and postscoring analysis. Finally, new inhibitors have been designed using the de novo approach of Ludi, and the activities of the most promising hits have been predicted with the CoMFA and CoRIA models.

Understanding interactions of gastric inhibitory polypeptide (GIP) with its G-protein coupled receptor through NMR and molecular modeling

Gastric inhibitory polypeptide (GIP, or glucose-dependent insulinotropic polypeptide) is a 42-amino acid incretin hormone moderating glucose-induced insulin secretion. Antidiabetic therapy based on GIP holds great promise because of the fact that its insulinotropic action is highly dependent on the level of glucose, overcoming the sideeffects of hypoglycemia associated with the current therapy of Type 2 diabetes. The truncated peptide, GIP(1-30)NH2, has the same activity as the full length native peptide. We have studied the structure of GIP(1-30)NH2 and built a model of its G-protein coupled receptor (GPCR). The structure of GIP(1-30)NH2 in DMSO-d6 and H2O has been studied using 2D NMR (total correlation spectroscopy (TOCSY), nuclear overhauser effect spectroscopy (NOESY), double quantum filtered-COSY (DQF-COSY), 13C-heteronuclear single quantum correlation (HSQC) experiments, and its conformation built by MD simulations with the NMR data as constraints. The peptide in DMSO-d6 exhibits an alpha-helix between residues Ile12 and Lys30 with a discontinuity at residues Gln19 and Gln20. In H2O, the alpha-helix starts at Ile7, breaks off at Gln19, and then continues right through to Lys30. GIP(1-30)NH2 has all the structural features of peptides belonging to family B1 GPCRs, which are characterized by a coil at the N-terminal and a long C-terminal alpha-helix with or without a break. A model of the seven transmembrane (TM) helices of the GIP receptor (GIPR) has been built on the principles of comparative protein modeling, using the crystal structure of bovine rhodopsin as a template. The N-terminal domain of GIPR has been constructed from the NMR structure of the N-terminal of corticoptropin releasing factor receptor (CRFR), a family B1 GCPR. The intra and extra cellular loops and the C-terminal have been modeled from fragments retrieved from the PDB. On the basis of the experimental data available for some members of family B1 GPCRs, four pairs of constraints between GIP(1-30)NH2 and its receptor were used in the FTDOCK program, to build the complete model of the GIP(1-30)NH2:GIPR complex. The model can rationalize the various experimental observations including the potency of the truncated GIP peptide. This work is the first complete model at the atomic level of GIP(1-30)NH2 and of the complex with its GPCR.

Pharmacophore modeling in drug discovery and development: an overview

Pharmacophore mapping is one of the major elements of drug design in the absence of structural data of the target receptor. The tool initially applied to discovery of lead molecules now extends to lead optimization. Pharmacophores can be used as queries for retrieving potential leads from structural databases (lead discovery), for designing molecules with specific desired attributes (lead optimization), and for assessing similarity and diversity of molecules using pharmacophore fingerprints. It can also be used to align molecules based on the 3D arrangement of chemical features or to develop predictive 3D QSAR models. This review begins with a brief historical overview of the pharmacophore evolution followed by a coverage of the developments in methodologies for pharmacophore identification over the period from inception of the pharmacophore concept to recent developments of the more sophisticated tools such as Catalyst, GASP, and DISCO. In addition, we present some very recent successes of the widely used pharmacophore generation methods in drug discovery.

Comparative residue interaction analysis (CoRIA): a 3D-QSAR approach to explore the binding contributions of active site residues with ligands

A novel approach termed comparative residue-interaction analysis (CoRIA), emphasizing the trends and principles of QSAR in a ligand-receptor environment has been developed to analyze and predict the binding affinity of enzyme inhibitors. To test this new approach, a training set of 36 COX-2 inhibitors belonging to nine families was selected. The putative binding (bioactive) conformations of inhibitors in the COX-2 active site were searched using the program DOCK. The docked configurations were further refined by a combination of Monte Carlo and simulated annealing methods with the Affinity program. The non-bonded interaction energies of the inhibitors with the individual amino acid residues in the active site were then computed. These interaction energies, plus specific terms describing the thermodynamics of ligand-enzyme binding, were correlated to the biological activity with G/PLS. The various QSAR models obtained were validated internally by cross validation and boot strapping, and externally using a test set of 13 molecules. The QSAR models developed on the CoRIA formalism were robust with good r (2), q (2) and r (pred) (2) values. The major highlights of the method are: adaptation of the QSAR formalism in a receptor setting to answer both the type (qualitative) and the extent (quantitative) of ligand-receptor binding, and use of descriptors that account for the complete thermodynamics of the ligand-receptor binding. The CoRIA approach can be used to identify crucial interactions of inhibitors with the enzyme at the residue level, which can be gainfully exploited in optimizing the inhibitory activity of ligands. Furthermore, it can be used with advantage to guide point mutation studies. As regards the COX-2 dataset, the CoRIA approach shows that improving Coulombic interaction with Pro528 and reducing van der Waals interaction with Tyr385 will improve the binding affinity of inhibitors.

Comparative protein modeling of methionine S-adenosyltransferase (MAT) enzyme from Mycobacterium tuberculosis: a potential target for antituberculosis drug discovery

Mycobacterium tuberculosis (Mtb) is a successful pathogen that overcomes the numerous challenges presented by the immune system of the host. In the last 40 years few anti-TB drugs have been developed, while the drug-resistance problem is increasing; there is thus a pressing need to develop new anti-TB drugs active against both the acute and chronic growth phases of the mycobacterium. Methionine S-adenosyltransferase (MAT) is an enzyme involved in the synthesis of S-adenosylmethionine (SAM), a methyl donor essential for mycolipid biosynthesis. As an anti-TB drug target, Mtb-MAT has been well validated. A homology model of MAT has been constructed using the X-ray structures of E. coli MAT (PDB code: 1MXA) and rat MAT (PDB code: 1QM4) as templates, by comparative protein modeling principles. The resulting model has the correct stereochemistry as gauged from the Ramachandran plot and good three-dimensional (3D) structure compatibility as assessed by the Profiles-3D score. The structurally and functionally important residues (active site) of Mtb-MAT have been identified using the E. coli and rat MAT crystal structures and the reported point mutation data. The homology model conserves the topological and active site features of the MAT family of proteins. The differences in the molecular electrostatic potentials (MEP) of Mtb and human MAT provide evidences that selective and specific Mtb-MAT inhibitors can be designed using the homology model, by the structure-based drug design approaches.