Transferable scoring function based on semiempirical quantum mechanical PM6-DH2 method: CDK2 with 15 structurally diverse inhibitors.

A semiempirical quantum mechanical PM6-DH2 method accurately covering the dispersion interaction and H-bonding was used to score fifteen structurally diverse CDK2 inhibitors. The geometries of all the complexes were taken from the X-ray structures and were reoptimised by the PM6-DH2 method in continuum water. The total scoring function was constructed as an estimate of the binding free energy, i.

Performance of Molecular Mechanics Force Fields for RNA Simulations: Stability of UUCG and GNRA Hairpins.

The RNA hairpin loops represent important RNA topologies with indispensable biological functions in RNA folding and tertiary interactions. 5'-UNCG-3' and 5'-GNRA-3' RNA tetraloops are the most important classes of RNA hairpin loops. Both tetraloops are highly structured with characteristic signature three-dimensional features and are recurrently seen in functional RNAs and ribonucleoprotein particles.

Graphene fluoride: a stable stoichiometric graphene derivative and its chemical conversion to graphene.

Stoichoimetric graphene fluoride monolayers are obtained in a single step by the liquid-phase exfoliation of graphite fluoride with sulfolane. Comparative quantum-mechanical calculations reveal that graphene fluoride is the most thermodynamically stable of five studied hypothetical graphene derivatives; graphane, graphene fluoride, bromide, chloride, and iodide. The graphene fluoride is transformed into graphene via graphene iodide, a spontaneously decomposing intermediate.

Explicit Water Models Affect the Specific Solvation and Dynamics of Unfolded Peptides While the Conformational Behavior and Flexibility of Folded Peptides Remain Intact.

Conventional molecular dynamics simulations on 50 ns to 1 μs time scales were used to study the effects of explicit solvent models on the conformational behavior and solvation of two oligopeptide solutes: α-helical EK-peptide (14 amino acids) and a β-hairpin chignolin (10 amino acids). The widely used AMBER force fields (ff99, ff99SB, and ff03) were combined with four of the most commonly used explicit solvent models (TIP3P, TIP4P, TIP5P, and SPC/E). Significant differences in the specific solvation of chignolin among the studied water models were identified.

Flexibility of human cytochrome P450 enzymes: molecular dynamics and spectroscopy reveal important function-related variations.

To gain more complete insight into flexibility and malleability of five forms of human liver cytochrome P450 enzymes, which play major roles in drug metabolism (CYPs 1A2, 2A6, 2C9, 2D6 and 3A4), we employed UV/VIS and resonance Raman spectroscopy in combination with all-atomic molecular dynamics simulations under normal and high pressure conditions (300 MPa). In general, the high pressure reduces the flexibility of CYPs, which become more dense and compact as their radii of gyration and temperature B-factors diminish. The flexibility of CYPs spans the regions, which are localized in solvent exposed loops.

Large-scale compensation of errors in pairwise-additive empirical force fields: comparison of AMBER intermolecular terms with rigorous DFT-SAPT calculations.

The intermolecular interaction energy components for several molecular complexes were calculated using force fields available in the AMBER suite of programs and compared with Density Functional Theory-Symmetry Adapted Perturbation Theory (DFT-SAPT) values. The extent to which such comparison is meaningful is discussed. The comparability is shown to depend strongly on the intermolecular distance, which means that comparisons made at one distance only are of limited value.

Protonation states of the key active site residues and structural dynamics of the glmS riboswitch as revealed by molecular dynamics.

The glmS catalytic riboswitch is part of the 5'-untranslated region of mRNAs encoding glucosamine-6-phosphate (GlcN6P) synthetase (glmS) in numerous gram-positive bacteria. Binding of the cofactor GlcN6P induces site-specific self-cleavage of the RNA. However, the detailed reaction mechanism as well as the protonation state of the glmS reactive form still remains elusive.

Extensive molecular dynamics simulations showing that canonical G8 and protonated A38H+ forms are most consistent with crystal structures of hairpin ribozyme.

The hairpin ribozyme is a prominent member of the group of small catalytic RNAs (RNA enzymes or ribozymes) because it does not require metal ions to achieve catalysis. Biochemical and structural data have implicated guanine 8 (G8) and adenine 38 (A38) as catalytic participants in cleavage and ligation catalyzed by the hairpin ribozyme, yet their exact role in catalysis remains disputed. To gain insight into dynamics in the active site of a minimal self-cleaving hairpin ribozyme, we have performed extensive classical, explicit-solvent molecular dynamics (MD) simulations on time scales of 50-150 ns.

A cancer-derived mutation in the PSTAIRE helix of cyclin-dependent kinase 2 alters the stability of cyclin binding.

Cyclin-dependent kinase 2 (cdk2) is a central regulator of the mammalian cell cycle. Here we describe the properties of a mutant form of cdk2 identified during large-scale sequencing of protein kinases from cancerous tissue. The mutation substituted a leucine for a proline in the PSTAIRE helix, the central motif in the interaction of the cdk with its regulatory cyclin subunit.

Dependence of A-RNA simulations on the choice of the force field and salt strength.

We present an extensive molecular dynamics study (0.6 micros in total) on three A-RNA duplexes. The dependence of the A-RNA geometry on force fields (Parm99 and Parmbsc0) and salt strength conditions (approximately 0.

Theoretical studies on the intermolecular interactions of potentially primordial base-pair analogues.

Recent experimental studies on the Watson-Crick type base pairing of triazine and aminopyrimidine derivatives suggest that acid/base properties of the constituent bases might be related to the duplex stabilities measured in solution. Herein we use high-level quantum chemical calculations and molecular dynamics simulations to evaluate the base pairing and stacking interactions of seven selected base pairs, which are common in that they are stabilized by two N-H..

Glycine-rich loop of mitochondrial processing peptidase alpha-subunit is responsible for substrate recognition by a mechanism analogous to mitochondrial receptor Tom20.

Tryptophan fluorescence measurements were used to characterize the local dynamics of the highly conserved glycine-rich loop (GRL) of the mitochondrial processing peptidase (MPP) alpha-subunit in the presence of the substrate precursor. Reporter tryptophan residue was introduced into the GRL of the yeast alpha-MPP (Y299W) or at a proximal site (Y303W). Time-resolved and steady-state fluorescence spectroscopy demonstrated that for Trp299, the primary contact with the yeast malate dehydrogenase precursor evokes a change of the local GRL mobility.

In silico structural and functional analysis of fragments of the ankyrin repeat protein p18(INK4c).

Ankyrin repeat proteins (ARPs) are ubiquitous proteins that play critical regulatory roles in organisms and consist of repeating motifs (ankyrin repeats) stacked in non-globular, almost linear, "quasi one-dimensional" configurations. They also have highly unusual mechanical properties, notably ARPs can behave as nano-springs. Both their essential cellular functions and distinctive nano-mechanical properties have aroused interest in ARPs for potential applications in medicine and nanotechnology.

Molecular dynamics and quantum mechanics of RNA: conformational and chemical change we can believe in.

Structure and dynamics are both critical to RNA's vital functions in biology. Numerous techniques can elucidate the structural dynamics of RNA, but computational approaches based on experimental data arguably hold the promise of providing the most detail. In this Account, we highlight areas wherein molecular dynamics (MD) and quantum mechanical (QM) techniques are applied to RNA, particularly in relation to complementary experimental studies.

Changes in electrostatic surface potential of Na+/K+-ATPase cytoplasmic headpiece induced by cytoplasmic ligand(s) binding.

A set of single-tryptophan mutants of the Na(+)/K(+)-ATPase isolated, large cytoplasmic loop connecting transmembrane helices M4 and M5 (C45) was prepared to monitor effects of the natural cytoplasmic ligands (i.e., Mg(2+) and/or ATP) binding.

Redesigning dehalogenase access tunnels as a strategy for degrading an anthropogenic substrate.

Engineering enzymes to degrade anthropogenic compounds efficiently is challenging. We obtained Rhodococcus rhodochrous haloalkane dehalogenase mutants with up to 32-fold higher activity than wild type toward the toxic, recalcitrant anthropogenic compound 1,2,3-trichloropropane (TCP) using a new strategy. We identified key residues in access tunnels connecting the buried active site with bulk solvent by rational design and randomized them by directed evolution.

Theoretical studies of RNA catalysis: hybrid QM/MM methods and their comparison with MD and QM.

Hybrid QM/MM methods combine the rigor of quantum mechanical (QM) calculations with the low computational cost of empirical molecular mechanical (MM) treatment allowing to capture dynamic properties to probe critical atomistic details of enzyme reactions. Catalysis by RNA enzymes (ribozymes) has only recently begun to be addressed with QM/MM approaches and is thus still a field under development. This review surveys methodology as well as recent advances in QM/MM applications to RNA mechanisms, including those of the HDV, hairpin, and hammerhead ribozymes, as well as the ribosome.

ATP and magnesium drive conformational changes of the Na+/K+-ATPase cytoplasmic headpiece.

Conformational changes of the Na(+)/K(+)-ATPase isolated large cytoplasmic segment connecting transmembrane helices M4 and M5 (C45) induced by the interaction with enzyme ligands (i.e. Mg(2+) and/or ATP) were investigated by means of the intrinsic tryptophan fluorescence measurement and molecular dynamic simulations.

Flexibility of human cytochromes P450: molecular dynamics reveals differences between CYPs 3A4, 2C9, and 2A6, which correlate with their substrate preferences.

Molecular dynamics (MD) simulations at normal and high temperature were used to study the flexibility and malleability of three microsomal cytochromes P450 (CYPs): CYP3A4, CYP2C9, and CYP2A6. Comparison of B-factors (describing the atomic fluctuations) between X-ray and MD data shows that the X-ray B-factors are significantly lower in the regions where the crystal contacts occur than for other regions. Consequently, the conclusions about CYP flexibility based solely on the X-ray data might be misleading.

Two C-terminal ankyrin repeats form the minimal stable unit of the ankyrin repeat protein p18INK4c.

Ankyrin repeat proteins (ARPs) appear to be abundant in organisms from all phyla, and play critical regulatory roles, mediating specific interactions with target biomolecules and thus ordering the sequence of events in diverse cellular processes. ARPs possess a non-globular scaffold consisting of repeating motifs named ankyrin (ANK) repeats, which stack on each other. The modular architecture of ARPs provides a new paradigm for understanding protein stability and folding mechanisms.

Regulatory phosphorylation of cyclin-dependent kinase 2: insights from molecular dynamics simulations.

The structures of fully active cyclin-dependent kinase-2 (CDK2) complexed with ATP and peptide substrate, CDK2 after the catalytic reaction, and CDK2 inhibited by phosphorylation at Thr14/Tyr15 were studied using molecular dynamics (MD) simulations. The structural details of the CDK2 catalytic site and CDK2 substrate binding box were described. Comparison of MD simulations of inhibited complexes of CDK2 was used to help understand the role of inhibitory phosphorylation at Thr14/Tyr15.

Functional flexibility of human cyclin-dependent kinase-2 and its evolutionary conservation.

Cyclin-dependent kinase 2 (CDK2) is the most thoroughly studied of the cyclin-dependent kinases that regulate essential cellular processes, including the cell cycle, and it has become a model for studies of regulatory mechanisms at the molecular level. This contribution identifies flexible and rigid regions of CDK2 based on temperature B-factors acquired from both X-ray data and molecular dynamics simulations. In addition, the biological relevance of the identified flexible regions and their motions is explored using information from the essential dynamics analysis related to conformational changes of CDK2 and knowledge of its biological function(s).

MOLE: a Voronoi diagram-based explorer of molecular channels, pores, and tunnels.

We have developed an algorithm, "MOLE," for the rapid, fully automated location and characterization of molecular channels, tunnels, and pores. This algorithm has been made freely available on the Internet (http://mole.chemi.