3D-QSAR studies on 4-Hydroxyphenylpyruvate Dioxygenase inhibitors by comparative molecular field analysis (CoMFA)

A comparative molecular field analysis (CoMFA) of alkanoic acid 3-oxo-cyclohex-1-enyl ester and 2-acylcyclohexane-1,3-dione derivatives of 4-hydroxyphenylpyruvate dioxygenase inhibitors has been performed to determine the factors required for the activity of these compounds. The substrate's conformation abstracted from dynamic modeling of the enzyme-substrate complex was used to build the initial structures of the inhibitors. Satisfactory results were obtained after an all-space searching procedure, performing a leave-one out (LOO) cross-validation study with cross-validation q(2) and conventional r(2) values of 0.779 and 0.989, respectively. The results provide the tools for predicting the affinity of related compounds, and for guiding the design and synthesis of new HPPD ligands with predetermined affinities.

Reduced-size polarized basis sets for calculations of molecular electric properties. III. Second-row atoms

Reduced-size polarized (ZmPolX) basis sets are developed for the second-row atoms X = Si, P, S, and Cl. The generation of these basis sets follows from a simple physical model of the polarization effect of the external electric field which leads to highly compact polarization functions to be added to the chosen initial basis set. The performance of the ZmPolX sets has been investigated in calculations of molecular dipole moments and polarizabilities. Only a small deterioration of the quality of the calculated molecular electric properties has been found. Simultaneously the size of the present reduced-size ZmPolX basis sets is about one-third smaller than that of the usual polarized (PolX) sets. This reduction considerably widens the range of applications of the ZmPolX sets in calculations of molecular dipole moments, dipole polarizabilities, and related properties.

On Entropy-Based Molecular Descriptors: Statistical Analysis of Real and Synthetic Chemical Structures

This paper presents an analysis of entropy-based molecular descriptors. Specifically, we use real chemical structures, as well as synthetic isomeric structures, and investigate properties of and among descriptors with respect to the used data set by a statistical analysis. Our numerical results provide evidence that synthetic chemical structures are notably different to real chemical structures and, hence, should not be used to investigate molecular descriptors. Instead, an analysis based on real chemical structures is favorable. Further, we find strong hints that molecular descriptors can be partitioned into distinct classes capturing complementary information.

Molecular dynamics and principal components analysis of human telomeric quadruplex multimers.

Guanine-rich DNA repeat sequences located at the terminal ends of chromosomal DNA can fold in a sequence-dependent manner into G-quadruplex structures, notably the terminal 150–200 nucleotides at the 3' end, which occur as a single-stranded DNA overhang. The crystal structures of quadruplexes with two and four human telomeric repeats show an all-parallel-stranded topology that is readily capable of forming extended stacks of such quadruplex structures, with external TTA loops positioned to potentially interact with other macromolecules. This study reports on possible arrangements for these quadruplex dimers and tetramers, which can be formed from 8 or 16 telomeric DNA repeats, and on a methodology for modeling their interactions with small molecules. A series of computational methods including molecular dynamics, free energy calculations, and principal components analysis have been used to characterize the properties of these higher-order G-quadruplex dimers and tetramers with parallel-stranded topology. The results confirm the stability of the central G-tetrads, the individual quadruplexes, and the resulting multimers. Principal components analysis has been carried out to highlight the dominant motions in these G-quadruplex dimer and multimer structures. The TTA loop is the most flexible part of the model and the overall multimer quadruplex becoming more stable with the addition of further G-tetrads. The addition of a ligand to the model confirms the hypothesis that flat planar chromophores stabilize G-quadruplex structures by making them less flexible.

Larval Transport Modeling of Deep-Sea Invertebrates Can Aid the Search for Undiscovered Populations

Background: Many deep-sea benthic animals occur in patchy distributions separated by thousands of kilometres, yet because deep-sea habitats are remote, little is known about their larval dispersal. Our novel method simulates dispersal by combining data from the Argo array of autonomous oceanographic probes, deep-sea ecological surveys, and comparative invertebrate physiology. The predicted particle tracks allow quantitative, testable predictions about the dispersal of benthic invertebrate larvae in the south-west Pacific. Principal Findings: In a test case presented here, using non-feeding, non-swimming (lecithotrophic trochophore) larvae of polyplacophoran molluscs (chitons), we show that the likely dispersal pathways in a single generation are significantly shorter than the distances between the three known population centres in our study region. The large-scale density of chiton populations throughout our study region is potentially much greater than present survey data suggest, with intermediate ‘stepping stone’ populations yet to be discovered. Conclusions/Significance: We present a new method that is broadly applicable to studies of the dispersal of deep-sea organisms. This test case demonstrates the power and potential applications of our new method, in generating quantitative, testable hypotheses at multiple levels to solve the mismatch between observed and expected distributions: probabilistic predictions of locations of intermediate populations, potential alternative dispersal mechanisms, and expected population genetic structure. The global Argo data have never previously been used to address benthic biology, and our method can be applied to any non-swimming larvae of the deep-sea, giving information upon dispersal corridors and population densities in habitats that remain intrinsically difficult to assess.

Molecular dynamics studies of the STAT3 homodimer:DNA complex: relationships between STAT3 mutations and protein-DNA recognition.

Signal Transducers and Activators of Transcription (STAT) proteins are a group of latent cytoplasmic transcription factors involved in cytokine signaling. STAT3 is a member of the STAT family and is expressed at elevated levels in a large number of diverse human cancers and is now a validated target for anticancer drug discovery.. Understanding the dynamics of the STAT3 dimer interface, accounting for both protein-DNA and protein-protein interactions, with respect to the dynamics of the latent unphosphorylated STAT3 monomer, is important for designing potential small-molecule inhibitors of the activated dimer. Molecular dynamics (MD) simulations have been used to study the activated STAT3 homodimer:DNA complex and the latent unphosphorylated STAT3 monomer in an explicit water environment. Analysis of the data obtained from MD simulations over a 50 ns time frame has suggested how the transcription factor interacts with DNA, the nature of the conformational changes, and ways in which function may be affected. Examination of the dimer interface, focusing on the protein-DNA interactions, including involvement of water molecules, has revealed the key residues contributing to the recognition events involved in STAT3 protein-DNA interactions. This has shown that the majority of mutations in the DNA-binding domain are found at the protein-DNA interface. These mutations have been mapped in detail and related to specific protein-DNA contacts. Their structural stability is described, together with an analysis of the model as a starting-point for the discovery of novel small-molecule STAT3 inhibitors.

P2Y(1) antagonists: Combining receptor-based modeling and QSAR for a quantitative prediction of the biological activity based on consensus scoring

P2Y(1) is an ADP-activated G protein-coupled receptor (GPCR). Its antagonists impede platelet aggregation in vivo and are potential antithrombotic agents. Combining ligand and structure-based modeling we generated a consensus model (LIST-CM) correlating antagonist structures with their potencies. We docked 45 antagonists into our rhodopsin-based human P2Y(1) homology model and calculated docking scores and free binding energies with the Linear Interaction Energy (LIE) method in continuum-solvent. The resulting alignment was also used to build QSAR based on CoMFA, CoMSIA, and molecular descriptors. To benefit from the strength of each technique and compensate for their limitations, we generated our LIST-CM with a PLS regression based on the predictions of each methodology. A test set featuring untested substituents was synthesized and assayed in inhibition of 2-MeSADP-stimulated PLC activity and in radioligand binding. LIST-CM outperformed internal and external predictivity of any individual model to predict accurately the potency of 75% of the test set.

Defining the molecular basis for the first potent and selective orthosteric agonists of the FFA2 free fatty acid receptor

FFA2 is a G protein-coupled receptor that responds to short chain fatty acids (SCFAs) and has generated interest as a therapeutic target for metabolic and inflammatory conditions. However, definition of its functions has been slowed by a dearth of selective ligands that can distinguish it from the closely related FFA3. At present, the only selective ligands described for FFA2 suffer from either poor potency, altered signaling due to allosteric modes of action, or a lack of function at non-human orthologs of the receptor. To address the need for novel selective ligands, we synthesized two compounds potentially having FFA2 activity and examined the molecular basis of their function. These compounds were confirmed to be potent and selective FFA2 agonists that interact with the orthosteric binding site. A combination of ligand structure-activity relationship, pharmacological analysis, homology modeling, species ortholog comparisons and mutagenesis studies were then employed to define the molecular basis of selectivity and function of these ligands. From this, we identified key residues within both extracellular loop 2 (ECL2) and the transmembrane domain (TM) regions of FFA2 critical for ligand function. One of these ligands was active with reasonable potency at rodent orthologs of FFA2 and demonstrated the role of FFA2 in the regulation of lipolysis in murine 3T3-L1 adipocytes. Together, these findings describe the first potent and selective FFA2 orthosteric agonists and demonstrate key aspects of ligand interaction within the orthosteric binding site of FFA2 that will be invaluable in future ligand development at this receptor.

Two-Ray Modeling of Recombination Effects in the Presence of a Phase Conjugation Array

This article describes by means of a simple model how signal recombination effects behave under the influence of phase conjugating retrodirective array (RDA) technology. A two-ray ground reflection model is used to predict the operational advantages of RDA technology in multipath rich environments. The simulation results show that advantageous signal recombination occurs due to automatic self-phasing. As the number of elements in the RDA increases, the fading effect normally observed due to out of phase multipath signal is mitigated to the extent that the system approaches that of one operating in a free space environment. © 2013 Wiley Periodicals, Inc. Microwave Opt Technol Lett 55:1987–1989, 2013

Variance-Constrained Capacity of the Molecular Timing Channel with Synchronization Error

Molecular communication is set to play an important role in the design of complex biological and chemical systems. An important class of molecular communication systems is based on the timing channel, where information is encoded in the delay of the transmitted molecule - a synchronous approach. At present, a widely used modeling assumption is the perfect synchronization between the transmitter and the receiver. Unfortunately, this assumption is unlikely to hold in most practical molecular systems. To remedy this, we introduce a clock into the model - leading to the molecular timing channel with synchronization error. To quantify the behavior of this new system, we derive upper and lower bounds on the variance-constrained capacity, which we view as the step between the mean-delay and the peak-delay constrained capacity. By numerically evaluating our bounds, we obtain a key practical insight: the drift velocity of the clock links does not need to be significantly larger than the drift velocity of the information link, in order to achieve the variance-constrained capacity with perfect synchronization.

Structure and thermal properties of salicylate-based-protic ionic liquids as new heat storage media. COSMO-RS structure characterization and modeling of heat capacities

During this research, we present a study on the thermal properties, such as the melting, cold crystallization, and glass transition temperatures as well as heat capacities from 293.15 K to 323.15 K of nine in-house synthesized protic ionic liquids based on the 3-(alkoxymethyl)-1H-imidazol-3-ium salicylate ([H-Im-C₁OC_n][Sal]) with n = 3–11. The 3D structures, surface charge distributions and COSMO volumes of all investigated ions are obtained by combining DFT calculations and the COSMO-RS methodology. The heat capacity data sets as a function of temperature of the 3-(alkoxymethyl)-1H-imidazol-3-ium salicylate are then predicted using the methodology originally proposed in the case of ionic liquids by Ge et al. 3-(Alkoxymethyl)-1H-imidazol-3-ium salicylate based ionic liquids present specific heat capacities higher in many cases than other ionic liquids that make them suitable as heat storage media and in heat transfer processes. It was found experimentally that the heat capacity increases linearly with increasing alkyl chain length of the alkoxymethyl group of 3-(alkoxymethyl)-1H-imidazol-3-ium salicylate as was expected and predicted using the Ge et al. method with an overall relative absolute deviation close to 3.2% for temperatures up to 323.15 K.

Modeling Klebsiella pneumoniae pathogenesis by infection of the wax moth Galleria mellonella

The implementation of infection models that approximate human disease is essential for understanding pathogenesis at the molecular level and for testing new therapies before they are entered into clinical stages. Insects are increasingly being used as surrogate hosts because they share, with mammals, essential aspects of the innate immune response to infections. We examined whether the larva of the wax moth Galleria mellonella could be used as a host model to conceptually approximate Klebsiella pneumoniae-triggered pneumonia. We report that the G. mellonella model is capable of distinguishing between pathogenic and nonpathogenic Klebsiella strains. Moreover, K. pneumoniae infection of G. mellonella models some of the known features of Klebsiella-induced pneumonia, i.e., cell death associated with bacterial replication, avoidance of phagocytosis by phagocytes, and the attenuation of host defense responses, chiefly the production of antimicrobial factors. Similar to the case for the mouse pneumonia model, activation of innate responses improved G. mellonella survival against subsequent Klebsiella challenge. Virulence factors necessary in the mouse pneumonia model were also implicated in the Galleria model. We found that mutants lacking capsule polysaccharide, lipid A decorations, or the outer membrane proteins OmpA and OmpK36 were attenuated in Galleria. All mutants activated G. mellonella defensive responses. The Galleria model also allowed us to monitor Klebsiella gene expression. The expression levels of cps and the loci implicated in lipid A remodeling peaked during the first hours postinfection, in a PhoPQ- and PmrAB-governed process. Taken together, these results support the utility of G. mellonella as a surrogate host for assessing infections with K. pneumoniae.

Magnitude differences in bioactive compounds, chemical functional groups, fatty acid profiles, nutrient degradation and digestion, molecular structure and metabolic characteristics of protein in newly developed yellow-seeded and black-seeded canola lines.

Recently, new lines of yellow-seeded (CS-Y) and black-seeded canola (CS-B) have been developed with chemical and structural alteration through modern breeding technology. However, no systematic study was found on the bioactive compounds, chemical functional groups, fatty acid profiles, inherent structure, nutrient degradation and absorption, or metabolic characteristics between the newly developed yellow- and black-seeded canola lines. This study aimed to systematically characterize chemical, structural, and nutritional features in these canola lines. The parameters accessed include bioactive compounds and antinutrition factors, chemical functional groups, detailed chemical and nutrient profiles, energy value, nutrient fractions, protein structure, degradation kinetics, intestinal digestion, true intestinal protein supply, and feed milk value. The results showed that the CS-Y line was lower (P ≤ 0.05) in neutral detergent fiber (122 vs 154 g/kg DM), acid detergent fiber (61 vs 99 g/kg DM), lignin (58 vs 77 g/kg DM), nonprotein nitrogen (56 vs 68 g/kg DM), and acid detergent insoluble protein (11 vs 35 g/kg DM) than the CS-B line. There was no difference in fatty acid profiles except C20:1 eicosenoic acid content (omega-9) which was in lower in the CS-Y line (P < 0.05) compared to the CS-B line. The glucosinolate compounds differed (P < 0.05) in terms of 4-pentenyl, phenylethyl, 3-CH3-indolyl, and 3-butenyl glucosinolates (2.9 vs 1.0 μmol/g) between the CS-Y and CS-B lines. For bioactive compounds, total polyphenols tended to be different (6.3 vs 7.2 g/kg DM), but there were no differences in erucic acid and condensed tannins with averages of 0.3 and 3.1 g/kg DM, respectively. When protein was portioned into five subfractions, significant differences were found in PA, PB1 (65 vs 79 g/kg CP), PB2, and PC fractions (10 vs 33 g/kg CP), indicating protein degradation and supply to small intestine differed between two new lines. In terms of protein structure spectral profile, there were no significant differences in functional groups of amides I and II, α helix, and β-sheet structure as well as their ratio between the two new lines, indicating no difference in protein structure makeup and conformation between the two lines. In terms of energy values, there were significant differences in total digestible nutrient (TDN; 149 vs 133 g/kg DM), metabolizable energy (ME; 58 vs 52 MJ/kg DM), and net energy for lactation (NEL; 42 vs 37 MJ/kg DM) between CS-Y and CS-B lines. For in situ rumen degradation kinetics, the two lines differed in soluble fraction (S; 284 vs 341 g/kg CP), potential degradation fraction (D; 672 vs 590 g/kg CP), and effective degraded organic matter (EDOM; 710 vs 684 g/kg OM), but no difference in degradation rate. CS-Y had higher digestibility of rumen bypass protein in the intestine than CS-B (566 vs 446 g/kg of RUP, P < 0.05). Modeling nutrient supply results showed that microbial protein synthesis (MCP; 148 vs 171 g/kg DM) and rumen protein degraded balance (DPB; 108 vs 127 g/kg DM) were lower in the CS-Y line, but there were no differences in total truly digested protein in small intestine (DVE) and feed milk value (FMV) between the two lines. In conclusion, the new yellow line had different nutritional, chemical, and structural features compared to the black line. CS-Y provided better nutrient utilization and availability.

Machine learning methods for quantitative structure-property relationship modeling

Tese de doutoramento, Informática (Bioinformática), Universidade de Lisboa, Faculdade de Ciências, 2014