Myricetin protects Galleria mellonella against Staphylococcus aureus infection and inhibits multiple virulence factors.

Staphylococcus aureus is an opportunistic pathogen related to a variety of life-threatening infections but for which antimicrobial resistance is liming the treatment options. We report here that myricetin, but not its glycosylated form, can remarkably decrease the production of several S. aureus virulence factors, including adhesion, biofilm formation, hemolysis and staphyloxanthin production, without interfering with growth. Myricetin affects both surface proteins and secreted proteins which indicate that its action is unrelated to inhibition of the agr quorum sensing system. Analysis of virulence related gene expression and computational simulations of pivotal proteins involved in pathogenesis demonstrate that myricetin downregulates the saeR global regulator and interacts with sortase A and α-hemolysin. Furthermore, Myr confers a significant degree of protection against staphylococcal infection in the Galleria mellonella model. The present findings reveal the potential of Myr as an alternative multi-target antivirulence candidate to control S. aureus pathogenicity.

Membrane negative curvature induced by a hybrid peptide from pediocin PA-1 and plantaricin 149 as revealed by atomistic molecular dynamics simulations.

Antimicrobial peptides (AMPs) are cationic peptides that kill bacteria with a broad spectrum of action, low toxicity to mammalian cells and exceptionally low rates of bacterial resistance. These features have led to considerable efforts in developing AMPs as an alternative antibacterial therapy. In vitro studies have shown that AMPs interfere with membrane bilayer integrity via several possible mechanisms, which are not entirely understood. We have performed the synthesis, membrane lysis measurements, and biophysical characterization of a novel hybrid peptide. These measurements show that PA-Pln149 does not form nanopores, but instead promotes membrane rupture. It causes fast rupture of the bacterial model membrane (POPG-rich) at concentrations 100-fold lower than that required for the disruption of mammalian model membranes (POPC-rich). Atomistic molecular dynamics (MD) simulations were performed for single and multiple copies of PA-Pln149 in the presence of mixed and pure POPC/POPG bilayers to investigate the concentration-dependent membrane disruption by the hybrid peptide. These simulations reproduced the experimental trend and provided a potential mechanism of action for PA-Pln149. It shows that the PA-Pln149 does not form nanopores, but instead promotes membrane destabilization through peptide aggregation and induction of membrane negative curvature with the collapse of the lamellar arrangement. The sequence of events depicted for PA-Pln149 may offer insights into the mechanism of action of AMPs previously shown to induce negative deformation of membrane curvature and often associated with peptide translocation via non-bilayer intermediate structures.

Effect of dimerization on the mechanism of action of aurein 1.2.

The mechanism of action of antimicrobial peptides depends on physicochemical properties such as structure, concentration, and oligomerization. Here, we focused on the effect of dimerization on the mechanism of action of aurein 1.2 (AU). We designed a lysine-linked AU dimer, (AU)2K, and its interaction with membrane mimetics was studied using four biophysical techniques and molecular dynamics simulations. Circular dichroism and molecular dynamics studies showed that AU displayed a typical spectrum for disordered structures in aqueous solution whereas (AU)2K exhibited the typical spectrum of α-helices in a coiled-coil conformation, wherein helices are wrapped around each other. With the addition of large unilamellar vesicles (LUVs), AU adopted an α-helix structure whereas the coiled-coil structure of (AU)2K assumed an extended conformation. Carboxyfluorescein release experiments with LUVs showed that both peptides were able to permeabilize vesicles although the leakage response to increases in peptide concentration differed. Optical microscopy experiments showed that both peptides induced pore opening and the dimer eventually caused the vesicles to burst. Finally, calorimetric traces determined by isothermal titration calorimetry on the LUVs also showed significant differences in peptide-membrane interactions. Together, the results of our study demonstrated that dimerization changes the mechanism of action of AU.

Characterization of the outer membrane protein OprF of Pseudomonas aeruginosa in a lipopolysaccharide membrane by computer simulation.

The N-terminal domain of outer membrane protein OprF of Pseudomonas aeruginosa forms a membrane spanning eight-stranded antiparallel beta-barrel domain that folds into a membrane channel with low conductance. The structure of this protein has been modeled after the crystal structure of the homologous protein OmpA of Escherichia coli. A number of molecular dynamics simulations have been carried out for the homology modeled structure of OprF in an explicit molecular model for the rough lipopolysaccharide (LPS) outer membrane of P. aeruginosa. The structural stability of the outer membrane model as a result of the strong electrostatic interactions compared with simple lipid bilayers is restricting both the conformational flexibility and the lateral diffusion of the porin in the membrane. Constricting side-chain interactions within the pore are similar to those found in reported simulations of the protein in a solvated lipid bilayer membrane. Because of the strong interactions between the loop regions of OprF and functional groups in the saccharide core of the LPS, the entrance to the channel from the extracellular space is widened compared with the lipid bilayer simulations in which the loops are extruding in the solvent. The specific electrostatic signature of the LPS membrane, which results in a net intrinsic dipole across the membrane, is found to be altered by the presence of OprF, resulting in a small electrically positive patch at the position of the channel.

Validation of the GROMOS force-field parameter set 45Alpha3 against nuclear magnetic resonance data of hen egg lysozyme.

The quality of molecular dynamics (MD) simulations of proteins depends critically on the biomolecular force field that is used. Such force fields are defined by force-field parameter sets, which are generally determined and improved through calibration of properties of small molecules against experimental or theoretical data. By application to large molecules such as proteins, a new force-field parameter set can be validated. We report two 3.5 ns molecular dynamics simulations of hen egg white lysozyme in water applying the widely used GROMOS force-field parameter set 43Alpha1 and a new set 45Alpha3. The two MD ensembles are evaluated against NMR spectroscopic data NOE atom-atom distance bounds, (3)J(NHalpha) and (3)J(alphabeta) coupling constants, and (15)N relaxation data. It is shown that the two sets reproduce structural properties about equally well. The 45Alpha3 ensemble fulfills the atom-atom distance bounds derived from NMR spectroscopy slightly less well than the 43Alpha1 ensemble, with most of the NOE distance violations in both ensembles involving residues located in loops or flexible regions of the protein. Convergence patterns are very similar in both simulations atom-positional root-mean-square differences (RMSD) with respect to the X-ray and NMR model structures and NOE inter-proton distances converge within 1.0-1.5 ns while backbone (3)J(HNalpha)-coupling constants and (1)H-(15)N order parameters take slightly longer, 1.0-2.0 ns. As expected, side-chain (3)J(alphabeta)-coupling constants and (1)H-(15)N order parameters do not reach full convergence for all residues in the time period simulated. This is particularly noticeable for side chains which display rare structural transitions. When comparing each simulation trajectory with an older and a newer set of experimental NOE data on lysozyme, it is found that the newer, larger, set of experimental data agrees as well with each of the simulations. In other words, the experimental data converged towards the theoretical result.

Internal dynamics and ionization states of the macrophage migration inhibitory factor: comparison between wild-type and mutant forms.

The macrophage migration inhibitory factor (MIF) is a cytokine that shares a common structural architecture and catalytic strategy with three isomerases: 4-oxalocrotonate tautomerase, 5-carboxymethyl-2-hydroxymuconate isomerase, and D-dopachrome tautomerase. A highly conserved N-terminal proline acts as a base-acid during the proton transfer reaction catalyzed by these enzymes. Such unusual catalytic strategy appears to be possible only due to the N-terminal proline pK(a) shifted to 5.0-6.0 units. Mutations of this residue result in a significant decrease of the catalytic activity of MIF. Two hypotheses have been proposed to explain the catalytic inefficiency of MIF: the lower basicity of primary amines with regard to secondary ones and the increased flexibility resulting from the replacement of a proline by residues like glycine. To investigate that, we have performed molecular dynamics simulations of MIF wild-type and its mutant P1G, as well as calculated the protonation properties of several mutant forms. It was found that the N-terminal glycine does not show larger fluctuations compared to proline, but the former residue is more exposed to the solvent throughout the simulations. The apparent pK(a) of these residues displays very little change (as expected from the structural rigidity of MIF) and is not significantly affected by the surrounding ionizable residues. Instead, the hydrophobic character of the active site seems to be the main factor in determining the pKa of the N-terminal residue and the catalytic efficiency of MIF.

Revisiting the structural flexibility of the complex p21(ras)-GTP: the catalytic conformation of the molecular switch II.

The hydrolysis of GTP in p21(ras) triggers conformational changes that regulate the ras/ERK signaling pathway. An important active site residue is Gln61, which has been found to be mutated in 30% of human tumors. The dynamics of the active site conformation is studied by using molecular dynamics simulation of two independent structures of the GTP-bound uncomplexed enzyme. Two distinct conformations of the enzyme are observed, in which the side-chain residue Gln61 is in different orientations. Essential dynamics analysis is used to describe the essential motions in the transition between the two conformations. Results are compared with earlier simulations of p21(ras) and its complex with GTPase activating protein p21-GAP.

Investigations on human immunodeficiency virus type 1 integrase/DNA binding interactions via molecular dynamics and electrostatics calculations.

The complete three-dimensional structure of the active site region of the human immunodeficiency virus type 1 (HIV-1) integrase (IN) is not unambiguously known. This region includes a flexible loop comprising residues 141-148 and the N-terminal portion of the helix alpha-4, which contains E152, the third catalytic residue, and Y143, which plays a secondary role in catalysis. Relatively high B-factors exist for most of the residues in the aforementioned region. The HIV-1 IN belongs to the polynucleotidyl transferase superfamily, whose members have been proposed to use two divalent metal ions for catalysis. Although only the position of the first metal ion has been determined crystallographically for the HIV-1 IN, we recently have proposed a binding site for the second metal ion. Based on this information, we have performed two 500-psec molecular dynamics simulations of the catalytic domain of the HIV-1 IN containing two Mg(2)+ ions. In one of the simulations, we included a dianionic phosphate group (HPO(4)(2)-) in the active site to mimic a portion of the DNA backbone of a substrate for the integration reaction. Electrostatics calculations and ionization state predictions were carried out on representative structures taken from the molecular dynamics simulations. Different conformational behaviors of the enzyme were observed, depending upon whether two Mg(2)+ ions were bound or two Mg(2)+ ions plus phosphate. The electrostatic calculations performed on the dynamical structures provide a further refinement about which regions of the catalytic domain of the HIV-1 IN may be involved in the DNA binding.

Solvent accessibility to aspartyl and succinimidyl residues at positions 7 and 23 in the amyloid beta 1-28 peptide.

The water accessibilities to aspartyl residues at positions 7 and 23 in the amyloid beta 1-28 peptide associated with Alzheimer's Disease have been calculated using different techniques. These accessibilities of water were compared to those of the succinimidyl residues (SUC) replacing the aspartyl ones (ASP). It has been possible to ascertain that these modifications (ASP--->SUC) lead to a significant increase in the water accessibility to the backbone and alpha-carbon atom of the SUC7 and SUC23 residues. It is suggested that the spontaneous transformation of the ASP--->SUC might lead to an increase of the racemization rates due to the higher accessibility of water at these sites. It is also proposed that the behavior of the adjacent residues in the selectivity of the racemization is to control the water accessibility at the reactive residue.

Docking of 4-oxalocrotonate tautomerase substrates: implications for the catalytic mechanism.

The enzyme 4-oxalocrotonate tautomerase catalyzes the ketonization of dienols, which after further processing become intermediates in the Krebs cycle. The enzyme uses a general acid-base mechanism for proton transfer: the amino-terminal proline has been shown to function as the catalytic base and Arg39 has been implicated as the catalytic acid. We report the results of molecular docking simulations of 4-oxalocrotonate tautomerase with two substrates, 2-hydroxymuconate and 5-carboxymethyl-2-hydroxymuconate. pKa calculations are also performed for the free enzyme. The predicted binding mode of 2-hydroxymuconate is in agreement with experimental data. A model for the binding mode of 5-carboxymethyl-2-hydroxymuconate is proposed which explains the lower catalytic efficiency of the enzyme toward this substrate. The pKa predictions and docking simulations support residue Arg39 as the general acid for the enzyme catalysis.

Plural origins of molecular homochirality in our biota Part II. The relative stabilities of homochiral and mixed oligoribotides and peptides.

By computer simulations--molecular mechanics and molecular dynamics with the amber force field (Weiner et al, (1986), J. Comp. Chem. 7, 230-252)--we have determined the stabilities of oligoribotide strands built with D- and L-riboses, and of peptide chains with D- and L-amino acid residues. In particular, complementary double-chains of oligoribotides were studied, since they are an important feature of the growing mechanism of modern nucleic acids. Peptide chains on the other hand, grow without need of a template. We found that mixed oligoribotides are less stable than homochiral ones, and that this chiral effect is less noticeable in peptide chains. The results support the interpretation that L-riboses act as terminators to the template-assisted growth of oligo-r-GD (enantiometric cross-inhibition; Joyce et al., (1987), Proc. Natl. Acad. Sci. USA 84, 4398-4402). Based on this effect, a chemical pathway is proposed which could, under assumed prebiotic conditions, bypass the hindrance of homochiral growth.

Biochemical correlate of depression in children.

The degree of depression in 88 abandoned children was analysed through a depression rating scale adapted to prepubertal children. The items were grouped into three dimensions: sociological-relational, psychological and biological. In 46 children from this sample it was dosed plasmatic cortisol and the urinary excretion of catecholamine, VMA, HVA and 5-HIAA. When analysing the principal components, the sociological and psychological dimensions were the most important ones in the sample, followed by age and catecholamine variables. The group of male depressed children presented a higher level of catecholamine urinary excretion and a lower peak of plasmatic cortisol than the non-depressed group. The variable age, in both sexes, was correlated with the biochemical variable catecholamine. Biochemical alterations are present in depressed children but it is difficult to show a correlation of dependence between them and the phenomenological aspects of depression.

Biochemical correlate of depression in children

O grau de depressäo em 88 crianças abandonadas foi analisado por escala de avaliaçäo de depressäo adaptada para crianças na pré-puberdade. Os itens da escala foram agrupados em três dimensöes: sociológico-relacional, psicológica e biológica. Em 46 crianças desta amostra foram dosados o cortisol plasmático e a excreçäo urinária de catecolomina, VMA, HVA e 5-HIAA. Pela análise dos principais componentes, mostraram-se mais importantes, na amostra, as dimensöes sociológica e psicológica, seguidas dos componentes idade e catecolomina. O grupo de crianças do sexo masculino com depressäo apresentava maior nível de excreçäo urinária de catecolomina e menor pico de cortisol plasmático que o grupo sem depressäo. A variável idade, em ambos os sexos, correlacionava-se à variável catecolomina. Alteraçöes bioquímicas estäo presentes em crianças com depressäo, mas é difícil demonstrar correlaçäo de dependência entre elas e aspectos fenomenológicos da depressäo