Positions 299 and 302 of the GerAA subunit are important for function of the GerA spore germination receptor in Bacillus subtilis.

Bacillus subtilis, as a model spore-forming Gram-positive bacterium, has been extensively used for spore germination research. Within this field, nutrient-dependent germination with specific germinant receptors (GerA, responding to L-alanine or L-valine; GerB and GerK, acting together to start spore germination process in response to AGFK) has been the most studied. There are three different variants of the GerAA subunit (299T/302S, 299A/302P, 299A/302S) of the GerA germination receptor present in B. subtilis subs. subtilis laboratory strains. According to our research, the 299A/302P one, unlike the others, interferes with the spore's ability to germinate in L-alanine as assessed by the measurement of DPA release upon stimulation with the germinant. Multiple genetic manipulations described in this work followed by spore germination tests, together with secondary structure predictions led us to the following conclusions. First, position 302 of GerAA protein is crucial in terms of GerA germination receptor functionality; a proline residue at this position renders the GerA receptor non-functional, most probably due to a change in the protein secondary structure. Second, the 302P GerAA variant has most probably an impaired affinity to other components of GerA receptor. Together, these may explain the loss of GerA receptor's function. Analysis of the GerAA protein should get us closer to understanding the mechanism of GerA receptor function.

Does C60 fullerene act as a transporter of small aromatic molecules?

C60 fullerene is reported to directly interact with biomolecules, such as aromatic mutagens or anticancer drugs. Therefore, it is extensively studied for its potential application in the fields of drug delivery and chemoprevention. Understanding the nature of fullerene-drugs interactions might contribute to optimization and modification of the existing chemotherapy systems. Possible interactions between ICR-191, a model acridine mutagen, with well-established biophysical properties and mutagenic activity, and C60 fullerene aqueous solution were investigated by broad range of biophysical methods, such as Dynamic Light Scattering, Isothermal Titration Calorimetry, and Atomic Force Microscopy. Additionally, to determine biological activity of ICR-191-C60 fullerene mixtures, Ames mutagenicity test was employed. It was demonstrated that C60 fullerene interacts non-covalently with ICR-191 and has strong affinity to bacterial membranes. The obtained results provide practical insight into C60 fullerene interactions with aromatic compounds.

Evidence for interactions between homocysteine and genistein: insights into stroke risk and potential treatment.

Elevated plasma homocysteine (2-amino-4-sulfanylbutanoic acid) level is a risk factor for stroke. Moreover, it has been suggested that high levels of homocysteine in the acute phase of an ischemic stroke can predict mortality, especially in stroke patients with the large-vessel atherosclerosis subtype. In clinical studies, supplementation with genistein (5, 7-dihydroxy-3- (4-hydroxyphenyl)-4H-1-benzopyran-4-one) decreased plasma homocysteine levels considerably. Therefore, genistein could be considered as a potential drug for prevention and/or treatment of stroke. However, the mechanism of the effect of genistein on homocysteine level remains to be elucidated. In this report, direct functional interactions between homocysteine and genistein are demonstrated in in vitro experimental systems for determination of methylenetetrahydrofolate reductase (MetF) and glutathione peroxidase (GPx) activities, reconstructed with purified compounds, and in a simple in vivo system, based on measurement of growth rate of Vibrio harveyi and Bacillus subtilis cultures. Results of molecular modelling indicated that homocysteine can directly interact with genistein. Therefore, genistein-mediated decrease in plasma levels of homocysteine, and alleviation of biochemical and physiological effects of one of these compounds by another, might be ascribed to formation of homocysteine-genistein complexes in which biological activities of these molecules are abolished or alleviated.

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.

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.

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.

Structure of Patt1 human proapoptotic histone acetyltransferase.

The results of modeling of a novel human histone acetyltransferase Patt1 are presented here. This protein belongs to the GNAT GCN5 family and shows proapoptotic activity in human hepatocellular carcinoma cells. Patt1 is an attractive therapeutic target. The sequence analysis, fold recognition predictions and homology modeling of Patt1 protein structure were performed. N- and C- termini of Patt1 were unstructured. Central part revealed classical GNAT fold-central 7-stranded beta sheet core surrounded by intervening 4 alpha helices. The model was assessed with the methods for protein structure validation PROQ and MetaMQAPII. The all-atom 12 ns molecular dynamics simulation of Patt1 model with TIP3P water model and counterions was conducted. All assessment methods implemented resulted in conviction that the model was of quality that could provide confident structural information to infer sequence-structure-function relationships of Patt1. Phe186 and Cys137 were identified as residues engaged in acetyltransfer reaction and the clues for the identification of reaction mechanism were proposed. The knowledge of detailed molecular architecture of Patt1 is not only the key to understanding its mechanistic functional properties but it also opens the possibility of rational drug and protein design experiments, leading to development of effective therapeutic methods.

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.

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.

ATP and its N6-substituted analogues: parameterization, molecular dynamics simulation and conformational analysis.

In this work we used a combination of classical molecular dynamics and simulated annealing techniques to shed more light on the conformational flexibility of 12 adenosine triphosphate (ATP) analogues in a water environment. We present simulations in AMBER force field for ATP and 12 published analogues [Shah et al. (1997) Proc Natl Acad Sci USA 94: 3565-3570]. The calculations were carried out using the generalized Born (GB) solvation model in the presence of the cation Mg(2+). The ion was placed at a close distance (2 Å) from the charged oxygen atoms of the beta and gamma phosphate groups of the -3 negatively charged ATP analogue molecules. Analysis of the results revealed the distribution of inter-proton distances H8-H1' and H8-H2' versus the torsion angle ψ (C4-N9-C1'-O4') for all conformations of ATP analogues. There are two gaps in the distribution of torsion angle ψ values: the first is between -30 and 30 degrees and is described by cis-conformation; and the second is between 90 and 175 degrees, which mostly covers a region of anti conformation. Our results compare favorably with results obtained in experimental assays [Jiang and Mao (2002) Polyhedron 21:435-438].

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 alphaC-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.

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.

The influence of modification at position 2 on the side-chain conformation in oxytocin analogs.

The nonapeptide oxytocin (OT) is used in medicine to help begin and/or continue childbirth. Its analogs can be also used to control bleeding following fetus delivery. The main function of oxytocin is to stimulate contraction of uterus smooth muscle and the smooth muscle of mammary glands, thus regulating lactation. This paper describes theoretical simulations of the distribution of the torsional angles chi1 in the non-standard methylated phenylalanine residues of three oxytocin analogs: [(Phe)(2)o-Me]OT, [(Phe)(2)m-Me]OT, [(Phe)(2)p-Me]OT. The conformations of the oxytocin analogs were studied both in vacuum and in solution. We found some correlations between the biological activity of the considered peptides and the side-chain conformations of amino-acid residues 2 and 8.

Molecular dynamics study of 4-OH-phenylacetyl- D-Y(Me)FQNRPR-NH2 selectivity to V1a receptor.

G protein-coupled receptors relay diverse extracellular signals into cells via a common mechanism, involving activation of cytosol G proteins. The mechanism underlies the actions of approximately 50% of all drugs. In this work, we focus on simulating three protein-ligand complexes of the neurohypophyseal hormone analog 4-OH-phenylacetyl- D-Y(Me)FQNRPR-NH2 (I) with the human V1a, V2 and oxytocin receptors. The peptideI is a potent selective V1a receptor antagonist. To obtain relaxed models of the complexes, the following techniques were used: docking ofI into the vasopressin V1a, V2 and oxytocin receptor models, optimization of the geometry of the resulting complexes and molecular dynamics in a fully hydrated 1-palmitoyl-2-oleoyl- sn-glycero-3-phosphatidylcholine lipid bilayer. The results of the simulations allow us to draw some conclusions about the ligand selectivity to V1aR.

Molecular modeling of the neurohypophyseal receptor/atosiban complexes.

The neurohypophyseal nonapeptide hormone oxytocin (OT) is the strongest uterotonic substance known and is responsible for the initiation of labor. Conversely, oxytocin antagonists blocking uterine OT receptor can suppress uterus contraction. In this paper we describe a computer simulated docking pertinent to affinity of an oxytocin antagonist atosiban towards OT receptor, versus vasopressin V1a and V2 receptors.

Addition of side chains to a known backbone with defined side-chain centroids.

An automatic procedure is proposed for adding side chains to a protein backbone; it is based on optimization of a simplified energy function for peptide side chains, given its backbone and positions of side-chain centroids. The energy is expressed as a sum of the energies of interaction between side chains, and a harmonic penalty function accounting for the preservation of the positions of the C(alpha) atoms and the side-chain centroids. The energy of side-chain interactions is calculated with the soft-sphere ECEPP/3 potential. A Monte Carlo search is carried out to explore all possible side-chain orientations within a fixed backbone and side-chain centroid positions. The initial, usually extended, side-chain conformations are taken directly from the ECEPP/3 database. The procedure was tested on six experimental (X-ray or NMR) structures: immunoglobulin binding protein (PDB code 1IGD, an alpha+beta-protein); transcription factor PML (PDB code 1BOR, a 49-104 fragment of the ring finger domain, predominantly beta-protein); bovine pancreatic trypsin inhibitor (crystal form II) (PDB code 1BPI, an alpha+beta-protein); the monomer of human deoxyhemoglobin (PDB code 1BZ0, an alpha-helical structure); chain A of alcohol dehydrogenase from Drosophila lebanonensis (PDB code 1A4U); as well as on the 10-55 portion of the B domain of staphylococcal protein A (PDB code 1BDD). In all cases except 1BPI, the data for the algorithm (i.e. the backbone or C(alpha) coordinates and the positions of side-chain centroids) were taken from the experimental structures. For protein A, the C(alpha) coordinates and positions of side-chain centroids were also taken from the 1.9-A-resolution model predicted by the UNRES force field. In all comparisons with experimental structures, complete side-chain geometry was reconstructed with a root-mean-square (RMS) deviation of approximately 0.6-0.9 A from the heavy atoms when complete backbone and side-chain-centroid coordinates were used in reconstruction, or approximately 1.0 A when the C(alpha) and centroid coordinates were used.

Energy-based reconstruction of a protein backbone from its alpha-carbon trace by a Monte-Carlo method.

An automatic procedure is proposed for reconstruction of a protein backbone from its C(alpha)-trace; it is based on optimization of a simplified energy function of a peptide backbone, given its alpha-carbon trace. The energy is expressed as a sum of the energies of interaction between backbone peptide groups that are not neighbors in the sequence, the energies of local interactions within all amino acid residues, and a harmonic penalty function accounting for the conservation of standard bond angles. The energy of peptide group interactions is calculated using the assumption that each peptide group acts as a point dipole. For local interaction energy, use is made of a two-dimensional Fourier series expansion of the energies of model terminally blocked amino acid residues, calculated with the Empirical Conformational Energy Program for Peptides (ECEPP/3) force field in the angles lambda((1)) and lambda((2)) defining the rotation of peptide groups adjacent to a C(alpha) carbon atom about the corresponding C(alpha) em leader C(alpha) virtual-bond axes. To explore all possible rotations of peptide groups within a fixed C(alpha)-trace, a Monte Carlo search is carried out. The initial lambda angles are calculated by aligning the dipoles of the peptide groups that are close in space, subject to the condition of favorable local interactions. After the Monte Carlo search is accomplished with the simplified energy function, the energy of the structure is minimized with the ECEPP/3 force field, with imposition of distance constraints corresponding to the initial C(alpha)-trace geometry. The procedure was tested on model alpha-helices and beta-sheets, as well as on the crystal structure of the immunoglobulin binding protein (PDB code: 1IGD, an alpha/beta protein). In all cases, complete backbone geometry was reconstructed with a root-mean-square (rms) deviation of 0.5 A from the all-atom target structure.