Velocity cross-correlations and atomic momentum transfer in simple liquids with different potential cores

Time correlation functions between the velocity of a tagged particle and velocities of particles within specified ranges of initial separations have been obtained by molecular dynamics simulation. These correlation functions have allowed us to analyze the momentum transfer between particles in different coordination shells. Two simple liquids at very different densities and two purely repulsive potentials with very different softnesses have been considered. The longitudinal correlations, which are the velocity cross-correlations along the initial direction defined by the centers of two given particles, have been calculated separately. It has been proven that these correlations should be attributed to particles both in front of and behind the central one. As with propagating longitudinal modes, they are strongly dependent on the softness of the potential core. Some characteristic features of the velocity correlation functions after the initial rise should be related to nonlongitudinal correlations. It has been shown that velocity cross-correlations between distinct particles cannot only be attributed to the direct interactions among particles, but also to the motions induced by the movement of a tagged particle on their neighbors.

Diffusive transport properties in monovalent and divalent metal-ion halide melts: A computer simulation study

Self- and cross-velocity correlation functions and related transport coefficients of molten salts are studied by molecular-dynamics simulation. Six representative systems are considered, i.e., NaCl and KCl alkali halides, CuCl and CuBr noble-metal halides, and SrCl2 and ZnCl2 divalent metal-ion halides. Computer simulation results are compared with experimental self-diffusion coefficients and electrical conductivities. Special attention is paid to dynamic cross correlations and their dependence on the Coulomb interactions as well as on the size and mass differences between anions and cations.

Cooperative kinking at distant sites in mechanically stressed DNA.

In cells, DNA is routinely subjected to significant levels of bending and twisting. In some cases, such as under physiological levels of supercoiling, DNA can be so highly strained, that it transitions into non-canonical structural conformations that are capable of relieving mechanical stress within the template. DNA minicircles offer a robust model system to study stress-induced DNA structures. Using DNA minicircles on the order of 100 bp in size, we have been able to control the bending and torsional stresses within a looped DNA construct. Through a combination of cryo-EM image reconstructions, Bal31 sensitivity assays and Brownian dynamics simulations, we have been able to analyze the effects of biologically relevant underwinding-induced kinks in DNA on the overall shape of DNA minicircles. Our results indicate that strongly underwound DNA minicircles, which mimic the physical behavior of small regulatory DNA loops, minimize their free energy by undergoing sequential, cooperative kinking at two sites that are located about 180° apart along the periphery of the minicircle. This novel form of structural cooperativity in DNA demonstrates that bending strain can localize hyperflexible kinks within the DNA template, which in turn reduces the energetic cost to tightly loop DNA.

Somatostatin Subtype‑2 Receptor-Targeted Metal-Based Anticancer Complexes

Conjugates of a dicarba analogue of octreotide, a potent somatostatin agonist whose receptors are overexpressed on tumor cells, with [PtCl2(dap)] (dap = 1-(carboxylic acid)-1,2-diaminoethane) (3), [(η6-bip)Os(4-CO2-pico)Cl] (bip = biphenyl, pico = picolinate) (4), [(η6-p-cym)RuCl(dap)]+ (p-cym = p-cymene) (5), and [(η6-p-cym)RuCl(imidazole-CO2H)(PPh3)]+ (6), were synthesized by using a solid-phase approach. Conjugates 35 readily underwent hydrolysis and DNA binding, whereas conjugate 6 was inert to ligand substitution. NMR spectroscopy and molecular dynamics calculations showed that conjugate formation does not perturb the overall peptide structure. Only 6 exhibited antiproliferative activity in human tumor cells (IC50 = 63 ± 2 μM in MCF-7 cells and IC50 = 26 ± 3 μM in DU-145 cells) with active participation of somatostatin receptors in cellular uptake. Similar cytotoxic activity was found in a normal cell line (IC50 = 45 ± 2.6 μM in CHO cells), which can be attributed to a similar level of expression of somatostatin subtype-2 receptor. These studies provide new insights into the effect of receptor-binding peptide conjugation on the activity of metal-based anticancer drugs, and demonstrate the potential of such hybrid compounds to target tumor cells specifically.

Matching stages of heavy-ion collision models

Heavy-ion reactions and other collective dynamical processes are frequently described by different theoretical approaches for the different stages of the process, like initial equilibration stage, intermediate locally equilibrated fluid dynamical stage, and final freeze-out stage. For the last stage, the best known is the Cooper-Frye description used to generate the phase space distribution of emitted, noninteracting particles from a fluid dynamical expansion or explosion, assuming a final ideal gas distribution, or (less frequently) an out-of-equilibrium distribution. In this work we do not want to replace the Cooper-Frye description, but rather clarify the ways of using it and how to choose the parameters of the distribution and, eventually, how to choose the form of the phase space distribution used in the Cooper-Frye formula. Moreover, the Cooper-Frye formula is used in connection with the freeze-out problem, while the discussion of transition between different stages of the collision is applicable to other transitions also. More recently, hadronization and molecular dynamics models have been matched to the end of a fluid dynamical stage to describe hadronization and freeze-out. The stages of the model description can be matched to each other on space-time hypersurfaces (just like through the frequently used freeze-out hypersurface). This work presents a generalized description of how to match the stages of the description of a reaction to each other, extending the methodology used at freeze-out, in simple covariant form which is easily applicable in its simplest version for most applications.

Influence of the pseudopotential on the properties of liquid rubidium at 315 K

The influence of the pseudopotential on both the structure and the self-diffusion of liquid rubidium at the melting point has been investigated by means of molecular-dynamics calculations. The model potential considered has been computed from the pseudopotential of Ashcroft, the dielectric function of Geldart and Vosko, and a Born-Mayer term. Four different values for the core radius which enters as input in the pseudopotential have been considered. In this way we have been able to observe and interpret the effect of this contribution on the properties of the liquid.

Oxygen adsorption on clean and oxygen precovered Cu(100)

Kuparipinnan hapettumisen alkuvaiheet ovat vielä nykyisin tutkijoille epäselviä. Kuitenkin, jotta hapettumisprosessia voitaisiin säädellä, on sangen tärkeää ymmärtää mistä varsinainen hapettuminen lähtee liikkeelle ja mitkä ovat hapettumisen seuraavat vaiheet. Tähän kysymykseen haetaan vastauksia tässä työssä käyttäen puhtaasti teoreettisia menetelmiä pinnan käsittelyssä. Aikaisempien teoreettisten ja kokeellisten tutkimusten välillä on pieni ristiriita liittyen hapen tarttumistodennäköisyyteen. Teoreettisten tutkimusten mukaan happi ei puhtaalle pinnalle tullessaan näe potentiaalivallia, mutta kokeelliset tutkimukset osoittavat sellaisen kuitenkin olevan. Tuohon ristiriitaan pureudutaan käyttäen aikaisemmista laskuista poikkeavaa kvanttimekaaniseen molekyylidynamiikkaan perustuvaa lähestymistapaa. Työssä havaitaan, että aikaisemmin yleisesti käytetty menetelmä hukkaa huomattavan määrän tietoa ja siten tutkijat eivät voi ainoastaan tyytyä tarkastelemaan kyseisellä menetelmällä saatuja tuloksia. Kuparipinnalle havaittiin, että korkeilla molekyylin kineettisen energian arvolla aikaisemmin suoritetut laskut hajottavista trajektoreista pitävät paikkansa, mutta matalilla kineettisen energian arvoilla molekyyli kohtaa erittäin voimakkaan ``steering'' vaikutuksen ja trajektorit joiden piti olla hajottavia johtavatkin molekulaariseen adsorptioon. Kun hapen konsentraatio pinnalla on suurempi kuin 0.5 ML, pinta rekonstruoituu. Myös rekonstruktion jälkeistä pintaa on tutkittu samanlaisilla menetelmillä kuin puhdasta pintaa. Rekonstruoituneelle pinnalle ei löydetty hajottavia trajektoreita ja havaittiin, että hapelle annetun kineettisen energian matalilla arvoilla myös tässä tapauksessa on erittäin voimakas ``steering'' vaikutus.

The structure of human tRNALys3 anticodon bound to HIV genome is stabilized by modified nucleosides and adjacent mismatch base pairs

Replication of human immunodeficiency virus (HIV) requires base pairing of the reverse transcriptase primer, human tRNA(Lys3), to the viral RNA. Although the major complementary base pairing occurs between the HIV primer binding sequence (PBS) and the tRNA's 3'-terminus, an important discriminatory, secondary contact occurs between the viral A-rich Loop I, 5'-adjacent to the PBS, and the modified, U-rich anticodon domain of tRNA(Lys3). The importance of individual and combined anticodon modifications to the tRNA/HIV-1 Loop I RNA's interaction was determined. The thermal stabilities of variously modified tRNA anticodon region sequences bound to the Loop I of viral sub(sero)types G and B were analyzed and the structure of one duplex containing two modified nucleosides was determined using NMR spectroscopy and restrained molecular dynamics. The modifications 2-thiouridine, s(2)U(34), and pseudouridine, Psi(39), appreciably stabilized the interaction of the anticodon region with the viral subtype G and B RNAs. The structure of the duplex results in two coaxially stacked A-form RNA stems separated by two mismatched base pairs, U(162)*Psi(39) and G(163)*A(38), that maintained a reasonable A-form helix diameter. The tRNA's s(2)U(34) stabilized the interaction between the A-rich HIV Loop I sequence and the U-rich anticodon, whereas the tRNA's Psi(39) stabilized the adjacent mismatched pairs.

Structure-Based, Rational Design of T Cell Receptors.

Adoptive cell transfer using engineered T cells is emerging as a promising treatment for metastatic melanoma. Such an approach allows one to introduce T cell receptor (TCR) modifications that, while maintaining the specificity for the targeted antigen, can enhance the binding and kinetic parameters for the interaction with peptides (p) bound to major histocompatibility complexes (MHC). Using the well-characterized 2C TCR/SIYR/H-2K(b) structure as a model system, we demonstrated that a binding free energy decomposition based on the MM-GBSA approach provides a detailed and reliable description of the TCR/pMHC interactions at the structural and thermodynamic levels. Starting from this result, we developed a new structure-based approach, to rationally design new TCR sequences, and applied it to the BC1 TCR targeting the HLA-A2 restricted NY-ESO-1157-165 cancer-testis epitope. Fifty-four percent of the designed sequence replacements exhibited improved pMHC binding as compared to the native TCR, with up to 150-fold increase in affinity, while preserving specificity. Genetically engineered CD8(+) T cells expressing these modified TCRs showed an improved functional activity compared to those expressing BC1 TCR. We measured maximum levels of activities for TCRs within the upper limit of natural affinity, K D = ∼1 - 5 μM. Beyond the affinity threshold at K D < 1 μM we observed an attenuation in cellular function, in line with the "half-life" model of T cell activation. Our computer-aided protein-engineering approach requires the 3D-structure of the TCR-pMHC complex of interest, which can be obtained from X-ray crystallography. We have also developed a homology modeling-based approach, TCRep 3D, to obtain accurate structural models of any TCR-pMHC complexes when experimental data is not available. Since the accuracy of the models depends on the prediction of the TCR orientation over pMHC, we have complemented the approach with a simplified rigid method to predict this orientation and successfully assessed it using all non-redundant TCR-pMHC crystal structures available. These methods potentially extend the use of our TCR engineering method to entire TCR repertoires for which no X-ray structure is available. We have also performed a steered molecular dynamics study of the unbinding of the TCR-pMHC complex to get a better understanding of how TCRs interact with pMHCs. This entire rational TCR design pipeline is now being used to produce rationally optimized TCRs for adoptive cell therapies of stage IV melanoma.

Técnicas para análise de desempenho de computadores

A series of open source benchmarks for computer performance analysis of personal computers with a focus on computational chemistry calculations is presented. The results returned by these tests are discussed and used to correlate with the actual performance of a set of computers available for research on two computing intensive fields of chemistry, quantum chemical and molecular simulation calculations.

Técnicas de análise do perfil de execução e otimização de programas em química computacional

In this paper we review the basic techniques of performance analysis within the UNIX environment that are relevant in computational chemistry, with particular emphasis on the execution profile using the gprof tool. Two case studies (in ab initio and molecular dynamics calculations) are presented in order to illustrate how execution profiling can be used to effectively identify bottlenecks and to guide source code optimization. Using these profiling and optimization techniques it was possible to obtain significant speedups (of up to 30%) in both cases.

Tetrahydrobenzo[h][1,6]naphthyridine-6-chlorotacrine hybrids as a new family of anti-Alzheimer agents targeting beta-amyloid, tau, and cholinesterase pathologies

Optimization of an essentially inactive 3,4-dihydro-2H-pyrano[3,2-c]quinoline carboxylic ester derivative as acetylcholinesterase (AChE) peripheral anionic site (PAS)-binding motif by double O → NH bioisosteric replacement, combined with molecular hybridization with the AChE catalytic anionic site (CAS) inhibitor 6-chlorotacrine and molecular dynamics-driven optimization of the length of the linker has resulted in the development of the trimethylene-linked 1,2,3,4-tetrahydrobenzo[h][1,6]naphthyridine6-chlorotacrine hybrid 5a as a picomolar inhibitor of human AChE (hAChE). The tetra-, penta-, and octamethylene-linked homologues 5bd have been also synthesized for comparison purposes, and found to retain the nanomolar hAChE inhibitory potency of the parent 6-chlorotacrine. Further biological profiling of hybrids 5ad has shown that they are also potent inhibitors of human butyrylcholinesterase and moderately potent Aβ42 and tau anti-aggregating agents, with IC50 values in the submicromolar and low micromolar range, respectively. Also, in vitro studies using an artificial membrane model have predicted a good brain permeability for hybrids 5ad, and hence, their ability to reach their targets in the central nervous system. The multitarget profile of the novel hybrids makes them promising leads for developing anti-Alzheimer drug candidates with more balanced biological activities.

Tetrahydrobenzo[h][1,6]naphthyridine-6-chlorotacrine hybrids as a new family of anti-Alzheimer agents targeting beta-amyloid, tau, and cholinesterase pathologies

Optimization of an essentially inactive 3,4-dihydro-2H-pyrano[3,2-c]quinoline carboxylic ester derivative as acetylcholinesterase (AChE) peripheral anionic site (PAS)-binding motif by double O → NH bioisosteric replacement, combined with molecular hybridization with the AChE catalytic anionic site (CAS) inhibitor 6-chlorotacrine and molecular dynamics-driven optimization of the length of the linker has resulted in the development of the trimethylene-linked 1,2,3,4-tetrahydrobenzo[h][1,6]naphthyridine6-chlorotacrine hybrid 5a as a picomolar inhibitor of human AChE (hAChE). The tetra-, penta-, and octamethylene-linked homologues 5bd have been also synthesized for comparison purposes, and found to retain the nanomolar hAChE inhibitory potency of the parent 6-chlorotacrine. Further biological profiling of hybrids 5ad has shown that they are also potent inhibitors of human butyrylcholinesterase and moderately potent Aβ42 and tau anti-aggregating agents, with IC50 values in the submicromolar and low micromolar range, respectively. Also, in vitro studies using an artificial membrane model have predicted a good brain permeability for hybrids 5ad, and hence, their ability to reach their targets in the central nervous system. The multitarget profile of the novel hybrids makes them promising leads for developing anti-Alzheimer drug candidates with more balanced biological activities.

Tetrahydrobenzo[h][1,6]naphthyridine-6-chlorotacrine hybrids as a new family of anti-Alzheimer agents targeting beta-amyloid, tau, and cholinesterase pathologies

Optimization of an essentially inactive 3,4-dihydro-2H-pyrano[3,2-c]quinoline carboxylic ester derivative as acetylcholinesterase (AChE) peripheral anionic site (PAS)-binding motif by double O → NH bioisosteric replacement, combined with molecular hybridization with the AChE catalytic anionic site (CAS) inhibitor 6-chlorotacrine and molecular dynamics-driven optimization of the length of the linker has resulted in the development of the trimethylene-linked 1,2,3,4-tetrahydrobenzo[h][1,6]naphthyridine6-chlorotacrine hybrid 5a as a picomolar inhibitor of human AChE (hAChE). The tetra-, penta-, and octamethylene-linked homologues 5bd have been also synthesized for comparison purposes, and found to retain the nanomolar hAChE inhibitory potency of the parent 6-chlorotacrine. Further biological profiling of hybrids 5ad has shown that they are also potent inhibitors of human butyrylcholinesterase and moderately potent Aβ42 and tau anti-aggregating agents, with IC50 values in the submicromolar and low micromolar range, respectively. Also, in vitro studies using an artificial membrane model have predicted a good brain permeability for hybrids 5ad, and hence, their ability to reach their targets in the central nervous system. The multitarget profile of the novel hybrids makes them promising leads for developing anti-Alzheimer drug candidates with more balanced biological activities.

Tetrahydrobenzo[h][1,6]naphthyridine-6-chlorotacrine hybrids as a new family of anti-Alzheimer agents targeting beta-amyloid, tau, and cholinesterase pathologies

Optimization of an essentially inactive 3,4-dihydro-2H-pyrano[3,2-c]quinoline carboxylic ester derivative as acetylcholinesterase (AChE) peripheral anionic site (PAS)-binding motif by double O → NH bioisosteric replacement, combined with molecular hybridization with the AChE catalytic anionic site (CAS) inhibitor 6-chlorotacrine and molecular dynamics-driven optimization of the length of the linker has resulted in the development of the trimethylene-linked 1,2,3,4-tetrahydrobenzo[h][1,6]naphthyridine6-chlorotacrine hybrid 5a as a picomolar inhibitor of human AChE (hAChE). The tetra-, penta-, and octamethylene-linked homologues 5bd have been also synthesized for comparison purposes, and found to retain the nanomolar hAChE inhibitory potency of the parent 6-chlorotacrine. Further biological profiling of hybrids 5ad has shown that they are also potent inhibitors of human butyrylcholinesterase and moderately potent Aβ42 and tau anti-aggregating agents, with IC50 values in the submicromolar and low micromolar range, respectively. Also, in vitro studies using an artificial membrane model have predicted a good brain permeability for hybrids 5ad, and hence, their ability to reach their targets in the central nervous system. The multitarget profile of the novel hybrids makes them promising leads for developing anti-Alzheimer drug candidates with more balanced biological activities.