Conformational preferences and pKa value of cysteine residue

The garlic compound ajoene covalently binds vimentin, disrupts the vimentin network and exerts anti-metastatic activity in cancer cells

BACKGROUND

Garlic has been used for centuries for its flavour and health promoting properties that include protection against cancer. The vinyl disulfide-sulfoxide ajoene is one of the phytochemicals found in crushed cloves, hypothesised to act by S-thiolating reactive cysteines in target proteins.

METHODS

Using our fluorescently labelled ajoene analogue called dansyl-ajoene, ajoene's protein targets in MDA-MB-231 breast cancer cells were tagged and separated by 2D electrophoresis. A predominant band was identified by MALDI-TOF MS/MS to be vimentin. Target validation experiments were performed using pure recombinant vimentin protein. Computational modelling of vimentin bound to ajoene was performed using Schrödinger and pKa calculations by Epik software. Cytotoxicity of ajoene in MDA-MB-231 and HeLa cells was measured by the MTT assay. The vimentin filament network was visualised in ajoene-treated and non-treated cells by immunofluorescence and vimentin protein expression was determined by immunoblot. The invasion and migration activity was measured by wound healing and transwell assays using wildtype cells and cells in which the vimentin protein had been transiently knocked down by siRNA or overexpressed.

RESULTS

The dominant protein tagged by dansyl-ajoene was identified to be the 57 kDa protein vimentin. The vimentin target was validated to reveal that ajoene and dansyl-ajoene covalently bind to recombinant vimentin via a disulfide linkage at Cys-328. Computational modelling showed Cys-328 to be exposed at the termini of the vimentin tetramer. Treatment of MDA-MB-231 or HeLa cells with a non-cytotoxic concentration of ajoene caused the vimentin filament network to condense; and to increase vimentin protein expression. Ajoene inhibited the invasion and migration of both cancer cell lines which was found to be dependent on the presence of vimentin. Vimentin overexpression caused cells to become more migratory, an effect that was completely rescued by ajoene.

CONCLUSIONS

The garlic-derived phytochemical ajoene targets and covalently modifies vimentin in cancer cells by S-thiolating Cys-328. This interaction results in the disruption of the vimentin filament network and contributes to the anti-metastatic activity of ajoene in cancer cells.

Cysteine-Ag Cluster Hydrogel Confirmed by Experimental and Numerical Studies

The native cysteine (Cys)-Ag3 cluster hydrogel is approved for the first time by both experimental and theoretical studies. From the detailed molecular structure and energy information, three factors are found to ensure the self-assembly of Cys and Ag3 , and result in the hydrogel. First, the Ag-S bonds make Cys and Ag3 form Cys-Ag3 -Cys monomer. Second, intermolecular hydrogen bonds between carboxyl groups of adjacent monomer push them self-assembled. Third, more monomer precisely self-assemble to produce the -[Cys-Ag3 -Cys]n multimer, e.g., a single molecular chain with the left-handed helix conformation, via a benign thermodynamic process. These multimers entangle together to form micro-network to trap water and produce hydorgel in situ. The hydrogen bonds of hydrogel are sensitive to thermal and proton stimuli, and the hydrogel presents lysosome targeting properties via fluorescent imaging with biocompatibility.

Binding selectivity of dibenzo-18-crown-6 for alkali metal cations in aqueous solution: A density functional theory study using a continuum solvation model

Background

Dibenzo-18-crown-6 (DB18C6) exhibits the binding selectivity for alkali metal cations in solution phase. In this study, we investigate the main forces that determine the binding selectivity of DB18C6 for the metal cations in aqueous solution using the density functional theory (DFT) and the conductor-like polarizable continuum model (CPCM).

Results

The bond dissociation free energies (BDFE) of DB18C6 complexes with alkali metal cations (M⁺-DB18C6, M = Li, Na, K, Rb, and Cs) in aqueous solution are calculated at the B3LYP/6-311++G(d,p)//B3LYP/6-31 + G(d) level using the CPCM. It is found that the theoretical BDFE is the largest for K⁺-DB18C6 and decreases as the size of the metal cation gets larger or smaller than that of K⁺, which agrees well with previous experimental results.

Conclusion

The solvation energy of M⁺-DB18C6 in aqueous solution plays a key role in determining the binding selectivity of DB18C6. In particular, the non-electrostatic dispersion interaction between the solute and solvent, which depends strongly on the complex structure, is largely responsible for the different solvation energies of M⁺-DB18C6. This study shows that the implicit solvation model like the CPCM works reasonably well in predicting the binding selectivity of DB18C6 in aqueous solution.

Conformational preferences and pK(a) value of selenocysteine residue

The conformational preferences of the L-selenocysteine (Sec) dipeptides with selenol and selenolate groups (Ac-Sec-NHMe and Ac-Sec(-) -NHMe, respectively) and the apparent (i.e., macroscopic) pK(a) value of the Sec residue have been studied using the dispersion-corrected density functionals M06-2X and B2PLYP-D with the implicit solvation method in the gas phase and in water. In the gas phase, the backbone-to-backbone and/or side chain-to-backbone hydrogen bonds are found to contribute in stabilizing the most preferred conformations for the Sec and Sec(-) residues, as seen for the Cys and Cys(-) residues. However, the polyproline II-like conformations prevail over the conformations with the backbone-to-backbone hydrogen bonds in water because of the weakened hydrogen bonds by the favorable direct interactions between the backbone CO and HN groups and water molecules. The Sec and Sec(-) residues are found to adopt more various conformations than the Cys and Cys(-) residues in water, although the most preferred conformations of the neutral and/or anionic forms of the two residues are similar each other in the gas phase and in water. Using the statistically weighted free energies of the Sec and Sec(-) dipeptides in the gas phase and their solvation free energies, the pK(a) value of the Sec residue is estimated to be 5.47 at 25°C, which is in good agreement with the experimental value of 5.43 ± 0.02. It is found that the lower pK(a) value of the selenol side chain for the Sec residue by ∼3 units than the thiol side chain for the Cys residue is ascribed to the higher gas-phase acidity of the Sec residue.

Molecular acidity: A quantitative conceptual density functional theory description

Accurate predictions of molecular acidity using ab initio and density functional approaches are still a daunting task. Using electronic and reactivity properties, one can quantitatively estimate pKa values of acids. In a recent paper [S. B. Liu and L. G. Pedersen, J. Phys. Chem. A 113, 3648 (2009)], we employed the molecular electrostatic potential (MEP) on the nucleus and the sum of valence natural atomic orbital (NAO) energies for the purpose. In this work, we reformulate these relationships on the basis of conceptual density functional theory and compare the results with those from the thermodynamic cycle method. We show that MEP and NAO properties of the dissociating proton of an acid should satisfy the same relationships with experimental pKa data. We employ 27 main groups and first to third row transition metal-water complexes as illustrative examples to numerically verify the validity of these strong linear correlations. Results also show that the accuracy of our approach and that of the conventional method through the thermodynamic cycle are statistically similar.

Conformational preferences and cis-trans isomerization of L-3,4-dehydroproline residue

The conformational study of N-acetyl-N'-methylamide of L-3,4-dehydroproline (Ac-Dhp-NHMe, the Dhp dipeptide) is carried out using hybrid density functional methods with the self-consistent reaction field method in the gas phase and in solution (chloroform and water). The incorporation of a double bond between C(beta) and C(gamma) into the prolyl ring results in the puckering, backbone population, and barriers to prolyl cis-trans isomerization different from those of the Pro dipeptide. For local minima of the Dhp dipeptide in the gas phase and in water, the C(beta)-C(gamma) bonds become shorter by approximately 0.2 A and the bond angles C(alpha)-C(beta)-C(gamma) and C(beta)-C(gamma)-C(delta) are widened by approximately 8 degrees than those of the Pro dipeptide, and the puckering amplitude is computed to be 0.01-0.07 A, indicating that the 3,4-dehydroprolyl ring is quite less puckered. However, polyproline-like conformations become preferred and the relative stability of the conformation tC with a C(7) intramolecular hydrogen bond decreases as the solvent polarity increases, as found for the Pro dipeptide. The barriers to cis-trans isomerization of the Ac-Dhp peptide bond increase with the increase of solvent polarity and the isomerization is likely to proceed through the clockwise rotation in water, as found for the prolyl peptide bond. The hydrogen bond between the prolyl nitrogen and the following amide N-H group seems to contribute in stabilizing the transition state structures.