Computational Design of Miniproteins as SARS-CoV-2 Therapeutic Inhibitors

A rational therapeutic strategy is urgently needed for combating SARS-CoV-2 infection. Viral infection initiates when the SARS-CoV-2 receptor-binding domain (RBD) binds to the ACE2 receptor, and thus, inhibiting RBD is a promising therapeutic for blocking viral entry. In this study, the structure of lead antiviral candidate binder (LCB1), which has three alpha-helices (H1, H2, and H3), is used as a template to design and simulate several miniprotein RBD inhibitors. LCB1 undergoes two modifications: structural modification by truncation of the H3 to reduce its size, followed by single and double amino acid substitutions to enhance its binding with RBD. We use molecular dynamics (MD) simulations supported by ab initio density functional theory (DFT) calculations. Complete binding profiles of all miniproteins with RBD have been determined. The MD investigations reveal that the H3 truncation results in a small inhibitor with a −1.5 kcal/mol tighter binding to RBD than original LCB1, while the best miniprotein with higher binding affinity involves D17R or E11V + D17R mutation. DFT calculations provide atomic-scale details on the role of hydrogen bonding and partial charge distribution in stabilizing the minibinder:RBD complex. This study provides insights into general principles for designing potential therapeutics for SARS-CoV-2.

Collinearity of alpha-helices or beta strands in membrane proteins causes a characteristic peak centred on 4.9 Å resolution in diffraction intensity profiles, inducing higher diffraction anisotropy

Diffraction anisotropy is a phenomenon that impacts more specifically membrane proteins, compared to soluble ones, but the reasons for this discrepancy remained unclear. Often, it is referred to a difference in resolution limits between highest and lowest diffraction limits as a signature for anisotropy. We show in this article that there is no simple correlation between anisotropy and difference in resolution limits, with notably a substantial number of structures displaying various anisotropy with no difference in resolution limits. We further investigated diffraction intensity profiles, and observed a peak centred on 4.9Å resolution more predominant in membrane proteins. Since this peak is in the region corresponding to secondary structures, we investigated the influence of secondary structure ratio. We showed that secondary structure content has little influence on this profile, while secondary structure collinearity in membrane proteins correlate with a stronger peak. Finally, we could further show that the presence of this peak is linked to higher diffraction anisotropy.SynopsisMembrane protein diffraction anisotropy originates from a peak at 4.9 Å resolution in intensity profiles, due to secondary structure collinearity.

Clostridium perfringens Beta2 toxin forms highly cation-selective channels in lipid bilayers

Clostridium perfringens is a potent producer of a variety of toxins. Well studied from these are five toxins (alpha, Beta (CPB), epsilon, iota and CPE) that are produced by seven toxinotype strains (A–G) of C. perfringens. Besides these toxins, C. perfringens produces also another toxin that causes necrotizing enterocolitis in piglets. This toxin termed consensus Beta2 toxin (cCPB2) has a molecular mass of 27,620 Da and shows only little homology to CPB and no one to the other toxins of C. perfringens. Its primary action on cells remained unknown to date. cCPB2 was heterogeneously expressed as fusion protein with GST in Escherichia coli and purified to homogeneity. Although cCPB2 does not exhibit the typical structure of beta-stranded pore-forming proteins and contains no indication for the presence of amphipathic alpha-helices we could demonstrate that cCPB2 is a pore-forming component with an extremely high activity in lipid bilayers. The channels have a single-channel conductance of about 700 pS in 1 M KCl and are highly cation-selective as judged from selectivity measurements in the presence of salt gradients. The high cation selectivity is caused by the presence of net negative charges in or near the channel that allowed an estimate of the channel size being about 1.4 nm wide. Our measurements suggest that the primary effect of cCPB2 is the formation of cation-selective channels followed by necrotic enteritis in humans and animals. We searched in databases for homologs of cCPB2 and constructed a cladogram representing the phylogenetic relationship to the next relatives of cCPB2.

Interaction of Spike protein and lipid membrane of SARS-CoV-2 with Ursodeoxycholic acid, an in-silico analysis

Numerous repositioned drugs have been sought to decrease the severity of SARS-CoV-2 infection. It is known that among its physicochemical properties, Ursodeoxycholic Acid (UDCA) has a reduction in surface tension and cholesterol solubilization, it has also been used to treat cholesterol gallstones and viral hepatitis. In this study, molecular docking was performed with the SARS-CoV-2 Spike protein and UDCA. In order to confirm this interaction, we used Molecular Dynamics (MD) in “SARS-CoV-2 Spike protein-UDCA”. Using another system, we also simulated MD with six UDCA residues around the Spike protein at random, naming this “SARS-CoV-2 Spike protein-6UDCA”. Finally, we evaluated the possible interaction between UDCA and different types of membranes, considering the possible membrane conformation of SARS-CoV-2, this was named “SARS-CoV-2 membrane-UDCA”. In the “SARS-CoV-2 Spike protein-UDCA”, we found that UDCA exhibits affinity towards the central region of the Spike protein structure of − 386.35 kcal/mol, in a region with 3 alpha helices, which comprises residues from K986 to C1032 of each monomer. MD confirmed that UDCA remains attached and occasionally forms hydrogen bonds with residues R995 and T998. In the presence of UDCA, we observed that the distances between residues atoms OG1 and CG2 of T998 in the monomers A, B, and C in the prefusion state do not change and remain at 5.93 ± 0.62 and 7.78 ± 0.51 Å, respectively, compared to the post-fusion state. Next, in “SARS-CoV-2 Spike protein-6UDCA”, the three UDCA showed affinity towards different regions of the Spike protein, but only one of them remained bound to the region between the region's heptad repeat 1 and heptad repeat 2 (HR1 and HR2) for 375 ps of the trajectory. The RMSD of monomer C was the smallest of the three monomers with a value of 2.89 ± 0.32, likewise, the smallest RMSF was also of the monomer C (2.25 ± 056). In addition, in the simulation of “SARS-CoV-2 membrane-UDCA”, UDCA had a higher affinity toward the virion-like membrane; where three of the four residues remained attached once they were close (5 Å, to the centre of mass) to the membrane by 30 ns. However, only one of them remained attached to the plasma-like membrane and this was in a cluster of cholesterol molecules. We have shown that UDCA interacts in two distinct regions of Spike protein sequences. In addition, UDCA tends to stay bound to the membrane, which could potentially reduce the internalization of SARS-CoV-2 in the host cell.

Bucky Ball Is a Novel Zebrafish Vasa ATPase Activator

Many multicellular organisms specify germ cells during early embryogenesis by the inheritance of ribonucleoprotein (RNP) granules known as germplasm. However, the role of complex interactions of RNP granules during germ cell specification remains elusive. This study characterizes the interaction of RNP granules, Buc, and zebrafish Vasa (zfVasa) during germ cell specification. We identify a novel zfVasa-binding motif (Buc-VBM) in Buc and a Buc-binding motif (zfVasa-BBM) in zfVasa. Moreover, we show that Buc and zfVasa directly bind in vitro and that this interaction is independent of the RNA. Our circular dichroism spectroscopy data reveal that the intrinsically disordered Buc-VBM peptide forms alpha-helices in the presence of the solvent trifluoroethanol. Intriguingly, we further demonstrate that Buc-VBM enhances zfVasa ATPase activity, thereby annotating the first biochemical function of Buc as a zfVasa ATPase activator. Collectively, these results propose a model in which the activity of zfVasa is a central regulator of primordial germ cell (PGC) formation and is tightly controlled by the germplasm organizer Buc.

Computational elucidation of phylogenetic, functional and structural features of Methioninase from Pseudomonas, Escherichia, Clostridium and Citrobacter strains

Background: L-Methioninase (EC 4.4.1.11; MGL) is a pyridoxal phosphate (PLP)-dependent enzyme that is produced by a variety of bacteria, fungi, and plants. L-methioninase, especially from Pseudomonas and Citrobacter sp. is considered as the efficient therapeutic enzymes, particularly in cancers such as glioblastomas, medulloblastoma, and neuroblastoma that are more sensitive to methionine starvation. Objective: The low stability is one of the main drawbacks of the enzyme; in this regard, in the current study, different features of the enzyme, including phylogenetic, functional, and structural from Pseudomonas, Escherichia, Clostridium, and Citrobacter strains were evaluated to find the best bacterial L-Methioninase. Methods: After the initial screening of L-Methioninase sequences from the above-mentioned bacterial strains, the three-dimensional structure of enzymes from Escherichia fergusonii, Pseudomonas fluorescens, and Clostridium homopropionicum were determined through homology modeling via GalaxyTBM server and refined by GalaxyRefine server. Results & Conclusion: Afterwards, PROCHECK, verify 3D, and ERRAT servers were used for verification of the obtained models. Moreover, antigenicity, allergenicity, and physico-chemical analysis of enzymes were also carried out. In order to get insight into interaction of enzyme with other proteins the STRING server was used. The secondary structure of enzyme is mainly composed of random coils and alpha-helices. However, these outcomes should further be validated by wet lab investigations.

Phospholipase A2 (PLA2) Sequences in Rattus norvegicus Genome

Phospholipase A2 is enzyme that hydrolyses phospholipids at sn-2 position. This class of enzymes are significant due to their ability to cleave membrane phospholipids and hence causing inflammation. The PLA2 enzymes present in Rattus norvegicus is extensively studied to predict its properties. The Protparam analysis was performed to predict the physical properties like number of amino acids, Theoretical pH, stability index value, aliphatic index value and GRAVY value. The SOPMA analysis predicted its structural properties like the number of alpha-helices and beta-strands. Hence the focus of the present study was to perform a preliminary in silico analysis to identify the PLA2 protein sequences in the genome of Rattus norvegicus.

Open-Bundle Structure as the Unfolding Intermediate of Cytochrome c’ Revealed by Small Angle Neutron Scattering

The open-bundle structure of cytochrome c’ as an unfolding intermediate was determined by small-angle neutron scattering experiment (SANS). The four-α-helix bundle structure of Cyt c’ at neutral pH was transited to an open-bundle structure (at pD ~13), which is a joint-clubs consisting of four clubs (α-helices) connected by short loops. The compactly folded structure of Cyt c’ (radius of gyration, Rg = 18 Å for the Cyt c’ dimer) at neutral or mildly alkaline pD transitioned to a remarkably larger “open-bundle” structure at pD ~13 (Rg = 25 Å for the Cyt c’ monomer). Cyt c’ adopts an unstructured random coil structure at pD = 1.7 (Rg = 25 Å for the Cyt c’ monomer). Numerical partial scattering function analysis (joint-clubs) and ab initio modelling gave structures similar to the “open-bundle”, which retains the α-helices but loses the bundle structure.

THAP9 Transposase Cleaves DNA via Conserved Acidic Residues in an RNaseH-Like Domain

The catalytic domain of most ‘cut and paste’ DNA transposases have the canonical RNase-H fold, which is also shared by other polynucleotidyl transferases such as the retroviral integrases and the RAG1 subunit of V(D)J recombinase. The RNase-H fold is a mixture of beta sheets and alpha helices with three acidic residues (Asp, Asp, Glu/Asp—DDE/D) that are involved in the metal-mediated cleavage and subsequent integration of DNA. Human THAP9 (hTHAP9), homologous to the well-studied Drosophila P-element transposase (DmTNP), is an active DNA transposase that, although domesticated, still retains the catalytic activity to mobilize transposons. In this study we have modeled the structure of hTHAP9 using the recently available cryo-EM structure of DmTNP as a template to identify an RNase-H like fold along with important acidic residues in its catalytic domain. Site-directed mutagenesis of the predicted catalytic residues followed by screening for DNA excision and integration activity has led to the identification of candidate Ds and Es in the RNaseH fold that may be a part of the catalytic triad in hTHAP9. This study has helped widen our knowledge about the catalytic activity of a functionally uncharacterized transposon-derived gene in the human genome.