Molecular dynamics of DNA quadruplex molecules containing inosine, 6-thioguanine and 6-thiopurine

The ability of the four-stranded guanine (G)-DNA motif to incorporate nonstandard guanine analogue bases 6-oxopurine (inosine, I), 6-thioguanine (tG), and 6-thiopurine (tI) has been investigated using large-scale molecular dynamics simulations. The simulations suggest that a G-DNA stem can incorporate inosines without any marked effect on its structure and dynamics. The all-inosine quadruplex stem d(IIII)(4) shows identical dynamical properties as d(GGGG)(4) on the nanosecond time scale, with both molecular assemblies being stabilized by monovalent cations residing in the channel of the stem. However, simulations carried out in the absence of these cations show dramatic differences in the behavior of d(GGGG)(4) and d(IIII)(4). Whereas vacant d(GGGG)(4) shows large fluctuations but does not disintegrate, vacant d(IIII)(4) is completely disrupted within the first nanosecond. This is a consequence of the lack of the H-bonds involving the N2 amino group that is not present in inosine. This indicates that formation of the inosine quadruplex could involve entirely different intermediate structures than formation of the guanosine quadruplex, and early association of cations in this process appears to be inevitable. In the simulations, the incorporation of 6-thioguanine and 6-thiopurine sharply destabilizes four-stranded G-DNA structures, in close agreement with experimental data. The main reason is the size of the thiogroup leading to considerable steric conflicts and expelling the cations out of the channel of the quadruplex stem. The G-DNA stem can accommodate a single thioguanine base with minor perturbations. Incorporation of a thioguanine quartet layer is associated with a large destabilization of the G-DNA stem whereas the all-thioguanine quadruplex immediately collapses.     

Activation and inhibition of cyclin-dependent kinase-2 by phosphorylation; a molecular dynamics study reveals the functional importance of the glycine-rich loop

Nanoseconds long molecular dynamics (MD) trajectories of differently active complexes of human cyclin-dependent kinase 2 (inactive CDK2/ATP, semiactive CDK2/Cyclin A/ATP, fully active pT160-CDK2/Cyclin A/ATP, inhibited pT14-; pY15-; and pT14,pY15,pT160-CDK2/Cyclin A/ATP) were compared. The MD simulations results of CDK2 inhibition by phosphorylation at T14 and/or Y15 sites provide insight into the structural aspects of CDK2 deactivation. The inhibitory sites are localized in the glycine-rich loop (G-loop) positioned opposite the activation T-loop. Phosphorylation of T14 and both inhibitory sites T14 and Y15 together causes ATP misalignment for phosphorylation and G-loop conformational change. This conformational change leads to the opening of the CDK2 substrate binding box. The phosphorylated Y15 residue negatively affects substrate binding or its correct alignment for ATP terminal phospho-group transfer to the CDK2 substrate. The MD simulations of the CDK2 activation process provide results in agreement with previous X-ray data.     

Dynamics and binding modes of free cdk2 and its two complexes with inhibitors studied by computer simulations

This article presents a molecular dynamics (MD) study of the cdk2 enzyme and its two complexes with the inhibitors isopentenyladenine and roscovitine using the Cornell et al. force field from the AMBER software package. The results show that inserting an inhibitor into the enzyme active site does not considerably change enzyme structure but it seemingly changes the distribution of internal motions. The inhibitor causes differences in the domain motions in free cdk2 and in its complexes. It was found out that repulsion of roscovitine N9 substituent causes conformational change on Lys 33 side chain. Isopentenyladenine forms with Lys 33 side chain terminal amino group a hydrogen bond. It implies that the cavity, where N9 substituent of roscovitine is buried, can adopt larger substituent due to Lys 33 side chain flexibility. The composition of electrostatic and van der Waals interactions between the inhibitor and the enzyme were also calculated along both cdk2/inhibitor MD trajectories together with MM-PB/GBSA analysis. These results show that isopentenyladenine-like inhibitors could be more effective after modifications leading to an increase in their van der Waals contact with the enzyme. We suggest that a way leading to better inhibitors occupying isopentenyladenine binding mode could be: to keep N9 and N7 purine positions free, to keep 3,3-dimethylallylamino group at C6 position, and to add, e.g., benzylamino group at C2 position. The results support the idea that the isopentenyladenine binding mode can be used for cdk2 inhibitors design and that all possibilities to improve this binding mode were not uncovered yet.     

Conformational flexibility of two RNA trimers explored by computational tools and database search

Two RNA sequences, AAA and AUG, were studied by the conformational search program CICADA and by molecular dynamics (MD) in the framework of the AMBER force field, and also via thorough PDB database search. CICADA was used to provide detailed information about conformers and conformational interconversions on the energy surfaces of the above molecules. Several conformational families were found for both sequences. Analysis of the results shows differences, especially between the energy of the single families, and also in flexibility and concerted conformational movement. Therefore, several MD trajectories (altogether 16 ns) were run to obtain more details about both the stability of conformers belonging to different conformational families and about the dynamics of the two systems. Results show that the trajectories strongly depend on the starting structure. When the MD start from the global minimum found by CICADA, they provide a stable run, while MD starting from another conformational family generates a trajectory where several different conformational families are visited. The results obtained by theoretical methods are compared with the thorough database search data. It is concluded that all except for the highest energy conformational families found in theoretical result also appear in experimental data. Registry numbers: adenylyl-(3' --> 5')-adenylyl-(3' --> 5')-adenosine [917-44-2] adenylyl-(3' --> 5')-uridylyl-(3' --> 5')-guanosine [3494-35-7].     

Lipid peroxidation and scavenging enzyme levels in the liver of streptozotocin-induced diabetic rats

In this study, alterations in the liver antioxidant enzymes status and lipid peroxidation in short-term (8-weeks) and long-term (24-weeks) diabetic rats were examined. Glutathione peroxidase (GSH-Px) activity and malondialdehyde (MDA) levels were significantly increased, but superoxide dismutase (SOD) activity was significantly reduced in 8-weeks diabetic rats, compared to control. Catalase (CAT) activity, however, was found unchanged. In 24-weeks diabetic rats, while GSH-Px activity was unchanged, but SOD and CAT activities and MDA levels were significantly increased, compared to control. These results suggest that diabetes-induced alterations in tissue antioxidant system may reflect a generalized increase in tissue oxidative stress. It can be concluded that lipid peroxidation and antioxidant enzyme levels are elevated in diabetic condition. Hence, diabetes mellitus, if left untreated, may increase degenerative processes due to accumulation of oxidative free radicals.     

TRITON: in silico construction of protein mutants and prediction of their activities

MOTIVATION: One of the objectives of protein engineering is to propose and construct modified proteins with improved activity for the substrate of interest. Systematic computational investigation of many protein variants requires the preparation and handling of a large number of data files. The type of the data generated during the modelling of protein variants and the estimation of their activities offers the possibility of process automatization. RESULTS: The graphical program TRITON has been developed for modelling protein mutants and assessment of their activities. Protein mutants are modelled from the wild type structure by homology modelling using the external program MODELLER. Chemical reactions taking place in the mutants active site are modelled using the semi-empirical quantum mechanic program MOPAC. Semi-quantitative predictions of mutants activities can be achieved by evaluating the changes in energies of the system and partial atomic charges of active site residues during the reaction. The program TRITON offers graphical tools for the preparation of the input data files, for calculation and for the analysis of the generated output data. AVAILABILITY: The program TRITON can run under operating systems IRIX, Linux and NetBSD. The software is available at http://www.chemi.muni.cz/lbsd/triton.ht ml.     

The effect of salt stress on lipid peroxidation, antioxidative enzymes and proline content of sesame cultivars

The effect of increasing NaCl concentrations was studied on two different cultivars (cv. Orhangazi and cv. Cumhuriyet) of Sesamum indicum. Seedlings were grown for 40 days in half strength Hoagland solution and after 40 days treated with different NaCl concentrations (0, 50 and 100 mM) for 21 days. Differences in growth parameters, lipid peroxidation, antioxidative enzyme activities and proline accumulation were tested in order to put forward the relative tolerance or sensitivity of the cultivars. Results indicated that both parameters differ according to the cultivar's ability in coping oxidative stress caused by salinity. Constitutive levels of antioxidative enzyme activities were almost the same between the cultivars; however, cv. Cumhuriyet was able to induce antioxidative enzyme activities more efficiently when subjected to salt stress. Growth parameters, lipid peroxidation and proline accumulation results are also in good correlation with supporting this cultivar's being relatively tolerant.     

Morphology, and muscle- and papilla-volume ratios, of the tongue of Laudakia stellio (Agamidae, Squamata): a histological and stereological study

We examined the histological structure of the tongue of Laudakia stellio, the starred agama lizard (Agamidae, Squamata), under light microscopy. We also investigated the muscle and papilla volume ratios, with volumes of each aspect of interest estimated according to the Cavalieri method. The macroscopically short, thick and muscle-rich front tip of the tongue of L. stellio does not show any bifurcation, and under light microscopy, the oval-shaped papilla-free front tip was seen to be covered by keratinized stratified epithelium. The dorsal and ventral parts were different, with the former partially covered by keratinized stratified epithelium and rich in secretory glands and secretory cells. The ventral part, which contained keratinized stratified cells, had a flat surface with no papillae. The dorsal surface of the anterior and posterior parts contained fungiform papillae, with the apical parts of these papillae containing minimal keratin; the interpapillar space was covered by keratin-free squamous stratified epithelium. The middle section of the tongue contained cylindrical-type papillae, with serous and mucous secretory glands and ducts at their base. Finally, the frontal and middle parts of the ventral and dorsal surfaces did not contain any taste buds, although there were some in the hind part of the dorsal surface. As morphometric estimates of volumes of the muscles and papillae, the mean volume ratios (relative to total tongue volume)+/-standard deviation were 0.66+/-0.03 and 0.33+/-0.03, with mean coefficients of error of 0.02 and 0.03, respectively.