Differential stabilization of adenine quartets by anions and cations

We have investigated the structures and stabilities of four different adenine quartets with alkali and halide ions in the gas phase and in water, using dispersion-corrected density functional theory at the BLYP-D/TZ2P level. First, we examine the empty quartets and how they interact with alkali cations and halide anions with formation of adenine quartet–ion complexes. Second, we examine the interaction in a stack, in which a planar adenine quartet interacts with a cation or anion in the periphery as well as in the center of the quartet. Interestingly, for the latter situation, we find that both cations and anions can stabilize a planar adenine quartet in a stack.     

Coordination numbers in hydrated Cu(II) ions

The potential energy surface of [Cu(H₂O)_n]²⁺ clusters with n?=?12, 16, and 18 was explored by using a modified version of the simulated annealing method. Such exploration was carried out by using the PM7 semiempirical method to obtain around 100,000 isomers, which provide candidates to be optimized with PBE0-D3, M06-2X, and BHLYP exchange-correlation functionals coupled with the 6–311++G** basis set. These methods based on the Kohn-Sham approach delivered isomers with coordination numbers of 4, 5, and 6. The analysis used to obtain coordination numbers was based on geometrical parameters and the quantum theory of atoms in molecules (QTAIM) approach. Our methodology found only one isomer with fourfold coordination and its probabilities to appear in these clusters are quite small for high temperatures. The procedure used in this article predicts important populations of fivefold and sixfold coordination clusters, in fact, the fivefold coordination dominates for PBE0-D3 and BHLYP methods, although the sixfold coordination starts to be important when the number of water molecules is increased. The nature of axial and equatorial contacts is discussed in the context of the QTAIM and the noncovalent interaction index (NCI), which gives a clear classification of such orientations. Also, these methods suggest a partial covalent interaction between the Cu²⁺ and water molecules in both positions; equatorial and axial.     

Calculating the geometry and Raman spectrum of physiological bis(<Emphasis Type="SmallCaps">l</Emphasis>-histidinato)copper(II): an assessment of DFT functionals for aqueous and isolated systems

Reliable density functional theory (DFT) calculations can be performed in conjuction with spectroscopic measurements to elucidate the structural properties of physiologically important bis(amino acidato)copper(II) compounds in solutions. They can provide insight into the influence of intermolecular interactions on the molecular geometry in the crystal lattice or solution when compared with a DFT gas-phase minimum. Our previous paper [Markovi? et al. (2014) Eur J Inorg Chem 198] reported the DFT-determined geometries and Raman spectra for different conformers of physiological bis(l-histidinato)copper(II) with 20 explicit water molecules, as calculated using the B3LYP functional. The present study examined the reliability of those B3LYP results by applying the M06 functional instead, as it should better account for noncovalent interactions. The water molecules were positioned more compactly around the complex by M06 than by B3LYP. The accuracies of the two functionals when compared to relevant experimental data showed that M06 was better at reproducing in-plane Cu?N bond lengths but B3LYP gave more accurate axial Cu?O distances. Both functionals reproduced the experimental Raman spectrum at pH 8 to similar levels of accuracy and provided precise information on the Cu(II) coordination mode and conformation in aqueous solution. Additionally, we assessed several DFT and DFT-D functionals (BP86, B3LYP, B3LYP-D, M06, M06 L, wB97XD, mPW2PLYPD) by using them to model the geometries of experimental bis(l-histidinato)copper(II) crystalline conformations as isolated systems, and then benchmarking the results against those from high-level second-order pertubation Møller–Plesset (MP2) calculations. Although this assessment resulted in an equivocal conclusion because the MP2 results for the isolated complex were inconsistent with the corresponding DFT outcomes, it does provide new information on future benchmark options.     

Potential of mean force for protein-protein interaction studies.

Calculating protein-protein interaction energies is crucial for understanding protein-protein associations. On the basis of the methodology of mean-field potential, we have developed an empirical approach to estimate binding free energy for protein-protein interactions. This knowledge-based approach has been used to derive distance-dependent free energies of protein complexes from a nonredundant training set in the Protein Data Bank (PDB), with a careful treatment of homology. We calculate atom pair potentials for 16 pair interactions, which can reflect the importance of hydrophobic interactions and specific hydrogen-bonding interactions. The derived potentials for hydrogen-bonding interactions show a valley of favorable interactions at a distance of approximately 3 A, corresponding to that of an established hydrogen bond. For the test set of 28 protein complexes, the calculated energies have a correlation coefficient of 0.75 compared with experimental binding free energies. The performance of the method in ranking the binding energies of different protein-protein complexes shows that the energy estimation can be applied to value binding free energies for protein-protein associations.     

Identification of a critical lysine residue in apolipoprotein B-100 that mediates noncovalent interaction with apolipoprotein(a)

We have previously shown that lipoprotein(a) (Lp(a)) assembly involves an initial noncovalent interaction between sequences within apolipoprotein(a) (apo(a)) kringle IV types 5-8 and the amino terminus of apolipoprotein B-100 (sequences between amino acids 680 and 781 in apoB-100), followed by formation of a disulfide bond. In the present study, citraconylation of lysine residues in apoB-100 abolished the ability of the modified low density lipoprotein to associate with apo(a), thereby demonstrating a direct role for lysine residues in apoB in the first step of Lp(a) assembly. To identify specific lysine residues in the amino terminus of apoB that are required for the noncovalent interaction, we initially used an affinity chromatography method in which recombinant forms of apo(a) (r-apo(a)) were immobilized on Sepharose beads. Assessment of the ability of carboxyl-terminal truncations of apoB-18 to bind to r-apo(a)-Sepharose revealed that a 25-amino acid sequence in apoB (amino acids 680-704) bound specifically to apo(a) in a lysine-dependent manner; citraconylation of the lysine residues in the apoB derivative encoding this sequence abolished the binding interaction. Using fluorescence spectrometry, we found that a synthetic peptide corresponding to this sequence bound directly to apo(a); the peptide also reduced covalent Lp(a) formation. Lysine residues present in this sequence (Lys(680) and Lys(690)) were mutated to alanine in the context of apoB-18. We found that the apoB-18 species containing the Lys(680) mutation was incapable of binding to r-apo(a)-Sepharose columns, whereas the apoB-18 species containing the Lys(690) mutation exhibited slightly reduced binding to these columns. Taken together, our data indicate that Lys(680) is critical for the noncovalent interaction of apo(a) and apoB-100 that precedes covalent Lp(a) formation.     

Establishment of the Embryo-derived Stem (ES) Cell Lines from Mouse Blastocysts: Effects of the Feeder Cell Layer.

The ES ceii lines are embryo-derived stem cell lines directly isolated from the inner cell mass of mouse blastocysts using feeder cell layer. We have established a number of ES cell lines from 129 or C57BL/6 strain mice by using the feeder layer of the STO cells (from ATCC) or the primary embryonic fibroblasts, which was obtained by trypsinizing the 16-day-old BALB/c mouse fetus. The ES cell lines established on the STO feeder layer showed differentiation into various tissues in solid tumors when injected into syngenic mice. Karyotype was, however, nearly tetraploid. The ES cell lines established on the primary fibroblasts exhibited differentiation into larger variety of tissues in solid tumors. Karyotype was almost diploid and majority of the cells kept normal set of chromosomes in G-banding. We conclude that the primary fibroblasts are better feeder layer than the STO cells for establishment and maintenance of the ES cell lines.     

Noncovalent chiral domino effect on one-handed helix of nonapeptide containing a midpoint L-residue

We here clarify whether noncovalent chiral domino effect characterized by the terminal interaction of a helical peptide with a chiral small molecule can alter the helical stability of N-deprotected peptides containing an L-residue covalently incorporated into the inner position. Two nonapeptides consisting of the midpoint L-leucine (1) or L-phenylalanine (2) and the achiral helix-forming residues were employed. NMR and IR spectroscopy and energy calculation indicated that both peptides adopt a 3(10)-helical conformation in chloroform. They strongly preferred a right-handed screw sense because of the presence of the midpoint L-residue. These original right-handed screw senses were retained on addition of chiral Boc-amino acid, but their helical stabilities clearly depended on its added chirality. Here, Boc-L-amino acid stabilizes the original right-handed helix, whereas the corresponding Boc-D-amino acid tends to less stabilize or destabilize it. This tendency was not observed for the corresponding N-Boc-protected peptides 1 and 2, strongly suggesting that the N-terminal amino group is required for controlling the stabilization of the original right-handed helix. Therefore, noncovalent chiral domino effect in peptides 1 and 2 can contribute even to the helical stability of a chiral peptide prevailing one-handed helix strongly through the midpoint L-residue. In addition, the N-terminal moiety of a 3(10)-helical peptide was found to generate chiral discrimination in complexation process with racemic additives.     

Predicting Binding Free Energy Change Caused by Point Mutations with Knowledge-Modified MM/PBSA Method

A new methodology termed Single Amino Acid Mutation based change in Binding free Energy (SAAMBE) was developed to predict the changes of the binding free energy caused by mutations. The method utilizes 3D structures of the corresponding protein-protein complexes and takes advantage of both approaches: sequence- and structure-based methods. The method has two components: a MM/PBSA-based component, and an additional set of statistical terms delivered from statistical investigation of physico-chemical properties of protein complexes. While the approach is rigid body approach and does not explicitly consider plausible conformational changes caused by the binding, the effect of conformational changes, including changes away from binding interface, on electrostatics are mimicked with amino acid specific dielectric constants. This provides significant improvement of SAAMBE predictions as indicated by better match against experimentally determined binding free energy changes over 1300 mutations in 43 proteins. The final benchmarking resulted in a very good agreement with experimental data (correlation coefficient 0.624) while the algorithm being fast enough to allow for large-scale calculations (the average time is less than a minute per mutation).     

Atomistic simulations of materials for optical chemical sensors: DFT-D calculations of molecular interactions between gas-phase analyte molecules and simple substrate models

The structures of complexes of some small molecules (formaldehyde, acetaldehyde, ammonia, methylamine, methanol, ethanol, acetone, benzene, acetonitrile, ethyl acetate, chloroform, and tetrahydrofuran, considered as possible analytes) with ethylbenzene and silanol (C₆H₅C₂H₅ and SiH₃OH, considered as models of polystyrene and silica gel substrates) and with acridine (C₁₃H₉N, considered as a model of an indicator dye molecule of the acridine series) and the corresponding interaction energies have been calculated using the DFT-D approximation. The PBE exchange-correlation potential was used in the calculations. The structures of complexes between the analyte and the substrate were determined by optimizing their ground-state geometry using the SVP split-valence double-zeta plus polarization basis set. The complex formation energies were refined by single-point calculations at the calculated equilibrium geometries using the sufficiently large triple-zeta TZVPP basis set. The calculated interaction energies are used to assess the possibility of using dyes of the acridine series adsorbed on a polystyrene or silica substrate for detecting the small molecules listed above.     

An improved generalized AMBER force field (GAFF) for urea

We describe an improved force field parameter set for the generalized AMBER force field (GAFF) for urea. Quantum chemical computations were used to obtain geometrical and energetic parameters of urea dimers and larger oligomers using AM1 semiempirical MO theory, density functional theory at the B3LYP/6-31G(d,p) level, MP2 and CCSD ab initio calculations with the 6-311++G(d,p), aug-cc-pVDZ, aug-cc-pVTZ, and aug-cc-pVQZ basis sets, and with the CBS-QB3 and CBS-APNO complete basis set methods. Seven different urea dimer structures were optimized at the MP2/aug-cc-pVDZ level to obtain accurate interaction energies. Atomic partial charges were calculated at the MP2/aug-cc-pVDZ level with the restrained electrostatic potential (RESP) fitting approach. The interaction energies computed with these new RESP charges in the force field are consistent with those obtained from CCSD and MP2 calculations. The linear dimer structure calculated using the force field with modified geometrical parameters and the new RESP charge set agrees well with available experimental data.     

Basis set validation for polyatomic cation-water interactions

A study of the effects of counterpoise (CP) corrections in polyatomic cation-water interactions is added to the systematic analysis performed in the past of the basis set superposition error (BSSE) in neutral and anionic adducts. The interaction with water of the ammonium cation and its methyl and -CHO derivatives is considered due to the need to model accurately this functional group, which is common in biological molecules. The basis sets employed are the STO-3G and MINI-1 minimal basis sets and the 3-21G, 4-31G and 6-31G** extended ones. In addition, the 6-311++G (2d,p) and 6-311++G (3d, 2p) basis sets have been used for the smallest system, at the SCF and MP2 levels, both without and with CP correction. These basis sets give an equilibrium distance slightly larger and an interaction energy less favourable than the 6-31G** basis set at the corresponding level, while the inclusion of correlation corrections produces a stronger H-bond at a shorter distance. The results confirm the previous hint of a lower incidence of BSSEs in medium size cationic systems, at a different extent for the various basis sets. While the STO-3G basis set is sharply affected by BSSEs in both the equilibrium distance and the interaction energy, the MINI-1 basis set shows a small BSSE, though its trend is not completely satisfactory because the charge transfer component has an anomalous behaviour with respect to our reference basis set. The 4-31G basis set is the only one able to hold comparison with the 6-31G**, even though the interaction energy produced is slightly overestimated. The 3-21G basis set, when corrected, almost parallels the 4-31G one in this set of compounds. The reliability of the CPED correction is checked and discussed.     

MKS5 and CEP290 Dependent Assembly Pathway of the Ciliary Transition Zone

Cilia have a unique diffusion barrier (“gate”) within their proximal region, termed transition zone (TZ), that compartmentalises signalling proteins within the organelle. The TZ is known to harbour two functional modules/complexes (Meckel syndrome [MKS] and Nephronophthisis [NPHP]) defined by genetic interaction, interdependent protein localisation (hierarchy), and proteomic studies. However, the composition and molecular organisation of these modules and their links to human ciliary disease are not completely understood. Here, we reveal Caenorhabditis elegans CEP-290 (mammalian Cep290/Mks4/Nphp6 orthologue) as a central assembly factor that is specific for established MKS module components and depends on the coiled coil region of MKS-5 (Rpgrip1L/Rpgrip1) for TZ localisation. Consistent with a critical role in ciliary gate function, CEP-290 prevents inappropriate entry of membrane-associated proteins into cilia and keeps ARL-13 (Arl13b) from leaking out of cilia via the TZ. We identify a novel MKS module component, TMEM-218 (Tmem218), that requires CEP-290 and other MKS module components for TZ localisation and functions together with the NPHP module to facilitate ciliogenesis. We show that TZ localisation of TMEM-138 (Tmem138) and CDKL-1 (Cdkl1/Cdkl2/Cdkl3/Cdlk4 related), not previously linked to a specific TZ module, similarly depends on CEP-290; surprisingly, neither TMEM-138 or CDKL-1 exhibit interdependent localisation or genetic interactions with core MKS or NPHP module components, suggesting they are part of a distinct, CEP-290-associated module. Lastly, we show that families presenting with Oral-Facial-Digital syndrome type 6 (OFD6) have likely pathogenic mutations in CEP-290-dependent TZ proteins, namely Tmem17, Tmem138, and Tmem231. Notably, patient fibroblasts harbouring mutated Tmem17, a protein not yet ciliopathy-associated, display ciliogenesis defects. Together, our findings expand the repertoire of MKS module-associated proteins—including the previously uncharacterised mammalian Tmem80—and suggest an MKS-5 and CEP-290-dependent assembly pathway for building a functional TZ.     

Revealing the nature of intermolecular interaction and configurational preference of the nonpolar molecular dimers (H2)2, (N2)2, and (H2)(N2)

Understanding the nature of noncovalent interactions between nonpolar small molecules is not only theoretically interesting but also important for practical purposes. The interaction mechanism of three prototype dimers (H₂)₂, (N₂)₂, and (H₂)(N₂) are investigated by state-of-the-art quantum chemistry calculations and energy decomposition analysis. It is shown that their configuration preferences are essentially controlled by the electrostatic component rather than the dispersion effect though the monomers have zero dipole moment. These configuration preferences can also be fairly well and conveniently interpreted by visually examining the electrostatic potential map.     

Synergistic cardioprotective effects of rAAV9-CyclinA2 combined with fibrin glue in rats after myocardial infarction

The present study aimed to investigate the protective effects of rAAV9-CyclinA2 combined with fibrin glue (FG) in vivo in rats after myocardial infarction (MI). Ninety male Sprague–Dawley rats were randomized into 6 groups (15 in each group): sham, MI, rAAV9-green fluorescent protein (GFP)?+?MI, rAAV9-CyclinA2?+?MI, FG?+?MI, and rAAV9-CyclinA2?+?FG?+?MI. Packed virus (5?×?10¹¹vg/ml) in 150 µl of normal saline or FG was injected into the infarcted myocardium at five locations in rAAV9-GFP?+?MI, rAAV9-CyclinA2?+?MI, and rAAV9-CyclinA2?+?FG?+?MI groups. The sham, MI, and FG?+?MI groups were injected with an equal volume of normal saline or FG at the same sites. Five weeks after injection, echocardiography was performed to evaluate the left ventricular function. The expressions of CyclinA2, proliferating cell nuclear antigen (PCNA), and phospho-histone-H3 (H3P), vascular density, and infarct area were assessed by Western blot, immunohistochemistry, immunofluorescence, and Masson staining. As a result, the combination of rAAV9-CyclinA2 and FG increased ejection fraction and fractional shortening compared with FG or rAAV9-CyclinA2 alone. The expression level of CyclinA2 was significantly higher in the rAAV9-CyclinA2?+?FG?+?MI group compared with the rAAV9-CyclinA2?+?MI and FG?+?MI groups (70.1?±?1.86% vs. 14.74?±?2.02%, P?<?0.01; or vs. 50.13?±?3.80%; P?<?0.01). A higher expression level of PCNA and H3P was found in the rAAV9-CyclinA2?+?FG?+?MI group compared with other groups. Comparing with other experiment groups, collagen deposition and the infarct size significantly decreased in rAAV9-CyclinA2?+?Fibrin?+?MI group. The vascular density was much higher in the rAAV9-CyclinA2?+?FG?+?MI group compared with the rAAV9-CyclinA2?+?MI group. We concluded that fibrin glue combined with rAAV9-CyclinA2 was found to be effective in cardiac remodeling and improving myocardial protection.     

A theoretical investigation on the proton transfer tautomerization mechanisms of 2-thioxanthine within microsolvent and long range solvent

A relative complete study on the mechanisms of the proton transfer reactions of 2-thioxanthine was carried out with density functional theory. The models were designed with monohydrated and dihydrated microsolvent catalyses either with or without the presence of water solvent considered with the polarized continuum model (PCM). A total number of 114 complexes and 67 transition states were found with the B3LYP/6-311+G** calculations. The energies were refined with both B3LYP/aug-cc-pVTZ and PCM-B3LYP/aug-cc-pVTZ methods. The activation energies were reported with respect to the Gibbs free energies obtained in conjunction with the standard statistical thermodynamics. Possible reaction pathways were confirmed with the intrinsic reaction coordinates. Pathways via C8 atom on the imidazole ring, via the bridged C4 and C5 atoms between pyrimidine and imidazole rings and via N, O and S atom on the pyrimidine ring were examined. The results show that the most feasible pathway is the proton transfers within the long range solvent surrounding via the N, O and S atoms in the pyrimidine ring with di-hydrated catalysis: N(7)H?+?2H₂O?→?IM89?→?IM90?→?P13?+?2H₂O?→?IM91?→?IM92?→?P6?+?2H₂O?→?IM71?→?IM72?→?P7?+?2H₂O?→?IM107?→?IM108?→?P18?+?2H₂O?→?IM111?→?IM112?→?P19?+?2H₂O?→?IM113?→?IM114?→?P17?+?2H₂O?→?IM105?→?IM106?→?N(9)H?+?2H₂O that has the highest energy barrier of 44.0 kJ mol^?1 in the transition of IM89 to IM90 via TS54. The small energy barrier is in good agreement with the experimental observation that 2-TX tautomerizes at room temperature in water. In the aqueous phase, the most stable intermediate is found to be IM21 [N(7)H?+?2H₂O] and the possible co-existing species are the monohydrated IM1, IM9, IM39 and IM46, and the di-hydrated IM5, IM8, IM13, IM16, IM81, IM89, IM90, IM91 and IM106 complexes that have a relative concentration larger than 10^?6 (1 ppm) with respect to IM21.

The role of amino-acid residues Q39 and E451 in the determination of substrate specificity of the Spodoptera frugiperda beta-glycosidase.

The design of beta-glycosidases with planed substrate specificity for biotechnological application has received little attention. This is mostly a consequence of the lack of data on the molecular basis of the beta-glycosidase specificity, namely data on the energy of the noncovalent interactions in the enzyme-transition state complex. In an attempt to fill this gap, site-directed mutagenesis and enzyme steady-state kinetic experiments with different substrates were conducted, using as model a digestive beta-glycosidase (glycoside hydrolase family 1) from Spodoptera frugiperda (Lepidoptera) (Sfbetagly50). The active site of this enzyme was modeled based on its sequence and on crystallographic data of similar enzymes. Energy of noncovalent interactions in transition state between Sfbetagly50 amino acids and glycone hydroxyls was determined. Sfbetagly50 residue E451 seems to be a key residue in determining beta-glycosidase preference for glucosides vs. galactosides based on the following data: (a) The energy of the noncovalent interaction between glycone equatorial hydroxyl 4 with E451 in the transition state is about 60% higher than its interaction with Q39. (b) The energy of the E451-hydroxyl 4 interaction decreases more than the Q39-hydroxyl 4 interaction when hydroxyl 4 is changed from equatorial to axial position. (c) A Sfbetagly50 mutant, where E451 was substituted by A, hydrolyzes galactosides faster than glucosides. It was also shown that glycone hydroxyl 6 interacts favorably with Q39, but not with E451, probably due to steric hindrance. These interactions result in the beta-glycosidase hydrolyzing fucosides (6-deoxygalactosides) faster than glucosides and galactosides.     

Molecular mechanical perspective on halogen bonding

The nature and strength of halogen bonding in halo molecule-Lewis base complexes were studied in terms of molecular mechanics using our recently developed positive extra-point (PEP) approach, in which the σ-hole on the halogen atom is represented by an extra point of positive charge. The contributions of the σ-hole (i.e., positively charged extra point) and the halogen atom to the strength of this noncovalent interaction were clarified using the atomic parameter contribution to the molecular interaction (APCtMI) approach. The molecular mechanical results revealed that the halogen bond is electrostatic and van der Waals in nature, and its strength depends on three types of interaction: (1) the attractive electrostatic interaction between the σ-hole and the Lewis base, (2) the repulsive electrostatic interaction between the negative halogen atom and the Lewis base, and (3) the repulsive/attractive van der Waals interactions between the halogen atom and the Lewis base. The strength of the halogen bond increases with increasing σ-hole size (i.e., magnitude of the extra-point charge) and increasing halogen atom size. The van der Waals interaction's contribution to the halogen bond strength is most favorable in chloro complexes, whereas the electrostatic interaction is dominant in iodo complexes. The idea that the chloromethane molecule can form a halogen bond with a Lewis base was revisited in terms of quantum mechanics and molecular mechanics. Although chloromethane does produce a positive region along the C-Cl axis, basis set superposition error corrected second-order M?ller-Plesset calculations showed that chloromethane-Lewis base complexes are unstable, producing halogen-Lewis base contacts longer than the sum of the van der Waals radii of the halogen and O/N atoms. Molecular mechanics using the APCtMI approach showed that electrostatic interactions between chloromethane and a Lewis base are unfavorable owing to the high negative charge on the chlorine atom, which overcomes the corresponding favorable van der Waals interactions.     

Quantitative analysis of the interaction strength and dynamics of human IgG4 half molecules by native mass spectrometry

Native mass spectrometry (MS) is a powerful technique for studying noncovalent protein-protein interactions. Here, native MS was employed to examine the noncovalent interactions involved in homodimerization of antibody half molecules (HL) in hinge-deleted human IgG4 (IgG4Δhinge). By analyzing the concentration dependence of the relative distribution of monomer HL and dimer (HL)(2) species, the apparent dissociation constant (K(D)) for this interaction was determined. In combination with site-directed mutagenesis, the relative contributions of residues at the CH3-CH3 interface to this interaction could be characterized and corresponding K(D) values quantified over a range of 10(-10)-10(-4) M. The critical importance of this noncovalent interaction in maintaining the intact dimeric structure was also proven for the full-length IgG4 backbone. Using time-resolved MS, the kinetics of the interaction could be measured, reflecting the dynamics of IgG4 HL exchange. Hence, native MS has provided a quantitative view of local structural features that define biological properties of IgG4.