Computational study of protein secondary structure elements: Ramachandran plots revisited

In Silico Targeting and Immunological Profiling of PpiA in Mycobacterium tuberculosis: A Computational Approach

Tuberculosis (TB) remains a leading cause of mortality, with drug resistance highlighting the need for new vaccine targets. Peptidyl-prolyl isomerase A (PpiA), a conserved Mycobacterium tuberculosis (Mtb) protein, plays a role in bacterial stress adaptation and immune evasion, making it a potential target for immunotherapy. This study uses computational methods to assess PpiA's antigenicity, structural integrity, and immunogenic potential. The PpiA sequence was retrieved from NCBI and analyzed for antigenicity and allergenicity using VaxiJen, AllerTOP, and AllergenFP. Physicochemical properties were evaluated using ProtParam, and structural models were generated through PSIPRED and SWISS-MODEL. Structural validation was performed with MolProbity, QMEANDisCo, and ProSA-Web. B-cell epitopes were predicted using BepiPred 2.0 and IEDB, while T-cell epitopes were mapped via IEDB's MHC-I and MHC-II tools. Epitope conservation across Mtb strains was confirmed using ConSurf. Results indicate PpiA is highly antigenic, non-allergenic, and stable, with several immunogenic epitopes identified for both B- and T-cells. This study supports PpiA as a promising immunogenic target for TB vaccine development.

Exploring stevioside binding affinity with various proteins and receptors actively involved in the signaling pathway and a future candidate for diabetic patients

Introduction and Background: Diabetes is a chronic metabolic disease characterized by elevated blood glucose levels and is one of the main global health concerns. Synthetic sugar substrate has many side effects such as leukemia, bladder cancer, hepatotoxicity, breast cancer, headache, and brain toxicity. The WHO and FDA has recently banned some of the synthetic sugar alternatives due to their carcinogenic effects. Objective and Methodology: Therefore, the main objective of the current study was to investigate the safety and binding affinity of Stevioside with Glucose Transpoter-4 (GLUT-4), Akt, Insulin Receptor (IR) and Insulin Receptor Substrate-1 (IRS-1) to confirmed that Stevioside is one the potent natural sweetener/drug for diabetes. This study delves into the molecular interaction between Stevioside and key diabetic proteins: GLUT-4, Akt, IR and IRS-1. A precise molecular docking approach was used to simulate the binding affinity of Stevioside to these proteins. The pharmacokinetic properties of the molecule should be taken into consideration as important variables throughout the virtual screening process. Results: The result of active site analysis of GLUT-4, Akt, IR and IRS-1 showed a zone of 2158.359 Ǻ2, 579.259 Ǻ2, 762.651 Ǻ2, and 152.167 Ǻ2 and a volume of 2765.094 Ǻ³, 355.567 Ǻ³, 686.806 Ǻ³, and 116.874 Ǻ³, respectively. Docking analysis of the Stevioside compound showed the highest docking energy with scores of -9.9 with GLUT-4, -6.7 with Akt, -8.0 with IR and -8.8 with IRS-1. Studies indicated that it remains undigested by stomach acids and enzymes and is not absorbed in the upper small intestine. Further, tests revealed no hepatotoxicity, AMES toxicity, or skin sensitivity, making it a promising candidate for safe consumption as drug metabolism. Conclusion and Recommendations: Instead of other sugar alternatives, Stevioside will help diabetic patients with a lower chance of infections, lowered blood pressure/blood sugar, and increased glucose uptake in diabetic muscles. Stevioside is a natural sweetener, and the current study recommends its usage in various dietary products for diabetic patients.

Conformational preference of dipeptide zwitterions in aqueous solvents

Proper description of solvent effects is challenging for theoretical methods, particularly if the solute is a zwitterion. Here, a series of theoretical procedures are used to determine the preferred solvated conformations of twelve hydrophobic dipeptides (Leu-Leu, Leu-Phe, Phe-Leu, Ile-Leu, Phe-Phe, Ala-Val, Val-Ala, Ala-Ile, Ile-Ala, Ile-Val, Val-Ile and Val-Val) in the zwitterionic state. First, the accuracy of density functional theory (DFT), combined with different implicit solvent models, for describing zwitterions in aqueous solvent is assessed by comparing the predicted against the experimental glycine tautomerization energy, i.e., the energetic difference between canonical and zwitterionic glycine in aqueous solvents. It is found that among the tested solvation schemes, the charge-asymmetric nonlocally determined local-electric solvation model (CANDLE) predicts an energetic difference in excellent agreement with the experimental value. Next, DFT-CANDLE is used to determine the most favorable solvated conformation for each of the investigated dipeptide zwitterions. The CANDLE-solvated structures are obtained by exploring the conformational space of each dipeptide zwitterion concatenating DFT calculations, in vacuum, with classical molecular dynamics simulations, in explicit solvents, and DFT calculations including explicit water molecules. It is found that the energetically most favorable conformations are similar to those of the dipeptide zwitterions in their respective crystal structures. Such structural agreement is indicative of the DFT-CANDLE accomplishment of the description of solvated zwitterions, and suggests that these biomolecules self-assemble as quasi-rigid objects.

Identification and In-Silico study of non-synonymous functional SNPs in the human SCN9A gene

Single nucleotide polymorphisms are the most common form of DNA alterations at the level of a single nucleotide in the genomic sequence. Genome-wide association studies (GWAS) were carried to identify potential risk genes or genomic regions by screening for SNPs associated with disease. Recent studies have shown that SCN9A comprises the NaV1.7 subunit, Na+ channels have a gene encoding of 1988 amino acids arranged into 4 domains, all with 6 transmembrane regions, and are mainly found in dorsal root ganglion (DRG) neurons and sympathetic ganglion neurons. Multiple forms of acute hypersensitivity conditions, such as primary erythermalgia, congenital analgesia, and paroxysmal pain syndrome have been linked to polymorphisms in the SCN9A gene. Under this study, we utilized a variety of computational tools to explore out nsSNPs that are potentially damaging to heath by modifying the structure or activity of the SCN9A protein. Over 14 potentially damaging and disease-causing nsSNPs (E1889D, L1802P, F1782V, D1778N, C1370Y, V1311M, Y1248H, F1237L, M936V, I929T, V877E, D743Y, C710W, D623H) were identified by a variety of algorithms, including SNPnexus, SNAP-2, PANTHER, PhD-SNP, SNP & GO, I-Mutant, and ConSurf. Homology modeling, structure validation, and protein-ligand interactions also were performed to confirm 5 notable substitutions (L1802P, F1782V, D1778N, V1311M, and M936V). Such nsSNPs may become the center of further studies into a variety of disorders brought by SCN9A dysfunction. Using in-silico strategies for assessing SCN9A genetic variations will aid in organizing large-scale investigations and developing targeted therapeutics for disorders linked to these variations.

Molecular modeling of ARF6 dysregulation caused by mutations in IQSEC2

IQSEC2 gene mutations are associated with epilepsy, autism, and intellectual disability. The primary function IQSEC2, mediated via its Sec 7 domain, is to act as a guanine nucleotide exchange factor for ARF6. We sought to develop a molecular model, which may explain the aberrant Sec 7 activity on ARF6 of different human IQSEC2 mutations. We integrated experimental data of IQSEC2 mutants with protein structure prediction by the RaptorX server combined with molecular modeling and molecular dynamics simulations. Normally, apocalmodulin (apoCM) binds to IQSEC2 resulting in its N-terminal fragment inhibiting access of its Sec 7 domain to ARF6. An increase in Ca2+ concentration destabilizes the interaction of IQSEC2 with apoCM and removes steric hindrance of Sec 7 binding with ARF6. Mutations at amino acid residue 350 of IQSEC2 result in loss of steric hindrance of Sec 7 binding with ARF6 leading to constitutive activation of ARF6 by Sec 7. On the other hand, a mutation at amino acid residue 359 of IQSEC2 results in constitutive hindrance of Sec 7 binding to ARF6 leading to the loss of the ability of IQSEC2 to activate ARF6. These studies provide a model for dysregulation of IQSEC2 Sec 7 activity by mutant IQSEC2 proteins.Communicated by Ramaswamy H. Sarma.

Computational study of the motor neuron protein KIF5A to identify nsSNPs, bioactive compounds, and its key regulators

Introduction: Kinesin family member 5A (KIF5A) is a motor neuron protein expressed in neurons and involved in anterograde transportation of organelles, proteins, and RNA. Variations in the KIF5A gene that interfere with axonal transport have emerged as a distinguishing feature in several neurodegenerative disorders, including hereditary spastic paraplegia (HSP10), Charcot-Marie-Tooth disease type 2 (CMT2), and Amyotrophic Lateral Sclerosis (ALS). Methods: In this study, we implemented a computational structural and systems biology approach to uncover the role of KIF5A in ALS. Using the computational structural biology method, we explored the role of non-synonymous Single Nucleotide Polymorphism (nsSNPs) in KIF5A. Further, to identify the potential inhibitory molecule against the highly destabilizing structure variant, we docked 24 plant-derived phytochemicals involved in ALS. Results: We found KIF5AS291F variant showed the most structure destabilizing behavior and the phytocompound "epigallocatechin gallate" showed the highest binding affinity (-9.0 Kcal/mol) as compared to wild KIF5A (-8.4 Kcal/mol). Further, with the systems biology approach, we constructed the KIF5A protein-protein interaction (PPI) network to identify the associated Kinesin Families (KIFs) proteins, modules, and their function. We also constructed a transcriptional and post-transcriptional regulatory network of KIF5A. With the network topological parameters of PPIN (Degree, Bottleneck, Closeness, and MNC) using CytoHubba and computational knock-out experiment using Network Analyzer, we found KIF1A, 5B, and 5C were the significant proteins. The functional modules were highly enriched with microtubule motor activity, chemical synaptic transmission in neurons, GTP binding, and GABA receptor activity. In regulatory network analysis, we found KIF5A post-transcriptionally down-regulated by miR-107 which is further transcriptionally up-regulated by four TFs (HIF1A, PPARA, SREBF1, and TP53) and down-regulated by three TFs (ZEB1, ZEB2, and LIN28A). Discussion: We concluded our study by finding a crucial variant of KIF5A and its potential therapeutic target (epigallocatechin gallate) and KIF5A associated significant genes with important regulators which could decrypt the novel therapeutics in ALS and other neurodegenerative diseases.

Oxa-376 and Oxa-530 variants of β-lactamase: computational study uncovers potential therapeutic targets of Acinetobacter baumannii

Antimicrobial resistance is a major global health crisis, resulting in thousands of deaths each year. Antibiotics' effectiveness against microorganisms deteriorates over time as multidrug resistance (MDR) develops, which is exacerbated by irregular antibiotic use, poor disease management, and the evasive nature of bacteria. The World Health Organization has recognized multidrug resistance as a critical public health concern, and Acinetobacter baumannii has been at the center of attention due to its ability to develop multidrug resistance (MDR). It generally produces carbapenem-hydrolyzing oxacillinase, which has been identified as the primary source of beta-lactam resistance in MDR bacteria. Recently, point mutations in A. baumannii have been identified as a key factor of multidrug resistance, making them a prime concern for researchers. The goal of the current work was to establish a unique way of finding multidrug-resistant variants and identify the most damaging mutations in the existing databases. We characterized the deleterious variants of oxacillinases using several computational tools. Following a thorough analysis, Oxa-376 and Oxa-530 were found to be more damaging when compared with the wild-type Oxa-51. The mutants' 3D structures were then prepared and refined with RaptorX, GalaxyRefine, and SAVES servers. Our research incorporates seven antimicrobial agents to illustrate the resistance capability of the variants of oxacillinase by evaluating binding affinity in Autodock-vina and Schrodinger software. RMSD, RMSF, Radius of gyration analysis, the solvent-accessible surface area (SASA), hydrogen bonding analysis and MM-GBSA from Molecular Dynamics Simulation revealed the dynamic nature and stability of wild-type and Oxa-376 and Oxa-530 variants. Our findings will benefit researchers looking for the deleterious mutations of Acinetobacter baumannii and new therapeutics to combat those variants. However, further studies are necessary to evaluate the mechanism of hydrolyzing activity and antibiotic resistance of these variants.

Determining novel therapeutic targets of Acinetobacter baumannii. Deleterious variants, causing antibiotic resistance, were identified by molecular docking and molecular dynamics simulation suggesting new therapeutic targets Oxa-376 and Oxa-530.

PSIQUE: Protein Secondary Structure Identification on the Basis of Quaternions and Electronic Structure Calculations

The secondary structure is important in protein structure analysis, classification, and modeling. We have developed a novel method for secondary structure assignment, termed PSIQUE, based on the potential energy surface (PES) of polyalanine obtained using an infinitely long chain model and density functional theory calculations. First, uniform protein segments are determined in terms of a difference of quaternions between neighboring amino acids along the protein backbone. Then, the identification of the secondary structure motifs is carried out based on the minima found in the PES. PSIQUE shows good agreement with other secondary structure assignment methods. However, it provides better discrimination of subtle secondary structures (e.g., helix types) and termini and produces more uniform segments while also accounting for local distortions. Overall, PSIQUE provides a precise and reliable assignment of secondary structures, so it should be helpful for the detailed characterization of the protein structure.

In silico drug design and molecular docking studies targeting Akt1 (RAC-alpha serine/threonine-protein kinase) and Akt2 (RAC-beta serine/threonine-protein kinase) proteins and investigation of CYP (cytochrome P450) inhibitors against MAOB (monoamine oxidase B) for OSCC (oral squamous cell carcinoma) treatment

The overexpression of Akt1 (RAC-alpha serine/threonine-protein Kinase) and Akt2 (RAC-beta serine/threonine-protein Kinase) is a hallmark of Oral Squamous Cell Carcinoma (OSCC). Because of the elevated frequency of OSCC occurrence in South Asian countries, novel therapeutic approaches are indispensable. Drugs that inhibit the overexpression of Akt1 and Akt2 proteins in Akt pathway and do not cause reduced expression of MAOB can be leads for OSCC treatment. In this study, Akt1, Akt2 and MAOB were targeted and 100 CYP inhibitors were screened through several in silico approaches and Galuteolin and Linarin were identified as potential leads for OSCC treatment as they inhibited Akt1 proteins with strong binding affinities of -12.3 and -11.5 kcal/mol respectively and also Akt2 proteins with strong binding affinities of -11.4 and -11.1 kcal/mol respectively, but they did not inhibit MAOB. Decreased expression of MAOB in tissues causes OSCC but overexpression is also responsible for other types of diseases and cancers. From the investigation of CYP inhibitors against MAOB, five CYP inhibitors- Diosmetin, Acacetin, Epicatechin, Eriodictyol and Capillin have expressed inhibitory action against MAOB without any interference with Akt1 and Akt2. This study mainly represents that Galuteolin and Linarin in the Akt pathway can be perceived for OSCC treatment and other five CYP inhibitors - Diosmetin, Acacetin, Epicatechin, Eriodictyol and Capillin for the treatment of other diseases and cancers caused by overexpression of MAOB. ADMET properties of CYP inhibitors obtained from admetSAR 2.0 and were compared with reference drugs for validation. Communicated by Ramaswamy H. Sarma.

Structure-based drug design, synthesis and screening of MmaA1 inhibitors as novel anti-TB agents

Tuberculosis is one of the leading causes of death across the world. The treatment regimens for tuberculosis are well established, but still the control of the disease faces many challenges such as lengthy treatment protocols, drug resistance and toxicity. In the present work, mycolic acid methyl transferase (MmaA1), a protein involved in the maturation of mycolic acids in the biochemical pathway of the Mycobacterium, was studied for novel drug discovery. The homology model of the MmaA1 protein was built and validated by using computational techniques. The MmaA1 protein has 286 amino acid residues consisting of 10 α-helices and 7 β-sheets. The active site of the MmaA1 protein was identified using CASTp, SiteMap and PatchDock. Virtual screening studies were performed with two small molecule ligand databases: Asinex synergy and Diverse_Elite_Gold_Platinum databases having a total of 43,446 molecules and generated 1,30,814 conformers against the predicted and validated active site of the MmaA1 protein. Binding analysis showed that the residues ASP 19, PHE 22, TRP 30, TYR 32, TRP 74 and ALA 77 of MmaA1 protein have consistent interactions with the ligands. The hit ligands were further filtered by in silico ADME properties to eliminate potentially toxic molecules. Of the top 10 molecules, 3-(2-morpholinoacetamido)-N-(1,4-dihydro-4-oxoquinazolin-6-yl) benzamide was synthesised and screened for in vitro anti-TB activity against Mtb H37Rv using MABA assay. The compound and its intermediates exhibited good in vitro anti-TB activity which can be taken up for future lead optimisation studies. Structure based virtual screening study was performed using a validated homology model against small molecules from two virtual compound libraries. Synthesised the lead compound 3-(2-morpholinoacetamido)-N-(1,4-dihydro-4-oxoquinazolin-6-yl)benzamide obtained from virtual screening. In vitro activity against Mtb H37Rv has given a promising result.

Identification and molecular characterization of a psychrophilic GH1 β-glucosidase from the subtropical soil microorganism Exiguobacterium sp. GXG2

The products of bacterial β-glucosidases with favorable cold-adapted properties have industrial applications. A psychrophilic β-glucosidase gene named bglG from subtropical soil microorganism Exiguobacterium sp. GXG2 was isolated and characterized by function-based screening strategy. Results of multiple alignments showed that the derived protein BglG shared 45.7% identities with reviewed β-glucosidases in the UniProtKB/Swiss-Prot database. Functional characterization of the β-glucosidase BglG indicated that BglG was a 468 aa protein with a molecular weight of 53.2 kDa. The BglG showed the highest activity in pH 7.0 at 35 °C and exhibited consistently high levels of activity within low temperatures ranging from 5 to 35 °C. The BglG appeared to be a psychrophilic enzyme. The values of Km, Vmax, kcat, and kcat/Km of recombinant BglG toward ρNPG were 1.1 mM, 1.4 µg/mL/min, 12.7 s-1, and 11.5 mM/s, respectively. The specific enzyme activity of BglG was 12.14 U/mg. The metal ion of Ca2+ and Fe3+ could stimulate the activity of BglG, whereas Mn2+ inhibited the activity. The cold-adapted β-glucosidase BglG displayed remarkable biochemical properties, making it a potential candidate for future industrial applications.

Newly Synthesized Oxygenated Xanthones as Potential P-Glycoprotein Activators: In Vitro, Ex Vivo, and In Silico Studies

P-glycoprotein (P-gp) plays a crucial role in the protection of susceptible organs, by significantly decreasing the absorption/distribution of harmful xenobiotics and, consequently, their toxicity. Therefore, P-gp has been proposed as a potential antidotal pathway, when activated and/or induced. Knowing that xanthones are known to interact with P-gp, the main goal was to study P-gp induction or/and activation by six new oxygenated xanthones (OX 1-6). Furthermore, the potential protection of Caco-2 cells against paraquat cytotoxicity was also assessed. The most promising compound was further tested for its ability to increase P-gp activity ex vivo, using everted intestinal sacs from adult Wistar-Han rats. The oxygenated xanthones interacted with P-gp in vitro, increasing P-gp expression and/or activity 24 h after exposure. Additionally, after a short-incubation period, several xanthones were identified as P-gp activators, as they immediately increased P-gp activity. Moreover, some xanthones decreased PQ cytotoxicity towards Caco-2 cells, an effect prevented under P-gp inhibition. Ex vivo, a significant increase in P-gp activity was observed in the presence of OX6, which was selectively blocked by a model P-gp inhibitor, zosuquidar, confirming the in vitro results. Docking simulations between a validated P-gp model and the tested xanthones predicted these interactions, and these compounds also fitted onto previously described P-gp induction and activation pharmacophores. In conclusion, the in vitro, ex vivo, and in silico results suggest the potential of some of the oxygenated xanthones in the modulation of P-gp, disclosing new perspectives in the therapeutics of intoxications by P-gp substrates.

Everyone Is a Protagonist: Residue Conformational Preferences in High-Resolution Protein Structures

In many structural bioinformatics problems, there is a broad range of unanswered questions about protein dynamics and amino acid properties. Proteins are not strictly static objects, but rather populate ensembles of conformations. One way to understand these particularities is to analyze the information available in experimental databases. The Ramachandran plot, despite being more than half a century old, remains an utterly useful tool in the study of protein conformation. Based on its assumptions, we inspected a large data set (11,130 protein structures, amounting to 5,255,768 residues) and discriminated the conformational preferences of each residue type regarding their secondary structure participation. These data were studied for phi [Formula: see text], psi [Formula: see text], and side chain chi [Formula: see text] angles, being presented in non-Ramachandranian plots. In the largest analysis of protein conformation made so far, we propose an original plot to depict conformational preferences in relation to different secondary structure elements. Despite confirming previous observations, our results strongly support a unique character for each residue type, whereas also reinforcing the observation that side chains have a major contribution to secondary structure and, by consequence, on protein conformation. This information can be further used in the development of more robust methods and computational strategies for structural bioinformatics problems.

Computational Modeling of complete HOXB13 protein for predicting the functional effect of SNPs and the associated role in hereditary prostate cancer

The human HOXB13 gene encodes a 284 amino acid transcription factor belonging to the homeobox gene family containing a homeobox and a HoxA13 N-terminal domain. It is highly linked to hereditary prostate cancer, the majority of which is manifested as a result of a Single Nucleotide Polymorphism (SNP). In silico analysis of 95 missense SNP's corresponding to the non-homeobox region of HOXB13 predicted 21 nsSNP's to be potentially deleterious. Among 123 UTR SNPs analysed by UTRScan, rs543028086, rs550968159, rs563065128 were found to affect the UNR_BS, GY-BOX and MBE UTR signals, respectively. Subsequent analysis by PolymiRTS revealed 23 UTR SNPs altering the miRNA binding site. The complete HOXB13_M26 protein structure was modelled using MODELLER v9.17. Computational analysis of the 21 nsSNP's mapped into the HOXB13_M26 protein revealed seven nsSNP's (rs761914407, rs8556, rs138213197, rs772962401, rs778843798, rs770620686 and rs587780165) seriously resulting in a damaging and deleterious effect on the protein. G84E, G135E, and A128V resulted in increased, while, R215C, C66R, Y80C and S122R resulted in decreased protein stability, ultimately predicted to result in the altered binding patterns of HOXB13. While the genotype-phenotype based effects of nsSNP's were assessed, the exact biological and biochemical mechanism driven by the above predicted SNPs still needs to be extensively evaluated by in vivo and GWAS studies.

Experimental conformational energy maps of proteins and peptides

We have presented an extensive analysis of the peptide backbone dihedral angles in the PDB structures and computed experimental Ramachandran plots for their distributions seen under a various constraints on X-ray resolution, representativeness at different sequence identity percentages, and hydrogen bonding distances. These experimental distributions have been converted into isoenergy contour plots using the approach employed previously by F. M. Pohl. This has led to the identification of energetically favored minima in the Ramachandran (ϕ, ψ) plots in which global minima are predominantly observed either in the right-handed α-helical or the polyproline II regions. Further, we have identified low energy pathways for transitions between various minima in the (ϕ,ψ) plots. We have compared and presented the experimental plots with published theoretical plots obtained from both molecular mechanics and quantum mechanical approaches. In addition, we have developed and employed a root mean square deviation (RMSD) metric for isoenergy contours in various ranges, as a measure (in kcal.mol^-1 ) to compare any two plots and determine the extent of correlation and similarity between their isoenergy contours. In general, we observe a greater degree of compatibility with experimental plots for energy maps obtained from molecular mechanics methods compared to most quantum mechanical methods. The experimental energy plots we have investigated could be helpful in refining protein structures obtained from X-ray, NMR, and electron microscopy and in refining force field parameters to enable simulations of peptide and protein structures that have higher degree of consistency with experiments. Proteins 2017; 85:979-1001. © 2017 Wiley Periodicals, Inc.

In Silico Template Selection of Short Antimicrobial Peptide Viscotoxin for Improving Its Antimicrobial Efficiency in Development of Potential Therapeutic Drugs

Rapid increase in antibiotic resistance has posed a worldwide threat, due to increased mortality, morbidity, and expenditure caused by antibiotic-resistant microbes. Recent development of the antimicrobial peptides like viscotoxin (Vt) has been successfully comprehended as a substitute for classical antibiotics. A structurally stable peptide, Vt can enhance antimicrobial property and can be used for various developmental purposes. Thus, structural stability among the antimicrobial peptides, Vt A1 (3C8P), A2 (1JMN), A3 (1ED0), B (1JMP), and C (1ORL) of Viscus album was computationally analyzed. In specific, the static confirmation of VtA3 showed high number of intramolecular interactions, along with an increase in hydrophobicity than others comparatively. Further, conformational sampling was used to analyze various geometrical parameters such as root mean square deviation, root mean square fluctuation, radius of gyration, and ovality which also revealed the structural stability of VtA3. Moreover, the statistically validated contours of surface area, lipophilicity, and distance constraints of disulfide bonds also supported the priority of VtA3 with respect to stability. Finally, the functional activity of peptides was accessed by computing their free energy of membrane association and membrane interactions, which defined VtA3 as functionally stable. Currently, peptide-based antibiotics and nanoparticles have attracted the pharmaceutical industries for their potential therapeutic applications. Thereby, it is proposed that viscotoxin A3 (1ED0) could be used as a preeminent template for scaffolding potentially efficient antimicrobial peptide-based drugs and nanomaterials in future.

Beyond basins: φ,ψ preferences of a residue depend heavily on the φ,ψ values of its neighbors

The Ramachandran plot distributions of nonglycine residues from experimentally determined structures are routinely described as grouping into one of six major basins: β, PII , α, αL , ξ and γ'. Recent work describing the most common conformations adopted by pairs of residues in folded proteins [i.e., (φ,ψ)2 -motifs] showed that commonly described major basins are not true single thermodynamic basins, but are composed of distinct subregions that are associated with various conformations of either the preceding or following neighbor residue. Here, as documentation of the extent to which the conformational preferences of a central residue are influenced by the conformations of its two neighbors, we present a set of φ,ψ-plots that are delimited simultaneously by the φ,ψ-angles of its neighboring residues on both sides. The level of influence seen here is typically greater than the influence associated with considering the identities of neighboring residues, implying that the use of this heretofore untapped information can improve the accuracy of structure prediction algorithms and low resolution protein structure refinement.

The dipeptide conformations of all twenty amino acid types in the context of biosynthesis

There have been many studies of dipeptide structure at a high level of accuracy using quantum chemical methods. Such calculations are resource-consuming (in terms of memory, CPU and other computational imperatives) which is the reason why most previous studies were restricted to the two simplest amino-acid residue types, glycine and alanine. We improve on this by extending the scope of residue types to include all 20 naturally occurring residue types. Our results reveal differences in secondary structure preferences for the all residue types. There are in most cases very deep energy troughs corresponding either to the polyproline II (collagen) helix and the α-helix or both. The β-strand was not strongly favoured energetically although the extent of this depression in the energy surface is, while not “deeper” (energetically), has a wider extent than the other two types of secondary structure. There is currently great interest in the question of cotranslational folding, the extent to which the nascent polypeptide begins to fold prior to emerging from the ribosome exit tunnel. Accordingly, while most previous quantum studies of dipeptides were carried out in the (simulated) gas or aqueous phase, we wished to consider the first step in polypeptide biosynthesis on the ribosome where neither gas nor aqueous conditions apply. We used a dielectric constant that would be compatible with the water-poor macromolecular (ribosome) environment.