Search

Ligand-based virtual screening (LBVS) methods are widely used to explore the vast chemical space in the search of novel compounds resorting to a variety of properties encoded in 1D, 2D or 3D descriptors. The success of 3D-LBVS is affected by the overlay of molecular pairs, thus making selection of the template compound, search of accessible conformational space and choice of the query conformation to be potential factors that modulate the successful retrieval of actives. This study examines the impact of adopting different choices for the query conformation of the template, paying also attention to the influence exerted by the structural similarity between templates and actives. The analysis is performed using PharmScreen, a 3D LBVS tool that relies on similarity measurements of the hydrophobic/philic pattern of molecules, and Phase Shape, which is based on the alignment of atom triplets followed by refinement of the volume overlap. The study is performed for the original DUD-E+ database and a Morgan Fingerprint filtered version (denoted DUD-E+-Diverse; available in ), which was prepared to minimize the 2D resemblance between template and actives. Although in most cases the query conformation exhibits a mild influence on the overall performance, a critical analysis is made to disclose factors, such as the content of structural features between template and actives and the induction of conformational strain in the template, that underlie the drastic impact of the query definition in the recovery of actives for certain targets. The findings of this research also provide valuable guidance for assisting the selection of the query definition in 3D LBVS campaigns. [ABSTRACT FROM AUTHOR]

In recent years the use of partition systems other than the widely used biphasic n-octanol/water has received increased attention to gain insight into the molecular features that dictate the lipophilicity of compounds. Thus, the difference between n-octanol/water and toluene/water partition coefficients has proven to be a valuable descriptor to study the propensity of molecules to form intramolecular hydrogen bonds and exhibit chameleon-like properties that modulate solubility and permeability. In this context, this study reports the experimental toluene/water partition coefficients (log Ptol/w) for a series of 16 drugs that were selected as an external test set in the framework of the Statistical Assessment of the Modeling of Proteins and Ligands (SAMPL) blind challenge. This external set has been used by the computational community to calibrate their methods in the current edition (SAMPL9) of this contest. Furthermore, the study also investigates the performance of two computational strategies for the prediction of log Ptol/w. The first relies on the development of two machine learning (ML) models, which are built up by combining the selection of 11 molecular descriptors in conjunction with either the multiple linear regression (MLR) or the random forest regression (RFR) model to target a dataset of 252 experimental log Ptol/w values. The second consists of the parametrization of the IEF-PCM/MST continuum solvation model from B3LYP/6-31G(d) calculations to predict the solvation free energies of 163 compounds in toluene and benzene. The performance of the ML and IEF-PCM/MST models has been calibrated against external test sets, including the compounds that define the SAMPL9 log Ptol/w challenge. The results are used to discuss the merits and weaknesses of the two computational approaches. [ABSTRACT FROM AUTHOR]

Fatty acid amide hydrolase (FAAH) plays a key role in the control of cannabinoid signaling and it represents a promising therapeutic strategy for the treatment of a wide range of diseases, including neuropathic pain and chronic inflammation. Starting from kinetics experiments carried out in our previous work for the most potent inhibitor 2-amino-3-chloropyridine amide (TPA14), we have investigated its non-competitive mechanism of action using molecular dynamics, thermodynamic integration and QM-MM/GBSA calculations. The computational studies highlighted the impact of mutations on the receptor binding pockets and elucidated the molecular basis of the non-competitive inhibition mechanism of TPA14, which prevents the endocannabinoid anandamide (AEA) from reaching its pro-active conformation. Our study provides a rationale for the design of non-competitive potent FAAH inhibitors for the treatment of neuropathic pain and chronic inflammation. [ABSTRACT FROM AUTHOR]

Within the scope of SAMPL7 challenge for predicting physical properties, the Integral Equation Formalism of the Miertus-Scrocco-Tomasi (IEFPCM/MST) continuum solvation model has been used for the blind prediction of n-octanol/water partition coefficients and acidity constants of a set of 22 and 20 sulfonamide-containing compounds, respectively. The log P and pKa were computed using the B3LPYP/6-31G(d) parametrized version of the IEFPCM/MST model. The performance of our method for partition coefficients yielded a root-mean square error of 1.03 (log P units), placing this method among the most accurate theoretical approaches in the comparison with both globally (rank 8th) and physical (rank 2nd) methods. On the other hand, the deviation between predicted and experimental pKa values was 1.32 log units, obtaining the second best-ranked submission. Though this highlights the reliability of the IEFPCM/MST model for predicting the partitioning and the acid dissociation constant of drug-like compounds compound, the results are discussed to identify potential weaknesses and improve the performance of the method. [ABSTRACT FROM AUTHOR]

Many (poly‐)phenolic natural products, for example, curcumin and taxifolin, have been studied for their activity against specific hallmarks of neurodegeneration, such as amyloid‐β 42 (Aβ42) aggregation and neuroinflammation. Due to their drawbacks, arising from poor pharmacokinetics, rapid metabolism, and even instability in aqueous medium, the biological activity of azobenzene compounds carrying a pharmacophoric catechol group, which have been designed as bioisoteres of curcumin has been examined. Molecular simulations reveal the ability of these compounds to form a hydrophobic cluster with Aβ42, which adopts different folds, affecting the propensity to populate fibril‐like conformations. Furthermore, the curcumin bioisosteres exceeded the parent compound in activity against Aβ42 aggregation inhibition, glutamate‐induced intracellular oxidative stress in HT22 cells, and neuroinflammation in microglial BV‐2 cells. The most active compound prevented apoptosis of HT22 cells at a concentration of 2.5 μm (83 % cell survival), whereas curcumin only showed very low protection at 10 μm (21 % cell survival). [ABSTRACT FROM AUTHOR]

Amyloids are characterized by their capacity to bind Congo red (CR), one of the most used amyloid‐specific dyes. The structural features of CR binding were unknown for years, mainly because of the lack of amyloid structures solved at high resolution. In the last few years, solid‐state NMR spectroscopy enabled the determination of the structural features of amyloids, such as the HET‐s prion forming domain (HET‐s PFD), which also has recently been used to determine the amyloid–CR interface at atomic resolution. Herein, we combine spectroscopic data with molecular docking, molecular dynamics, and excitonic quantum/molecular mechanics calculations to examine and rationalize CR binding to amyloids. In contrast to a previous assumption on the binding mode, our results suggest that CR binding to the HET‐s PFD involves a cooperative process entailing the formation of a complex with 1:1 stoichiometry. This provides a molecular basis to explain the bathochromic shift in the maximal absorbance wavelength when CR is bound to amyloids. [ABSTRACT FROM AUTHOR]

Amyloids are characterized by their capacity to bind Congo red (CR), one of the most used amyloid‐specific dyes. The structural features of CR binding were unknown for years, mainly because of the lack of amyloid structures solved at high resolution. In the last few years, solid‐state NMR spectroscopy enabled the determination of the structural features of amyloids, such as the HET‐s prion forming domain (HET‐s PFD), which also has recently been used to determine the amyloid–CR interface at atomic resolution. Herein, we combine spectroscopic data with molecular docking, molecular dynamics, and excitonic quantum/molecular mechanics calculations to examine and rationalize CR binding to amyloids. In contrast to a previous assumption on the binding mode, our results suggest that CR binding to the HET‐s PFD involves a cooperative process entailing the formation of a complex with 1:1 stoichiometry. This provides a molecular basis to explain the bathochromic shift in the maximal absorbance wavelength when CR is bound to amyloids. [ABSTRACT FROM AUTHOR]

The IEFPCM/MST continuum solvation model is used for the blind prediction of n-octanol/water partition of a set of 11 fragment-like small molecules within the SAMPL6 Part II Partition Coefficient Challenge. The partition coefficient of the neutral species (log P) was determined using an extended parametrization of the B3LYP/6-31G(d) version of the Miertus–Scrocco–Tomasi continuum solvation model in n-octanol. Comparison with the experimental data provided for partition coefficients yielded a root-mean square error (rmse) of 0.78 (log P units), which agrees with the accuracy reported for our method (rmse = 0.80) for nitrogen-containing heterocyclic compounds. Out of the 91 sets of log P values submitted by the participants, our submission is within those with an rmse < 1 and among the four best ranked physical methods. The largest errors involve three compounds: two with the largest positive deviations (SM13 and SM08), and one with the largest negative deviations (SM15). Here we report the potentiometric determination of the log P for SM13, leading to a value of 3.62 ± 0.02, which is in better agreement with most empirical predictions than the experimental value reported in SAMPL6. In addition, further inclusion of several conformations for SM08 significantly improved our results. Inclusion of these refinements led to an overall error of 0.51 (log P units), which supports the reliability of the IEFPCM/MST model for predicting the partitioning of neutral compounds. [ABSTRACT FROM AUTHOR]

A first-order perturbation treatment of the polarization contribution to the free energy of hydration is presented. Very simple expressions for the computation of the total electrostatic free energy of solvation, the polarization contribution, and its components (distortion and stabilization) are derived. These equations can be used with either continuum (quantum and classical) or discrete approaches. The reliability of these equations is examined by comparison with rigorous expressions derived previously within the framework of the self-consistent reaction field theory. Indeed, the suitability of the classical expressions for the distortion and stabilization terms used in the context of discrete strategies has been explored by comparison between self-consistent reaction field and molecular dynamics–free-energy perturbation results. The excellent agreement found between the two techniques allows us to envisage a procedure to account for the polarization in force field-derived methods. © 1995 American Institute of Physics. [ABSTRACT FROM AUTHOR]

The magnitude of the polarization contribution to the free energy of hydration is analyzed using several approaches. Monte Carlo quantum mechanical-molecular mechanical (MC-QM/MM), self-consistent reaction field (SCRF), and classical free energy perturbation calculations were used to obtain independent estimates of the polarization contribution to the free energy of hydration. The results derived from these techniques are in reasonable agreement, which supports the reliability of the calculated polarization values. The polarization contribution to the free energy of hydration was partitioned into components to elucidate the nature of the polarization in the solvation process. Strategies to capture the polarization effects in force-field simulations are discussed. © 1995 American Institute of Physics. [ABSTRACT FROM AUTHOR]

The effect of solvation, as studied by the self-consistent reaction field procedure, on the shapes, sizes, and anisotropies of eight small prototypical anions (OH-, NH-2, CH-3, CN-, SH-, N-3, NO-2, and NO-3) has been studied at the ab initio 6–31++G(d) level. Both the general molecular electrostatic potential distribution and the molecular electrostatic potential topography have been used for examining the effect of solvation. The results demonstrate both the complexity and specificity of the hydration effect on the solute charge distribution. It is observed that in general anions tend to shrink upon solvation. The molecular electrostatic potential, in general, becomes more negative and the electron density is increased upon solvation. The effect of the solvation on the chemical reactivity and anisotropy of the anions in aqueous environment is discussed. [ABSTRACT FROM AUTHOR]

A Monte Carlo quantum mechanical and molecular mechanical (QM/MM) simulation method using the combined Austin model 1 (AM1) and TIP3P potential is proposed to evaluate the average molecular electrostatic potentials (MEP) in aqueous solution. It was demonstrated that the average MEP in water can provide a quantitative account of the solvent polarization effect on molecular dipole moments and atomic charges. The computed induced dipole moments from the Monte Carlo AM1/TIP3P simulation method were found to be in good agreement with those predicted by continuum self-consistent reaction field theory using the AM1 and 6-31G(d) wave functions. The findings also indicate that the minimal basis set adopted by the AM1 method can capture the solvent effect effectively. Finally, good correlations in the MEP-fitted charges were obtained between different methods. These results should be useful for developing polarizable intermolecular potential functions (PIPF). [ABSTRACT FROM AUTHOR]

Copyright of Retos: Nuevas Perspectivas de Educación Física, Deporte y Recreación is the property of Federacion Espanola de Asociaciones de Docentes de Educacion Fisica and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)

Abstract: Multiple multicomponent reactions rapidly assemble complex structures. Despite being very productive, the lack of selectivity and the reduced number of viable transformations restrict their general application in synthesis. Hereby, we describe a rationale for a selective version of these processes based in the preferential generation of intermediates which are less reactive than the initial substrates. In this way, applying the Groebke–Blackburn–Bienaymé reaction on a range of α‐polyamino‐polyazines, we prepared a family compact heterocyclic scaffolds with relevant applications in medicinal and biological chemistry (live cell imaging probes, selective binders for DNA quadruplexes, and antiviral agents against human adenoviruses). The approach has general character and yields complex molecular targets in a selective, tunable and direct manner. [ABSTRACT FROM AUTHOR]

The light-harvesting mechanisms of cryptophyte antenna complexes have attracted considerable attention due to their ability to exhibit maximal photosynthetic activity under very low-light conditions and to display several colors, as well as the observation of vibronic coherent features in their two-dimensional electronic spectra. However, detailed investigations on the interplay between the protein environment and their light-harvesting properties are hampered by the uncertainty related to the protonation state of the underlying bilin pigments. Here we study the protonation preferences of four types of bilin pigments including 15,16-dihydrobiliverdin (DBV), phycoerythrobilin (PEB), phycocyanobilin (PCB) and mesobiliverdin (MBV), which are found in phycoerythrin PE545 and phycocyanin PC577, PC612, PC630 and PC645 complexes. We apply quantum chemical calculations coupled to continuum solvation calculations to predict the intrinsic acidity of bilins in aqueous solution, and then combine molecular dynamics simulations with empirical pKa estimates to investigate the impact of the local protein environment on the acidity of the pigments. We also report measurements of the absorption spectra of the five complexes in a wide range of pH in order to validate our simulations and investigate possible changes in the light harvesting properties of the complexes in the range of physiological pH found in the lumen (pH ∼ 5–7). The results suggest a pKa ＞ 7 for DBV and MBV pigments in the α polypeptide chains of PE545 and PC630/PC645 complexes, which are not coordinated to a negatively charged amino acid. For the other PEB, DBV and PCB pigments, which interact with a Glu or Asp side chain, higher pKa values (pKa ＞ 8) are estimated. Overall, the results support a preferential population of the fully protonated state for bilins in cryptophyte complexes under physiological conditions regardless of the specific type of pigment and local protein environment. [ABSTRACT FROM AUTHOR]

Abstract: Inhibition of abnormal protein self‐aggregation is an attractive strategy against amyloidogenic diseases, but has found limited success due to the complexity of protein self‐assembly, the absence of fully reproducible aggregation assays, and the scarce knowledge of the inhibition mechanisms by small molecules. In this context, catechol‐containing compounds may lead to covalent adducts with amyloid fibrils that interfere with the aggregation process. In particular, the covalent adduct formed between the oxidized form of (+)‐taxifolin and an β‐amyloid (Aβ42) suggests the involvement of a specific recognition motif that enables the chemical reaction with Aβ42. In this study, we have examined the mechanisms implicated in the aza‐Michael addition of the o‐quinone species of (+)‐taxifolin with Aβ42 fibrils. The results support the binding of (+)‐taxifolin to the hydrophobic groove delimited by the edges defined by Lys16 and Glu22 residues in the fibril. The chemical reaction proceeds through the nucleophilic attack of the deprotonated amino group of a Lys16 residue in a process activated by the interaction between the o‐quinone ring with a vicinal Lys16 residue, as well as by a water‐assisted proton transfer, which is the rate‐limiting step of the reaction. This specific inhibition mechanism, which may explain the enhanced anti‐aggregating activity of oxidized flavonoids compared to fresh compounds, holds promise for developing disease‐modifying therapies. [ABSTRACT FROM AUTHOR]

The article presents research briefs related to medicinal chemistry, as of January 2018. Topics include a study on an allosteric inhibitor with the most complex mode of action to combat tuberculosis; use of selective indazole ether-based glucocorticoid receptor modulators for treating asthma and chronic obstructive pulmonary disease; and developing potent proteolysis-targeting chimeras (protacs) in anticancer drug discovery.

Methyl-substituted Nazarov reagent 4 reacts stereoselectively with Nind-Boc indoloquinolizidine lactams to give the expected 3-H/15-H cis pentacyclic yohimbine-type adducts when using 1,8-diazabicyclo[5.4.0]undec-7-ene (DBU) as the base. However, a dramatic change of the facial selectivity was observed when the reaction was performed in the presence of Cs2CO3, leading to the corresponding trans adducts. This annulation is the key step of a stereocontrolled synthesis of the 17acarba analogue of the heteroyohimbine alkaloid akuammigine. Theoretical calculations were used to rationalize the facial selectivity of these double Michael addition reactions. [ABSTRACT FROM AUTHOR]

The critical role of BACE-1 in the formation of neurotoxic ß-amyloid peptides in the brain makes it an attractive target for an efficacious treatment of Alzheimer’s disease. However, the development of clinically useful BACE-1 inhibitors has proven to be extremely challenging. In this study we examine the binding mode of a novel potent inhibitor (compound 1, with IC50 80 nM) designed by synergistic combination of two fragments—huprine and rhein—that individually are endowed with very low activity against BACE-1. Examination of crystal structures reveals no appropriate binding site large enough to accommodate 1. Therefore we have examined the conformational flexibility of BACE-1 through extended molecular dynamics simulations, paying attention to the highly flexible region shaped by loops 8–14, 154–169 and 307–318. The analysis of the protein dynamics, together with studies of pocket druggability, has allowed us to detect the transient formation of a secondary binding site, which contains Arg307 as a key residue for the interaction with small molecules, at the edge of the catalytic cleft. The formation of this druggable “floppy” pocket would enable the binding of multisite inhibitors targeting both catalytic and secondary sites. Molecular dynamics simulations of BACE-1 bound to huprine-rhein hybrid compounds support the feasibility of this hypothesis. The results provide a basis to explain the high inhibitory potency of the two enantiomeric forms of 1, together with the large dependence on the length of the oligomethylenic linker. Furthermore, the multisite hypothesis has allowed us to rationalize the inhibitory potency of a series of tacrine-chromene hybrid compounds, specifically regarding the apparent lack of sensitivity of the inhibition constant to the chemical modifications introduced in the chromene unit. Overall, these findings pave the way for the exploration of novel functionalities in the design of optimized BACE-1 multisite inhibitors. [ABSTRACT FROM AUTHOR]

Cation–π interactions of aromatic rings and positively charged groups are among the most important interactions in structural biology. The role and energetic characteristics of these interactions are well established. However, the occurrence of cation–π–cation interactions is an unexpected motif, which raises intriguing questions about its functional role in proteins. We present a statistical analysis of the occurrence, composition and geometrical preferences of cation–π–cation interactions identified in a set of non-redundant protein structures taken from the Protein Data Bank. Our results demonstrate that this structural motif is observed at a small, albeit non-negligible frequency in proteins, and suggest a preference to establish cation–π–cation motifs with Trp, followed by Tyr and Phe. Furthermore, we have found that cation–π–cation interactions tend to be highly conserved, which supports their structural or functional role. Finally, we have performed an energetic analysis of a representative subset of cation–π–cation complexes combining quantum-chemical and continuum solvation calculations. Our results point out that the protein environment can strongly screen the cation–cation repulsion, leading to an attractive interaction in 64% of the complexes analyzed. Together with the high degree of conservation observed, these results suggest a potential stabilizing role in the protein fold, as demonstrated recently for a miniature protein (Craven et al., J. Am. Chem. Soc. 2016, 138, 1543). From a computational point of view, the significant contribution of non-additive three-body terms challenges the suitability of standard additive force fields for describing cation–π–cation motifs in molecular simulations. [ABSTRACT FROM AUTHOR]

Aggregation is a major problem for the anti-microbial photodynamic applications of hydrophobic photosensitizers since it strongly reduces the amount of singlet oxygen generated in aqueous solutions. Binding of hypericin (Hyp) to the milk whey protein β-lactoglobulin (βLG), occurring at the two hydrophobic cavities located at the interface of the protein homodimer, can be exploited to confer water-solubility and biocompatibility to the photosensitizer. The introduction of a small amount of the organic cosolvent dimethyl sulfoxide (DMSO) leads to a remarkable improvement of the photophysical properties of the complex Hyp–βLG by increasing its fluorescence emission and singlet oxygen photosensitization quantum yields. Surprisingly, the ability of the complex to photo-inactivate bacteria of the strain Staphylococcus aureus is strongly reduced in the presence of DMSO, despite the higher yield of photosensitization. The reasons for this apparently contradictory behavior are investigated, providing new insights into the use of carrier systems for hydrophobic photosensitizers. [ABSTRACT FROM AUTHOR]

Trimethylsilyl chloride is an efficient activating agent for azines in isocyanide-based reactions, which then proceed through a key insertion of the isocyanide into a N−Si bond. The reaction is initiated by N activation of the azine, followed by nucleophilic attack of an isocyanide in a Reissert-type process. Finally, a second equivalent of the same or a different isocyanide inserts into the N−Si bond leading to the final adduct. The use of distinct nucleophiles leads to a variety of α-substituted dihydroazines after a selective cascade process. Based on computational studies, a mechanistic hypothesis for the course of these reactions was proposed. The resulting products exhibit significant activity against Trypanosoma brucei and T. cruzi, featuring favorable drug-like properties and safety profiles. [ABSTRACT FROM AUTHOR]

Trimethylsilyl chloride is an efficient activating agent for azines in isocyanide-based reactions, which then proceed through a key insertion of the isocyanide into a N−Si bond. The reaction is initiated by N activation of the azine, followed by nucleophilic attack of an isocyanide in a Reissert-type process. Finally, a second equivalent of the same or a different isocyanide inserts into the N−Si bond leading to the final adduct. The use of distinct nucleophiles leads to a variety of α-substituted dihydroazines after a selective cascade process. Based on computational studies, a mechanistic hypothesis for the course of these reactions was proposed. The resulting products exhibit significant activity against Trypanosoma brucei and T. cruzi, featuring favorable drug-like properties and safety profiles. [ABSTRACT FROM AUTHOR]

Since the development of structure-activity relationships about 50 years ago, 3D-QSAR methods belong to the most refined ligand-based in silico techniques for prediction of biological data using physicochemical molecular fields. In this scenario, this study reports the development and validation of quantum mechanical (QM)-based hydrophobic descriptors derived from the parametrized MST continuum solvation model to be used in 3D-QSAR studies within the framework of the Hydrophobic Pharmacophore (HyPhar) method. To this end, five sets of compounds reported in the literature (dopamine D2/D4 antagonists, antifungal 2-aryl-4-chromanones, and inhibitors of GSK-3, cruzain and thermolysin) have been revisited. The results derived from the QM/MST-based hydrophobic descriptors have been compared with previous CoMFA and CoMSIA studies, and examined in light of the available X-ray crystallographic structures of the targets. The analysis reveals that the combination of electrostatic and nonelectrostatic components of the octanol/water partition coefficient yields pharmacophoric models fully comparable with the predictive potential of standard 3D-QSAR techniques. Moreover, the graphical representation of the hydrophobic maps provides a direct linkage with the pattern of interactions found in crystallographic structures. Overall, the introduction of the QM/MST-based descriptors, which could be easily adapted to other continuum solvation formalisms, paves the way to novel computational strategies for disclosing structure-activity relationships in drug design. © 2016 Wiley Periodicals, Inc. [ABSTRACT FROM AUTHOR]

Background: Combined fatty acid amide hydrolase (FAAH) and cyclooxygenase (COX) inhibition is a promising approach for pain-relief. The Flu-AM1 and Ibu-AM5 derivatives of flurbiprofen and ibuprofen retain similar COX-inhibitory properties and are more potent inhibitors of FAAH than the parent compounds. However, little is known as to the nature of their interaction with FAAH, or to the importance of their chirality. This has been explored here. Methodology/Principal Findings: FAAH inhibitory activity was measured in rat brain homogenates and in lysates expressing either wild-type or FAAHT488A-mutated enzyme. Molecular modelling was undertaken using both docking and molecular dynamics. The (R)- and (S)-enantiomers of Flu-AM1 inhibited rat FAAH with similar potencies (IC50 values of 0.74 and 0.99 μM, respectively), whereas the (S)-enantiomer of Ibu-AM5 (IC50 0.59 μM) was more potent than the (R)-enantiomer (IC50 5.7 μM). Multiple inhibition experiments indicated that both (R)-Flu-AM1 and (S)-Ibu-AM5 inhibited FAAH in a manner mutually exclusive to carprofen. Computational studies indicated that the binding site for the Flu-AM1 and Ibu-AM5 enantiomers was located between the acyl chain binding channel and the membrane access channel, in a site overlapping the carprofen binding site, and showed a binding mode in line with that proposed for carprofen and other non-covalent ligands. The potency of (R)-Flu-AM1 was lower towards lysates expressing FAAH mutated at the proposed carprofen binding area than in lysates expressing wild-type FAAH. Conclusions/Significance: The study provides kinetic and structural evidence that the enantiomers of Flu-AM1 and Ibu-AM5 bind in the substrate channel of FAAH. This information will be useful in aiding the design of novel dual-action FAAH: COX inhibitors. [ABSTRACT FROM AUTHOR]

Recent advances in polarizable force fields have revealed that major reparameterization is necessary when the polarization energy is treated explicitly. This study is focused on the torsional parameters, which are crucial for the accurate description of conformational equilibria in biomolecules. In particular, attention is paid to the influence of polarization on the (i) transferability of dihedral terms between molecules, (ii) transferability between different environments, and (iii) additivity of dihedral energies. To this end, three polarizable force fields based on the induced point dipole model designed for use in AMBER are tested, including two recent ff02 reparameterizations. Attention is paid to the contributions due to short range interactions (1-2, 1-3, and 1-4) within the four atoms defining the dihedral angle. The results show that when short range 1-2 and 1-3 polarization interactions are omitted, as for instance in ff02, the 1-4 polarization contribution is rather small and unlikely to improve the description of the torsional energy. Conversely, when screened 1-2 and 1-3 interactions are included, the polarization contribution is sizeable and shows potential to improve the transferability of parameters between different molecules and environments as well as the additivity of dihedral terms. However, to reproduce intramolecular polarization effects accurately, further fine-tuning of the short range damping of polarization is necessary. © 2015 Wiley Periodicals, Inc. [ABSTRACT FROM AUTHOR]

The facial selectivity of double Michael addition reactions of the silylated Nazarov reagent 4 to unsaturated indolo[2,3- a]quinolizidine lactams 3 has been studied. Pentacyclic 3-H/15-H trans adducts 5 are generated from Nind-unsubstituted lactams, but the corresponding cis isomers 6 are formed when the indole nitrogen has a tert-butyloxycarbonyl (Boc) substituent. This reversal in the facial selectivity of the annulation has been rationalized by means of theoretical calculations, which indicate that the initial nucleophilic attack under stereoelectronic control is hampered by the presence of the bulky Boc group. The synthetic usefulness of the pentacyclic Nazarov-derived adducts is demonstrated by their conversion into allo and epiallo yohimbine-type targets. [ABSTRACT FROM AUTHOR]

UbcH10 is a component of the Ubiquitin Conjugation Enzymes (Ubc; E2) involved in the ubiquitination cascade controlling the cell cycle progression, whereby ubiquitin, activated by E1, is transferred through E2 to the target protein with the involvement of E3 enzymes. In this work we propose the first three dimensional model of the tetrameric complex formed by the human UbA1 (E1), two ubiquitin molecules and UbcH10 (E2), leading to the transthiolation reaction. The 3D model was built up by using an experimentally guided incremental docking strategy that combined homology modeling, protein-protein docking and refinement by means of molecular dynamics simulations. The structural features of the in silico model allowed us to identify the regions that mediate the recognition between the interacting proteins, revealing the active role of the ubiquitin crosslinked to E1 in the complex formation. Finally, the role of these regions involved in the E1–E2 binding was validated by designing short peptides that specifically interfere with the binding of UbcH10, thus supporting the reliability of the proposed model and representing valuable scaffolds for the design of peptidomimetic compounds that can bind selectively to Ubcs and inhibit the ubiquitylation process in pathological disorders. [ABSTRACT FROM AUTHOR]

Graphite deposits result from the metamorphism of sedimentary rocks rich in carbonaceous matter or from precipitation from carbon-bearing fluids (or melts). The latter process forms vein deposits which are structurally controlled and usually occur in granulites or igneous rocks. The origin of carbon, the mechanisms of transport, and the factors controlling graphite deposition are discussed in relation to their geological settings. Carbon in granulite-hosted graphite veins derives from sublithospheric sources or from decarbonation reactions of carbonate-bearing lithologies, and it is transported mainly in CO-rich fluids from which it can precipitate. Graphite precipitation can occur by cooling, water removal by retrograde hydration reactions, or reduction when the CO-rich fluid passes through relatively low-fO rocks. In igneous settings, carbon is derived from assimilation of crustal materials rich in organic matter, which causes immiscibility and the formation of carbon-rich fluids or melts. Carbon in these igneous-hosted deposits is transported as CO and/or CH and eventually precipitates as graphite by cooling and/or by hydration reactions affecting the host rock. Independently of the geological setting, vein graphite is characterized by its high purity and crystallinity, which are required for applications in advanced technologies. In addition, recent discovery of highly crystalline graphite precipitation from carbon-bearing fluids at moderate temperatures in vein deposits might provide an alternative method for the manufacture of synthetic graphite suitable for these new applications. [ABSTRACT FROM AUTHOR]

We present here an exhaustive characterization of the structure and properties of 6-selenoguanine, an isoster of guanine, and the impact of its introduction in DNA. This study reports the results of state-of-the-art quantum mechanical calculations and atomistic molecular dynamics simulations carried out to shed light on the impact of the replacement of guanine (G) by 6-selenoguanine (SeG) in different forms of DNA. The results point out that the G → SeG substitution leads to stable DNA duplex, antiparallel triplex and G-quadruplex structures, though local distortions are also found. These structural changes affect the thermodynamic stability of the mutation leading to a clear destabilization for all studied systems. Interestingly, the lowest effect has been found when the mutation was placed in the triplex-forming oligonucleotide strand in a reverse Hoogsteen orientation, which favours the antiparallel triplex formation regarding the G-tetraplex formation. Detailed QM studies strongly suggest that SeG impacts the HOMO–LUMO gap and accordingly the transfer properties of DNA, opening the way to modulate the conductivity properties of non-natural DNAs. [ABSTRACT FROM AUTHOR]

The evolution of a ternary molecular system (imine, diene and nitrile) is analyzed to disclose the pathways leading to a divergent synthetic outcome. The Lewis acid catalyzed reaction between cyclohexadiene, 2-phenyl-indol-3-one and acetonitrile yields the imino-Diels-Alder adduct as the major product together with minor amounts of the Mannich-Ritter-amidine product. The experimental and computational data show that the relative orientation of the initial reactants dictates the synthetic outcome. The exo approach between imine and diene leads to the Diels-Alder adduct in a concerted process, whereas the endo mode leads to a polarized intermediate, which is trapped by acetonitrile to yield the multicomponent adduct. [ABSTRACT FROM AUTHOR]

The presence of cavities and tunnels in the interior of proteins, in conjunction with the structural plasticity arising from the coupling to the thermal fluctuations of the protein scaffold, has profound consequences on the pathways followed by ligands moving through the protein matrix. In this perspective we discuss how quantitative analysis of experimental rebinding kinetics from laser flash photolysis, trapping of unstable conformational states by embedding proteins within the nanopores of silica gels, and molecular simulations can synergistically converge to gain insight into the migration mechanism of ligands. We show how the evaluation of the free energy landscape for ligand diffusion based on the outcome of computational techniques can assist the definition of sound reaction schemes, leading to a comprehensive understanding of the broad range of chemical events and time scales that encompass the transport of small ligands in hemeproteins. [ABSTRACT FROM AUTHOR]

This study reports the application of the quantum mechanical self-consistent reaction field MST method to compute the solvation profile in water of the twenty natural amino acids. The aim is to derive intrinsic fractional contributions to the hydration free energy and to examine their transferability to peptides. To this end, IEF-MST calculations have been performed at the B3LYP/6-31G(d) level for the series of acetyl amino acid amides, which were chosen as model compounds. In order to account for the flexibility of both the backbone and the side chain in deriving the hydration fractional contributions, calculations have been performed for representative conformers taken from the Dunbrack's backbone-dependent conformational library. The results allow us to dissect the hydration free energy into backbone and side chain contributions and examine the conformational dependence of the fragmental contributions to hydration. For the backbone, different hydration contributions are found for α-helical and β-sheet conformations, which mainly reflect differences in the electrostatic contribution to hydration of the carbonyl group. In contrast, the conformational flexibility of the side chain is found to have little impact on the fractional contribution to hydration. These findings should be valuable to refine semiempirical methods for predicting solvation properties of peptides and proteins in large-scale genomic studies. [ABSTRACT FROM AUTHOR]

Cytoglobin (Cygb) was recently discovered in the human genome and localized in different tissues. It was suggested to play tissue-specific protective roles, spanning from scavenging of reactive oxygen species in neurons to supplying oxygen to enzymes in fibroblasts. To shed light on the functioning of such versatile machinery, we have studied the processes supporting transport of gaseous heme ligands in Cygb. Carbon monoxide rebinding shows a complex kinetic pattern with several distinct reaction intermediates, reflecting rebinding from temporary docking sites, second order recombination, and formation (and dissociation) of a bis-histidyl heme hexacoordinated reaction intermediate. Ligand exit to the solvent occurs through distinct pathways, some of which exploit temporary docking sites. The remarkable change in energetic barriers, linked to heme bis-histidyl hexacoordination by HisE7, may be responsible for active regulation of the flux of reactants and products to and from the reaction site on the distal side of the heme. A substantial change in both protein dynamics and inner cavities is observed upon transition from the CO-liganded to the pentacoordinated and bis-histidyl hexacoordinated species, which could be exploited as a signalling state. These findings are consistent with the expected versatility of the molecular activity of this protein. [ABSTRACT FROM AUTHOR]

The truncated hemoglobin N, HbN, of Mycobacterium tuberculosis is endowed with a potent nitric oxide dioxygenase (NOD) activity that allows it to relieve nitrosative stress and enhance in vivo survival of its host. Despite its small size, the protein matrix of HbN hosts a two-branched tunnel, consisting of orthogonal short and long channels, that connects the heme active site to the protein surface. A novel dual-path mechanism has been suggested to drive migration of O2 and NO to the distal heme cavity. While oxygen migrates mainly by the short path, a ligand-induced conformational change regulates opening of the long tunnel branch for NO, via a phenylalanine (PheE15) residue that acts as a gate. Site-directed mutagenesis and molecular simulations have been used to examine the gating role played by PheE15 in modulating the NOD function of HbN. Mutants carrying replacement of PheE15 with alanine, isoleucine, tyrosine and tryptophan have similar O2/CO association kinetics, but display significant reduction in their NOD function. Molecular simulations substantiated that mutation at the PheE15 gate confers significant changes in the long tunnel, and therefore may affect the migration of ligands. These results support the pivotal role of PheE15 gate in modulating the diffusion of NO via the long tunnel branch in the oxygenated protein, and hence the NOD function of HbN. [ABSTRACT FROM AUTHOR]

Over the last decades, molecular simulations have spread through the drug discovery arena. This trend is expected to continue in the foreseeable future thanks to increased performance and the positive impact they can exert on productivity. In this article we highlight three aspects of molecular modelling for which we expect significant improvements over the next 25 years. Increased computational resources, faster algorithms and novel methods to sample rare events will provide a better handle on target flexibility and its relation with ligand binding. More accurate target druggability predictions will improve the success, but also broaden the scope of target-based drug discovery strategies. Finally, the use of higher levels of theory will increase the accuracy of protein-ligand binding affinity predictions, resulting in better hit identification success rates as well as more efficient lead optimization processes. [ABSTRACT FROM AUTHOR]

Motivation: A variety of pocket detection algorithms are now freely or commercially available to the scientific community for the analysis of static protein structures. However, since proteins are dynamic entities, enhancing the capabilities of these programs for the straightforward detection and characterization of cavities taking into account protein conformational ensembles should be valuable for capturing the plasticity of pockets, and therefore allow gaining insight into structure–function relationships.Results: This article describes a new method, called MDpocket, providing a fast, free and open-source tool for tracking small molecule binding sites and gas migration pathways on molecular dynamics (MDs) trajectories or other conformational ensembles. MDpocket is based on the fpocket cavity detection algorithm and a valuable contribution to existing analysis tools. The capabilities of MDpocket are illustrated for three relevant cases: (i) the detection of transient subpockets using an ensemble of crystal structures of HSP90; (ii) the detection of known xenon binding sites and migration pathways in myoglobin; and (iii) the identification of suitable pockets for molecular docking in P38 Map kinase.Availability: MDpocket is free and open-source software and can be downloaded at http://fpocket.sourceforge.net.Contact: pschmidtke@ub.eduSupplementary Information: Supplementary data are available at Bioinformatics online. [ABSTRACT FROM AUTHOR]