The Effects of Sequence Variation on Genome-wide NRF2 Binding--New Target Genes and Regulatory SNPs.

Transcription factor binding specificity is crucial for proper target gene regulation. Motif discovery algorithms identify the main features of the binding patterns, but the accuracy on the lower affinity sites is often poor. Nuclear factor E2-related factor 2 (NRF2) is a ubiquitous redox-activated transcription factor having a key protective role against endogenous and exogenous oxidant and electrophile stress. Herein, we decipher the effects of sequence variation on the DNA binding sequence of NRF2, in order to identify both genome-wide binding sites for NRF2 and disease-associated regulatory SNPs (rSNPs) with drastic effects on NRF2 binding. Interactions between NRF2 and DNA were studied using molecular modelling, and NRF2 chromatin immunoprecipitation-sequence datasets together with protein binding microarray measurements were utilized to study binding sequence variation in detail. The binding model thus generated was used to identify genome-wide binding sites for NRF2, and genomic binding sites with rSNPs that have strong effects on NRF2 binding and reside on active regulatory elements in human cells. As a proof of concept, miR-126-3p and -5p were identified as NRF2 target microRNAs, and a rSNP (rs113067944) residing on NRF2 target gene (Ferritin, light polypeptide, FTL) promoter was experimentally verified to decrease NRF2 binding and result in decreased transcriptional activity.

Evaluation of the antiviral efficacy of bis[1,2]dithiolo[1,4]thiazines and bis[1,2]dithiolopyrrole derivatives against the nucelocapsid protein of the Feline Immunodeficiency Virus (FIV) as a model for HIV infection.

A diverse library of bis[1,2]dithiolo[1,4]thiazines and bis[1,2]dithiolopyrrole derivatives were prepared for evaluation of activity against the nucleocapsid protein of the Feline Immunodeficiency Virus (FIV) as a model for HIV, using an in vitro cell culture approach, yielding nanomolar active compounds with low toxicity.

Combining NMR ensembles and molecular dynamics simulations provides more realistic models of protein structures in solution and leads to better chemical shift prediction.

While chemical shifts are invaluable for obtaining structural information from proteins, they also offer one of the rare ways to obtain information about protein dynamics. A necessary tool in transforming chemical shifts into structural and dynamic information is chemical shift prediction. In our previous work we developed a method for 4D prediction of protein (1)H chemical shifts in which molecular motions, the 4th dimension, were modeled using molecular dynamics (MD) simulations. Although the approach clearly improved the prediction, the X-ray structures and single NMR conformers used in the model cannot be considered fully realistic models of protein in solution. In this work, NMR ensembles (NMRE) were used to expand the conformational space of proteins (e.g. side chains, flexible loops, termini), followed by MD simulations for each conformer to map the local fluctuations. Compared with the non-dynamic model, the NMRE+MD model gave 6-17% lower root-mean-square (RMS) errors for different backbone nuclei. The improved prediction indicates that NMR ensembles with MD simulations can be used to obtain a more realistic picture of protein structures in solutions and moreover underlines the importance of short and long time-scale dynamics for the prediction. The RMS errors of the NMRE+MD model were 0.24, 0.43, 0.98, 1.03, 1.16 and 2.39 ppm for (1)Hα, (1)HN, (13)Cα, (13)Cß, (13)CO and backbone (15)N chemical shifts, respectively. The model is implemented in the prediction program 4DSPOT, available at http://www.uef.fi/4dspot.

Overexpression of SUMO perturbs the growth and development of Caenorhabditis elegans.

Small ubiquitin-related modifiers (SUMOs) are important regulator proteins. Caenorhabditis elegans contains a single SUMO ortholog, SMO-1, necessary for the reproduction of C. elegans. In this study, we constructed transgenic C. elegans strains expressing human SUMO-1 under the control of pan-neuronal (aex-3) or pan-muscular (myo-4) promoter and SUMO-2 under the control of myo-4 promoter. Interestingly, muscular overexpression of SUMO-1 or -2 resulted in morphological changes of the posterior part of the nematode. Movement, reproduction and aging of C. elegans were perturbed by the overexpression of SUMO-1 or -2. Genome-wide expression analyses revealed that several genes encoding components of SUMOylation pathway and ubiquitin-proteasome system were upregulated in SUMO-overexpressing nematodes. Since muscular overexpression of SMO-1 also brought up reproductive and mobility perturbations, our results imply that the phenotypes were largely due to an excess of SUMO, suggesting that a tight control of SUMO levels is important for the normal development of multicellular organisms.

Molecular mechanism of allosteric communication in the human PPARalpha-RXRalpha heterodimer.

The peroxisome proliferator-activated receptor alpha (PPARalpha) is a nuclear receptor (NR) that forms a heterodimeric transcription factor complex with the retinoid X receptor alpha (RXRalpha). The phenomenon that the heterodimer can be activated by both PPARalpha and RXRalpha ligands, while both ligands have a synergistic effect on its activity suggests that there is an allosteric communication within the heterodimer. In this study, the molecular mechanism of this allosteric signaling was studied by molecular dynamics (MD) simulations and some of the residues involved in this communication were tested experimentally. Multiple MD simulations were done for the PPARalpha-RXRalpha heterodimer ligand-binding domains (LBDs) without ligands, with agonistic ligand bound to RXRalpha or PPARalpha, and ligand bound to both receptors. Fluctuation calculations and structural clustering analysis of the heterodimer MD simulations showed that ligand binding to RXRalpha decreases fluctuations of large parts of PPARalpha, most notably helices 3 and 4 at the coactivator binding site, which presumably stabilizes the coactivator binding to heterodimer complex. The dynamics of helix 8-9 loop and helix 10/11 located at the heterodimeric interface were affected by RXRalpha ligand binding, suggesting that these parts of the dimer are involved in allosteric communication. Experimental data complemented this view by showing that a large set of residues at the heterodimerization surface has a role in the communication. These results provided evidence that RXRalpha ligand binding-induced stabilization of PPARalpha coactivator binding site has a role in the permissive and synergistic activation of the PPARalpha-RXRalpha heterodimer. Proteins 2010. (c) 2009 Wiley-Liss, Inc.

Identification of novel peptide inhibitors for human trypsins.

Human trypsin isoenzymes share extensive sequence similarity, but certain differences in their activity and susceptibility to inhibitors have been observed. Using phage display technology, we identified seven different peptides that bind to and inhibit the activity of trypsin-3, a minor trypsin isoform expressed in pancreas and brain. All of the peptides contain at least two of the amino acids tryptophan, alanine and arginine, whereas proline was found closer to the N-terminus in all but one peptide. All peptides contain two or more cysteines, suggesting a cyclic structure. However, we were able to make synthetic linear variants of these peptides without losing bioactivity. Alanine replacement experiments for one of the peptides suggest that the IPXXWFR motif is important for activity. By molecular modeling the same amino acids were found to interact with trypsin-3. The peptides also inhibit trypsin-1, but only weakly, if at all, trypsin-2 and -C. As trypsin is a highly active enzyme which can activate protease-activated receptors and enzymes that participate in proteolytic cascades involved in tumor invasion and metastasis, these peptides might be useful lead molecules for the development of drugs for diseases associated with increased trypsin activity.

Synthesis and biological evaluation of phenolic 4,5-dihydroisoxazoles and 3-hydroxy ketones as estrogen receptor alpha and beta agonists.

In this work, 52 diphenyl-4,5-dihydroisoxazoles and -3-hydroxy ketones were prepared and their estrogen receptor alpha (ERalpha) and estrogen receptor beta (ERbeta) activities were explored in order to systematize and maximize their biological activity. The biological activity was firstly screened by using ERE reporter assay to find out how aromatic hydroxylation and methylation of the chiral centers of the compounds affect the ability of ER to mediate biological responses. For selected 19 compounds, the relative binding affinities (RBA, relative to 3,17beta-estradiol) and ability to induce transcription of primary E2 target gene pS2 in human MCF-7 breast cancer cells were determined. In the reporter assay, many compounds showed even stronger activity than E2 and some of them showed RBA larger than 1%. The highest RBAs were determined for the enantiomers of 1-hydroxy-6-(4-hydroxy-phenyl)-1-phenyl-hexan-3-one (50a and 50b). Isomer 50a showed high binding affinity both to ERalpha (with RBA approximately 200%) and ERbeta (with RBA approximately 60%), while the RBAs of 50b were ca. 40% of those. Some of the other compounds (with RBA approximately 1-16%) showed also notable ERalpha binding selectivity. When four most promising ligands (50a, 50b, 45a, and 45b) were studied with respect to their ability to induce the transcription of primary E2 target gene pS2, the compounds acted as agonists or partial agonists. Computer modeling was used to predict receptor binding conformations and to rationalize the RBA differences of the compounds.

4D prediction of protein (1)H chemical shifts.

A 4D approach for protein (1)H chemical shift prediction was explored. The 4th dimension is the molecular flexibility, mapped using molecular dynamics simulations. The chemical shifts were predicted with a principal component model based on atom coordinates from a database of 40 protein structures. When compared to the corresponding non-dynamic (3D) model, the 4th dimension improved prediction by 6-7%. The prediction method achieved RMS errors of 0.29 and 0.50 ppm for Halpha and HN shifts, respectively. However, for individual proteins the RMS errors were 0.17-0.34 and 0.34-0.65 ppm for the Halpha and HN shifts, respectively. X-ray structures gave better predictions than the corresponding NMR structures, indicating that chemical shifts contain invaluable information about local structures. The (1)H chemical shift prediction tool 4DSPOT is available from http://www.uku.fi/kemia/4dspot .

Ligand unbinding pathways from the vitamin D receptor studied by molecular dynamics simulations.

Molecular dynamics simulation techniques have been used to study the unbinding pathways of 1alpha,25-dihydroxyvitamin D(3) from the ligand-binding pocket of the vitamin D receptor (VDR). The pathways observed in a large number of relatively short (< 200 ps) random acceleration molecular dynamics (RAMD) trajectories were found to be in fair agreement, both in terms of pathway locations and deduced relative preferences, compared to targeted molecular dynamics (TMD) and streered molecular dynamics simulations (SMD). However, the high-velocity ligand expulsions of RAMD tend to favor straight expulsion trajectories and the observed relative frequencies of different pathways were biased towards the probability of entering a particular exit channel. Simulations indicated that for VDR the unbinding pathway between the H1-H2 loop and the beta-sheet between H5 and H6 is more favorable than the pathway located between the H1-H2 loop and H3. The latter pathway has been suggested to be the most likely unbinding path for thyroid hormone receptors (TRs) and a likely path for retinoic acid receptor. Ligand entry/exit through these two pathways would not require displacement of H12 from its agonistic position. Differences in the packing of the H1, H2, H3 and beta-sheet region explain the changed relative preference of the two unbinding pathways in VDR and TRs. Based on the crystal structures of the ligand binding domains of class 2 nuclear receptors, whose members are VDR and TRs, this receptor class can be divided in two groups according to the packing of the H1, H2, H3 and beta-sheet region.

Synthesis and evaluation of estrogen agonism of diaryl 4,5-dihydroisoxazoles, 3-hydroxyketones, 3-methoxyketones, and 1,3-diketones: a compound set forming a 4D molecular library.

In this paper, the preparation and systematic evaluation of estrogen receptor alpha (ER alpha) and estrogen receptor beta (ER beta) activities of some diaryl-1,3-diones and their synthetic intermediates, diaryl-4,5-dihydroisoxazoles, diaryl-3-hydroxyketones, diaryl-3-methoxyketones, and diaryl-2-(dimethyl-lambda 4-sulfanylidene)-1,3-diones, is described. The set of 72 compounds constitutes a general schematic structure aryl1-linker1-spacer-linker2-aryl2, where the linker1-spacer-linker2 length varies between 4 and 8 carbons. The set of compounds was applied here to map and explore the active sites of subtypes ER alpha and ER beta. The highest activities were obtained with dihydroisoxazole and hydroxyketone spacers, but even the most flexible diones with unsubstituted aryl groups showed some agonism. Most compounds were found to be ER alpha selective or to activate both receptors, but in some cases we saw also clearly stronger ER beta activation.

Gene regulatory potential of nonsteroidal vitamin D receptor ligands.

The seco-steroid 1alpha,25-dihydroxyvitamin D3 [1alpha,25(OH)2D3] is a promising drug candidate due to its pleiotropic function including the regulation of calcium homeostasis, bone mineralization and cellular proliferation, differentiation, and apoptosis. We report here a novel class of nonsteroidal compounds, represented by the bis-aromatic molecules CD4409, CD4420, and CD4528, as ligands of the 1alpha,25(OH)2D3 receptor (VDR). Taking the known diphenylmethane derivative LG190178 as a reference, this study provides molecular evaluation of the interaction of nonsteroidal ligands with the VDR. All four nonsteroidal compounds were shown to induce VDR-retinoid X receptor heterodimer complex formation on a 1alpha,25(OH)2D3 response element, stabilize the agonistic conformation of the VDR ligand-binding domain, enable the interaction of VDR with coactivator proteins and contact with their three hydroxyl groups the same residues within the ligand-binding pocket of the VDR as 1alpha,25(OH)2D3. Molecular dynamics simulations demonstrated that all four nonsteroidal ligands take a shape within the ligand-binding pocket of the VDR that is very similar to that of the natural ligand. CD4528 is mimicking the natural hormone best and was found to be in vitro at least five times more potent than LG190178. In living cells, CD4528 was only two times less potent than 1alpha,25(OH)2D3 and induced mRNA expression of the VDR target gene CYP24 in a comparable fashion. At a noncalcemic dose of 150 microg/kg, CD4528 showed in vivo a clear induction of CYP24 expression and therefore may be used as a lead compound for the development of therapeutics against psoriasis, osteoporosis, and cancer.

Antagonist- and inverse agonist-driven interactions of the vitamin D receptor and the constitutive androstane receptor with corepressor protein.

Ligand-dependent signal transduction by nuclear receptors (NRs) includes dynamic exchanges of coactivator (CoA) and corepressor (CoR) proteins. Here we focused on the structural determinants of the antagonist- and inverse agonist-enhanced interaction of the endocrine NR vitamin D receptor (VDR) and the adopted orphan NR constitutive androstane receptor (CAR) from two species with the CoR NR corepressor. We found that the pure VDR antagonist ZK168281 and the human CAR inverse agonist clotrimazole are both effective inhibitors of the CoA interaction of their respective receptors, whereas ZK168281 resembled more the mouse CAR inverse agonist androstanol in its ability to recruit CoR proteins. Molecular dynamics simulations resulted in comparable models for the CoR receptor interaction domain peptide bound to VDR/antagonist or CAR/inverse agonist complexes. A salt bridge between the CoR and a conserved lysine in helix 4 of the NR is central to this interaction, but also helix 12 was stabilized by direct contacts with residues of the CoR. Fixation of helix 12 in the antagonistic/inverse agonistic conformation prevents an energetically unfavorable free floatation of the C terminus. The comparable molecular mechanisms that explain the similar functional profile of antagonist and inverse agonists are likely to be extended from VDR and CAR to other members of the NR superfamily and may lead to the design of even more effective ligands.

Critical role of helix 12 of the vitamin D(3) receptor for the partial agonism of carboxylic ester antagonists.

The carboxy-terminal alpha-helix of a nuclear receptor ligand-binding domain (LBD), helix 12, contains a critical, ligand-modulated interface for the interaction with coactivator proteins. In this study, using the example of the vitamin D receptor (VDR) and the partial antagonist ZK159222, the role of helix 12 (residues 417-427) for both antagonistic and agonistic receptor actions was investigated. Amino acid residue G423 was demonstrated to be critical for partial agonism of ZK159222, but not for the activity of the natural VDR agonist, 1alpha,25-dihydroxyvitamin D(3) (1alpha,25(OH)(2)D(3)). The amount of partial agonism of ZK159222 increased when helix 12 was truncated by the last four amino acid residues (Delta424-27) and augmented even more, when in addition helix 12 of VDR's dimerization partner, retinoid X receptor (RXR), was truncated. In contrast, the low agonism of a structural derivative of ZK159222, ZK168281, was not affected comparably, whereas other close structural relatives of ZK159222 even demonstrated the same agonistic activity as that of 1alpha,25(OH)(2)D(3). The amount of agonism of ZK159222 and ZK168281 at different variations of helix 12 correlated well with VDR's ability to complex with coactivator proteins and inversely correlated with the strength of the compound's antagonistic action on 1alpha,25(OH)(2)D(3) signalling. Molecular dynamics simulations of the LBD complexed with the two antagonists could explain their different action by demonstrating a more drastic displacement of helix 12 through ZK168281 than through ZK159222. Moreover, the modelling could indicate a kink of helix 12 at amino acid residue G423, which provides the last four amino acid residues of helix 12 with a modulatory role for the partial agonism of some VDR antagonists, such as ZK159222. In conclusion, partial agonism of a VDR antagonist is lower the more it disturbs helix 12 in taking the optimal position for coactivator interaction.

A structural basis for the species-specific antagonism of 26,23-lactones on vitamin D signaling.

The 26,23-lactone derivative of 1alpha,25-dihydroxyvitamin D3, TEI-9647, is a partial antagonist of the of human vitamin D receptor (VDR). However, we found that TEI-9647 in rat cells behaves as a weak VDR agonist. This behavior could be mimicked in human cells by the double mutagenesis of human VDR (specifically C403S and C410N). The increased agonistic action of TEI-9647 correlates to a gain in the interaction of the VDR with coactivator protein and a decreased stabilization of the antagonistic conformation of the receptor. Molecular dynamics simulations indicated that TEI-9647 acts as antagonist of human VDR by reducing the stability of helix 12 of the ligand binding domain. In contrast, N410 of the rat VDR stabilized, via backbone contacts, the interaction between helices 11 and 12. This results in TEI-9647 becoming a weak agonist in this organism.

Birch PR-10c interacts with several biologically important ligands.

PR-10c is a unique member of PR-10 proteins in birch, since it is the only one known to be post-translationally modified by glutathione and is not constitutively expressed in pollen. Both reduced and S-glutathiolated forms of PR-10c show low ribonuclease activity. However, the major function of the protein is apparently not yet resolved. Our protein-ligand interaction studies with saturation transfer difference (STD) NMR revealed that PR-10c interacts with several biologically important molecules, including cytokinin, flavonoid glycosides, sterols and emodin. Competition study with deoxycholate and kinetin revealed no statistically significant binding interference, indicating that these ligands have different binding sites in PR-10c. Ligand docking studies with a molecular model of PR-10c support the STD NMR results of ligand binding and binding epitopes, suggesting that there are three potential binding sites in PR-10c: two in the hydrophobic cavity and one in the glycine-rich loop. Our docking calculations suggested that only kinetin interacts with the glycine-rich loop, the binding occurring through its adenine moiety. Clear ligand specificity could be observed in the binding of nucleotide derivatives. S-glutathiolation of PR-10c did not affect kinetin binding. The present results suggest that birch PR-10c is a multifunctional protein, which has diverse roles in plant stress responses.

Corepressor excess shifts the two-side chain vitamin D analog Gemini from an agonist to an inverse agonist of the vitamin D receptor.

The vitamin D receptor (VDR) is an endocrine nuclear receptor that binds with high affinity its natural ligand 1alpha,25-dihydroxyvitamin D3. Gemini is a 1alpha,25-dihydroxyvitamin D3 analog with two identical side chains that, despite its significantly increased volume, binds to the VDR and can function as a potent agonist. This study demonstrates that, at excess corepressor (CoR) levels, Gemini shifts from an agonist to an inverse agonist that actively recruits CoR proteins to the VDR and mediates superrepression. Under these conditions Gemini stabilizes the VDR into a silent conformation, in which helix 12 of the ligand-binding domain is repositioned and thus unable to contribute to coactivator interaction. Amino acid F422 has been described as the lock of helix 12 and seems to be the most critical VDR residue in the inverse agonistic action of Gemini. Molecular dynamics simulations of the Gemini-VDR complex support these observation by indicating that the second side chain of Gemini induces tension to the receptor structure that can be released by a shift of helix 12. Taken together, Gemini is the first described (conditional) inverse agonist to an endocrine nuclear receptor and may function as a sensor for the cell-specific coactivator/CoR ratio.

Analysis of the binding energies of testosterone, 5alpha-dihydrotestosterone, androstenedione and dehydroepiandrosterone sulfate with an antitestosterone antibody.

Molecular dynamics simulations and molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) free energy calculations were used to study the binding of testosterone (TES), 5alpha-dihydrotestosterone (5ADHT), androstenedione (AND), and dehydroepiandrosterone sulfate (DHEAS) to the monoclonal antitestosterone antibody 3-C(4)F(5). The relative binding free energy of TES and AND was also calculated with free energy perturbation (FEP) simulations. The antibody 3-C(4)F(5) has a relatively high affinity (3 x 10(8) M(-1)) and on overall good binding profile for testosterone but its cross-reactivity with DHEAS has been the main reason for the failure to use this antibody in clinical immunoassays. The relative binding free energies obtained with the MM-PBSA method were 1.5 kcal/mol for 5ADHT, 3.8 kcal/mol for AND, and 4.3 kcal/mol for DHEAS, as compared to TES. When a water molecule of the ligand binding site, observed in the antibody-TES crystal structure, was explicitly included in MM-PBSA calculations, the relative binding energies were 3.4, 4.9, and 5.4 kcal/mol for 5ADHT, AND, and DHEAS, respectively. The calculated numbers are in correct order but larger than the corresponding experimental energies of 1.3, 1.5, and 2.6 kcal/mol, respectively. The fact that the MM-PBSA method reproduced the relative binding free energies of DHEAS, a steroid having a negatively charged sulfate group, and the neutrally charged TES, 5ADHT, and AND in satisfactory agreement with experiment shows the robustness of the method in predicting relative binding affinities. The 800-ps FEP simulations predicted that the antibody 3-C(4)F(5) binds TES 1.3 kcal/mol tighter than AND. Computational mutagenesis of selected amino acid residues of the ligand binding site revealed that the lower affinities of AND and DHEAS as compared to TES are due to a combined effect of several residues, each contributing a small fraction to the tighter binding of TES. An exception to this is Tyr99H, whose mutation to Ala lowered the binding of DHEAS 0.7 kcal/mol more than the binding of TES. This is probably due to the hydrogen bonding interaction formed between the OH group of Tyr99H and the sulfate group of DHEAS. Computational mutagensis data also showed that the affinity of the steroids to the antitestosterone antibody 3-C(4)F(5) would be enhanced if Trp47H were repositioned so that it would make more extensive contacts with the bound ligands. In addition, the binding of steroids to antitestosterone, antiprogesterone, and antiestradiol antibodies is discussed.

Free energy simulations and MM-PBSA analyses on the affinity and specificity of steroid binding to antiestradiol antibody.

Antiestradiol antibody 57-2 binds 17beta-estradiol (E2) with moderately high affinity (K(a) = 5 x 10(8) M(-1)). The structurally related natural estrogens estrone and estriol as well synthetic 17-deoxy-estradiol and 17alpha-estradiol are bound to the antibody with 3.7-4.9 kcal mol(-1) lower binding free energies than E2. Free energy perturbation (FEP) simulations and the molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) method were applied to investigate the factors responsible for the relatively low cross-reactivity of the antibody with these four steroids, differing from E2 by the substituents of the steroid D-ring. In addition, computational alanine scanning of the binding site residues was carried out with the MM-PBSA method. Both the FEP and MM-PBSA methods reproduced the experimental relative affinities of the five steroids in good agreement with experiment. On the basis of FEP simulations, the number of hydrogen bonds formed between the antibody and steroids, which varied from 0 to 3 in the steroids studied, determined directly the magnitude of the steroid-antibody interaction free energies. One hydrogen bond was calculated to contribute about 3 kcal mol(-1) to the interaction energy. Because the relative binding free energies of estrone (two antibody-steroid hydrogen bonds), estriol (three hydrogen bonds), 17-deoxy-estradiol (no hydrogen bonds), and 17alpha-estradiol (two hydrogen bonds) are close to each other and clearly lower than that of E2 (three hydrogen bonds), the water-steroid interactions lost upon binding to the antibody make an important contribution to the binding free energies. The MM-PBSA calculations showed that the binding of steroids to the antiestradiol antibody is driven by van der Waals interactions, whereas specificity is solely due to electrostatic interactions. In addition, binding of steroids to the antiestradiol antibody 57-2 was compared to the binding to the antiprogesterone antibody DB3 and antitestosterone antibody 3-C4F5, studied earlier with the MM-PBSA method.

Ab Initio Quantum Mechanical Study on the Origin of the pK(a) Differences of the Proton Sponges 1,8-Bis(dimethylamino)naphthalene, 1,8-Bis(dimethylamino)-2,7-dimethoxynaphthalene, 1,6-Dimethyl-1,6-diazacyclodecane, and 1,6-Diazabicyclo[4.4.4]tetradecane.

Ab initio quantum mechanical calculations were used in studying the origin of the exceptionally high basicities of four diamines (13-16) with pK(a1) values ranging from 12.1 to 25. The computational approach involved the calculation of the gas-phase proton affinities of the molecules studied at the MP2/6-31G//HF/6-31G level and the solvation energies with the polarizable continuum model at the HF/6-31G level. The calculated gas-phase and aqueous-phase proton affinities of a structurally diverse series of amines were compared with the corresponding experimental gas-phase proton affinities and pK(a1) values. The calculated values were found to be in reasonable agreement with the experimental ones. The high basicities of the studied diamines were found to originate from the nitrogen lone-pair repulsion, solvation effects, and strong intramolecular hydrogen bonds. Each of these factors were found to be able to increase the pK(a1) values of the high-basicity diamines by 2-6 pK(a) units. The relative contributions of the factors varied between the compounds. The nitrogen lone-pair repulsion was estimated to be the most important factor in increasing the pK(a1) values. In addition, barriers for proton transfers between the nitrogens of selected diamines were calculated, and comparison was made between the barrier heights and the geometries of the diamines.

Correlative motions and memory effects in molecular dynamics simulations of molecules: principal components and rescaled range analysis suggest that the motions of native BPTI are more correlated than those of its mutants.

In this work MD simulations of the native bovine pancreatic trypsin inhibitor (BPTI) and 16 mutants were done in vacuum in order to study memory effects in the mutants using principal component analysis (PCA) and the rescaled range analysis (Hurst exponents). Both PCA and the rescaled range analysis support our previous proposition, based on PCA of lysozyme, that the motions of a native protein are more correlated than those of mutants. The methods are compared, the nature and applications of the rule and the role of the long-range correlations in MD time series (i.e. memory) are discussed in the context of collective motions.