Complex macrocycle exploration: parallel, heuristic, and constraint-based conformer generation using ForceGen.

ForceGen is a template-free, non-stochastic approach for 2D to 3D structure generation and conformational elaboration for small molecules, including both non-macrocycles and macrocycles. For conformational search of non-macrocycles, ForceGen is both faster and more accurate than the best of all tested methods on a very large, independently curated benchmark of 2859 PDB ligands. In this study, the primary results are on macrocycles, including results for 431 unique examples from four separate benchmarks. These include complex peptide and peptide-like cases that can form networks of internal hydrogen bonds. By making use of new physical movements ("flips" of near-linear sub-cycles and explicit formation of hydrogen bonds), ForceGen exhibited statistically significantly better performance for overall RMS deviation from experimental coordinates than all other approaches. The algorithmic approach offers natural parallelization across multiple computing-cores. On a modest multi-core workstation, for all but the most complex macrocycles, median wall-clock times were generally under a minute in fast search mode and under 2 min using thorough search. On the most complex cases (roughly cyclic decapeptides and larger) explicit exploration of likely hydrogen bonding networks yielded marked improvements, but with calculation times increasing to several minutes and in some cases to roughly an hour for fast search. In complex cases, utilization of NMR data to constrain conformational search produces accurate conformational ensembles representative of solution state macrocycle behavior. On macrocycles of typical complexity (up to 21 rotatable macrocyclic and exocyclic bonds), design-focused macrocycle optimization can be practically supported by computational chemistry at interactive time-scales, with conformational ensemble accuracy equaling what is seen with non-macrocyclic ligands. For more complex macrocycles, inclusion of sparse biophysical data is a helpful adjunct to computation

Random search conformational analysis of piperazine and piperadine analogs of GBR12909 : implicit aqueous solvation effects

The object of this work was to study the effect of solvent on the conformational potential energy surface (PES) of GBR12909 analogs. Local minima on the PES\u27s were found by the Random Search algorithm using the Sybyl molecular modeling package from Tripos, Inc., and an implicit solvent model. Two force-field/charge models were employed in the analysis: the Tripos force field with Gasteiger-Huckel charges and the MMFF94 force field with MMFF94 charges. The effect of solvent on the location of minima in multi-dimensional torsional angle space was studied by comparison to the vacuum phase results. Minima were plotted in torsional angle space using successive pairs of torsional angles. The results showed that, at least for the simple implicit solvent model used here, solvation does not significantly affect the location of the conformational minima for either of the two analogs investigated. With the MMIFF94 model, some of the conformational energy minima were found in a narrower range of torsional angle space in the solvent compared to the vacuum phase, while there were no consistent differences with the Tripos model. One notable exception was the Tripos solvent phase results for the piperadine analog where a cluster of minima was found in a region of torsional angle space where no minima were present in the original vacuum phase results. This is believed to be an anomaly arising from incomplete searching of the conformational PES of the molecule. This was supported by the results of a second vacuum phase random search, which were similar to those of the solvation case. Further study remains to be performed employing other more sophisticated solvation models and conformational search techniques

Modeling of flexible drug-like molecules : qsar of GBR 12909 analog dat/sert selectivity

The dopamine reuptake inhibitor GBR 12909 and related dialkyl piperazine and piperidine analogs have been studied as agonist substitution therapies acting on the dopamine transporter (DAT) to treat cocaine addiction. Undesirable binding to the serotonin transporter (SERT) can vary greatly depending on the specific substituents on the molecule. This study uses Comparative Molecular Field Analysis (CoMFA) and Comparative Molecular Similarity Indices (CoMSIA) techniques to determine a stable and predictive model for DAT/SERT selectivity for a set of flexible GBR 12909 analogs. Families of analogs were constructed from six pairs of naphthyl-substituted piperazine and piperidine templates identified by hierarchical clustering as representative conformers. Three-dimensional quantitative structure-activity relationship (3D-QSAR) studies led to focused models that were stable to y-value scrambling. Test set correlation validation led to one acceptable model (q2=0.508, two components, r2=0.685, average residual = 0.00 for the training set, 0.22 for the extended test set). DAT/SERT selectivities higher than that of the most active compound in the QSAR series were predicted for nine novel compounds. This is the first CoMFA/CoMSIA study of the highly flexible GBR 12909 class of dopamine reuptake inhibitors. Previously, molecular modeling was based on more rigid dopamine reuptake inhibitors, and often only on global energy minimum (GEM) structures. Flexible molecules like GBR 12909 have multiple possible binding conformations, distributed across the potential energy surface in key torsional angle space, which can vary from the GEM by as much as 20 kcal/mol or more. The significance of this study lies in the combining of a clustering technique for identifying representative conformers from a set of low-energy (less than 20 kcal/mol from the GEM) conformers with an extensive 3D-QSAR analysis based on each representative conformer and analogs in a similar potential bioactive conformation

On evaluating moleculardocking methods for pose prediction and enrichment factors.

Four of the most well-known, commercially available docking programs, FlexX, GOLD, GLIDE, and ICM, have been examined for their ligand-docking and virtual-screening capabilities. The relative performance of the programs in reproducing the native ligand conformation from starting SMILES strings for 164 highresolution protein-ligand complexes is presented and compared. Applying only the native scoring functions, the latest versions of these four docking programs were also used to conduct virtual screening for 12 protein targets of therapeutic interest, involving both publicly available structures and AstraZeneca in-house structures. The capability of the four programs to correctly rank-order target-specific active compounds over alternative binders and nonbinders (decoys plus randomly selected compounds) and thereby enrich a small subset of a screening library is compared. Enrichments from the virtual-screening experiments are contrasted with those obtained with alternative 3D shape-matching and 2D similarity database-search methods

Examination of Molecular Recognition in Protein-Ligand Interactions

This dissertation is a compilation of two main projects that were investigated during my thesis research. The first project was a prospective study which identified and characterized drug-like inhibitors of a prototype of bacterial two-component signal transduction response regulator using computational and experimental methods. The second project was the development and validation of a scoring function, PHOENIX, derived using high-resolution structures and calorimetry measurements to predict binding affinities of protein-ligand interactions. Collectively, my thesis research aimed to better understand the underlying driving forces and principles which govern molecular recognition and molecular design. A prospective study coupled computational predictions with experimental validation resulted in the discovery of first-in-class inhibitors targeting a signal transduction module important for bacterial virulence. Development and validation of the PHOENIX scoring function for binding affinity prediction derived using high-resolution structures and calorimetry measurements should guide future molecular recognition studies and endeavors in computer-aided molecular design. To request for an electronic copy of this dissertation, please email the author: yattang at gmail dot com)

In silico studies of the effect of phenolic compounds from grape seed extracts on the activity of phosphoinositide 3-kinase (PI3K) and the farnesoid x receptor (FXR)

In silico studies of the effect of phenolic compounds from grape seed extracts on the activity of phosphoinositide 3-kinase (PI3K) and farnesoid X receptor (FXR)Montserrat Vaqué Marquès En aquesta tesis es pretén aplicar metodologies computacionals (generació de farmacòfors i docking proteïna lligand) en l'àmbit de la nutigenòmica (ciència que pretén entendre, a nivell molecular, com els nutrients afecten la salut). S'aplicaran metodologies in silico per entendre a nivell molecular com productes naturals com els compostos fenòlics presents en la nostra dieta, poden modular la funció d'una diana comportant un efect en la salut. Aquest efecte es creu que podria ser degut a la seva interacció directa amb proteïnes de vies de senyalització molecular o bé a la modificació indirecta de l'expressió gènica. Donat que el coneixement de l'estructura del complex lligand-receptor és bàsic per entendre el mecanisme d'acció d'aquests lligands s'aplica la metodologia docking per predir l'estructura tridimensional del complex. En aquest sentit, un dels programes de docking és AutoGrid/AutoDock (un dels més citats). No obstant, l'automatització d'AutoGrid/AutoDock no és trivial tan per (a) la cerca virtual en una llibreria de lligands contra un grup de possibles receptors, (b) l'ús de flexibilitat, i (c) realitzar un docking a cegues utilitzant tota la superfície del receptor. Per aquest motiu, es dissenya una interfície gràfica de fàcil ús per utilitzar AutoGrid/AutoDock. Blind Docking Tester (BDT) és una aplicació gràfica que s'executa sobre quatre programes escrits en Fortran i que controla les condicions de les execucions d'AutoGrid i AutoDock. BDT pot ser utilitzat per equips d'investigadors en el camp de la química i de ciències de la vida interessats en dur a terme aquest tipus d'experiments però que no tenen suficient habilitats en programació. En la modulació del metabolisme de la glucosa, treballs in vivio i in vitro en el nostre grup de recerca s'han atribuït els efectes beneficiosos de l'extracte de pinyol de raïm en induir captació de glucosa (punt crític pel manteniment de l'homeostasis de la glucosa). No obstant alguns compostos fenòlics no tenen efecte en la captació de la glucosa, d'altres l'inhibeixen reversiblement. En alguns casos aquesta inhibició és el resultat de la competició dels compostos fenòlics amb ATP pel lloc d'unió de l'ATP de la subunitat catalítica de la fosfatidil inositol 3-kinasa (PI3K). Estudis recents amb inhibidors específics d'isoforma han identificat la p110α (la subunitat catalítica de PI3Kα) com la isoforma crucial per la captació de glucosa estimulada per insulina en algunes línies cel·lulars. Els programes computacionals han estat aplicats per tal de correlacionar l'activitat biològica dels compostos fenòlics amb informació estructural per obtenir una relació quantitativa estructura-activitat (3D-QSAR) i obtenir informació dels requeriments estructura-lligand per augmentar l'afinitat i/o selectivitat amb la diana (proteïna). Tot hi haver-se demostrat que l'adició d'extractes de compostos fenòlics en l'aliment pot tenir en general un benefici per la salut, s'ha de tenir en compte que l'estudi 3D-QSAR (construït a partir d'inhibidors sintètics de p110α) prediu que algunes d'aquestes molècules poden agreujar la resistència a la insulina en individus susceptibles dificultant la capatació de glucosa en múscul i teixit adipós i, per tant, produir un efecte secundari indesitjat. Resultats en el nostre grup de recerca han demostrat que compostos fenòlics presents en extractes de llavor de raïm incrementen l'activitat del receptor "farnesoid x receptor" (FXR) de manera dosi depenent quan el lligand natural de FXR (CDCA) és present. Les metodologies in silico, docking i 3D-QSAR, han estat aplicades juntament amb dades biològiques d'agonistes no esteroidals de FXR que s'uneixen a un lloc d'unió proper però diferent al lligand esteroidal 6CDCA. Els resultats han mostrat que els compostos fenòlics no són capaços d'activar FXR per ells mateixos però poden afegir noves interaccions que estabilitzarien la conformació activa de FXR en presència del lligand natural CDCA. Els compostos fenòlics podrien induir canvis conformacionals específics que augmentarien l'activitat de FXR. In silico studies of the effect of phenolic compounds from grape seed extracts on the activity of phosphoinositide 3-kinase (PI3K) and farnesoid X receptor (FXR)Montserrat Vaqué Marquès This thesis was written with the aim of applying computational methods that have already been developed for molecular design and simulation (i.e. pharmacophore generation and protein-ligand docking) to nutrigenomics. So, in silico tools that are routinely used by the pharmaceutical industry to develop drugs have been used to understand, at the molecular level, how natural products such as phenolic compounds (i.e. molecules that are commonly found in fruits and vegetables) can improve health and prevent diseases. Therefore, we first focused on predicting the structure of protein-ligand complexes. The docking algorithms can use the individual structures from receptor and ligand to predict (1) whether they can form a complex and (2) if so, the structure of the resulting complex. This prediction can be made, for instance, with AutoGrid/AutoDock, the most cited docking software in the literature. The automation of AutoGrid/AutoDock is not trivial for tasks such as (1) the virtual screening of a library of ligands against a set of possible receptors; (2) the use of receptor flexibility and (3) making a blind-docking experiment with the whole receptor surface. Therefore, in order to circumvent these limitations, we have designed BDT (i.e. blind-docking tester; http://www.quimica.urv.cat/~pujadas/BDT), an easy-to-use graphic interface for using AutoGrid/AutoDock. BDT is a Tcl/Tk graphic front-end application that runs on top of four Fortran programs and which controls the conditions of the AutoGrid and AutoDock runs. As far as the modulation of the glucose metabolism is concerned, several in vivo and in vitro results obtained by our group have shown that grape seed procyanidin extracts (GSPE) stimulate glucose uptake in 3T3-L1 adipocytes and thus help to maintain their glucose homeostasis. In contrast, it is also well known that although some phenolic compounds do not affect glucose uptake, others reversibly inhibit it in several cell lines. Moreover, for at least some of these phenolic compounds, this inhibition is the result of their competition with ATP for the ATP-binding site in p110α (i.e. the α isoform of the catalytic subunit of phosphoinositide 3-kinase or PI3Kα). Furthermore, recent studies with isoform-specific inhibitors have identified p110α as the crucial isoform for insulin-stimulated glucose-uptake in some cell lines. Therefore, although it has been proved that the addition of phenolic compound extracts to food can have an overall benefit on health, it should be taken into account that some of these molecules may exacerbate insulin resistance in susceptible individuals via impaired glucose uptake in muscle and adipose tissues and, therefore, produce an undesirable side effect. In this context, we have applied computational approaches (i.e. protein-ligand docking and 3D-QSAR) to predict the IC50 (i.e. the concentration that reduces the p110α activity to 50%). Our results agree with previous experimental results and predict that some compounds are potential inhibitors of this enzyme. Recent results in our research group have demonstrated that the phenolic compounds in GSPE increase the activity of the farnesoid X receptor (i.e. FXR) in a dose-dependent way when the natural ligand of FXR (i.e. CDCA) is also present. The phenolic compounds might induce specific conformational changes that increase FXR activity and then contribute to cardioprotection through mechanisms that are independent of their intrinsic antioxidant capacities but that involve direct interaction with FXR to modulate gene expression. Taking into account this hypothesis a 3D-QSAR analysis was made in an attempt to understand how phenolic compounds activate FXR. So, our results explain why phenolic compounds cannot activate FXR by themselves and how they can add new interactions to stabilize the active conformation of FXR when its natural ligand (i.e. CDCA) is present. Therefore, we proposed a mechanism of FXR activation by dietary phenolic compounds in which they may enhance bile acid-bound FXR activity

The molecular modelling of calixarene inclusion complexes

Chapter 1: This consists of a literature review and an introduction into the variety of computational techniques and methodologies available. Chapter 2: The molecular modelling of alkali metal complexes with calixarenes. We have published three papers in this area1. The HyperChem molecular modelling software package has been demonstrated to be a straightforward and computationally undemanding technique that is surprisingly good at reproducing the structures of a range of metallocalixarenes. Chapter 3: The molecular modelling of the enantioselective complexation of chiral amines by chiral calix[4]arenes. With reference to the enantioselective complexation studies we have shown that by using molecular modelling software, an insight can be gained into the interactions that occur on complexation of chiral amines by chiral calixarenes. By observing the hydrogen bonding patterns formed between the host and guest we can postulate why one enantiomer interacts more strongly than the other and, hence, provide a plausible explanation for the observed enantiomeric selectivity observed in the fluorescence quenching experiments. Chapter 4: The catalysis of the Diels Alder reaction by ionic liquids. The first example of the usage of dialkylimidazolium salts (R2 lm+X‘) as heterogeneous and homogeneous Lewis acids catalysts in the Diels-Alder reaction is detailed in a published paper

Los receptores para el reconocimiento de patrones moleculares: aportaciones de la química computacional para el diseño de fármacos y la modulación de la inmunidad innata

Tesis inédita de la Universidad Complutense de Madrid, Facultad de Farmacia, Departamento de Química Orgánica y Farmacéutica, leída el 18/11/2019In this Thesis we have aimed the study of the molecular recognition processes of receptors involved in the innate immunity. More concretely, we have focused in two different types of lectins, Galectins and DC-SIGN, and in Toll-like receptor 4. We have made use of computational techniques, including docking and virtual screening, molecular dynamics simulations, conformational analysis and quantum mechanics calculations. The work has been organized into several chapters that are summarized as follows: Chapter 1 corresponds to the current knowledge and perspectives about receptors related to immunity, in particular: galectins, DC-SIGN, and Toll-like receptor 4, corresponding to the molecular recognition events and modulation by small molecules. Chapter 2 describes the state-of-the-art methods in molecular modeling and computational chemistry applied to the study of molecular recognition processes and drug design...En esta tesis hemos estudiado los procesos reconocimiento molecular de receptores involucrados en la inmunidad innata. Más concretamente, nos hemos centrado en dos tipos diferentes de lectinas, Galectinas y DC-SIGN, y en el receptor Toll-like 4 (TLR4). Hemos utilizado técnicas computacionales, incluyendo docking y cribado virtual, simulaciones de dinámica molecular, análisis conformacional y cálculos de mecánica cuántica. El trabajo se ha organizado en diferentes capítulos que se resumen como sigue: El capítulo 1 corresponde al estado del arte y las perspectivas relacionadas con los estudios de reconocimiento molecular proteína-carbohidrato y diseño de nuevos moduladores con actividad biológica en receptores de la inmunidad, en particular galectinas, DC-SIGN y el receptor Toll-like 4. El capítulo 2 describe el estado actual de los métodos en modelado molecular y química computacional aplicados al estudio de los procesos de reconocimiento molecular y diseño de fármacos...Fac. de FarmaciaTRUEunpu