Search

Multitopic supramolecular guests with finely tuned affinities toward widely explored cucurbit[n]urils (CBs) and cyclodextrins (CDs) have been recently designed and tested as functional components of advanced supramolecular systems. We employed various spacers between the adamantane cage and a cationic moiety as a tool for tuning the binding strength toward CB7 to prepare a set of model guests with KCB7 and Kβ-CD values of (0.6-5.0) × 1010 M-1 and (0.6-2.6) × 106 M-1, respectively. These accessible adamantylphenyl-based binding motifs open a way toward supramolecular components with an outstanding affinity toward β-cyclodextrin. 1H NMR experiments performed in 30% CaCl2/D2O at 273 K along with molecular dynamics simulations allowed us to identify two arrangements of the guest@β-CD complexes. The approach, joining experimental and theoretical methods, provided a better understanding of the structure of cyclodextrin complexes and related molecular recognition, which is highly important for the rational design of drug delivery systems, molecular sensors and switches.

Mismatch repair is a highly conserved cellular pathway responsible for repairing mismatched dsDNA. Errors are detected by the MutS enzyme, which most likely senses altered mechanical property of damaged dsDNA rather than a specific molecular pattern. While the curved shape of dsDNA in crystallographic MutS/DNA structures suggests the role of DNA bending, the theoretical support is not fully convincing. Here, we present a computational study focused on a base-pair opening into the minor groove, a specific base-pair motion observed upon interaction with MutS. Propensities for the opening were evaluated in terms of two base-pair parameters: Opening and Shear. We tested all possible base pairs in anti/anti, anti/syn and syn/anti orientations and found clear discrimination between mismatches and canonical base-pairs only for the opening into the minor groove. Besides, the discrimination gap was also confirmed in hotspot and coldspot sequences, indicating that the opening could play a more significant role in the mismatch recognition than previously recognized. Our findings can be helpful for a better understanding of sequence-dependent mutability. Further, detailed structural characterization of mismatches can serve for designing anti-cancer drugs targeting mismatched base pairs.

A novel application of strip and pixel silicon radiation detectors for study and characterization of run-away electron events in tokamaks is presented. Main goal was to monitor runaway electrons both directly and indirectly. The strip detector was placed inside the tokamak vacuum chamber in order to monitor the run-away electrons directly. Whereas the pixel detector was placed outside the tokamak chamber behind a pin hole for monitoring the run-away electrons indirectly via radiation produce by interaction of the electrons with the plasma facing material. Results obtained using the silicon detectors are compared with already existing diagnostic methods consisting of scintillation devices detecting X-rays and photo-neutrons, providing the same results in the observable comparisons. Tests with the pixel detector proved that the pinhole camera is able to extract spatial information of interaction point (a place where the runaway electrons hit on the facing material) and the strip detectors indicate presence of additional signal from throughout the discharge. The performed experiments are innovative, illustrating possible development of new and easy to use diagnostic method.

Stilbene derivatives are well-recognised substructures of molecular switches based on photochemically and/or thermally induced (E)/(Z) isomerisation. We combined a stilbene motif with two benzimidazolium arms to prepare new sorts of supramolecular building blocks and examined their binding properties towards cucurbit[n]urils (n=7, 8) and cyclodextrins (β-CD, γ-CD) in water. Based on the 1 H NMR data and molecular dynamics simulations, we found that two distinct complexes with different stoichiometry, i. e., guest@β-CD and guest@β-CD2 , coexist in equilibrium in a water solution of the (Z)-stilbene-based guests. We also demonstrated that the bis(benzimidazolio)stilbene guests can be transformed from the (E) into the (Z) form via UV irradiation and back via thermal treatment in DMSO.

We report the synthesis of the conformationally mobile S-shaped glycoluril pentamer building block 3a and two new acyclic CB[n]-type receptors P1 and P2. P1 (9 mM) and P2 (11 mM) have moderate aqueous solubility but their host•guest complexes are poorly soluble. Host P1 does not undergo intermolecular self-association whereas P2 does (K(s) = 189±27 M(−1)). (1)H NMR titrations show that P1 and P2 are poor hosts toward hydrophobic (di)cations 6 – 11 (P1: K(a) = 375–1400 M(−1); P2: K(a) = 1950–19800 M(−1)) compared to Tet1 and Tet2 (Tet1: K(a) = 3.09 × 10(6) to 4.69 × 10(8) M(−1); Tet2: K(a) = 4.59 × 10(8) to 1.30 × 10(10) M(−1)). Molecular modelling shows that P1 and P2 exist as a mixture of three different conformers due to the two S-shaped methylene bridged glycoluril dimer subunits that each possess two different conformations. The lowest energy conformers of P1 and P2 do not feature a well-defined central cavity. In the presence of guests, P2 adapts its conformation to form 1:1 P2•guest complexes; the binding free energy pays the energetic price of conformer selection. This energetically unfavorable conformer selection results in significantly decreased K(a) values of P1 and P2 compared to Tet1 and Tet2.

This paper presents two scenarios used for generation of a runaway electron (RE) beam in the COMPASS tokamak with a focus on the decay phase and control of the beam. The first scenario consists of massive gas injection of argon into the current ramp-up phase, leading to a disruption accompanied by runaway plateau generation. In the second scenario, injection of a smaller amount of gas is used in order to isolate the RE beam from high-temperature plasma. The performances of current control and radial and vertical position feedback control in the second scenario were experimentally studied and analysed. The role of RE energy in the radial position stability of the RE beam seems to be crucial. A comparison of the decay phase of the RE beam in various amounts of Ar or Ne was studied using absolute extreme ultraviolet (AXUV) tomography and hard x-ray (HXR) intensity measurement. Argon clearly leads to higher HXR fluxes for the same current decay rate than neon, while radiated power based on AXUV measurements is larger for Ne in the same set of discharges.

We studied sequence-dependent retention properties of synthetic 5′-terminal phosphate absent trinucleotides containing adenine, guanine and thymine through reversed-phase liquid chromatography (RPLC) and QSRR modelling. We investigated the influence of separation conditions, namely mobile phase composition (ion interaction agent content, pH and organic constituent content), on sequence-dependent separation by means of ion-interaction RPLC (II-RPLC) using two types of models: experimental design–artificial neural networks (ED-ANN), and linear regression based on molecular dynamics data. The aim was to determine those properties of the above-mentioned analytes responsible for the retention dependence of the sequence. Our results show that there is a deterministic relation between sequence and II-RPLC retention properties of the studied trinucleotides. Further, we can conclude that the higher the content of ion-interaction agent in the mobile phase, the more prominent these properties are. We also show that if we approximate the polar component of solvation energy in QSRR by the electrostatic work in transferring molecules from vacuum to water, and the non-polar component by the solvent accessible surface area, these parameters best describe the retention properties of trinucleotides. There are some exceptions to this finding, namely sequences 5′-NAN-3′, 5′-ANN-3′, 5′-TGN-3′, 5′-NTA-3′and 5′-NGA-3′ (N stands for generic nucleotide). Their role is still unknown, but since linear regression including these specific constellations showed a higher observable variance coverage than the model with only the basic descriptors, we may assume that solvent-analyte interactions are responsible for the exceptional behaviour of 5′-NAN-3′ & 5′-ANN-3′ trinucleotides and some intramolecular interactions of neighbouring nucleobases for 5′-TGN-3′, 5′-NTA-3′and 5′-NGA-3′ trinucleotides.

Chitin is one of the most abundant structural biomolecules in nature, where it occurs in the form of nanofibrils that are the smallest building blocks for many biological structural materials, such as the exoskeleton of Arthropoda. Despite this fact, little is known about the structural properties of these nanofibrils. Here, we present a theoretical study of a single chitin molecule and 10 α-chitin nanofibrils with different numbers of chains in an aqueous environment that mimics the conditions in natural systems during self-assembly. Our extensive analysis of the molecular dynamics trajectories, including free energy calculations for every model system, reveals not only the structural properties of the nanofibrils, but also provides insight into the principles of nanofibril formation. We identified the fundamental phenomena occurring in the chitin nanofibrils such as their hydrogen bonding pattern and resulting helical shape. With increasing size, the nanofibrils become increasingly stable and their structural properties approach those of crystalline α-chitin if they consist of more than 20 chains. Interestingly, this coincides with the typical size of chitin nanofibrils observed in natural systems, suggesting that their evolutionary success was at least partially driven by these specific structure–property relations.

A novel macrocycle, decamethylpressocucurbit[5]uril (Me10prCB[5]), was synthesized by acid-catalyzed condensation of propanediurea and paraformaldehyde. This macrocycle binds methane with higher affinity than cucurbit[5]uril and its permethylated derivative.

BsoBI is a type II restriction endonuclease belonging to the EcoRI family. There is only one previously published X-ray structure for this endonuclease: it shows a homodimer of BsoBI completely encircling DNA in a tunnel. In this work, molecular dynamics simulations were employed to elucidate possible ways in which DNA is loaded into this complex prior to its cleavage. We found that the dimer does not open spontaneously when DNA is removed from the complex on the timescale of our simulations (~ 0.5 μs). A biased simulation had to be used to facilitate the opening, which revealed a possible way for the two catalytic domains to separate. The α-helices connecting the catalytic and helical domains were found to act as a hinge during the separation. In addition, we found that the opening of the BsoBI dimer was influenced by the type of counterions present in the environment. A reference simulation of the BsoBI/DNA complex further showed spontaneous reorganization of the active sites due to the binding of solvent ions, which led to an active-site structure consistent with other experimental structures of type II restriction endonucleases determined in the presence of metal ions.

Methyl viologen hexafluorophosphate (MV2+·2PF6–) and dodecamethylbambus[6]uril (BU6) form crystals in which the layers of viologen dications alternate with those of a 1:2 supramolecular complex of BU6 and PF6–. This arrangement allows for a one-electron reduction of MV2+ ions upon UV irradiation to form MV+• radical cations within the crystal structure with half-lives of several hours in air. The mechanism of this photoinduced electron transfer in the solid state and the origin of the long-lived charge-separated state were studied by steady-state and transient spectroscopies, cyclic voltammetry, and electron paramagnetic resonance spectroscopy. Our experiments are supported by quantum-chemical calculations showing that BU6 acts as a reductant. In addition, analogous photochemical behavior is also demonstrated on other MV2+/BU6 crystals containing either BF4– or Br– counterions.

Guanine DNA quadruplexes are interesting and important biological objects because they represent potential targets for regulatory drugs. Their use as building blocks for biomaterial applications is also being investigated. This contribution reports the in silico design of artificial building blocks derived from xanthine. Methods of quantum chemistry were used to evaluate the properties of xanthine structures relative to those containing guanine, the natural reference used. Tailoring the xanthine core showed that the base stacking and the ion coordination were significantly enhanced in the designed systems, while the ion-transport properties were not affected. Our study suggests that the 9-deaza-8-haloxanthine bases (where the halogen is fluorine or chlorine) are highly promising candidates for the development of artificial quadruplexes and quadruplex-active ligands.

Imidazolium-based guests containing two distinct binding epitopes are capable of binding β-cyclodextrin and cucurbit[6/7]uril (CB) simultaneously to form heteroternary 1:1:1 inclusion complexes. In the final configuration, the hosts occupy binding sites disfavored in the binary complexes because of the chemically induced reorganization of the intermediate 1:1 aggregate. In addition, the reported guests are capable of binding two CBs to form either 1:2 or 1:1:1 ternary assemblies despite consisting of a single cationic moiety. Whereas the adamantane site binds CB solely via hydrophobic interactions, the CB unit at the butyl site is stabilized by a combination of hydrophobic and ion–dipole interactions.

Cucurbiturils are the most potent artificial receptors known for many organic molecules in water. However, little is known about their supramolecular chemistry in organic solvents. Here we present a new cucurbituril derivative, 1, and investigate its supramolecular properties in methanol. The macrocycle resembles a five-membered cucurbituril in which four glycoluril units are replaced with propanediurea. Macrocycle 1 can bind to one cation such as potassium or anilinium via each of its opposed portals. The stability of these complexes in methanol at nanomolar concentrations exceeds that of complexes between metal cations and crown ethers. Moreover, macrocycle 1 forms a self-assembled tetrameric aggregate in the solid state and in methanol. The tetramer is stabilized by the addition of up to 1 equiv of a cation but is fully disassembled in the presence of 2 equiv of the cation. Cations can thus be used to tune the aggregation of 1 in solution.

The isotropic (129)Xe nuclear magnetic resonance (NMR) chemical shift (CS) in Xe@C60 dissolved in liquid benzene was calculated by piecewise approximation to faithfully simulate the experimental conditions and to evaluate the role of different physical factors influencing the (129)Xe NMR CS. The (129)Xe shielding constant was obtained by averaging the (129)Xe nuclear magnetic shieldings calculated for snapshots obtained from the molecular dynamics trajectory of the Xe@C60 system embedded in a periodic box of benzene molecules. Relativistic corrections were added at the Breit-Pauli perturbation theory (BPPT) level, included the solvent, and were dynamically averaged. It is demonstrated that the contribution of internal dynamics of the Xe@C60 system represents about 8% of the total nonrelativistic NMR CS, whereas the effects of dynamical solvent add another 8%. The dynamically averaged relativistic effects contribute by 9% to the total calculated (129)Xe NMR CS. The final theoretical value of 172.7 ppm corresponds well to the experimental (129)Xe CS of 179.2 ppm and lies within the estimated errors of the model. The presented computational protocol serves as a prototype for calculations of (129)Xe NMR parameters in different Xe atom guest-host systems.

Adamantylated bisimidazolium cations exhibit a distinct fragmentation pathway in contrast to their cucurbit[7]uril (CB7) complexes. The observed alternative fragmentation of the guest molecule in a complex clearly correlates to the supposed sterically hindered or allowed slippage of the macrocycle over the axel molecule.

Folded RNA molecules are shaped by an astonishing variety of highly conserved noncanonical molecular interactions and backbone topologies. The dinucleotide platform is a widespread recurrent RNA modular building submotif formed by the side-by-side pairing of bases from two consecutive nucleotides within a single strand, with highly specific sequence preferences. This unique arrangement of bases is cemented by an intricate network of noncanonical hydrogen bonds and facilitated by a distinctive backbone topology. The present study investigates the gas-phase intrinsic stabilities of the three most common RNA dinucleotide platforms — 5′-GpU-3′, ApA, and UpC — via state-of-the-art quantum-chemical (QM) techniques. The mean stability of base-base interactions decreases with sequence in the order GpU > ApA > UpC. Bader’s atoms-in-molecules analysis reveals that the N2(G)…O4(U) hydrogen bond of the GpU platform is stronger than the corresponding hydrogen bonds in the other two platforms. The mixed-pucker sugar-phosphate backbone conformation found in most GpU platforms, in which the 5′-ribose sugar (G) is in the C2′-endo form and the 3′-sugar (U) in the C3′-endo form, is intrinsically more stable than the standard A-RNA backbone arrangement, partially as a result of a favorable O2′…O2P intra-platform interaction. Our results thus validate the hypothesis of Lu et al. (Lu Xiang-Jun, et al. Nucleic Acids Res. 2010, 38, 4868-4876), that the superior stability of GpU platforms is partially mediated by the strong O2′…O2P hydrogen bond. In contrast, ApA and especially UpC platform-compatible backbone conformations are rather diverse and do not display any characteristic structural features. The average stabilities of ApA and UpC derived backbone conformers are also lower than those of GpU platforms. Thus, the observed structural and evolutionary patterns of the dinucleotide platforms can be accounted for, to a large extent, by their intrinsic properties as described by modern QM calculations. In contrast, we show that the dinucleotide platform is not properly described in the course of atomistic explicit-solvent simulations. Our work also gives methodological insights into QM calculations of experimental RNA backbone geometries. Such calculations are inherently complicated by rather large data and refinement uncertainties in the available RNA experimental structures, which often preclude reliable energy computations.

The isotropic 129Xe NMR chemical shift of atomic Xe dissolved in liquid benzene was simulated by combining classical molecular dynamics and quantum chemical calculations of 129Xe nuclear magnetic shielding. Snapshots from the molecular dynamics trajectory of xenon atom in a periodic box of benzene molecules were used for the quantum chemical calculations of isotropic 129Xe chemical shift using nonrelativistic density functional theory as well as relativistic Breit–Pauli perturbation corrections. Thus, the correlation and relativistic effects as well as the temperature and dynamics effects could be included in the calculations. Theoretical results are in a very good agreement with the experimental data. The most of the experimentally observed isotropic 129Xe shift was recovered in the nonrelativistic dynamical region, while the relativistic effects explain of about 8% of the total 129Xe chemical shift.

In this study, we have investigated the supramolecular interaction between series of 1-alkyl-3-methylimidazolium guests with variable alkyl substituent lengths and cucurbit[6]uril (CB6) in the solution and the solid state. Correct interpretation of (1)H NMR spectra was a key issue for determining the binding modes of the complexes in solution. Unusual chemical shifts of some protons in the (1)H NMR spectra were explained by the polarization of the imidazolium aromatic ring upon the complexation with the host. The formation of 1:1 complex between 1-ethyl-3-methylimidazolium and CB6 is in disagreement with previously reported findings describing an inclusion of two guest molecules in the CB6 cavity.

We present the results of a series of 10-ns molecular dynamics simulations on Pseudomonas aeruginosa lectin-II (PA-IIL) and its complexes with four different monosaccharides. We compare the saccharide-free, saccharideoccupied, and saccharide- and ion-free forms of the lectin. The results are coupled with analysis of the water density map and calcium coordination. The water density pattern around the binding site in the free lectin molecular dynamics was fitted with that in the X-ray and with the hydroxyl groups of the monosaccharide within the lectin/monosaccharide complexes and the best ligand was predicted based on the best fit. Interestingly, the water density pattern around the binding site in the uncomplexed lectin exactly fitted the O2, O3, and O4 hydroxyl groups of the fucose complex with the lectin. This observation could lead to a hypothesis that the replacement of these three water molecules from the binding site by the monosaccharide decreases the entropy of the complex and increases the entropy of the water molecules, which favors the binding. It suggests that the high density peaks of the solvent around the binding site in the free protein could be the tool to predict hydroxyl group orientation of the sugar in the protein/sugar complexes. The high affinity of PA-IIL binding site is also attributed to the presence of two calcium ions, each of them making five to six coordinations with the protein part and two coordinations with either water or the monosaccharide. When the calcium ions are removed from the simulated system, they are replaced by sodium ions from the solvent. These observations rationalize the high binding affinity of PA-IIL towards fucose.

Possible rearrangement mechanisms of hydrogen-bond arrays formed at the lower rim of tetrahydroxycalix[4]arene and tetrahydroxythiacalix[4]arene were studied by means of density functional theory and the resolution identity approximation modification of Møller-Plesset perturbation theory (RI-MP2). Influence of solvent to height of energy barriers was quantified by use of the conductorlike screening model (COSMO) of implicit solvent (chloroform). Generally, two types of mechanisms were investigated. The first is represented by a synchronous single-step jump of all four hydroxyl protons. Pathways of the second mechanism include the rotation of one or more hydroxyl groups around the CAr-O bond. Theoretical results, in agreement with recently published experimental data (Lang et al. J. Chem. Phys. 2005, 122, 044056), prefer a jump mechanism for the methylene-bridged calix[4]arene. Concerning the thiacalix[4]arene, results obtained by COSMO as well as RI-MP2 calculations show that the rotational mechanism is very competitive and it could even be more favorable.

Molecular dynamics simulations and adaptive biasing force analysis of the quadruplex DNA dynamics in an explicit solvent reveal fundamentally different mechanisms of Na(+) transport in xanthine- and guanine-based DNA systems. The barrier to the transport of K(+) through the xanthine-based quadruplex is significantly lower than those reported for the guanine-based analogs.

Molecular dynamics simulations have been performed on fully solvated alpha-(1-->4)-galactosyltransferase LgtC from Neisseria meningitidis with and without the donor substrate UDP-Gal and in the presence of the manganese ion. The analysis of the trajectories revealed a limited movement in the loop X (residues 75-80) and a larger conformational change in the loop Y (residues 246-251) in the simulation, when UDP-Gal was not present. In this case, the loops X and Y open by almost 10A, exposing the active site to the solvent. The 'hinge region' responsible for the opening is composed of residues 246-247. We have also analyzed the behavior of the manganese ion in the simulations. The coordination number is 6 when UDP-Gal is present and it increases to 7 when it is absent. In the latter case, three water molecules become coordinated to the ion. In both cases, the coordination is very stable implying that the manganese ion is tightly bound in the active site of the enzyme even if UDP-Gal is not present. Further analysis of the structural water molecules location confirmed that the mobility of water molecules in the active site and the accessibility of this site for solvent are higher in the absence of the substrate.

Two different reaction mechanisms of a mixed criss-cross cycloaddition with opposite sequence of reaction steps, intra-intermolecular and inter-intramolecular, were explored by quantum-chemical calculations at the MP2 and B3LYP levels of theory and with cc-pVDZ basis set. It was found that the rate-determining step in the both mechanisms is the first step regardless of the mechanism. The Gibbs activation barrier of the intramolecular step of the intra-intermolecular sequence is by 10.4 kcal mol-1 lower than that of the intermolecular step of the inter-intramolecular sequence at the MP2 level of theory and almost the same at the B3LYP level of theory. This together with known experimental data confirms that the intra-intermolecular sequence is the most probable mechanism.

A study of proton transfer in models of a single peptide unit (N-methylacetamide) and diamide (N2-acetyl-N1-methylglycinamide) as well as the influence of a single water molecule on proton transfer is presented here. Three proton pathways in protonated N-methylacetamide are considered: isomerization, inversion, and 1,3-proton shift. The isomerization step exhibits the lowest energy barrier. When a single water molecule was added, no significant influence on proton isomerization was observed. In the diamide model, the isomerization−jump mechanism of proton transfer along diamide carbonyl oxygens was inspected, and the proton isomerization steps were found to be the most energy-demanding processes (∼17 kcal mol-1). The presence of a single water molecule leads to a different, lower-energy-barrier proton-transfer mechanism with proton exchange. The highest energy barrier is only 7.6 kcal mol-1. Possible competing pathways are also discussed.

Mutations in human genes can be responsible for inherited genetic disorders and cancer. Mutations can arise due to environmental factors or spontaneously. It has been shown that certain DNA sequences are more prone to mutate. These sites are termed hotspots and exhibit a higher mutation frequency than expected by chance. In contrast, DNA sequences with lower mutation frequencies than expected by chance are termed coldspots. Mutation hotspots are usually derived from a mutation spectrum, which reflects particular population where an effect of a common ancestor plays a role. To detect coldspots/hotspots unaffected by population bias, we analysed the presence of germline mutations obtained from HGMD database in the 5-nucleotide segments repeatedly occurring in genes associated with common inherited disorders, in particular, the PAH, LDLR, CFTR, F8, and F9 genes. Statistically significant sequences (mutational motifs) rarely associated with mutations (coldspots) and frequently associated with mutations (hotspots) exhibited characteristic sequence patterns, e.g. coldspots contained purine tract while hotspots showed alternating purine-pyrimidine bases, often with the presence of CpG dinucleotide. Using molecular dynamics simulations and free energy calculations, we analysed the global bending properties of two selected coldspots and two hotspots with a G/T mismatch. We observed that the coldspots were inherently more flexible than the hotspots. We assume that this property might be critical for effective mismatch repair as DNA with a mutation recognized by MutSα protein is noticeably bent.

Thermally induced intra-intermolecular criss-cross cycloaddition of nonsymmetrical azines in the presence of phenyl isocyanate is described. This one-pot reaction gives heterocyclic compounds containing three fused five-membered rings. The molecular and crystal structure of the products is reported and discussed.

Vectorial proton transfer among carbonyl oxygen atoms was studied in two models of tripeptide via quantum chemical calculations using the hybrid B3LYP functional and the 6-31++G** basis set. Two principal proton transfer pathways were found: a first path involving isomerization of the proton around the double bond of the carbonyl group, and a second based on the large conformational flexibility of the tripeptide model where all carbonyl oxygen atoms cooperate. The latter pathway has a rate-determining step energy barrier that is only around half of that for the first pathway. As conformational flexibility plays a crucial role in second pathway, the effect of attaching methyl groups to the alpha carbon atoms was studied. The results obtained are presented for all four possible stereochemical configurations.

Free energy calculations tightly bind together experimental observations and computer-aided design through corresponding simulations at atomic level. Nowadays, free energy calculations represent a cornerstone for obtaining deeper insight into any biomolecular system, its molecular structure and especially dynamic behaviour. Despite this fact, several challenges have to be considered before routine utilization of this type of calculations. At first one has to deal properly with so-called sampling problem. Secondly, the appropriate way how to treat their extreme time demands has to be applied for obtaining converged and reliable data comparable with the experimental ones. Here we present the experience obtained through design, implementation and subsequent application of Multiple Walkers Approach (MWA) connected with Adaptive Biasing Force (ABF) and corresponding advantages and added values of MWA ABF run within distributed, heterogeneous grid environment under the EUAsia virtual organization (VO).

We have prepared organic guest molecules in which two pyridinium rings are connected through an aromatic/aliphatic bridge bearing a carboxyl group. The supramolecular interactions between these guests and macrocyclic hosts cucurbit[7]uril (CB7) and cucurbit[8]uril (CB8) has been studied. We have demonstrated that the binding modes of the complexes depend on the type of central bridge present in the guest molecules and the size of the macrocycle. We have also showed that the binding mode between cucurbiturils and guests with aromatic bridges is pH independent. On the other hand, a guest containing an aliphatic bridge and CB7 formed a pseudorotaxane, which behaved as a pH-driven molecular switch.

The virtual organizations form a key concept for seamless utilization of all advanced features of available worldwide grid environments, the EGEE Grid in particular. However, apart from large, well organized user communities there is a substantial, non-trivial activation barrier in adopting grid concept and creating a new VO for small research groups. To minimize this barrier, a so-called “catch-all” virtual organization approach has been implemented. Here we present two successful use cases of catch-all virtual organizations (Virtual Organization for Central Europe and for EUAsia), demonstrating the strengths of this approach, discussing their setups, encountered challenges and suggesting corresponding solutions.

Functional RNA molecules such as ribosomal RNAs (rRNAs) frequently contain highly conserved internal loops with a 5'-UAA/5'-GAN (UAA/GAN) consensus sequence. The UAA/GAN internal loops adopt a distinctive structure inconsistent with secondary structure predictions. The structure has a narrow major groove and forms a trans Hoogsteen/Sugar edge (tHS) A/G base pair followed by an unpaired stacked adenine, a trans Watson-Crick/Hoogsteen (tWH) U/A base pair, and finally a bulged nucleotide (N). The structure is further stabilized by a three-adenine stack and base-phosphate interaction. In the ribosome, the UAA/GAN internal loops are involved in extensive tertiary contacts, mainly as donors of A-minor interactions. Further, this sequence can adopt an alternative 2D/3D pattern stabilized by a four-adenine stack involved in a smaller number of tertiary interactions. The solution structure of an isolated UAA/GAA internal loop shows substantially rearranged base pairing with three consecutive non-Watson-Crick base pairs. Its A/U base pair adopts an incomplete cis Watson-Crick/Sugar edge (cWS) A/U conformation instead of the expected Watson-Crick arrangement. We performed 3.1 μs of explicit solvent molecular dynamics (MD) simulations of the X-ray and NMR UAA/GAN structures, supplemented by molecular mechanics, Poisson-Boltzmann, and surface area free energy calculations; locally enhanced sampling (LES) runs; targeted MD (TMD); and nudged elastic band (NEB) analysis. We compared parm99 and parmbsc0 force fields and net-neutralizing Na(+) versus excess salt KCl ion environments. Both force fields provide a similar description of the simulated structures, with the parmbsc0 leading to modest narrowing of the major groove. The excess salt simulations also cause a similar effect. While the NMR structure is entirely stable in simulations, the simulated X-ray structure shows considerable widening of the major groove, a loss of base-phosphate interaction, and other instabilities. The alternative X-ray geometry even undergoes a conformational transition toward the solution 2D structure. Free energy calculations confirm that the X-ray arrangement is less stable than the solution structure. LES, TMD, and NEB provide a rather consistent pathway for interconversion between the X-ray and NMR structures. In simulations, the incomplete cWS A/U base pair of the NMR structure is water-mediated and alternates with the canonical A-U base pair, which is not indicated by the NMR data. Completion of the full cWS A/U base pair is prevented by the overall internal loop arrangement. In summary, the simulations confirm that the UAA/GAN internal loop is a molecular switch RNA module that adopts its functional geometry upon specific tertiary contexts.

Eight mutants of the DhaA haloalkane dehalogenase carrying mutations at the residues lining two tunnels, previously observed by protein X-ray crystallography, were constructed and biochemically characterized. The mutants showed distinct catalytic efficiencies with the halogenated substrate 1,2,3-trichloropropane. Release pathways for the two dehalogenation products, 2,3-dichloropropane-1-ol and the chloride ion, and exchange pathways for water molecules, were studied using classical and random acceleration molecular dynamics simulations. Five different pathways, denoted p1, p2a, p2b, p2c, and p3, were identified. The individual pathways showed differing selectivity for the products: the chloride ion releases solely through p1, whereas the alcohol releases through all five pathways. Water molecules play a crucial role for release of both products by breakage of their hydrogen-bonding interactions with the active-site residues and shielding the charged chloride ion during its passage through a hydrophobic tunnel. Exchange of the chloride ions, the alcohol product, and the waters between the buried active site and the bulk solvent can be realized by three different mechanisms: (i) passage through a permanent tunnel, (ii) passage through a transient tunnel, and (iii) migration through a protein matrix. We demonstrate that the accessibility of the pathways and the mechanisms of ligand exchange were modified by mutations. Insertion of bulky aromatic residues in the tunnel corresponding to pathway p1 leads to reduced accessibility to the ligands and a change in mechanism of opening from permanent to transient. We propose that engineering the accessibility of tunnels and the mechanisms of ligand exchange is a powerful strategy for modification of the functional properties of enzymes with buried active sites.

Hybrid quantum mechanics/molecular mechanics (QM/MM) simulations provide a powerful tool for studying chemical reactions, especially in complex biochemical systems. In most works to date, the quantum region is kept fixed throughout the simulation and is defined in an ad hoc way based on chemical intuition and available computational resources. The simulation errors associated with a given choice of the quantum region are, however, rarely assessed in a systematic manner. Here we study the dependence of two relevant quantities on the QM region size: the force error at the center of the QM region and the free energy of a proton transfer reaction. Taking lysozyme as our model system, we find that in an apolar region the average force error rapidly decreases with increasing QM region size. In contrast, the average force error at the polar active site is considerably higher, exhibits large oscillations and decreases more slowly, and may not fall below acceptable limits even for a quantum region radius of 9.0 A. Although computation of free energies could only be afforded until 6.0 A, results were found to change considerably within these limits. These errors demonstrate that the results of QM/MM calculations are heavily affected by the definition of the QM region (not only its size), and a convergence test is proposed to be a part of setting up QM/MM simulations.

The selection of a proper reaction coordinate is a major bottleneck in simulations of chemical reactions in complex systems. Increasing the number of variables that are used to bias the reaction largely affects the convergence and leads to an unbearable increase in computational price. This problem can be overcome by employing a complex reaction coordinate that depends on many geometrical variables of the system, such as the energy gap (EGAP) in the empirical valence bond (EVB) method. EGAP depends on all of the coordinates of the system, and its robustness has been demonstrated for a variety of enzymatic reactions. In this work, we demonstrate that EGAP, derived from a classical representation, can be used as a reaction coordinate in systems described with any quantum chemistry Hamiltonian. Benefits of using EGAP as a reaction coordinate as compared to a traditional geometrical variable are illustrated in the case of a symmetric nucleophilic substitution reaction in water solution. EGAP is shown to provide a significantly more efficient sampling and allows a better localization of the transition state as compared to a geometrical reaction coordinate.

We present the results of a series of 10-ns molecular dynamics simulations on Pseudomonas aeruginosa lectin-II (PA-IIL) and its complexes with four different monosaccharides. We compare the saccharide-free, saccharide-occupied, and saccharide- and ion-free forms of the lectin. The results are coupled with analysis of the water density map and calcium coordination. The water density pattern around the binding site in the free lectin molecular dynamics was fitted with that in the X-ray and with the hydroxyl groups of the monosaccharide within the lectin/monosaccharide complexes and the best ligand was predicted based on the best fit. Interestingly, the water density pattern around the binding site in the uncomplexed lectin exactly fitted the O2, O3, and O4 hydroxyl groups of the fucose complex with the lectin. This observation could lead to a hypothesis that the replacement of these three water molecules from the binding site by the monosaccharide decreases the entropy of the complex and increases the entropy of the water molecules, which favors the binding. It suggests that the high density peaks of the solvent around the binding site in the free protein could be the tool to predict hydroxyl group orientation of the sugar in the protein/sugar complexes. The high affinity of PA-IIL binding site is also attributed to the presence of two calcium ions, each of them making five to six coordinations with the protein part and two coordinations with either water or the monosaccharide. When the calcium ions are removed from the simulated system, they are replaced by sodium ions from the solvent. These observations rationalize the high binding affinity of PA-IIL towards fucose.

The choreography of restriction endonuclease catalysis is a long-standing paradigm in molecular biology. Bivalent metal ions are required almost for all PD..D/ExK type enzymes, but the number of cofactors essential for the DNA backbone scission remained ambiguous. On the basis of crystal structures and biochemical data for various restriction enzymes, three models have been developed that assign critical roles for one, two, or three metal ions during the phosphodiester hydrolysis. To resolve this apparent controversy, we investigated the mechanism of BamHI catalysis using quantum mechanical/molecular mechanical simulation techniques and determined the activation barriers of three possible pathways that involve a Glu-113 or a neighboring water molecule as a general base or an external nucleophile that penetrated from bulk solution. The extrinsic mechanism was found to be the most favorable with an activation free energy of 23.4 kcal/mol, in reasonable agreement with the experimental data. On the basis of the effect of the individual metal ions on the activation barrier, metal ion A was concluded to be pivotal for the reaction, while the enzyme lacking metal ion B still has moderate efficiency. Thus, we propose that the catalytic scheme of BamHI does not involve a general base for nucleophile generation and requires one obligatory metal ion for catalysis that stabilizes the attacking nucleophile and coordinates it throughout the nucleophilic attack. Such a model may also explain the variation in the number of metal ions in the crystal structures and thus could serve as a framework for a unified catalytic scheme of type II restriction endonucleases.

International audience; Explicit solvent molecular dynamics simulations (in total almost 800 ns including locally enhanced sampling runs) were applied with different ion conditions and with two force fields (AMBER and CHARMM) to characterize typical geometries adopted by the flanking bases in the RNA kissing-loop complexes. We focus on flanking base positions in multiple x-ray and NMR structures of HIV-1 DIS kissing complexes and kissing complex from the large ribosomal subunit of Haloarcula marismortui. An initial x-ray open conformation of bulged-out bases in HIV-1 DIS complexes, affected by crystal packing, tends to convert to a closed conformation formed by consecutive stretch of four stacked purine bases. This is in agreement with those recent crystals where the packing is essentially avoided. We also observed variants of the closed conformation with three stacked bases, while nonnegligible populations of stacked geometries with bulged-in bases were detected, too. The simulation results reconcile differences in positions of the flanking bases observed in x-ray and NMR studies. Our results suggest that bulged-out geometries are somewhat more preferred, which is in accord with recent experiments showing that they may mediate tertiary contacts in biomolecular assemblies or allow binding of aminoglycoside antibiotics.

Large RNA molecules exhibit an astonishing variability of base-pairing patterns, while many of the RNA base-pairing families have no counterparts in DNA. The cis Watson-Crick/sugar edge (cis WC/SE) RNA base pairing is investigated by ab initio quantum chemical calculations. A detailed structural and energetic characterization of all 13 crystallographically detected members of this family is provided by means of B3LYP/ 6-31G** and RIMP2/aug-cc-pVDZ calculations. Further, a prediction is made for the remaining 3 cis WC/ SE base pairs which are yet to be seen in the experiments. The interaction energy calculations point at the key role of the 2 '-OH group in stabilizing the sugar-base contact and predict all 16 cis WC/SE base-pairing patterns to be nearly isoenergetic. The perfect correlation of the main geometrical parameters in the gas-phase optimized and X-ray structures shows that the principle of isosteric substitutions in RNA is rooted from the intrinsic structural similarity of the isolated base pairs. The present quantum chemical calculations for the first time analyze base pairs involving the ribose 2 '-OH group and unambiguously correlate the structural information known from experiments with the energetics of interactions. The calculations further show that the relative importance and absolute value of the dispersion energy in the cis WC/SE base pairs are enhanced compared to the standard base pairs. This may by an important factor contributing to the strength of such interactions when RNA folds in its polar environment. The calculations further demonstrate that the Cornell et al. force field commonly used in molecular modeling and simulations provides satisfactory performance for this type of RNA interactions.

The pinched cone-pinched cone interconversion (PCI) of tetraethoxy substituted calix[4]arene and thiacalix[4]arene was studied by means of molecular dynamics using the GAFF force field and quantum mechanics. Influence of solvent was expressed by changes in cavity geometries as well as in free energy barriers that were calculated using constrained MD simulations. Water and chloroform were found to reduce the energy barriers in comparison with vacuum simulations. Only qualitative agreement between energy barriers calculated by MD and experimental values was achieved. Also, quantum chemical calculations were performed with the RI-MP2 method, and a better description of PCI was achieved than that previously obtained by HF or DFT methods. The combination of qualitative MD data on solvated systems with RI-MP2 gas-phase data gave us very good agreement with experimental results.

The study of the temperature-sensitive photochemical release of a carboxylic acid from 2,5-dimethylphenacyl ester is reported. Quantum yields of the benzoate ester degradation in benzene increased from 0.22 at room temperature to 0.28 at 50 degrees C whereas a more significant increase (nearly by a factor of 3) was observed in methanol and ethanol, reaching a high reaction efficiency (0.25) typically found in non-polar solvents. The reaction proceeds predominantly via the triplet pathway and the E-photoenol in the whole temperature range in methanol solution. A higher quantum efficiency in heated methanol is explained by enhancing the E-photoenol population. This picture was partially confirmed by the quantum chemical calculations. The 2,5-dimethylphenacyl chromophore is proposed as an efficient photoremovable protecting group for carboxylic acids in solutions under conventional or microwave-assisted heating for applications in organic synthesis, such as the solid-phase synthesis.

A temperature-sensitive photochemical nucleophilic aromatic substitution on 4-nitroanisole by a hydroxide ion in homogeneous solutions, in a two-phase system under phase-transfer catalysis conditions, and in the microwave field is reported. It was found that reaction regioselectivity dramatically changes with temperature in the region of -20 to 196 degrees C. The quantum yield of the 4-methoxyphenol formation was found to be temperature independent, in contrast to that of the 4-nitrophenol formation, suggesting that there is a temperature dependent process occurring after the partitioning between replacement of the nitro group and the methoxy group has taken place. The reaction was also investigated by using quantum chemical calculations. A technique for microwave-assisted photochemical synthesis is proposed as an efficient and practical tool for organic synthesis.

Very important part of theoretical simulations is a process of structure preparation, geometry editing and detailed examination of properties. Article introduces advanced molecular modelling package called Nemesis.

The study presents quantum-chemical calculations on proton transfer in protonated N-acetylglycyl-N1-methylglycinamide (AGA) as a short oligopeptide model. All calculations employ the B3LYP functional and the 6-31++G** basis set. Two different mechanisms of proton transfer are discussed. The rate-determining step of the first mechanism exhibits an energy barrier of about 17.7 kcal mol-1, and it is represented by an isomerization of the proton around the double bond of the carbonyl group. The second mechanism is based on the large conformational flexibility of AGA, where all carbonyl oxygens cooperate. The rate-determining step of this mechanism exhibits an energy barrier of only 8.3 kcal mol-1.

Since the discovery of the cucurbit[6]uril (CB6) macrocycle, the binding modes of its supramolecular complexes with organic cations have been elucidated using NMR spectroscopy. The observed changes in chemical shifts of guests induced by this complexation follow generally accepted rules. Here we report some exceptions to this usual behavior for complexes of CB6 with pyridinium and pyrazinium guests, as determined by using 1H, 13C, and 15N NMR spectroscopy. We use quantum chemical calculations to rationalize our observations by the concept of redistribution of the electron density in the aromatic part of the guest caused by complexation with the macrocycle. Furthermore, we used a simple model employing point charges that mimic the CB6host to analyze the electron deformation density induced by complexation. It indicates that the changes in the nuclear magnetic shielding observed for the aromatic part of the guest result from the redistribution of electrons induced mainly by the electrostatic repulsion of the oxygen atoms of the CB6 portal.