What limits the primary sequence space of natural proteins?

The determinants of the rarity of nucleic and peptide short sequences in nature

The prevalence of nucleic and peptide short sequences across organismal genomes and proteomes has not been thoroughly investigated. We examined 45 785 reference genomes and 21 871 reference proteomes, spanning archaea, bacteria, eukaryotes and viruses to calculate the rarity of short sequences in them. To capture this, we developed a metric of the rarity of each sequence in nature, the rarity index. We find that the frequency of certain dipeptides in rare oligopeptide sequences is hundreds of times lower than expected, which is not the case for any dinucleotides. We also generate predictive regression models that infer the rarity of nucleic and proteomic sequences across nature or within each domain of life and viruses separately. When examining each of the three domains of life and viruses separately, the R² performance of the model predicting rarity for 5-mer peptides from mono- and dipeptides ranged between 0.814 and 0.932. A separate model predicting rarity for 10-mer oligonucleotides from mono- and dinucleotides achieved R² performance between 0.408 and 0.606. Our results indicate that the mono- and dinucleotide composition of nucleic sequences and the mono- and dipeptide composition of peptide sequences can explain a significant proportion of the variance in their frequencies in nature.

In-silico Investigations for the Identification of Novel Inhibitors Targeting Hepatitis C Virus RNA-dependent RNA Polymerase

BACKGROUND

Hepatitis C is an inflammatory condition of the liver caused by the hepatitis C virus, exhibiting acute and chronic manifestations with severity ranging from mild to severe and lifelong illnesses leading to liver cirrhosis and cancer. According to the World Health Organization's global estimates, a population of about 58 million have chronic hepatitis C virus infection, with around 1.5 million new infections occurring every year.

OBJECTIVE

The present study aimed to identify novel molecules targeting the Hepatitis C viral RNA Dependent RNA polymerases, which play a crucial role in genome replication, mRNA synthesis, etc. Methods: Structure-based virtual screening of chemical libraries of small molecules was done using AutoDock/Vina. The top-ranking pose for every ligand was complexed with the protein and used for further protein-ligand interaction analysis using the Protein-ligand interaction Profiler. Molecules from virtual screening were further assessed using the pkCSM web server. The proteinligand interactions were further subjected to molecular dynamics simulation studies to establish dynamic stability.

RESULTS

Molecular docking-based virtual screening of the database of small molecules, followed by screening based on pharmacokinetic and toxicity parameters, yielded eight probable RNA Dependent RNA polymerase inhibitors. The docking scores for the proposed candidates ranged from - 8.04 to -9.10 kcal/mol. The potential stability of the ligands bound to the target protein was demonstrated by molecular dynamics simulation studies.

CONCLUSION

Data from exhaustive computational studies proposed eight molecules as potential anti-viral candidates, targeting Hepatitis C viral RNA Dependent RNA polymerases, which can be further evaluated for their biological potential.

Secretory Conservation in Insulin Producing Cells: Is There a System-Level Law of Mass Action in Biology?

Altered secretion of insulin from pancreatic β-cells can manifest into disorders. For example, a lack of endogenously produced and/or secreted insulin results in Type 1 diabetes (and other associated subtypes). Pancreatic β-cells are the endocrine secretory cells that promote insulin secretion in response to glucose stimulation. Secretion in response to extracellular triggers is an interplay among various signaling pathways, transcription factors, and molecular mechanisms. The Mouse Insulinoma 6 (MIN6) cell line serves as a model system for gaining mechanistic insights into pancreatic β-cell functions. It is obvious that higher glucose consumption and increased insulin secretion are correlated. However, it has been reported that intracellular ATP levels remain ∼ constant beyond the extracellular glucose (EG) concentration of 10 mM. Therefore, any cause-effect relationship between glucose consumption (GC) and enhanced insulin secretion (eIS) remains unclear. We also found that total cellular protein, as well as total protein content in the culture "supernatant," remains constant regardless of varying EG concentrations. This indicated that eIS may be at the cost of (a) intracellular synthesis of other proteins and (b) secretion of other secretory proteins, or both (a) and (b), somehow coupled with GC by cells. To gain insights into the above, we carried out a transcriptome study of MIN6 cells exposed to hypoglycemic (HoG = 2.8 mM EG) and hyperglycemic (HyG = 25 mM EG) conditions. Expression of transcripts was analyzed in terms of Fragments Per Kilobase of transcript per Million mapped reads and Transcripts Per Million (FPKM and TPM) as well as values obtained by normalizing w.r.t. "∑(FPKM)" and "∑(TPM)." We report that HyG extracellular conditions lead to an ∼2-fold increase in insulin secretion compared to HoG measured by the enzyme-linked immunosorbent assay (ELISA) and transcripts of secreted proteins as well as their isoforms decreased in HyG conditions compared to HoG. Our results show for the first time that eIS in HyG conditions is at the cost of reduced transcription of other secreted proteins and is coupled with higher GC. The higher GC at increased extracellular glucose also indicates a yet undiscovered role of glucose molecules enhancing insulin secretion, since ATP levels resulting from glucose metabolism have been reported to be constant above an EG concentration of 10 mM. While extrapolation of our results to clinical implications is ambitious at best, this work reports novel cellular level aspects that seem relevant in some clinical observations pertaining to Type 1 diabetes. In addition, the conservatory nature of cellular secretions in insulin-secreting cells, discovered here, may be a general feature in cell biology.

Pseudopeptosomes: non-lipidated vesicular assemblies from bispidine-appended pseudopeptides

We report a novel molecular topology-based approach for creating reproducible vesicular assemblies in different solvent environments (including aqueous) using specifically designed pseudopeptides. Deviating from the classical "polar head group and hydrophobic tail" model of amphiphiles, we showed (reversible) self-assembly of synthesized pseudopeptides into vesicles. Naming these new type/class of vesicles "pseudopetosomes", we characterized them by high-resolution microscopy (scanning electron, transmission electron, atomic force, epifluorescence and confocal) along with dynamic light scattering. While accounting for hydropathy index of the constituent amino acids (side chains) of pseudopeptides, we probed molecular interactions, resulting in assembly of pseudopeptosomes by spectroscopy (fourier-transform infrared and fluorescence). Molecular characterization by X-ray crystallography and circular dichroism revealed "tryptophan (Trp)-Zip" arrangements and/or hydrogen-bonded one-dimensional assembly depending on specific pseudopeptides and solvent environments. Our data indicated that pseudopeptosomes are formed in solutions by self-assembly of bispidine pseudopeptides (of Trp, leucine and alanine amino-acid constituents) into sheets that transform into vesicular structures. Thus, we showed that assembly of pseudopeptosomes utilizes the full spectrum of all four weak interactions essential in biological systems. Our findings have direct implications in chemical and synthetic biology, but may also provide a new avenue of investigations on origins of life via pseudopeptosome-like assemblies. We also showed that these designer peptides can act as carriers for cellular transport.

Aspects of Biological Replication and Evolution Independent of the Central Dogma: Insights from Protein-Free Vesicular Transformations and Protein-Mediated Membrane Remodeling

Biological membrane remodeling is central to living systems. In spite of serving as “containers” of whole-living systems and functioning as dynamic compartments within living systems, biological membranes still find a “blue collar” treatment compared to the “white collar” nucleic acids and proteins in biology. This may be attributable to the fact that scientific literature on biological membrane remodeling is only 50 years old compared to ~ 150 years of literature on proteins and a little less than 100 years on nucleic acids. However, recently, evidence for symbiotic origins of eukaryotic cells from data only on biological membranes was reported. This, coupled with appreciation of reproducible amphiphilic self-assemblies in aqueous environments (mimicking replication), has already initiated discussions on origins of life beyond nucleic acids and proteins. This work presents a comprehensive compilation and meta-analyses of data on self-assembly and vesicular transformations in biological membranes—starting from model membranes to establishment of Influenza Hemagglutinin-mediated membrane fusion as a prototypical remodeling system to a thorough comparison between enveloped mammalian viruses and cellular vesicles. We show that viral membrane fusion proteins, in addition to obeying “stoichiometry-driven protein folding”, have tighter compositional constraints on their amino acid occurrences than general-structured proteins, regardless of type/class. From the perspective of vesicular assemblies and biological membrane remodeling (with and without proteins) we find that cellular vesicles are quite different from viruses. Finally, we propose that in addition to pre-existing thermodynamic frameworks, kinetic considerations in de novo formation of metastable membrane structures with available “third-party” constituents (including proteins) were not only crucial for origins of life but also continue to offer morphological replication and/or functional mechanisms in modern life forms, independent of the central dogma.

Design and synthesis of amino acid derivatives of substituted benzimidazoles and pyrazoles as Sirt1 inhibitors

Owing to its presence in several biological processes, Sirt1 acts as a potential therapeutic target for many diseases. Here, we report the structure-based designing and synthesis of two distinct series of novel Sirt1 inhibitors, benzimidazole mono-peptides and amino-acid derived 5-pyrazolyl methylidene rhodanine carboxylic acid. The compounds were evaluated for in vitro enzyme-based and cell-based Sirt1 inhibition assay, and cytotoxic-activity in both liver and breast cancer cells. The tryptophan conjugates i.e.13h (IC50 = 0.66 μM, ΔG bind = -1.1 kcal mol-1) and 7d (IC50 = 0.77 μM, ΔG bind = -4.4 kcal mol-1) demonstrated the maximum efficacy to inhibit Sirt1. The MD simulation unveiled that electrostatic complementarity at the substrate-binding-site through a novel motif "SLxVxP(V/F)A" could be a cause of increased Sirt1 inhibition by 13h and 13l over Sirt2 in cell-based assay, as compared to the control Ex527 and 7d. Finally, this study highlights novel molecules 7d and 13h, along with a new key hot-spot in Sirt1, which could be used as a starting lead to design more potent and selective sirtuin inhibitors as a potential anticancer molecule.

Virtual screening-driven drug discovery of SARS-CoV2 enzyme inhibitors targeting viral attachment, replication, post-translational modification and host immunity evasion infection mechanisms

The novel coronavirus SARS-CoV2, the causative agent of the pandemic disease COVID-19, emerged in December 2019 forcing lockdown of communities in many countries. The absence of specific drugs and vaccines, the rapid transmission of the virus, and the increasing number of deaths worldwide necessitated the discovery of new substances for anti-COVID-19 drug development. With the aid of bioinformatics and computational modelling, ninety seven antiviral secondary metabolites from fungi were docked onto five SARS-CoV2 enzymes involved in viral attachment, replication, post-translational modification, and host immunity evasion infection mechanisms followed by molecular dynamics simulation and in silico ADMET prediction (absorption, distribution, metabolism, excretion and toxicity) of the hit compounds. Thus, three fumiquinazoline alkaloids scedapin C (15), quinadoline B (19) and norquinadoline A (20), the polyketide isochaetochromin D1 (8), and the terpenoid 11a-dehydroxyisoterreulactone A (11) exhibited high binding affinities on the target proteins, papain-like protease (PLpro), chymotrypsin-like protease (3CLpro), RNA-directed RNA polymerase (RdRp), non-structural protein 15 (nsp15), and the spike binding domain to GRP78. Molecular dynamics simulation was performed to optimize the interaction and investigate the stability of the top-scoring ligands in complex with the five target proteins. All tested complexes were found to have dynamic stability. Of the five top-scoring metabolites, quinadoline B (19) was predicted to confer favorable ADMET values, high gastrointestinal absorptive probability and poor blood-brain barrier crossing capacities.Communicated by Ramaswamy H. Sarma.

Identification of a novel dual-target scaffold for 3CLpro and RdRp proteins of SARS-CoV-2 using 3D-similarity search, molecular docking, molecular dynamics and ADMET evaluation

The new SARS-CoV-2 coronavirus is the causative agent of the COVID-19 pandemic outbreak that affected whole the world with more than 6 million confirmed cases and over 370,000 deaths. At present, there are no effective treatments or vaccine for this disease, which constitutes a serious global health crisis. As the pandemic still spreading around the globe, it is of interest to use computational methods to identify potential inhibitors for the virus. The crystallographic structures of 3CLpro (PDB: 6LU7) and RdRp (PDB 6ML7) were used in virtual screening of 50000 chemical compounds obtained from the CAS Antiviral COVID19 database using 3D-similarity search and standard molecular docking followed by ranking and selection of compounds based on their binding affinity, computational techniques for the sake of details on the binding interactions, absorption, distribution, metabolism, excretion, and toxicity prediction; we report three 4-(morpholin-4-yl)-1,3,5-triazin-2-amine derivatives; two compounds (2001083-68-5 and 2001083-69-6) with optimal binding features to the active site of the main protease and one compound (833463-19-7) with optimal binding features to the active site of the polymerase for further consideration to fight COVID-19. The structural stability and dynamics of lead compounds at the active site of 3CLpro and RdRp were examined using molecular dynamics (MD) simulation. Essential dynamics demonstrated that the three complexes remain stable during simulation of 20 ns, which may be suitable candidates for further experimental analysis. As the identified leads share the same scaffold, they may serve as promising leads in the development of dual 3CLpro and RdRp inhibitors against SARS-CoV-2.Communicated by Ramaswamy H. Sarma.

Structural disorder originates beyond narrow stoichiometric margins of amino acids in naturally occurring folded proteins

Rigorous analyses of Euclidean distances between non-peptide bonded residues in structures of several thousand naturally occurring folded proteins yielded a surprising "margin of life" for percentage occurrence of individual amino acids in naturally occurring folded proteins. On one hand, the concept of "margin of life", referring to lower than expected variances in average stoichiometric occurrences of individual amino acids in folded proteins, remains unchallenged since its discovery a decade ago. On the other hand, within this past decade there has been a strong emergence of a gradual paradigm shift in biology, from sequence-structure-function in proteins to sequence-disorder-function, fuelled by discoveries on functional implications of intrinsically disordered proteins (primary sequences that do not form stable structures). Thus the applicability of "margin of life" to peptide-bonded residues in all known natural proteins, adopting stable structures vis-à-vis intrinsically disordered needs to be explored. Therefore in this work, we analyze compositions of the complete naturally occurring primary sequence space (over 560000 sequences) after dividing it into mutually exclusive subsets of structured and intrinsically disordered proteins along with a subset without any structural information. While finding that occurrence of different peptides (up to pentapeptides) is a direct consequence of the relative occurrences of their constituting residues in folded proteins, we report that structural disorder in natural proteins originates beyond the narrow stoichiometric margins of amino acids found in structured proteins.Communicated by Ramaswamy H. Sarma.