PEP-FOLD: an online resource for de novo peptide structure prediction

Rational peptide design and large-scale prediction of peptide structure from sequence remain a challenge for chemical biologists. We present PEP-FOLD, an online service, aimed at de novo modelling of 3D conformations for peptides between 9 and 25 amino acids in aqueous solution. Using a hidden Markov model-derived structural alphabet (SA) of 27 four-residue letters, PEP-FOLD first predicts the SA letter profiles from the amino acid sequence and then assembles the predicted fragments by a greedy procedure driven by a modified version of the OPEP coarse-grained force field. Starting from an amino acid sequence, PEP-FOLD performs series of 50 simulations and returns the most representative conformations identified in terms of energy and population. Using a benchmark of 25 peptides with 9–23 amino acids, and considering the reproducibility of the runs, we find that, on average, PEP-FOLD locates lowest energy conformations differing by 2.6 Å Cα root mean square deviation from the full NMR structures. PEP-FOLD can be accessed at http://bioserv.rpbs.univ-paris-diderot.fr/PEP-FOL

Vasopressin V2R-Targeting Peptide Carrier Mediates siRNA Delivery into Collecting Duct Cells

Internalization of receptor proteins after interacting with specific ligands has been proposed to facilitate siRNA delivery into the target cells via receptor-mediated siRNA transduction. In this study, we demonstrated a novel method of vasopressin V2 receptor (V2R)-mediated siRNA delivery against AQP2 in primary cultured inner medullary collecting duct (IMCD) cells of rat kidney. We synthesized the dDAVP conjugated with nine D-arginines (dDAVP-9r) as a peptide carrier for siRNA delivery. The structure of synthetic peptide carrier showed two regions (i.e., ligand domain to V2R (dDAVP) and siRNA carrying domain (nine D-arginine)) bisected with a spacer of four glycines. The results revealed that 1) synthesized dDAVP-9r peptides formed a stable polyplex with siRNA; 2) siRNA/dDAVP-9r polyplex could bind to the V2R of IMCD cells and induced AQP2 phosphorylation (Ser 256); 3) siRNA/dDAVP-9r polyplex was stable in response to the wide range of different osmolalities, pH levels, or to the RNases; 4) fluorescein-labeled siRNA was delivered into V2R-expressing MDCK and LLC-PK1 cells by siRNA/dDAVP-9r polyplex, but not into the V2R-negative Cos-7 cells; and 5) AQP2-siRNA/dDAVP-9r polyplex effectively delivered siRNA into the IMCD cells, resulting in the significant decrease of protein abundance of AQP2, but not AQP4. Therefore, for the first time to our knowledge, we demonstrated that V2R-mediated siRNA delivery could be exploited to deliver specific siRNA to regulate abnormal expression of target proteins in V2R-expressing kidney cells. The methods could be potentially used in vivo to regulate abnormal expression of proteins associated with disease conditions in the V2R-expressing kidney cells

An innovative ethosuximide granule formulation designed for pediatric use: Comparative pharmacokinetics, safety, tolerability, and palatability profile versus reference syrup.

Ethosuximide, the first-line therapy for childhood absence epilepsy, is currently formulated as a syrup (Zarontin®, Pfizer) with a bitter taste and high sugar content, poorly adapted to children, and a ketogenic diet. The collaborative European FP7 project KIEKIDS aimed at developing an innovative sugar-free, tasteless formulation convenient for pediatric use. This dual Phase-I study evaluated two granule formulations based on lipid multiparticulate (LMP) technology. Two panels of 6 healthy adult volunteers underwent a randomized, placebo-controlled, partly blinded, 3-way cross-over trial, comparing ethosuximide granules A or B with placebo granules and syrup at single 10 mg/kg doses. Corresponding plasma pharmacokinetic profiles of ethosuximide were compared, along with palatability, safety, and tolerability. The LMP granule A proved suboptimal due to bitterness and adherence to beaker walls, while the optimized granule B revealed excellent palatability, similar to placebo granules, and low adherence to glass. The relative bioavailability of granules A versus syrup, based on dose-normalized C <sub>max</sub> and AUC <sub>0-∞</sub> was 93.7% [90% CI: 76.3-115.1] and 96.1% [91.0-101.5], respectively. For granules B it was 87.6% [81.6-94.0] and 92.5% [88.5-96.6], respectively, with slightly delayed t <sub>max</sub> of 0.75 h [0.5-4.05] compared to syrup 0.5 h [0.3-0.8]. Tolerability visual analog scales revealed a trend for statistically non-significant improvement versus syrup at peak (30 min) for transient dizziness (both granules), fatigue (granules A), and anxiety (granules B). The innovative ethosuximide granule formulation B achieves a suitable profile for pediatric use, being sugar-free, tasteless, bioequivalent, and well-tolerated while enabling precise adjustment to body weight

fpocket: online tools for protein ensemble pocket detection and tracking

Computational small-molecule binding site detection has several important applications in the biomedical field. Notable interests are the identification of cavities for structure-based drug discovery or functional annotation of structures. fpocket is a small-molecule pocket detection program, relying on the geometric α-sphere theory. The fpocket web server allows: (i) candidate pocket detection—fpocket; (ii) pocket tracking during molecular dynamics, in order to provide insights into pocket dynamics—mdpocket; and (iii) a transposition of mdpocket to the combined analysis of homologous structures—hpocket. These complementary online tools allow to tackle various questions related to the identification and annotation of functional and allosteric sites, transient pockets and pocket preservation within evolution of structural families. The server and documentation are freely available at http://bioserv.rpbs.univ-paris-diderot.fr/fpocket

The Crystal Structure of PPIL1 Bound to Cyclosporine A Suggests a Binding Mode for a Linear Epitope of the SKIP Protein

BACKGROUND: The removal of introns from pre-mRNA is carried out by a large macromolecular machine called the spliceosome. The peptidyl-prolyl cis/trans isomerase PPIL1 is a component of the human spliceosome and binds to the spliceosomal SKIP protein via a binding site distinct from its active site. PRINCIPAL FINDINGS: Here, we have studied the PPIL1 protein and its interaction with SKIP biochemically and by X-ray crystallography. A minimal linear binding epitope derived from the SKIP protein could be determined using a peptide array. A 36-residue region of SKIP centred on an eight-residue epitope suffices to bind PPIL1 in pull-down experiments. The crystal structure of PPIL1 in complex with the inhibitor cyclosporine A (CsA) was obtained at a resolution of 1.15 A and exhibited two bound Cd(2+) ions that enabled SAD phasing. PPIL1 residues that have previously been implicated in binding of SKIP are involved in the coordination of Cd(2+) ions in the present crystal structure. Employing the present crystal structure, the determined minimal binding epitope and previously published NMR data, a molecular docking study was performed. In the docked model of the PPIL1.SKIP interaction, a proline residue of SKIP is buried in a hydrophobic pocket of PPIL1. This hydrophobic contact is encircled by several hydrogen bonds between the SKIP peptide and PPIL1. CONCLUSION: We characterized a short, linear epitope of SKIP that is sufficient to bind the PPIL1 protein. Our data indicate that this SKIP peptide could function in recruiting PPIL1 into the core of the spliceosome. We present a molecular model for the binding mode of SKIP to PPIL1 which emphasizes the versatility of cyclophilin-type PPIases to engage in additional interactions with other proteins apart from active site contacts despite their limited surface area

Multiscale Coarse-Graining of the Protein Energy Landscape

A variety of coarse-grained (CG) models exists for simulation of proteins. An outstanding problem is the construction of a CG model with physically accurate conformational energetics rivaling all-atom force fields. In the present work, atomistic simulations of peptide folding and aggregation equilibria are force-matched using multiscale coarse-graining to develop and test a CG interaction potential of general utility for the simulation of proteins of arbitrary sequence. The reduced representation relies on multiple interaction sites to maintain the anisotropic packing and polarity of individual sidechains. CG energy landscapes computed from replica exchange simulations of the folding of Trpzip, Trp-cage and adenylate kinase resemble those of other reduced representations; non-native structures are observed with energies similar to those of the native state. The artifactual stabilization of misfolded states implies that non-native interactions play a deciding role in deviations from ideal funnel-like cooperative folding. The role of surface tension, backbone hydrogen bonding and the smooth pairwise CG landscape is discussed. Ab initio folding aside, the improved treatment of sidechain rotamers results in stability of the native state in constant temperature simulations of Trpzip, Trp-cage, and the open to closed conformational transition of adenylate kinase, illustrating the potential value of the CG force field for simulating protein complexes and transitions between well-defined structural states

Lichen response to ammonia deposition defines the footprint of a penguin rookery

Ammonia volatilized from penguin rookeries is a major nitrogen source in Antarctic coastal terrestrial ecosystems. However, the spatial extent of ammonia dispersion from rookeries and its impacts have not been quantified previously. We measured ammonia concentration in air and lichen ecophysiological response variables proximate to an Adèlie penguin rookery at Cape Hallett, northern Victoria Land. Ammonia emitted from the rookery was 15N-enriched (δ15N value +6.9) and concentrations in air ranged from 36–75 µg m−3 at the rookery centre to 0.05 µg m−3 at a distance of 15.3 km. δ15N values and rates of phosphomonoesterase (PME) activity in the lichens Usnea sphacelata and Umbilicaria decussata were strongly negatively related to distance from the rookery and PME activity was positively related to thallus N:P mass ratio. In contrast, the lichen Xanthomendoza borealis, which is largely restricted to within an area 0.5 km from the rookery perimeter, had high N, P and 15N concentrations but low PME activity suggesting that nutrient scavenging capacity is suppressed in highly eutrophicated sites. An ammonia dispersion model indicates that ammonia concentrations sufficient to significantly elevate PME activity and δ15N values (≥0.1 µg NH3 m−3) occurred over c. 40–300 km2 surrounding the rookery suggesting that penguin rookeries potentially can generate large spatial impact zones. In a general linear model NH3 concentration and lichen species identity were found to account for 72 % of variation in the putative proportion of lichen thallus N originating from penguin derived NH3. The results provide evidence of large scale impact of N transfer from a marine to an N-limited terrestrial ecosystem

In Silico Prediction and Analysis of Caenorhabditis EF-hand Containing Proteins

Calcium (Ca+2) is a ubiquitous messenger in eukaryotes including Caenorhabditis. Ca+2-mediated signalling processes are usually carried out through well characterized proteins like calmodulin (CaM) and other Ca+2 binding proteins (CaBP). These proteins interact with different targets and activate it by bringing conformational changes. Majority of the EF-hand proteins in Caenorhabditis contain Ca+2 binding motifs. Here, we have performed homology modelling of CaM-like proteins using the crystal structure of Drosophila melanogaster CaM as a template. Molecular docking was applied to explore the binding mechanism of CaM-like proteins and IQ1 motif which is a ∼25 residues and conform to the consensus sequence (I, L, V)QXXXRXXXX(R,K) to serve as a binding site for different EF hand proteins. We made an attempt to identify all the EF-hand (a helix-loop-helix structure characterized by a 12 residues loop sequence involved in metal coordination) containing proteins and their Ca+2 binding affinity in Caenorhabditis by analysing the complete genome sequence. Docking studies revealed that F165, F169, L29, E33, F44, L57, M61, M96, M97, M108, G65, V115, F93, N104, E144 of CaM-like protein is involved in the interaction with IQ1 motif. A maximum of 170 EF-hand proteins and 39 non-EF-hand proteins with Ca+2/metal binding motif were identified. Diverse proteins including enzyme, transcription, translation and large number of unknown proteins have one or more putative EF-hands. Phylogenetic analysis revealed seven major classes/groups that contain some families of proteins. Various domains that we identified in the EF-hand proteins (uncharacterized) would help in elucidating their functions. It is the first report of its kind where calcium binding loop sequences of EF-hand proteins were analyzed to decipher their calcium affinities. Variation in Ca+2-binding affinity of EF-hand CaBP could be further used to study the behaviour of these proteins. Our analyses postulated that Ca+2 is likely to be key player in Caenorhabditis cell signalling

Diffusion, Crowding & Protein Stability in a Dynamic Molecular Model of the Bacterial Cytoplasm

A longstanding question in molecular biology is the extent to which the behavior of macromolecules observed in vitro accurately reflects their behavior in vivo. A number of sophisticated experimental techniques now allow the behavior of individual types of macromolecule to be studied directly in vivo; none, however, allow a wide range of molecule types to be observed simultaneously. In order to tackle this issue we have adopted a computational perspective, and, having selected the model prokaryote Escherichia coli as a test system, have assembled an atomically detailed model of its cytoplasmic environment that includes 50 of the most abundant types of macromolecules at experimentally measured concentrations. Brownian dynamics (BD) simulations of the cytoplasm model have been calibrated to reproduce the translational diffusion coefficients of Green Fluorescent Protein (GFP) observed in vivo, and “snapshots” of the simulation trajectories have been used to compute the cytoplasm's effects on the thermodynamics of protein folding, association and aggregation events. The simulation model successfully describes the relative thermodynamic stabilities of proteins measured in E. coli, and shows that effects additional to the commonly cited “crowding” effect must be included in attempts to understand macromolecular behavior in vivo