Structure based computational assessment of channel properties of assembled ORF‐8a from SARS‐CoV

ORF 8a is a short 39 amino acid bitopic membrane protein encoded by severe acute respiratory syndrome causing corona virus (SARS‐CoV). It has been identified to increase permeability of the lipid membrane for cations. Permeability is suggested to occur due to the assembly of helical bundles. Computational models of a pentameric assembly of 8a peptides are generated using the first 22 amino acids, which include the transmembrane domain. Low energy structures reveal a hydrophilic pore mantled by residues Thr‐8, and ?18, Ser‐11, Cys‐13, and Arg‐22. Potential of mean force (PMF) profiles for mono (Na⁺, K⁺, Cl^?) and divalent (Ca²⁺) ions along the pore are calculated. The data support experimental findings of a weak cation selectivity of the channel. Calculations on 8a are compared to data derived for a pentameric bundle consisting of the M2 helices of the bacterial pentameric ligand gated ion channel GLIC (3EHZ). PMF curves of both, bundles 8a and M2, show sigmoidal shaped profiles. In comparison to the data for the M2‐GLIC model, data of the 8a bundle show lower amplitude of the PMF values between maximum and minimum and less discrimination amongst ions. Proteins 2015; 83:300–308. © 2014 Wiley Periodicals, Inc.     

Application of biasing‐potential replica‐exchange simulations for loop modeling and refinement of proteins in explicit solvent

Comparative or homology modeling of a target protein based on sequence similarity to a protein with known structure is widely used to provide structural models of proteins. Depending on the target‐template similarity these model structures may contain regions of limited structural accuracy. In principle, molecular dynamics (MD) simulations can be used to refine protein model structures and also to model loop regions that connect structurally conserved regions but it is limited by the currently accessible simulation time scales. A recently developed biasing potential replica exchange (BP‐REMD) method was used to refine loops and complete decoy protein structures at atomic resolution including explicit solvent. In standard REMD simulations several replicas of a system are run in parallel at different temperatures allowing exchanges at preset time intervals. In a BP‐REMD simulation replicas are controlled by various levels of a biasing potential to reduce the energy barriers associated with peptide backbone dihedral transitions. The method requires much fewer replicas for efficient sampling compared with T‐REMD. Application of the approach to several protein loops indicated improved conformational sampling of backbone dihedral angle of loop residues compared to conventional MD simulations. BP‐REMD refinement simulations on several test cases starting from decoy structures deviating significantly from the native structure resulted in final structures in much closer agreement with experiment compared to conventional MD simulations. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Structural characterization of human tryptophan hydroxylase 2 reveals that L-Phe is superior to L-Trp as the regulatory domain ligand

1. Download : Download high-res image (253KB)
2. Download : Download full-size image

     

Structural analysis of the N‐terminal fragment of the antiangiogenic protein endostatin: A molecular dynamics study

Endostatin is a potent antiangiogenic protein derived from the noncollagenous domain 1 (NC1) of collagen XVIII. The mechanism by which endostatin exerts its antiangiogenic effect is still incompletely understood. It has been shown that the 27 amino acid N‐terminal fragment of murine endostatin has antitumor, antimigration, and antipermeability activities comparable to the full soluble protein. To understand how this peptide can exert such elaborate function, we performed structural analysis using molecular dynamics to evaluate the behavior of this fragment in aqueous environment. Here, we show that the N‐terminal peptide of murine endostatin is able to assume a well‐defined structure, folding into a zinc‐dependent β‐hairpin conformation. Analyzing the folding mechanism, we were able to understand why the N‐terminal peptide of human endostatin with the same length failed to acquire a stable conformation. Conversely, we were able to predict the successful folding of the R4Q mutant and of a shorter form of the human peptide with 25 residues. Finally, we show that the β‐hairpin conformation assumed by the zinc‐bound peptide of murine endostatin has a high structural similarity with fragments of another family of angiogenesis inhibitors: the integrin‐binding portion of the NC1 domain of collagen IV. Indeed, our docking simulations show that arresten, canstatin, and the endostatin peptide bind to the same spot of αVβ3 integrin, suggesting similar interactions via a common binding site on this receptor. Proteins 2011;. © 2011 Wiley‐Liss, Inc.     

Structural and functional characterization of peptides derived from the carboxy‐terminal region of a defensin from the tick Ornithodoros savignyi

Tick defensins may serve as templates for the development of multifunctional peptides. The purpose of this study was to evaluate shorter peptides derived from tick defensin isoform 2 (OsDef2) in terms of their antibacterial, antioxidant, and cytotoxic activities. We compared the structural and functional properties of a synthetic peptide derived from the carboxy‐terminal of the parent peptide (Os) to that of an analogue in which the three cysteine residues were omitted (Os–C). Here, we report that both peptides were bactericidal (MBC values ranging from 0.94–15 µg/ml) to both Gram‐positive and Gram‐negative bacteria, whereas the parent peptide only exhibited Gram‐positive antibacterial activity. The Os peptide was found to be two‐fold more active than Os–C against three of the four tested bacteria but equally active against Staphylococcus aureus. Os showed rapid killing kinetics against both Escherichia coli and Bacillus subtilis, whereas Os–C took longer, suggesting different modes of action. Scanning electron microscopy showed that in contrast to melittin for which blebbing of bacterial surfaces was observed, cells exposed to either peptide appeared flattened and empty. Circular dichroism data indicated that in a membrane‐mimicking environment, the cysteine‐containing peptide has a higher α‐helical content. Both peptides were found to be non‐toxic to mammalian cells. Moreover, the peptides displayed potent antioxidant activity and were 12 times more active than melittin. Multifunctional peptides hold potential for a wide range of clinical applications and further investigation into their mode of antibacterial and antioxidant properties is therefore warranted. Copyright © 2013 European Peptide Society and John Wiley & Sons, Ltd.     

Structural and biochemical characterization of a metagenome‐derived esterase with a long N‐terminal extension

The genes encoding six novel esterolytic/lipolytic enzymes, termed LC‐Est1～6, were isolated from a fosmid library of a leaf‐branch compost metagenome by functional screening using tributyrin agar plates. These enzymes greatly vary in size and amino acid sequence. The highest identity between the amino acid sequence of each enzyme and that available from the database varies from 44 to 73%. Of these metagenome‐derived enzymes, LC‐Est1 is characterized by the presence of a long N‐terminal extension (LNTE, residues 26–283) between a putative signal peptide (residues 1–25) and a C‐terminal esterase domain (residues 284–510). A putative esterase from Candidatus Solibacter usitatus (CSu‐Est) is the only protein, which shows the significant amino acid sequence identity (46%) to the entire region of LC‐Est1. To examine whether LC‐Est1 exhibits activity and its LNTE is important for activity and stability of the esterase domain, LC‐Est1 (residues 26–510), LC‐Est1C (residues 284–510), and LC‐Est1C* (residues 304–510) were overproduced in E. coli, purified, and characterized. LC‐Est1C* was only used for structural analysis. The crystal structure of LC‐Est1C* highly resembles that of the catalytic domain of Thermotoga maritima esterase, suggesting that LNTE is not required for folding of the esterase domain. The enzymatic activity of LC‐Est1C was lower than that of LC‐Est1 by 60%, although its substrate specificity was similar to that of LC‐Est1. LC‐Est1C was less stable than LC‐Est1 by 3.3°C. These results suggest that LNTE of LC‐Est1 rather exists as an independent domain but is required for maximal activity and stability of the esterase domain.     

Gαs protein C‐terminal α‐helix at the interface: does the plasma membrane play a critical role in the Gαs protein functionality?

The heterotrimeric guanine nucleotide-binding regulatory proteins (G proteins, Galphabetagamma) mediate the signalling process of a large number of receptors, known as G protein-coupled receptors. The C-terminal domain of the heterotrimeric G protein alpha-subunit plays a key role in the selective activation of G proteins by their cognate receptors. The interaction of this domain can take place at the end of a cascade including several successive conformational modifications. Galpha(s)(350-394) is the 45-mer peptide corresponding to the C-terminal region of the Galpha(s) subunit. In the crystal structure of the Galpha(s) subunit it encompasses the alpha4/beta6 loop, the beta6 beta-sheet segment and the alpha5 helix region. Following a previous study based on the synthesis, biological activity and conformational analysis of shorter peptides belonging to the same Galpha(s) region, Galpha(s)(350-394) was synthesized and investigated. The present study outlines the central role played by the residues involved in the alpha4/beta6 loop and beta6/alpha5 loops in the stabilization of the C-terminal Galpha(s)alpha-helix. H(2)O/(2)H(2)O exchange experiments, and NMR diffusion experiments show interesting evidence concerning the interaction between the SDS micelles and the polypeptide. These data prompt intriguing speculations on the role of the intracellular environment/cellular membrane interface in the stabilization and functionality of the C-terminal Galpha(s) region.     

DRESS: a database of REfined solution NMR structures

Several studies have shown that biomolecular NMR structures are often of lower quality when compared to crystal structures, and consequently they are often excluded from structural analyses. We present a publicly available database of re-refined NMR structures, exhibiting significantly improved quality. This database (available at http://www.cmbi.kun.nl/dress/) presents a uniformly refined and validated set of structural models that improves the value of these NMR structures as input for experimental and theoretical studies in many fields of research.     

Improved anticancer potency by head‐to‐tail cyclization of short cationic anticancer peptides containing a lipophilic β2,2‐amino acid

We have recently reported a series of synthetic anticancer heptapeptides (H‐KKWβ^2,2WKK‐NH₂) containing a central achiral and lipophilic β^2,2‐amino acid that display low toxicity against non‐malignant cells and high proteolytic stability. In the present study, we have further investigated the effects of increasing the rigidity and amphipathicity of two of our lead heptapeptides by preparing a series of seven to five residue cyclic peptides containing the two most promising β^2,2‐amino acid derivatives as part of the central lipophilic core. The peptides were tested for anticancer activity against human Burkitt's lymphoma (Ramos cells), haemolytic activity against human red blood cells (RBC) and cytotoxicity against healthy human lung fibroblast cells (MRC‐5). The results demonstrated a considerable increase in anticancer potency following head‐to‐tail peptide cyclization, especially for the shortest derivatives lacking a tryptophan residue. High‐resolution NMR studies and molecular dynamics simulations together with an annexin‐V‐FITC and propidium iodide fluorescent assay showed that the peptides had a membrane disruptive mode of action and that the more potent peptides penetrated deeper into the lipid bilayer. The need for new anticancer drugs with novel modes of action is demanding, and development of short cyclic anticancer peptides with an overall rigidified and amphipathic structure is a promising approach to new anticancer agents. Copyright © 2012 European Peptide Society and John Wiley & Sons, Ltd.     

Structural and functional consequences of the presence of a fourth disulfide bridge in the scorpion short toxins: solution structure of the potassium channel inhibitor HsTX1

We have determined the three-dimensional structure of the potassium channel inhibitor HsTX1, using nuclear magnetic resonance and molecular modeling. This protein belongs to the scorpion short toxin family, which essentially contains potassium channel blockers of 29 to 39 amino acids and three disulfide bridges. It is highly active on voltage-gated Kv1.3 potassium channels. Furthermore, it has the particularity to possess a fourth disulfide bridge. We show that HsTX1 has a fold similar to that of the three-disulfide-bridged toxins and conserves the hydrophobic core found in the scorpion short toxins. Thus, the fourth bridge has no influence on the global conformation of HsTX1. Most residues spatially analogous to those interacting with voltage-gated potassium channels in the three-disulfide-bridged toxins are conserved in HsTX1. Thus, we propose that Tyr21, Lys23, Met25, and Asn26 are involved in the biological activity of HsTX1. As an additional positively charged residue is always spatially close to the aromatic residue in toxins blocking the voltage-gated potassium channels, and as previous mutagenesis experiments have shown the critical role played by the C-terminus in HsTX1, we suggest that Arg33 is also important for the activity of the four disulfide-bridged toxin. Docking calculations confirm that, if Lys23 and Met25 interact with the GYGDMH motif of Kv1.3, Arg33 can contact Asp386 and, thus, play the role of the additional positively charged residue of the toxin functional site. This original configuration of the binding site of HsTX1 for Kv1.3, if confirmed experimentally, offers new structural possibilities for the construction of a molecule blocking the voltage-gated potassium channels.     

β‐Hairpin folding and stability: molecular dynamics simulations of designed peptides in aqueous solution

The structural properties of a 10‐residue and a 15‐residue peptide in aqueous solution were investigated by molecular dynamics simulation. The two designed peptides, SYINSDGTWT and SESYINSDGTWTVTE, had been studied previously by NMR at 278 K and the resulting model structures were classified as 3:5 β‐hairpins with a type I + G1 β‐bulge turn. In simulations at 278 K, starting from the NMR model structure, the 3:5 β‐hairpin conformers proved to be stable over the time period evaluated (30 ns). Starting from an extended conformation, simulations of the decapeptide at 278 K, 323 K and 353 K were also performed to study folding. Over the relatively short time scales explored (30 ns at 278 K and 323 K, 56 ns at 353 K), folding to the 3:5 β‐hairpin could only be observed at 353 K. At this temperature, the collapse to β‐hairpin‐like conformations is very fast. The conformational space accessible to the peptide is entirely dominated by loop structures with different degrees of β‐hairpin character. The transitions between different types of ordered loops and β‐hairpins occur through two unstructured loop conformations stabilized by a single side‐chain interaction between Tyr2 and Trp9, which facilitates the changes of the hydrogen‐bond register. In agreement with previous experimental results, β‐hairpin formation is initially driven by the bending propensity of the turn segment. Nevertheless, the fine organization of the turn region appears to be a late event in the folding process. Copyright © 2004 European Peptide Society and John Wiley & Sons, Ltd.     

Exploring the dynamic information content of a protein NMR structure: Comparison of a molecular dynamics simulation with the NMR and X‐ray structures of Escherichia coli ribonuclease HI

The multiconformer nature of solution nuclear magnetic resonance (NMR) structures of proteins results from the effects of intramolecular dynamics, spin diffusion and an uneven distribution of structural restraints throughout the molecule. A delineation of the former from the latter two contributions is attempted in this work for an ensemble of 15 NMR structures of the protein Escherichia coli ribonuclease HI (RNase HI). Exploration of the dynamic information content of the NMR ensemble is carried out through correlation with data from two crystal structures and a 1.7‐ns molecular dynamics (MD) trajectory of RNase HI in explicit solvent. Assessment of the consistency of the crystal and mean MD structures with nuclear Overhauser effect (NOE) data showed that the NMR ensemble is overall more compatible with the high‐resolution (1.48 Å) crystal structure than with either the lower‐resolution (2.05 Å) crystal structure or the MD simulation. Furthermore, the NMR ensemble is found to span more conformational space than the MD simulation for both the backbone and the sidechains of RNase HI. Nonetheless, the backbone conformational variability of both the NMR ensemble and the simulation is especially consistent with NMR relaxation measurements of two loop regions that are putative sites of substrate recognition. Plausible side‐chain dynamic information is extracted from the NMR ensemble on the basis of (i) rotamericity and syn‐pentane character of variable torsion angles, (ii) comparison of the magnitude of atomic mean‐square fluctuations (msf) with those deduced from crystallographic thermal factors, and (iii) comparison of torsion angle conformational behavior in the NMR ensemble and the simulation. Several heterogeneous torsion angles, while adopting non‐rotameric/syn‐pentane conformations in the NMR ensemble, exist in a unique conformation in the simulation and display low X‐ray thermal factors. These torsions are identified as sites whose variability is likely to be an artifact of the NMR structure determination procedure. A number of other torsions show a close correspondence between the conformations sampled in the NMR and MD ensembles, as well as significant correlations among crystallographic thermal factors and atomic msf calculated from the NMR ensemble and the simulation. These results indicate that a significant amount of dynamic information is contained in the NMR ensemble. The relevance of the present findings for the biological function of RNase HI, protein recognition studies, and previous investigations of the motional content of protein NMR structures are discussed. Proteins 1999;36:87–110. © 1999 Wiley‐Liss, Inc.     

High‐resolution structures of transient receptor potential vanilloid channels: Unveiling a functionally diverse group of ion channels

Transient receptor potential vanilloid (TRPV) channels are part of the superfamily of TRP ion channels and play important roles in widespread physiological processes including both neuronal and non‐neuronal pathways. Various diseases such as skeletal abnormalities, chronic pain, and cancer are associated with dysfunction of a TRPV channel. In order to obtain full understanding of disease pathogenesis and create opportunities for therapeutic intervention, it is essential to unravel how these channels function at a molecular level. In the past decade, incredible progress has been made in biochemical sample preparation of large membrane proteins and structural biology techniques, including cryo‐electron microscopy. This has resulted in high resolution structures of all TRPV channels, which has provided novel insights into the molecular mechanisms of channel gating and regulation that will be summarized in this review.     

Molecular simulations study of novel 1,4‐dihydropyridines derivatives with a high selectivity for Cav3.1 calcium channel

1,4‐Dihydropyridines (DHPs) have been developed to treat hypertension, angina, and nerve system disease. They are thought to mainly target the L‐type calcium channels, but low selectivity prompts them to block Cav1.2 and Cav3.1 channels simultaneously. Recently, some novel DHPs with different hydrophobic groups have been synthesized and among them M12 has a higher selectivity for Cav3.1. However, the structural information about Cav3.1‐DHPs complexes is not available in the experiment. Thus, we combined homology modeling, molecular docking, molecular dynamics simulations, and binding free energy calculations to quantitatively elucidate the inhibition mechanism of DHPs. The calculated results indicate that our model is in excellent agreement with experimental results. On the basis of conformational analysis, we identify the main interactions between DHPs and calcium channels and further elaborate on the different selectivity of ligands from the micro perspective. In conjunction with energy distribution, we propose that the binding sites of Cav3.1‐DHPs is characterized by several interspersed hydrophobic amino acid residues on the IIIS6 and IVS6 segments. We also speculate the favorable function groups on prospective DHPs. Besides, our model provides important information for further mutagenesis experiments.     

Coupling of activation and inactivation gate in a K+‐channel: potassium and ligand sensitivity

Potassium (K⁺)‐channel gating is choreographed by a complex interplay between external stimuli, K⁺ concentration and lipidic environment. We combined solid‐state NMR and electrophysiological experiments on a chimeric KcsA–Kv1.3 channel to delineate K⁺, pH and blocker effects on channel structure and function in a membrane setting. Our data show that pH‐induced activation is correlated with protonation of glutamate residues at or near the activation gate. Moreover, K⁺ and channel blockers distinctly affect the open probability of both the inactivation gate comprising the selectivity filter of the channel and the activation gate. The results indicate that the two gates are coupled and that effects of the permeant K⁺ ion on the inactivation gate modulate activation‐gate opening. Our data suggest a mechanism for controlling coordinated and sequential opening and closing of activation and inactivation gates in the K⁺‐channel pore.     

Stereochemistry of Schellman motifs in peptides: Crystal structure of a hexapeptide with a C‐terminus 6 → 1 hydrogen bond

The Schellman motif is a widely observed helix terminating structural motif in proteins, which is generated when the C‐terminus residue adopts a left‐handed helical (α_L) conformation. The resulting hydrogen‐bonding pattern involves the formation of an intramolecular 6 → 1 interaction. This helix terminating motif is readily mimicked in synthetic helical peptides by placing an achiral residue at the penultimate position of the sequence. Thus far, the Schellman motif has been characterized crystallographically only in peptide helices of length 7 residues or greater. The structure of the hexapeptide Boc–Pro–Aib–Gly–Leu–Aib–Leu–OMe in crystals reveal a short helical stretch terminated by a Schellman motif, with the formation of 6 → 1 C‐terminus hydrogen bond. The crystals are in the space group P2₁2₁2₁ with a = 18.155(3) Å, b = 18.864(8) Å, c = 11.834(4) Å, and Z = 4 . The final R₁ and wR₂ values are 7.68 and 14.6%, respectively , for 1524 observed reflections [F_o ≥ 3ς(F_o)]. A 6 → 1 hydrogen bond between Pro(1)CO · · · Leu(6)NH and a 5 → 2 hydrogen bond between Aib(2)CO · · · Aib(5)NH are observed. An analysis of the available oligopeptides having an achiral Aib residue at the penultimate position suggests that chain length and sequence effects may be the other determining factors in formation of Schellman motifs. © 1999 John Wiley & Sons, Inc. Biopoly 50: 13–22, 1999