Álvaro Cortés-Cabrera

MM-ISMSA: An Ultrafast and Accurate Scoring Function for Protein-Protein Docking.

An ultrafast and accurate scoring function for protein-protein docking is presented. It includes (1) a molecular mechanics (MM) part based on a 12-6 Lennard-Jones potential; (2) an electrostatic component based on an implicit solvent model (ISM) with individual desolvation penalties for each partner in the protein-protein complex plus a hydrogen bonding term; and (3) a surface area (SA) contribution to account for the loss of water contacts upon protein-protein complex formation. The accuracy and performance of the scoring function, termed MM-ISMSA, have been assessed by (1) comparing the total binding energies, the electrostatic term, and its components (charge-charge and individual desolvation energies), as well as the per residue contributions, to results obtained with well-established methods such as APBSA or MM-PB(GB)SA for a set of 1242 decoy protein-protein complexes and (2) testing its ability to recognize the docking solution closest to the experimental structure as that providing the most favorable total binding energy. For this purpose, a test set consisting of 15 protein-protein complexes with known 3D structure mixed with 10 decoys for each complex was used. The correlation between the values afforded by MM-ISMSA and those from the other methods is quite remarkable (r(2) ∼ 0.9), and only 0.2-5.0 s (depending on the number of residues) are spent on a single calculation including an all vs all pairwise energy decomposition. On the other hand, MM-ISMSA correctly identifies the best docking solution as that closest to the experimental structure in 80% of the cases. Finally, MM-ISMSA can process molecular dynamics trajectories and reports the results as averaged values with their standard deviations. MM-ISMSA has been implemented as a plugin to the widely used molecular graphics program PyMOL, although it can also be executed in command-line mode. MM-ISMSA is distributed free of charge to nonprofit organizations.

Molecular recognition of epothilones by microtubules and tubulin dimers revealed by biochemical and NMR approaches.

The binding of epothilones to dimeric tubulin and to microtubules has been studied by means of biochemical and NMR techniques. We have determined the binding constants of epothilone A (EpoA) and B (EpoB) to dimeric tubulin, which are 4 orders of magnitude lower than those for microtubules, and we have elucidated the conformation and binding epitopes of EpoA and EpoB when bound to tubulin dimers and microtubules in solution. The determined conformation of epothilones when bound to dimeric tubulin is similar to that found by X-ray crystallographic techniques for the binding of EpoA to the Tubulin/RB3/TTL complex; it is markedly different from that reported for EpoA bound to zinc-induced sheets obtained by electron crystallography. Likewise, only the X-ray structure of EpoA bound to the Tubulin/RB3/TTL complex at the luminal site, but not the electron crystallography structure, is compatible with the results obtained by STD on the binding epitope of EpoA bound to dimeric tubulin, thus confirming that the allosteric change (structuring of the M-loop) is the biochemical mechanism of induction of tubulin assembly by epothilones. TR-NOESY signals of EpoA bound to microtubules have been obtained, supporting the interaction with a transient binding site with a fast exchange rate (pore site), consistent with the notion that epothilones access the luminal site through the pore site, as has also been observed for taxanes. Finally, the differences in the tubulin binding affinities of a series of epothilone analogues has been quantitatively explained using the newly determined binding pose and the COMBINE methodology.

Esterase LpEst1 from Lactobacillus plantarum: a novel and atypical member of the αβ hydrolase superfamily of enzymes

The genome of the lactic acid bacterium Lactobacillus plantarum WCFS1 reveals the presence of a rich repertoire of esterases and lipases highlighting their important role in cellular metabolism. Among them is the carboxylesterase LpEst1 a bacterial enzyme related to the mammalian hormone-sensitive lipase, which is known to play a central role in energy homeostasis. In this study, the crystal structure of LpEst1 has been determined at 2.05 Å resolution; it exhibits an αβ-hydrolase fold, consisting of a central β-sheet surrounded by α-helices, endowed with novel topological features. The structure reveals a dimeric assembly not comparable with any other enzyme from the bacterial hormone-sensitive lipase family, probably echoing the specific structural features of the participating subunits. Biophysical studies including analytical gel filtration and ultracentrifugation support the dimeric nature of LpEst1. Structural and mutational analyses of the substrate-binding pocket and active site together with biochemical studies provided insights for understanding the substrate profile of LpEst1 and suggested for the first time the conserved Asp173, which is adjacent to the nucleophile, as a key element in the stabilization of the loop where the oxyanion hole resides.

Comparison of ultra-fast 2D and 3D ligand and target descriptors for side effect prediction and network analysis in polypharmacology.

Some existing computational methods are used to infer protein targets of small molecules and can therefore be used to find new targets for existing drugs, with the goals of re‐directing the molecule towards a different therapeutic purpose or explaining off‐target effects due to multiple targeting. Inherent limitations, however, arise from the fact that chemical analogy is calculated on the basis of common frameworks or scaffolds and also because target information is neglected. The method we present addresses these issues by taking into account 3D information from both the ligand and the target.

ALFA: Automatic Ligand Flexibility Assignment

ALFA is a fast computational tool for the conformational analysis of small molecules that uses a custom-made iterative algorithm to provide a set of representative conformers in an attempt to reproduce the diversity of states in which small molecules can exist, either isolated in solution or bound to a target. The results shown in this work prove that ALFA is fast enough to be integrated into massive high-throughput virtual screening protocols with the aim of incorporating ligand flexibility and also that ALFA reproduces crystallographic X-ray structures of bound ligands with great accuracy. Furthermore, the application includes a graphical user interface that allows its use through the popular molecular graphics program PyMOL to make it accessible to nonexpert users. ALFA is distributed free of charge upon request from the authors.

AtlasCBS: A web server to map and explore chemico-biological space

New approaches are needed that can help decrease the unsustainable failure in small-molecule drug discovery. Ligand Efficiency Indices (LEI) are making a great impact on early-stage compound selection and prioritization. Given a target-ligand database with chemical structures and associated biological affinities/activities for a target, the AtlasCBS server generates two-dimensional, dynamical representations of its contents in terms of LEI. These variables allow an effective decoupling of the chemical (angular) and biological (radial) components. BindingDB, PDBBind and ChEMBL databases are currently implemented. Proprietary datasets can also be uploaded and compared. The utility of this atlas-like representation in the future of drug design is highlighted with some examples. The web server can be accessed at http://ub.cbm.uam.es/atlascbs and https://www.ebi.ac.uk/chembl/atlascbs.

Structural rationale for the cross-resistance of tumor cells bearing the A399V variant of elongation factor eEF1A1 to the structurally unrelated didemnin B, ternatin, nannocystin A and ansatrienin B

At least four classes of structurally distinct natural products with potent antiproliferative activities target the translation elongation factor eEF1A1, which is best known as the G-protein that delivers amino acyl transfer RNAs (aa-tRNAs) to ribosomes during mRNA translation. We present molecular models in atomic detail that provide a common structural basis for the high-affinity binding of didemnin B, ternatin, ansatrienin B and nannocystin A to eEF1A1, as well as a rationale based on molecular dynamics results that accounts for the deleterious effect of replacing alanine 399 with valine. The proposed binding site, at the interface between domains I and III, is eminently hydrophobic and exists only in the GTP-bound conformation. Drug binding at this site is expected to disrupt neither loading of aa-tRNAs nor GTP hydrolysis but would give rise to stabilization of this particular conformational state, in consonance with reported experimental findings. The experimental solution of the three-dimensional structure of mammalian eEF1A1 has proved elusive so far and the highly homologous eEF1A2 from rabbit muscle has been crystallized and solved only as a homodimer in a GDP-bound conformation. Interestingly, in this dimeric structure the large interdomain cavity where the drugs studied are proposed to bind is occupied by a mostly hydrophobic α-helix from domain I of the same monomer. Since binding of this α-helix and any of these drugs to domain III of eEF1A(1/2) is, therefore, mutually exclusive and involves two distinct protein conformations, one intriguing possibility that emerges from our study is that the potent antiproliferative effect of these natural products may arise not only from inhibition of protein synthesis, which is the current dogma, but also from interference with some other non-canonical functions. From this standpoint, this type of drugs could be considered antagonists of eEF1A1/2 oligomerization, a hypothesis that opens up novel areas of research.

Stepwise Simulation of 3,5-Dihydro-5-methylidene-4H-imidazol-4-one (MIO) Biogenesis in Histidine Ammonia-lyase

A 3,5-dihydro-5-methylidene-4H-imidazol-4-one (MIO) electrophilic moiety is post-translationally and autocatalytically generated in homotetrameric histidine ammonia-lyase (HAL) and other enzymes containing the tripeptide Ala-Ser-Gly in a suitably positioned loop. The backbone cyclization step is identical to that taking place during fluorophore formation in green fluorescent protein from the tripeptide Ser-Tyr-Gly, but dehydration, rather than dehydrogenation by molecular oxygen, is the reaction that gives rise to the mature MIO ring system. To gain additional insight into this unique process and shed light on some still unresolved issues, we have made use of extensive molecular dynamics simulations and hybrid quantum mechanics/molecular mechanics calculations implementing the self-consistent charge density functional tight-binding method on a fully solvated tetramer of Pseudomonas putida HAL. Our results strongly support the idea that mechanical compression of the reacting loop by neighboring protein residues in the precursor state is absolutely required to prevent formation of inhibitory main-chain hydrogen bonds and to enforce proper alignment of donor and acceptor orbitals for bond creation. The consideration of the protein environment in our computations shows that water molecules, which have been mostly neglected in previous theoretical work, play a highly relevant role in the reaction mechanism and, more importantly, that backbone cyclization resulting from the nucleophilic attack of the Gly amide lone pair on the π* orbital of the Ala carbonyl precedes side-chain dehydration of the central serine.

Enantioselective oxidation of galactitol 1-phosphate by galactitol-1-phosphate 5-dehydrogenase from Escherichia coli

Galactitol-1-phosphate 5-dehydrogenase (GPDH) is a polyol dehydrogenase that belongs to the medium-chain dehydrogenase/reductase (MDR) superfamily. It catalyses the Zn(2+)- and NAD(+)-dependent stereoselective dehydrogenation of L-galactitol 1-phosphate to D-tagatose 6-phosphate. Here, three crystal structures of GPDH from Escherichia coli are reported: that of the open state of GPDH with Zn(2+) in the catalytic site and those of the closed state in complex with the polyols Tris and glycerol, respectively. The closed state of GPDH reveals no bound cofactor, which is at variance with the conformational transition of the prototypical mammalian liver alcohol dehydrogenase. The main intersubunit-contacting interface within the GPDH homodimer presents a large internal cavity that probably facilitates the relative movement between the subunits. The substrate analogue glycerol bound within the active site partially mimics the catalytically relevant backbone of galactitol 1-phosphate. The glycerol binding mode reveals, for the first time in the polyol dehydrogenases, a pentacoordinated zinc ion in complex with a polyol and also a strong hydrogen bond between the primary hydroxyl group and the conserved Glu144, an interaction originally proposed more than thirty years ago that supports a catalytic role for this acidic residue.

A reverse combination of structure-based and ligand-based strategies for virtual screening

A new approach is presented that combines structure- and ligand-based virtual screening in a reverse way. Opposite to the majority of the methods, a docking protocol is first employed to prioritize small ligands (“fragments”) that are subsequently used as queries to search for similar larger ligands in a database. For a given chemical library, a three-step strategy is followed consisting of (1) contraction into a representative, non-redundant, set of fragments, (2) selection of the three best-scoring fragments docking into a given macromolecular target site, and (3) expansion of the fragments’ structures back into ligands by using them as queries to search the library by means of fingerprint descriptions and similarity criteria. We tested the performance of this approach on a collection of fragments and ligands found in the ZINC database and the directory of useful decoys, and compared the results with those obtained using a standard docking protocol. The new method provided better overall results and was several times faster. We also studied the chemical diversity that both methods cover using an in-house compound library and concluded that the novel approach performs similarly but at a much smaller computational cost.