Hampering the early aggregation of PrP-E200K protein by charge-based inhibitors: a computational study

AbstractA multilayered computational workflow was designed to identify a druggable binding site on the surface of the E200K pathogenic mutant of the human prion protein, and to investigate the effect of the binding of small molecules in the inhibition of the early aggregation of this protein. At this purpose, we developed an efficient computational tool to scan the molecular interaction properties of a whole MD trajectory, thus leading to the characterization of plausible binding regions on the surface of PrP-E200K. These structural data were then employed to drive structure-based virtual screening and fragment-based approaches to the seeking of small molecular binders of the PrP-E200K. Six promising compounds were identified, and their binding stabilities were assessed by MD simulations. Therefore, analyses of the molecular electrostatic potential similarity between the bound complexes and unbound protein evidenced their potential activity as charged-based inhibitors of the PrP-E200K early aggregation.

Download Full-text

Biophysical characterization of the SARS-CoV-2 E protein

10.1101/2021.05.28.446179 ◽

2021 ◽

Author(s):

Aujan Mehregan ◽

Sergio Perez-Conesa ◽

Yuxuan Zhuang ◽

Ahmad Elbahnsi ◽

Diletta Pasini ◽

...

Keyword(s):

Recombinant Expression ◽

Structural Data ◽

Md Simulations ◽

Full Length ◽

Coarse Grained ◽

Biophysical Characterization ◽

Viral Budding ◽

Nmr Structures ◽

E Protein

SARS-CoV-2 is the virus responsible for the COVID-19 pandemic which continues to wreak havoc across the world, over a year and a half after its effects were first reported in the general media. Current fundamental research efforts largely focus on the SARS-CoV-2 Spike protein. Since successful antiviral therapies are likely to target multiple viral components, there is considerable interest in understanding the biophysical role of its other proteins, in particular structural membrane proteins. Here, we have focused our efforts on the characterization of the full-length E protein from SARS-CoV-2, combining experimental and computational approaches. Recombinant expression of the full-length E protein from SARS-CoV-2 reveals that this membrane protein is capable of independent multimerization, possibly as a tetrameric or smaller species. Fluorescence microscopy shows that the protein localizes intracellularly, and coarse-grained MD simulations indicate it causes bending of the surrounding lipid bilayer, corroborating a potential role for the E protein in viral budding. Although we did not find robust electrophysiological evidence of ion-channel activity, cells transfected with the E protein exhibited reduced intracellular Ca2+, which may further promote viral replication. However, our atomistic MD simulations revealed that previous NMR structures are relatively unstable, and result in models incapable of ion conduction. Our study highlights the importance of using high-resolution structural data obtained from a full-length protein to gain detailed molecular insights, and eventually permitting virtual drug screening.

Download Full-text

ConfID: an analytical method for conformational characterization of small molecules using molecular dynamics trajectories

Bioinformatics ◽

10.1093/bioinformatics/btaa130 ◽

2020 ◽

Vol 36 (11) ◽

pp. 3576-3577 ◽

Cited By ~ 1

Author(s):

Marcelo D Polêto ◽

Bruno I Grisci ◽

Marcio Dorn ◽

Hugo Verli

Keyword(s):

Molecular Dynamics ◽

Small Molecules ◽

Conformational Transition ◽

Clustering Algorithms ◽

Md Simulations ◽

Conformational Space ◽

Supplementary Information ◽

Conformational Sampling ◽

Conformational Population

Abstract Motivation The conformational space of small molecules can be vast and difficult to assess. Molecular dynamics (MD) simulations of free ligands in solution have been applied to predict conformational populations, but their characterization is often based on clustering algorithms or manual efforts. Results Here, we introduce ConfID, an analytical tool for conformational characterization of small molecules using MD trajectories. The evolution of conformational sampling and population frequencies throughout trajectories is calculated to check for sampling convergence while allowing to map relevant conformational transitions. The tool is designed to track conformational transition events and calculate time-dependent properties for each conformational population detected. Availability and implementation Toolkit and documentation are freely available at http://sbcb.inf.ufrgs.br/confid Contact [email protected] or [email protected] Supplementary information Supplementary data are available at Bioinformatics online.

Download Full-text

The glycomes of Caenorhabditis elegans and other model organisms

Biochemical Society Symposium ◽

10.1042/bss0690117 ◽

2002 ◽

Vol 69 ◽

pp. 117-134 ◽

Cited By ~ 47

Author(s):

Stuart M. Haslam ◽

David Gems ◽

Howard R. Morris ◽

Anne Dell

Keyword(s):

Caenorhabditis Elegans ◽

Current Knowledge ◽

Structural Data ◽

Model Organisms ◽

Entire Genome ◽

C Elegans ◽

The Face ◽

Area Of Concern ◽

Nematode Worm

There is no doubt that the immense amount of information that is being generated by the initial sequencing and secondary interrogation of various genomes will change the face of glycobiological research. However, a major area of concern is that detailed structural knowledge of the ultimate products of genes that are identified as being involved in glycoconjugate biosynthesis is still limited. This is illustrated clearly by the nematode worm Caenorhabditis elegans, which was the first multicellular organism to have its entire genome sequenced. To date, only limited structural data on the glycosylated molecules of this organism have been reported. Our laboratory is addressing this problem by performing detailed MS structural characterization of the N-linked glycans of C. elegans; high-mannose structures dominate, with only minor amounts of complex-type structures. Novel, highly fucosylated truncated structures are also present which are difucosylated on the proximal N-acetylglucosamine of the chitobiose core as well as containing unusual Fucα1–2Gal1–2Man as peripheral structures. The implications of these results in terms of the identification of ligands for genomically predicted lectins and potential glycosyltransferases are discussed in this chapter. Current knowledge on the glycomes of other model organisms such as Dictyostelium discoideum, Saccharomyces cerevisiae and Drosophila melanogaster is also discussed briefly.

Download Full-text

Effect of Ionic Strength on the Aggregation Propensity of Aβ1-42 Peptide: An In-silico Study

Current Chemical Biology ◽

10.2174/2212796814999200818103157 ◽

2020 ◽

Vol 14 (3) ◽

pp. 216-226

Author(s):

Priyanka Borah ◽

Venkata S.K. Mattaparthi

Keyword(s):

Diffusion Coefficient ◽

Ionic Strength ◽

Marked Effect ◽

Computational Study ◽

Conformational Dynamics ◽

Rapid Change ◽

Md Simulations ◽

Radius Of Gyration ◽

Aggregation Propensity ◽

In Silico Study

Background: Aggregation of misfolded proteins under stress conditions in the cell might lead to several neurodegenerative disorders. Amyloid-beta (Aβ1-42) peptide, the causative agent of Alzheimer’s disease, has the propensity to fold into β-sheets under stress, forming aggregated amyloid plaques. This is influenced by factors such as pH, temperature, metal ions, mutation of residues, and ionic strength of the solution. There are several studies that have highlighted the importance of ionic strength in affecting the folding and aggregation propensity of Aβ1-42 peptide. Objective: To understand the effect of ionic strength of the solution on the aggregation propensity of Aβ1-42 peptide, using computational approaches. Materials and Methods: In this study, Molecular Dynamics (MD) simulations were performed on Aβ1-42 peptide monomer placed in (i) 0 M, (ii) 0.15 M, and (iii) 0.30 M concentration of NaCl solution. To prepare the input files for the MD simulations, we have used the Amberff99SB force field. The conformational dynamics of Aβ1-42 peptide monomer in different ionic strengths of the solutions were illustrated from the analysis of the corresponding MD trajectory using the CPPtraj tool. Results: From the MD trajectory analysis, we observe that with an increase in the ionic strength of the solution, Aβ1-42 peptide monomer shows a lesser tendency to undergo aggregation. From RMSD and SASA analysis, we noticed that Aβ1-42 peptide monomer undergoes a rapid change in conformation with an increase in the ionic strength of the solution. In addition, from the radius of gyration (Rg) analysis, we observed Aβ1-42 peptide monomer to be more compact at moderate ionic strength of the solution. Aβ1-42 peptide was also found to hold its helical secondary structure at moderate and higher ionic strengths of the solution. The diffusion coefficient of Aβ1-42 peptide monomer was also found to vary with the ionic strength of the solution. We observed a relatively higher diffusion coefficient value for Aβ1-42 peptide at moderate ionic strength of the solution. Conclusion: Our findings from this computational study highlight the marked effect of ionic strength of the solution on the conformational dynamics and aggregation propensity of Aβ1-42 peptide monomer.

Download Full-text

Water and Deuterium Oxide Permeability through Aquaporin 1: MD Predictions and Experimental Verification

The Journal of General Physiology ◽

10.1085/jgp.200709810 ◽

2007 ◽

Vol 130 (1) ◽

pp. 111-116 ◽

Cited By ~ 31

Author(s):

Artem B. Mamonov ◽

Rob D. Coalson ◽

Mark L. Zeidel ◽

John C. Mathai

Keyword(s):

Deuterium Oxide ◽

Diffusion Constant ◽

Structural Data ◽

Md Simulations ◽

Water Model ◽

Aquaporin 1 ◽

Self Diffusion ◽

Osmotic Permeability ◽

Protein Channels ◽

And Diffusion

Determining the mechanisms of flux through protein channels requires a combination of structural data, permeability measurement, and molecular dynamics (MD) simulations. To further clarify the mechanism of flux through aquaporin 1 (AQP1), osmotic pf (cm3/s/pore) and diffusion pd (cm3/s/pore) permeability coefficients per pore of H2O and D2O in AQP1 were calculated using MD simulations. We then compared the simulation results with experimental measurements of the osmotic AQP1 permeabilities of H2O and D2O. In this manner we evaluated the ability of MD simulations to predict actual flux results. For the MD simulations, the force field parameters of the D2O model were reparameterized from the TIP3P water model to reproduce the experimentally observed difference in the bulk self diffusion constants of H2O vs. D2O. Two MD systems (one for each solvent) were constructed, each containing explicit palmitoyl-oleoyl-phosphatidyl-ethanolamine (POPE) phospholipid molecules, solvent, and AQP1. It was found that the calculated value of pf for D2O is ∼15% smaller than for H2O. Bovine AQP1 was reconstituted into palmitoyl-oleoyl-phosphatidylcholine (POPC) liposomes, and it was found that the measured macroscopic osmotic permeability coefficient Pf (cm/s) of D2O is ∼21% lower than for H2O. The combined computational and experimental results suggest that deuterium oxide permeability through AQP1 is similar to that of water. The slightly lower observed osmotic permeability of D2O compared to H2O in AQP1 is most likely due to the lower self diffusion constant of D2O.

Download Full-text

Database of ab initio L-edge X-ray absorption near edge structure

Scientific Data ◽

10.1038/s41597-021-00936-5 ◽

2021 ◽

Vol 8 (1) ◽

Author(s):

Yiming Chen ◽

Chi Chen ◽

Chen Zheng ◽

Shyam Dwaraknath ◽

Matthew K. Horton ◽

...

Keyword(s):

High Throughput ◽

Scattering Theory ◽

Research Community ◽

Metal Compounds ◽

Edge Structure ◽

X Ray ◽

Initial Release ◽

Computational Workflow ◽

X Ray Absorption

AbstractThe L-edge X-ray Absorption Near Edge Structure (XANES) is widely used in the characterization of transition metal compounds. Here, we report the development of a database of computed L-edge XANES using the multiple scattering theory-based FEFF9 code. The initial release of the database contains more than 140,000 L-edge spectra for more than 22,000 structures generated using a high-throughput computational workflow. The data is disseminated through the Materials Project and addresses a critical need for L-edge XANES spectra among the research community.

Download Full-text

Computational study for suppression of CD25/IL-2 interaction

Biological Chemistry ◽

10.1515/hsz-2020-0326 ◽

2020 ◽

Vol 0 (0) ◽

Author(s):

Moein Dehbashi ◽

Zohreh Hojati ◽

Majid Motovali-bashi ◽

Mazdak Ganjalikhani-Hakemi ◽

Akihiro Shimosaka ◽

...

Keyword(s):

Small Molecules ◽

Cancer Progression ◽

Immune Escape ◽

De Novo ◽

Computational Study ◽

Treg Cells ◽

Homology Search ◽

Small Interfering Rnas

AbstractCancer recurrence presents a huge challenge in cancer patient management. Immune escape is a key mechanism of cancer progression and metastatic dissemination. CD25 is expressed in regulatory T (Treg) cells including tumor-infiltrating Treg cells (TI-Tregs). These cells specially activate and reinforce immune escape mechanism of cancers. The suppression of CD25/IL-2 interaction would be useful against Treg cells activation and ultimately immune escape of cancer. Here, software, web servers and databases were used, at which in silico designed small interfering RNAs (siRNAs), de novo designed peptides and virtual screened small molecules against CD25 were introduced for the prospect of eliminating cancer immune escape and obtaining successful treatment. We obtained siRNAs with low off-target effects. Further, small molecules based on the binding homology search in ligand and receptor similarity were introduced. Finally, the critical amino acids on CD25 were targeted by a de novo designed peptide with disulfide bond. Hence we introduced computational-based antagonists to lay a foundation for further in vitro and in vivo studies.

Download Full-text