Biomedical subjects

Jaroslav Koca

Publications and source records attributed to Jaroslav Koca.

At least 19 recordsLinked to original sources

PubMed2006

Structure, dynamics, and elasticity of free 16s rRNA helix 44 studied by molecular dynamics simulations.

Molecular dynamics (MD) simulations were employed to investigate the structure, dynamics, and local base-pair step deformability of the free 16S ribosomal helix 44 from Thermus thermophilus and of a canonical A-RNA double helix. While helix 44 is bent in the crystal structure of the small ribosomal subunit, the simulated helix 44 is intrinsically straight. It shows, however, substantial instantaneous bends that are isotropic. The spontaneous motions seen in simulations achieve large degrees of bending seen in the X-ray structure and would be entirely sufficient to allow the dynamics of the upper part of helix 44 evidenced by cryo-electron microscopic studies. Analysis of local base-pair step deformability reveals a patch of flexible steps in the upper part of helix 44 and in the area proximal to the bulge bases, suggesting that the upper part of helix 44 has enhanced flexibility. The simulations identify two conformational substates of the second bulge area (bottom part of the helix) with distinct base pairing. In agreement with nuclear magnetic resonance (NMR) and X-ray studies, a flipped out conformational substate of conserved 1492A is seen in the first bulge area. Molecular dynamics (MD) simulations reveal a number of reversible alpha-gamma backbone flips that correspond to transitions between two known A-RNA backbone families. The flipped substates do not cumulate along the trajectory and lead to a modest transient reduction of helical twist with no significant influence on the overall geometry of the duplexes. Despite their considerable flexibility, the simulated structures are very stable with no indication of substantial force field inaccuracies.

Base Pairing↗

PubMed2006

CAVER: a new tool to explore routes from protein clefts, pockets and cavities.

BACKGROUND: The main aim of this study was to develop and implement an algorithm for the rapid, accurate and automated identification of paths leading from buried protein clefts, pockets and cavities in dynamic and static protein structures to the outside solvent. RESULTS: The algorithm to perform a skeleton search was based on a reciprocal distance function grid that was developed and implemented for the CAVER program. The program identifies and visualizes routes from the interior of the protein to the bulk solvent. CAVER was primarily developed for proteins, but the algorithm is sufficiently robust to allow the analysis of any molecular system, including nucleic acids or inorganic material. Calculations can be performed using discrete structures from crystallographic analysis and NMR experiments as well as with trajectories from molecular dynamics simulations. The fully functional program is available as a stand-alone version and as plug-in for the molecular modeling program PyMol. Additionally, selected functions are accessible in an online version. CONCLUSION: The algorithm developed automatically finds the path from a starting point located within the interior of a protein. The algorithm is sufficiently rapid and robust to enable routine analysis of molecular dynamics trajectories containing thousands of snapshots. The algorithm is based on reciprocal metrics and provides an easy method to find a centerline, i.e. the spine, of complicated objects such as a protein tunnel. It can also be applied to many other molecules. CAVER is freely available from the web site http://loschmidt.chemi.muni.cz/caver/.

Algorithms↗

PubMed2006

Effect of solvation on pinched cone-pinched cone interconversion of tetraethoxycalix[4]arene and tetraethoxythiacalix[4]arene.

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.

Calixarenes↗

PubMed2006

Different mechanisms of CDK5 and CDK2 activation as revealed by CDK5/p25 and CDK2/cyclin A dynamics.

A detailed analysis is presented of the dynamics of human CDK5 in complexes with the protein activator p25 and the purine-like inhibitor roscovitine. These and other findings related to the activation of CDK5 are critically reviewed from a molecular perspective. In addition, the results obtained on the behavior of CDK5 are compared with data on CDK2 to assess the differences and similarities between the two kinases in terms of (i) roscovitine binding, (ii) regulatory subunit association, (iii) conformational changes in the T-loop following CDK/regulatory subunit complex formation, and (iv) specificity in CDK/regulatory subunit recognition. An energy decomposition analysis, used for these purposes, revealed why the binding of p25 alone is sufficient to stabilize the extended active T-loop conformation of CDK5, whereas the equivalent conformational change in CDK2 requires both the binding of cyclin A and phosphorylation of the Thr(160) residue. The interaction energy of the CDK5 T-loop with p25 is about 26 kcal.mol(-1) greater than that of the CDK2 T-loop with cyclin A. The binding pattern between CDK5 and p25 was compared with that of CDK2/cyclin A to find specific regions involved in CDK/regulatory subunit recognition. The analyses performed revealed that the alphaNT-helix of cyclin A interacts with the alpha6-alpha7 loop and the alpha7 helix of CDK2, but these regions do not interact in the CDK5/p25 complex. Further differences between the CDK5/p25 and CDK2/cyclin A systems studied are discussed with respect to their specific functionality.

Cell Cycle Proteins↗

PubMed2006

RNA kink-turns as molecular elbows: hydration, cation binding, and large-scale dynamics.

The presence of Kink-turns (Kt) at key functional sites in the ribosome (e.g., A-site finger and L7/L12 stalk) suggests that some Kink-turns can confer flexibility on RNA protuberances that regulate the traversal of tRNAs during translocation. Explicit solvent molecular dynamics demonstrates that Kink-turns can act as flexible molecular elbows. Kink-turns are associated with a unique network of long-residency static and dynamical hydration sites that is intimately involved in modulating their conformational dynamics. An implicit solvent conformational search confirms the flexibility of Kink-turns around their X-ray geometries and identifies a second low-energy region with open structures that could correspond to Kink-turn geometries seen in solution experiments. An extended simulation of Kt-42 with the factor binding site (helices 43 and 44) shows that the local Kt-42 elbow-like motion fully propagates beyond the Kink-turn, and that there is no other comparably flexible site in this rRNA region. Kink-turns could mediate large-scale adjustments of distant RNA segments.

Binding Sites↗

PubMed2006

Molecular modeling of glycosyltransferases.

Glycosyltransferases, the enzymes that build oligosaccharides and glycoconjugates, have received much interest in recent years owing to their biological functions and their potential uses in biotechnology. The analysis of the wealth of sequences that are now available in databases allowed the classification in different families characterized by conserved peptide motifs. Nevertheless, only a limited number of crystal structures is available and molecular modeling appears to be an inescapable tool for rationalization of binding data, engineering of enzyme properties, and design of inhibitors that would be of interest as therapeutic compounds. Because of sequence diversity and limited experimental data, molecular modeling of these enzymes is not straightforward and utilizes the most recent tools, such as fold recognition programs.

Amino Acid Sequence↗

PubMed2005

Structures and mechanisms of glycosyltransferases.

Glycosyltransferases (GTs) catalyze the transfer of a sugar moiety from an activated donor sugar onto saccharide and nonsaccharide acceptors. A sequence-based classification spreads GTs in many families thus reflecting the variety of molecules that can be used as acceptors. In contrast, this enzyme family is characterized by a more conserved three-dimensional architecture. Until recently, only two different folds (GT-A and GT-B) have been identified for solved crystal structures. The recent report of a structure for a bacterial sialyltransferase allows the definition of a new fold family. Progress in the elucidation of the structures and mechanisms of GTs are discussed in this review. To accommodate the growing number of crystal structures, we created the 3D-Glycosyltransferase database to gather structural information concerning this class of enzymes.

Animals↗

PubMed2005

Hinge-like motions in RNA kink-turns: the role of the second a-minor motif and nominally unpaired bases.

Kink-turn (K-turn) motifs are asymmetric internal loops found at conserved positions in diverse RNAs, with sharp bends in phosphodiester backbones producing V-shaped structures. Explicit-solvent molecular dynamics simulations were carried out for three K-turns from 23S rRNA, i.e., Kt-38 located at the base of the A-site finger, Kt-42 located at the base of the L7/L12 stalk, and Kt-58 located in domain III, and for the K-turn of human U4 snRNA. The simulations reveal hinge-like K-turn motions on the nanosecond timescale. The first conserved A-minor interaction between the K-turn stems is entirely stable in all simulations. The angle between the helical arms of Kt-38 and Kt-42 is regulated by local variations of the second A-minor (type I) interaction between the stems. Its variability ranges from closed geometries to open ones stabilized by insertion of long-residency waters between adenine and cytosine. The simulated A-minor geometries fully agree with x-ray data. Kt-58 and Kt-U4 exhibit similar elbow-like motions caused by conformational change of the adenosine from the nominally unpaired region. Despite the observed substantial dynamics of K-turns, key tertiary interactions are stable and no sign of unfolding is seen. We suggest that some K-turns are flexible elements mediating large-scale ribosomal motions during the protein synthesis cycle.

Adenine↗

PubMed2005

High affinity fucose binding of Pseudomonas aeruginosa lectin PA-IIL: 1.0 A resolution crystal structure of the complex combined with thermodynamics and computational chemistry approaches.

PA-IIL is a fucose-binding lectin from Pseudomonas aeruginosa that is closely related to the virulence factors of the bacterium. Previous structural studies have revealed a new carbohydrate-binding mode with direct involvement of two calcium ions (Mitchell E, Houles C, Sudakevitz D, Wimmerova M, Gautier C, Perez S, Wu AM, Gilboa-Garber N, Imberty A. Structural basis for selective recognition of oligosaccharides from cystic fibrosis patients by the lectin PA-IIL of Pseudomonas aeruginosa. Nat Struct Biol 2002;9:918-921). A combination of thermodynamic, structural, and computational methods has been used to study the basis of the high affinity for the monosaccharide ligand. A titration microcalorimetry study indicated that the high affinity is enthalpy driven. The crystal structure of the tetrameric PA-IIL in complex with fucose and calcium was refined to 1.0 A resolution and, in combination with modeling, allowed a proposal to be made for the hydrogen-bond network in the binding site. Calculations of partial charges using ab initio computational chemistry methods indicated that extensive delocalization of charges between the calcium ions, the side chains of the protein-binding site and the carbohydrate ligand is responsible for the high enthalpy of binding and therefore for the unusually high affinity observed for this unique mode of carbohydrate recognition.

Adhesins, Bacterial↗

PubMed2005

The mechanism of inhibition of the cyclin-dependent kinase-2 as revealed by the molecular dynamics study on the complex CDK2 with the peptide substrate HHASPRK.

Molecular dynamics (MD) simulations were used to explain structural details of cyclin-dependent kinase-2 (CDK2) inhibition by phosphorylation at T14 and/or Y15 located in the glycine-rich loop (G-loop). Ten-nanosecond-long simulations of fully active CDK2 in a complex with a short peptide (HHASPRK) substrate and of CDK2 inhibited by phosphorylation of T14 and/or Y15 were produced. The inhibitory phosphorylations at T14 and/or Y15 show namely an ATP misalignment and a G-loop shift (~5 A) causing the opening of the substrate binding box. The biological functions of the G-loop and GxGxxG motif evolutionary conservation in protein kinases are discussed. The position of the ATP gamma-phosphate relative to the phosphorylation site (S/T) of the peptide substrate in the active CDK2 is described and compared with inhibited forms of CDK2. The MD results clearly provide an explanation previously not known as to why a basic residue (R/K) is preferred at the P(2) position in phosphorylated S/T peptide substrates.

Adenosine Triphosphate↗

PubMed2005

Long-range electrostatic interactions in molecular dynamics: an endothelin-1 case study.

An extensive conformational search in explicit solvent was performed in order to compare the influence of different long-range electrostatic interaction treatments in molecular dynamics. The short peptide endothelin-1 was selected as the subject of molecular dynamics studies that started from both X-ray and NMR obtained structures. Electrostatic interactions were treated using two of the most common methods--residue-based cutoff and particle mesh Ewald (PME). Analyses of free energy calculations (MM-PBSA method used), secondary structure elements and hydrogen bonds were performed, and there suggested that there is no unambiguous conclusion about which of the two methods of long-range electrostatics treatment should be used in MD simulations in this case. The most reliable data was provided by a trajectory that started with the NMR structure and used the cutoff method to treat electrostatic interactions. This leads to a recommendation that the choice of electrostatics treatment should be made carefully and not automatically by choosing the PME method simply because it is the most widely used.

Endothelin-1↗

PubMed2004

Long-residency hydration, cation binding, and dynamics of loop E/helix IV rRNA-L25 protein complex.

Molecular dynamics simulations of RNA-protein complex between Escherichia coli loop E/helix IV (LE/HeIV) rRNA and L25 protein reveal a qualitative agreement between the experimental and simulated structures. The major groove of LE is a prominent rRNA cation-binding site. Divalent cations rigidify the LE major groove geometry whereas in the absence of divalent cations LE extensively interacts with monovalent cations via inner-shell binding. The HeIV region shows bistability of its major groove explaining the observed differences between x-ray and NMR structures. In agreement with the experiments, the simulations suggest that helix-alpha1 of L25 is the least stable part of the protein. Inclusion of Mg2+ cations into the simulations causes perturbation of basepairing at the LE/HeIV junction, which does not, however, affect the protein binding. The rRNA-protein complex is mediated by a number of highly specific hydration sites with long-residing water molecules and two of them are bound throughout the entire 24-ns simulation. Long-residing water molecules are seen also outside the RNA-protein contact areas with water-binding times substantially enhanced compared to simulations of free RNA. Long-residency hydration sites thus represent important elements of the three-dimensional structure of rRNA.

Binding Sites↗

PubMed2004

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.

CDC2-CDC28 Kinases↗

PubMed2004

Analysis of CDK2 active-site hydration: a method to design new inhibitors.

The interactions between the protein and the solvent were analyzed, and protein regions with a high density of water molecules, as well as structural water molecules, were determined by using molecular dynamics (MD) simulations. A number of water molecules that were in contact with the protein for the whole trajectory were determined. Their interaction energies and hydrogen bonds with protein residues were analyzed. Altogether, 39, 27, 49, and 32 water molecules bound to the protein were found for trajectories of the free CDK2, CDK2/ATP, CDK2/roscovitine, and CDK2/isopentenyladenine complexes, respectively. Positions of observed water molecules were compared with X-ray crystallography data. Special attention was paid to water molecules in the active site of the enzyme, and especially to the deep pocket, where the N9 roscovitine side-chain is buried. Exchange of active-site water molecules with bulk water through the tunnel from the pocket was observed. In the CDK2/isopentenyladenine complex simulation, two water molecules that arrange interaction between the inhibitor and the enzyme via an H-bond were observed. Two stable water molecules in the trajectory of the free CDK2 were found that occupy the same position as the nitrogens N3 and N9 of the isopentenyladenine or N1 and N6 nitrogens of the adenosine triphosphate (ATP). The positions of structural water molecules were compared with the positions of substrate polar groups and crystallographic water molecules found in the Brookhaven Protein Data Bank for various CDK2 complexes. It was concluded that tracing tightly bound water molecules may substantially help in designing new inhibitors.

Adenine↗

PubMed2004

Conformational analysis of complex oligosaccharides: the CICADA approach to the uromodulin O-glycans.

Uromodulin is the pregnancy-associated Tamm-Horsfall glycoprotein, with the enhanced ability to inhibit T-cell proliferation. Pregnancy-associated structural changes mainly occur in the O-glycosylation of this glycoprotein. These include up to 12 glycan structures, made up of an unusual core type 2 sequence terminated with one, two, or three sialyl Lewis(x) sequences; this type of O-glycans could serve as E- and P-selectin ligands. The present work focuses on the most complex one; a tetradecamer made up of a type 2 core carrying three sialyl Lewis(x) branches. Five different monosaccharides are assembled by 14 glycosidic linkages. The conformational behavior of the constituting disaccharide segments was evaluated using the flexible residue procedure of the MM3 molecular mechanics procedure. For each disaccharide, the adiabatic energy surface, along with the local energy minima were established. All these results were used for the generation, prior to complete optimization of the tetradecamer. This was followed by a complete exploration of conformational hyperspace throughout the use of the single coordinate method as implemented in the CICADA program. Despite the potential flexibility of the tetradecasaccharide, only four conformational families occur, accounting for more than 95% of the total low energy conformations. For each family, the molecular properties (electrostatic, lipophilicity, and hydrogen potential) were studied. The shape of the tetradecasaccharide is best described as a flat ribbon, flanked by three branches having terminal sialyl residues. Two of the branches interact through nonbonded interactions, bringing further energy stabilization, and limiting the conformational flexibility of the sialyl residues. Only one branch maintains the original conformational features of sialyl Lewis(x). This O-glycan can be seen as a fascinating example of 'dendrimeric' structure, where the spatial arrangement of three S-Le(x) epitopes may favor its complementary 'presentations' for the interactions with E- and P-selectins.

Carbohydrate Conformation↗

PubMed2004

Molecular dynamics simulations of glycosyltransferase LgtC.

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.

Amino Acid Sequence↗

PubMed2004

Enantioselective chromatography and molecular modeling of novel aryloxyaminopropan-2-ols with the alkyl carbamate function.

A series of different racemic aryloxyaminopropan-2-ol derivatives 1a-d-3a-d with potential beta-adrenergic blocking effects related to propanolol 4 and atenolol 5 was resolved by HPLC using Chiralcel OD-H and Chiralpak AD as chiral stationary phases. Mobile phases consisted of a hexane/alcohol (propan-2-ol or ethanol) mixture doped with a modifier (DEA or TFA). The retention behavior of the compounds depended on the position of the carbamate attached to the aryloxy moiety and on the length of the alkyl residue in the carbamate. Enantiomers of the title compounds were baseline separated with the separation factors alpha and resolutions R(s) varying in the range of 1.34-4.55 and 1.50-10.65, respectively. The chromatographic systems developed can be used for the determination of the enantiomeric purity of the title compounds. Molecular modelling using empirical molecular mechanics and ab initio quantum chemistry methods provided low-energy structures in which sites of potential interactions responsible for retention behavior and chiral recognition could be identified.

Journal Article↗

PubMed2004

Ribosomal RNA kink-turn motif--a flexible molecular hinge.

Ribosomal RNA K-turn motifs are asymmetric internal loops characterized by a sharp bend in the phosphodiester backbone resulting in "V" shaped structures, recurrently observed in ribosomes and showing a high degree of sequence conservation. We have carried out extended explicit solvent molecular dynamics simulations of selected K-turns, in order to investigate their intrinsic structural and dynamical properties. The simulations reveal an unprecedented dynamical flexibility of the K-turns around their X-ray geometries. The K-turns sample, on the nanosecond timescale, different conformational substates. The overall behavior of the simulations suggests that the sampled geometries are essentially isoenergetic and separated by minimal energy barriers. The nanosecond dynamics of isolated K-turns can be qualitatively considered as motion of two rigid helix stems controlled by a very flexible internal loop which then leads to substantial hinge-like motions between the two stems. This internal dynamics of K-turns is strikingly different for example from the bacterial 5S rRNA Loop E motif or BWYV frameshifting pseudoknot which appear to be rigid in the same type of simulations. Bistability and flexibility of K-turns was also suggested by several recent biochemical studies. Although the results of MD simulations should be considered as a qualitative picture of the K-turn dynamics due to force field and sampling limitations, the main advantage of the MD technique is its ability to investigate the region close to K-turn ribosomal-like geometries. This part of the conformational space is not well characterized by the solution experiments due to large-scale conformational changes seen in the experiments. We suggest that K-turns are well suited to act as flexible structural elements of ribosomal RNA. They can for example be involved in mediation of large-scale motions or they can allow a smooth assembling of the other parts of the ribosome.

Base Sequence↗