Rajmund Kaźmierkiewicz

Recent improvements in prediction of protein structure by global optimization of a potential energy function.

Recent improvements of a hierarchical ab initio or de novo approach for predicting both α and β structures of proteins are described. The united-residue energy function used in this procedure includes multibody interactions from a cumulant expansion of the free energy of polypeptide chains, with their relative weights determined by Z-score optimization. The critical initial stage of the hierarchical procedure involves a search of conformational space by the conformational space annealing (CSA) method, followed by optimization of an all-atom model. The procedure was assessed in a recent blind test of protein structure prediction (CASP4). The resulting lowest-energy structures of the target proteins (ranging in size from 70 to 244 residues) agreed with the experimental structures in many respects. The entire experimental structure of a cyclic α-helical protein of 70 residues was predicted to within 4.3 Å α-carbon (Cα) rms deviation (rmsd) whereas, for other α-helical proteins, fragments of roughly 60 residues were predicted to within 6.0 Å Cα rmsd. Whereas β structures can now be predicted with the new procedure, the success rate for α/β- and β-proteins is lower than that for α-proteins at present. For the β portions of α/β structures, the Cα rmsds are less than 6.0 Å for contiguous fragments of 30–40 residues; for one target, three fragments (of length 10, 23, and 28 residues, respectively) formed a compact part of the tertiary structure with a Cα rmsd less than 6.0 Å. Overall, these results constitute an important step toward the ab initio prediction of protein structure solely from the amino acid sequence.

Molecular modeling of the human vasopressin V2 receptor/agonist complex.

The V2 vasopressin renal receptor (V2R), which controls antidiuresis in mammals, is a member of the large family of heptahelical transmembrane (7TM) G protein-coupled receptors (GPCRs). Using the automated GPCR modeling facility available via Internet (http://expasy.hcuge.ch/swissmod/SWISS-MODEL.html) for construction of the 7TM domain in accord with the bovine rhodopsin (RD) footprint, and the SYBYL software for addition of the intra- and extracellular domains, the human V2R was modeled. The structure was further refined and its conformational variability tested by the use of a version of the Constrained Simulated Annealing (CSA) protocol developed in this laboratory. An inspection of the resulting structure reveals that the V2R (likewise any GPCR modeled this way) is much thicker and accordingly forms a more spacious TM cavity than most of the hitherto modeled GPCR constructs do, typically based on the structure of bacteriorhodopsin (BRD). Moreover, in this model the 7TM helices are arranged differently than they are in any BRD-based model. Thus, the topology and geometry of the TM cavity, potentially capable of receiving ligands, is in this model quite different than it is in the earlier models. In the subsequent step, two ligands, the native [arginine8]vasopressin (AVP) and the selective agonist [d-arginine8]vasopressin (DAVP) were inserted, each in two topologically non-equivalent ways, into the TM cavity and the resulting structures were equilibrated and their conformational variabilities tested using CSA as above. The best docking was selected and justified upon consideration of ligand-receptor interactions and structure-activity data. Finally, the amino acid residues were indicated, mainly in TM helices 3-7, as potentially important in both AVP and DAVP docking. Among those Cys112, Val115-Lys116, Gln119, Met123 in helix 3; Glu174 in helix 4; Val206, Ala210, Val213-Phe214 in helix 5; Trp284, Phe287-Phe288, Gln291 in helix 6; and Phe307, Leu310, Ala314 and Asn317 in helix 7 appeared to be the most important ones. Many of these residues are invariant for either the GPCR superfamily or the neurophyseal (vasopressin V2R, V1aR and V1bR and oxytocin OR) subfamily of receptors. Moreover, some of the equivalent residues in V1aR have already been found critical for the ligand affinity [Mouillac et al., J. Biol. Chem, 270 (1995) 25771].

Molecular Dynamics Simulation by GROMACS Using GUI Plugin for PyMOL

The molecular models stored as PDB formatted files are static, but most of the biomolecular systems display a dynamic behavior, in other words their conformations depend on time. To get the dynamic model from the static one, one needs to perform the molecular dynamics (MD) simulation using tools like GROMACS. This paper describes functionality of the newly created plugin for PyMOL (the popular and easy to use program for displaying and manipulating molecule models). This plugin enables the easy use of molecular dynamics simulations using GROMACS through a graphic interface. It transfers the results of those calculations and displays them back in PyMOL. All the components of the stack are open source and are available free of charge. This strategy gives researchers easy access to the molecular dynamics PYMOL plugin and creates an opportunity to modify its source when needed.

Conantokin-Br from Conus brettinghami and selectivity determinants for the NR2D subunit of the NMDA receptor.

Conantokins are venom peptides from marine cone snails that are NMDA receptor antagonists. Here, we report the characterization of a 24 AA conantokin from Conus brettinghami Coomans , H. E. , Moolenbeek , R. G. and Wils , E. ( 1982 ) Basteria 46 ( 1/4 ), 3 - 67 , conantokin-Br (con-Br), the first conantokin that does not have the conserved glutamate residue at position 2. Molecular modeling studies suggest that con-Br has a helical structure between residues 2-13. In contrast to other characterized conantokins, con-Br has a high potency for NMDA receptors with NR2D subunits. To identify determinants for NR2D potency, we synthesized chimeras of con-Br and conantokin-R (con-R); the latter has a approximately 30-fold lower potency for the NR2D subtype. The characterization of two reciprocal chimeras (con-Br/R and con-R/Br), comprising the first 9-10 N-terminal AAs of each conantokin followed by the corresponding C-terminal AAs of the other conantokin demonstrates that determinants for NR2D selectivity are at the N-terminal region. Additional analogues comprising 1-3 amino acid substitutions from each peptide into the homologous region of the other led to the identification of a key determinant; a Tyr residue in position 5 increases potency for NR2D, while Val at this locus causes a decrease. The systematic definition of key determinants in the conantokin peptides for NMDA receptor subtype selectivity is an essential component in the development of conantokin peptides that are highly selective for each specific NMDA receptor subtype.

Caffeine and other methylxanthines as interceptors of food-borne aromatic mutagens: inhibition of Trp-P-1 and Trp-P-2 mutagenic activity.

Caffeine is one of the most important biologically active food components. In this article, we demonstrate that caffeine and other methylxanthines significantly reduce the mutagenic activity of two food-derived heterocyclic aromatic amines, Trp-P-1 and Trp-P-2 in the Salmonella typhimurium TA98 strain. Moreover, protection against Trp-P-1-induced mutagenicity was independent of liver S9 enzymatic fraction, suggesting that mechanisms other than modulation of mutagen bioactivation can contribute to the observed protective effects. UV-vis spectroscopy and computational studies revealed that methylxanthines intercept Trp-P-1 and Trp-P-2 in noncovalent molecular complexes, with association constants (KAC) in the 10(2) M(-1) range. Enthalpy values (ΔH about -30 kJ·mol(-1)) of mutagen-methylxanthine heterocomplexation obtained microcalorimetrically correspond to stacking (π-π) interactions. Finally, we demonstrated that the biological activity of Trp-P-1 and Trp-P-2 is strictly dependent on the presence of the mutagen in a free (unbound with methylxanthine) form, suggesting that mutagen sequestration in stacking heterocomplexes with methylxanthines can decrease its bioavailability and diminish its biological effects.

Phosphorylation and ATP-binding induced conformational changes in the PrkC, Ser/Thr kinase from B. subtilis.

Recent studies on the PrkC, serine-threonine kinase show that that the enzyme is located at the inner membrane of endospores and is responsible for triggering spore germination. The activity of the protein increases considerably after phosphorylation of four threonine residues placed on the activation loop and one serine placed in the C-terminal lobe of the PrkC. The molecular relationship between phosphorylation of these residues and enzyme activity is not known. In this work molecular dynamics simulation is performed on four forms of the protein kinase PrkC from B. subtilis—phosphorylated or unphosphorylated; with or without ATP bound—in order to gain insight into phosphorylation and ATP binding on the conformational changes and functions of the protein kinase. Our results show how phosphorylation, as well as the presence of ATP, is important for the activity of the enzyme through its molecular interaction with the catalytic core residues. Three of four threonine residues were found to be involved in the interactions with conservative motifs important for the enzyme activity. Two of the threonine residues (T167 and T165) are involved in ionic interactions with an arginine cluster from αC-helix. The third residue (T163) plays a crucial role, interacting with His-Arg-Asp triad (HRD). Last of the threonine residues (T162), as well as the serine (S214), were indicated to play a role in the substrate recognition or dimerization of the enzyme. The presence of ATP in the unphosphorylated model induced conformational instability of the activation loop and Asp-Phe-Gly motif (DFG). Based on our calculations we put forward a hypothesis suggesting that the ATP binds after phosphorylation of the activation loop to create a fully active conformation in the closed position.

Improvements in GROMACS plugin for PyMOL including implicit solvent simulations and displaying results of PCA analysis

In order to get the dynamic molecule model from the static one, the molecular dynamics (MD) simulation needs to be performed. Some software sets such as GROMACS are used for that purpose. Unfortunately they lack GUI. The Dynamics PyMOL plugin allows researcher to perform MD simulations directly from the PyMOL software by GUI-based interface of GROMACS tools. This paper describes many improvements introduced into the Dynamics PyMOL plugin 2.0 including: an integration with ProDy library, possibility to use the implicit solvents, an ability to interpret the MD simulations, and implementation of some more GROMACS functionality.

Factors Determining Staphylococcus aureus Susceptibility to Photoantimicrobial Chemotherapy: RsbU Activity, Staphyloxanthin Level, and Membrane Fluidity

Photoantimicrobial chemotherapy (PACT) constitutes a particular type of stress condition, in which bacterial cells induce a pleiotropic and as yet unexplored effect. In light of this, the key master regulators are of putative significance to the overall phototoxic outcome. In Staphylococcus aureus, the alternative sigma factor σB controls the expression of genes involved in the response to environmental stress. We show that aberration of any sigB operon genes in S. aureus USA300 isogenic mutants causes a pronounced sensitization (>5 log10 reduction in CFU drop) to PACT with selected photosensitizers, namely protoporphyrin diarginate, zinc phthalocyanine and rose bengal. This effect is partly due to aberration-coupled staphyloxanthin synthesis inhibition. We identified frequent mutations in RsbU, a σB activator, in PACT-vulnerable clinical isolates of S. aureus, resulting in σB activity impairment. Locations of significant changes in protein structure (IS256 insertion, early STOP codon occurrence, substitutions A230T and A276D) were shown in a theoretical model of S. aureus RsbU. As a phenotypic hallmark of PACT-vulnerable S. aureus strains, we observed an increased fluidity of bacterial cell membrane, which is a result of staphyloxanthin content and other yet unidentified factors. Our research indicates σB as a promising target of adjunctive antimicrobial therapy and suggests that enhanced cell membrane fluidity may be an adjuvant strategy in PACT.

Genistein inhibits activities of methylenetetrahydrofolate reductase and lactate dehydrogenase, enzymes which use NADH as a substrate

Genistein (5, 7-dihydroxy-3- (4-hydroxyphenyl)-4H-1-benzopyran-4-one) is a natural isoflavone revealing many biological activities. Thus, it is considered as a therapeutic compound in as various disorders as cancer, infections and genetic diseases. Here, we demonstrate for the first time that genistein inhibits activities of bacterial methylenetetrahydrofolate reductase (MetF) and lactate dehydrogenase (LDH). Both enzymes use NADH as a substrate, and results of biochemical as well as molecular modeling studies with MetF suggest that genistein may interfere with binding of this dinucleotide to the enzyme. These results have implications for our understanding of biological functions of genistein and its effects on cellular metabolism.

Vasopressin V2 receptor/bioligand interactions

We predict some essential interactions between the V2 vasopressin renal receptor (V2R) and its agonists [Arg8]vasopressin (AVP) and [D-Arg8]vasopressin (DAVP), and the non-peptide antagonist OPC-31260. V2R controls antidiuresis and belongs to the superfamily of heptahelical transmembrane (7TM) G-protein-coupled receptors (GPCRs). The receptor was built, the ligands were docked and the structures relaxed using advanced molecular modeling techniques. Docked agonists and antagonists appear to prefer similar V2R compartments. A number of receptor amino acid residues are indicated, mainly in the TM3-TM7 helices, as potentially important in ligand binding. Many of these residues are invariant for either the GPCR superfamily or the subfamily of related (vasopressin V2R, V1aR and V1bR and oxytocin OR) receptors. Moreover, some of the equivalent residues in V1aR have already been found critical for ligand affinity [Mouillac et al., J. Biol. Chem., 270 (1995) 25771].