Vladimir Frecer

De novo design of potent antimicrobial peptides.

ABSTRACT Lipopolysaccharide (LPS), shed by gram-negative bacteria during infection and antimicrobial therapy, may lead to lethal endotoxic shock syndrome. A rational design strategy based on the presumed mechanism of antibacterial effect was adopted to design cationic antimicrobial peptides capable of binding to LPS through tandemly repeated sequences of alternating cationic and nonpolar residues. The peptides were designed to achieve enhanced antimicrobial potency due to initial bacterial membrane binding with a reduced risk of endotoxic shock. The peptides designed displayed binding affinities to LPS and lipid A (LA) in the low micromolar range and by molecular modeling were predicted to form amphipathic β-hairpin-like structures when they bind to LPS or LA. They also exhibited strong effects against gram-negative bacteria, with MICs in the nanomolar range, and low cytotoxic and hemolytic activities at concentrations significantly exceeding their MICs. Quantitative structure-activity relationship (QSAR) analysis of peptide sequences and their antimicrobial, cytotoxic, and hemolytic activities revealed that site-directed substitutions of residues in the hydrophobic face of the amphipathic peptides with less lipophilic residues selectively decrease the hemolytic effect without significantly affecting the antimicrobial or cytotoxic activity. On the other hand, the antimicrobial effect can be enhanced by substitutions in the polar face with more polar residues, which increase the amphipathicity of the peptide. On the basis of the QSARs, new analogs that have strong antimicrobial effects but that lack hemolytic activity can be proposed. The findings highlight the importance of peptide amphipathicity and allow a rational method that can be used to dissociate the antimicrobial and hemolytic effects of cationic peptides, which have potent antimicrobial properties, to be proposed.

Dynamic Behavior of Avian Influenza A Virus Neuraminidase Subtype H5N1 in Complex with Oseltamivir, Zanamivir, Peramivir, and Their Phosphonate Analogues

The outbreak of avian influenza A subtype H5N1 virus has raised a global concern for both animal as well as human health. Recently, drug resistance in H5N1 infections has been widely reported due to neuraminidase mutations. Consequently, the understanding of inhibitor-neuraminidase interactions at the molecular level represents the main goal of our study. Molecular dynamics simulations were carried out for the neuraminidase N1 in complex with six inhibitors--oseltamivir, zanamivir, peramivir, and their phosphonate analogues. Molecular dynamics trajectories were extensively analyzed in terms of important interactions between inhibitors and the enzyme target. Results show that open and closed forms (defined by the relative position of the flexible 150-loop) of neuraminidase N1 interchange during the course of 20 ns molecular dynamics simulation of the protein-inhibitor complexes. Reported free energies of closing indicate that the carboxylate inhibitors prefer the closed form more than their phosphonate analogues. This can be understood in view of the negative total charge (-1 e0) of the phosphonate inhibitors, which repels the Asp151 residue of the loop away from the inhibitor and drives the complex into the open form. Obtained results constitute new valuable information to assist further drug development of inhibitors against the H5N1 avian influenza A virus and could also inspire similar studies for other systems of the influenza family such as the 2009 influenza A (H1N1) virus.

Structure-based design of inhibitors of NS3 serine protease of hepatitis C virus

We have designed small focused combinatorial library of hexapeptide inhibitors of NS3 serine protease of the hepatitis C virus (HCV) by structure-based molecular design complemented by combinatorial optimisation of the individual residues. Rational residue substitutions were guided by the structure and properties of the binding pockets of the enzymes active site. The inhibitors were derived from peptides known to inhibit the NS3 serine protease by using unusual amino acids and alpha-ketocysteine or difluoroaminobutyric acid, which are known to bind to the S1 pocket of the catalytic site. Inhibition constants (Ki) of the designed library of inhibitors were predicted from a QSAR model that correlated experimental Ki of known peptidic inhibitors of NS3 with the enthalpies of enzyme-inhibitor interaction computed via molecular mechanics and the solvent effect contribution to the binding affinity derived from the continuum model of solvation. The library of the optimised inhibitors contains promising drug candidates-water-soluble anionic hexapeptides with predicted Ki* in the picomolar range.

Temperature dependence of estrogen binding: importance of a subzone in the ligand binding domain of a novel piscine estrogen receptor.

The full length estrogen receptor from Oreochromis aureus (OaER) was cloned and expressed in vitro and in vivo as a functional transcription factor. Amino acid residues involved in the thermal stability of the receptor are located at/near subzones beta1 and beta3, which are highly conserved in other non-piscine species but not in OaER. Hormone binding studies, however, indicate that OaER is thermally stable but exhibited a approximately 3-fold reduced affinity for estrogen at elevated temperatures. Transfection of OaER into various cell lines cultured at different temperatures displayed a significant estrogen dose-response shift compared with that of chicken ER (cER). At 37 degrees C, OaER requires approximately 80-fold more estrogen to achieve half-maximal stimulation of CAT. Lowering of the incubation temperature from 37 degrees C to 25 degrees C or 20 degrees C resulted in a 4-fold increase in its affinity for estrogen. The thermally deficient transactivation of OaER at temperatures above 25 degrees C was fully prevented by high levels of estrogen. Thus, compared to cER, the OaER exhibits reduced affinity for estrogen at elevated temperature as reflected in its deficient transactivation capability. Amino acid replacements of OaER beta3 subzones with corresponding amino acids from cER could partially rescue this temperature sensitivity. The three-dimensional structure of the OaER ligand binding domain (LBD) was modelled based on conformational similarity and sequence homology with human RXRalpha apo, RARgamma holo and ERalpha LBDs. Unliganded and 17beta-estradiol-liganded OaER LBD retained the overall folding pattern of the nuclear receptor LBDs. The residues at/near the subzone beta3 of the LBD constitute the central core of OaER structure. Thus, amino acid alteration at this region potentially alters the structure and consequently its temperature-dependent ligand binding properties.

Design of oseltamivir analogs inhibiting neuraminidase of avian influenza virus H5N1

Neuraminidase is an important target for design of antiviral agents in the prophylaxis and treatment of avian influenza virus infections. We have shown the applicability of computer-assisted combinatorial techniques in the design, focusing and in silico screening of a virtual library of analogs of oseltamivir (Tamiflu) with the goal to find potent inhibitors of influenza A neuraminidase N1 that fill the cavity found adjacent to the active site. Crystal structure of oseltamivir-N1 complex was used in the structure-based focusing and virtual screening of the designed library. A target-specific Piecewise Linear Potential type 1 scoring function fitted for a training set of 14 carbocyclic inhibitors and validated for three other inhibitors was used to select virtual hits with predicted inhibitory activities in the subnanomolar range. The results of this computational study are useful as a rational guide for synthetic and medicinal chemists who are developing new drugs against the avian influenza virus H5N1.

A SALL4 zinc finger missense mutation predicted to result in increased DNA binding affinity is associated with cranial midline defects and mild features of Okihiro syndrome

Truncating mutations of the gene SALL4 on chromosome 20q13.13–13.2 cause Okihiro and acro-renal-ocular syndromes. Pathogenic missense mutations within the SALL4 or SALL1 genes have not yet been reported, raising the question which phenotypic features would be associated with them. Here we describe the first missense mutation within the SALL4 gene. The mutation results in an exchange of a highly conserved zinc-coordinating Histidine crucial for zinc finger (ZF) structure within a C2H2 double ZF domain to an Arginine. Molecular modeling predicts that this exchange does not result in a loss of zinc ion binding but leads to an increased DNA-binding affinity of the domain. The index patient shows mild features of Okihiro syndrome, but in addition cranial midline defects (pituitary hypoplasia and single central incisor). This finding illustrates that the phenotypic and functional effects of SALL4 missense mutations are difficult to predict, and that other SALL4 missense mutations might lead to phenotypes not overlapping with Okihiro syndrome.

Molecular dynamics study on lipid A from Escherichia coli: insights into its mechanism of biological action.

Structural properties of the Escherichia coli lipid A moiety were analysed by means of molecular mechanics and molecular dynamics simulations and compared to synthetic monophospho and dephospho analogues with different biological activities in the Limulus assay. The conformation of glucosamine disaccharide headgroup, order and packing of fatty acid chains, solvation of phosphate groups, coordination by water molecules, sodium counterions and models of cationic amino acid side chains were described in terms of mean values, mean residence times, radial distribution functions, coordination numbers, solvation and interaction energies. Solvation and polar interactions of the phosphate groups were correlated to known biological activities the lipid A variants. The observed relationship between the biological effect and the number and position of the phosphate groups were explained with the help of simple mechanistic models of lipid A action. The possible mechanism of action involving specific binding of lipid A disaccharide headgroup to cationic residues of a receptor model was compared with an alternative mechanism, which assumes a relationship between the ability to adopt non-lamellar supramolecular structures and the biological activity. Conclusions are drawn about the probable mode of lipid A action. Implications for rational drug design of endotoxin-neutralising agents are discussed.

Design, structure-based focusing and in silico screening of combinatorial library of peptidomimetic inhibitors of Dengue virus NS2B-NS3 protease.

Serine protease activity of the NS3 protein of Dengue virus is an important target of antiviral agents that interfere with the viral polyprotein precursor processing catalyzed by the NS3 protease (NS3pro), which is important for the viral replication and maturation. Recent studies showed that substrate-based peptidomimetics carrying an electrophilic warhead inhibit the NS2B-NS3pro cofactor-protease complex with inhibition constants in the low micromolar concentration range when basic amino acid residues occupy P1 and P2 positions of the inhibitor, and an aldehyde warhead is attached to the P1. We have used computer-assisted combinatorial techniques to design, focus using the NS2B-NS3pro receptor 3D structure, and in silico screen a virtual library of more than 9,200 peptidomimetic analogs targeted around the template inhibitor Bz-Nle-Lys-Arg-Arg-H (Bz—benzoyl) that are composed mainly of unusual amino acid residues in all positions P1–P4. The most promising virtual hits were analyzed in terms of computed enzyme-inhibitor interactions and Adsorption, Distribution, Metabolism and Excretion (ADME) related physico-chemical properties. Our study can direct the interest of medicinal chemists working on a next generation of antiviral chemotherapeutics against the Dengue Fever towards the explored subset of the chemical space that is predicted to contain peptide aldehydes with NS3pro inhibition potencies in nanomolar range which display ADME-related properties comparable to the training set inhibitors.Graphical abstract

The extended polarizable continuum model for calculation of solvent effects

Abstract An extended version of the polarizable continuum model of salvation is presented. Principles of the particular models for electrostatic, dispersion, repulsion, and cavitation terms evaluation are shown and the recent state in methodology is described. Modifications suitable for large biomolecules, fields of application and possible further development are discussed. An illustrative example showing the significance of individual salvation Gibbs free energy contributions is given.

Molecular Interfaces of the Galactose-binding Protein Tectonin Domains in Host-Pathogen Interaction

β-Propeller proteins function in catalysis, protein-protein interaction, cell cycle regulation, and innate immunity. The galactose-binding protein (GBP) from the plasma of the horseshoe crab, Carcinoscorpius rotundicauda, is a β-propeller protein that functions in antimicrobial defense. Studies have shown that upon binding to Gram-negative bacterial lipopolysaccharide (LPS), GBP interacts with C-reactive protein (CRP) to form a pathogen-recognition complex, which helps to eliminate invading microbes. However, the molecular basis of interactions between GBP and LPS and how it interplays with CRP remain largely unknown. By homology modeling, we showed that GBP contains six β-propeller/Tectonin domains. Ligand docking indicated that Tectonin domains 6 to 1 likely contain the LPS binding sites. Protein-protein interaction studies demonstrated that Tectonin domain 4 interacts most strongly with CRP. Hydrogen-deuterium exchange mass spectrometry mapped distinct sites of GBP that interact with LPS and with CRP, consistent with in silico predictions. Furthermore, infection condition (lowered Ca2+ level) increases GBP-CRP affinity by 1000-fold. Resupplementing the system with a physiological level of Ca2+ did not reverse the protein-protein affinity to the basal state, suggesting that the infection-induced complex had undergone irreversible conformational change. We propose that GBP serves as a bridging molecule, participating in molecular interactions, GBP-LPS and GBP-CRP, to form a stable pathogen-recognition complex. The interaction interfaces in these two partners suggest that Tectonin domains can differentiate self/nonself, crucial to frontline defense against infection. In addition, GBP shares architectural and functional homologies to a human protein, hTectonin, suggesting its evolutionarily conservation for ∼500 million years, from horseshoe crab to human.