microRNA 1307 Is a Potential Target for SARS-CoV-2 Infection: An in Vitro Model.

microRNAs (miRs) are proposed as critical molecular targets in SARS-CoV-2 infection. Our recent in silico studies identified seven SARS-CoV-2 specific miR-like sequences, which are highly conserved with humans, including miR-1307-3p, with critical roles in COVID-19. In this current study, Vero cells were infected with SARS-CoV-2, and miR expression profiles were thereafter confirmed by qRT-PCR. miR-1307-3p was the most highly expressed miR in the infected cells; we, therefore, transiently inhibited its expression in both infected and uninfected cells. The 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl tetrazolium bromide (MTT) cell proliferation assay assessed cell viability following SARS-CoV-2 infection, identifying that miR-1307 expression is inversely correlated with cell viability. Lastly, changes in miR-1307-dependent pathways were analyzed through a detailed miRNOME and associated in silico analysis. In addition to our previously identified miRs, including miR-1307-3p, the upregulation of miR-193a-5p, miR-5100, and miR-23a-5p and downregulation of miR-130b-5p, miR34a-5p, miR-505-3p, miR181a-2-3p, miR-1271-5p, miR-598-3p, miR-34c-3p, and miR-129-5p were also established in Vero cells related to general lung disease-related genes following SARS-CoV-2 infection. Targeted anti-miR-1307-3p treatment rescued cell viability in infection when compared to SARS CoV-2 mediated cell cytotoxicity only. We furthermore identified by in silico analysis that miR-1307-3p is conserved in all SARS-CoV-2 sequences/strains, except in the BA.2 variant, possibly contributing to the lower disease severity of this variant, which warrants further investigation. Small RNA seq analysis was next used to evaluate alterations in the miRNOME, following miR-1307-3p manipulation, identifying critical pathobiological pathways linked to SARS-CoV-2 infection-mediated upregulation of this miR. On the basis of our findings, miRNAs like miR-1307-3p play a critical role in SARS-CoV-2 infection, including via effects on disease progression and severity.

The Role of CDK4 in the Pathogenesis of Pancreatic Cancer.

Pancreatic cancer (PC) continues to have the lowest overall survival and the lack of effective early diagnosis. Cyclin-dependent kinase 4 (CDK4) plays a fundamental role in the orderly progression of the cell cycle, binding to cyclin D to promote the progression through the G1/2 transition. The inhibition of CDK4/6 has therefore gained substantial interest in the hope of new and effective therapeutics in multiple cancers, such as advanced metastatic breast cancer. While the use of these agents is encouraging, their potential is yet to be fully explored. In this study we used the GLOBOCAN database to understand the most recent epidemiology of PC, Human Protein Atlas and KEGG to highlight the role, prevalence, and significance on patient survival of CDK4 in PC. We found that CDK4 cannot be used as prognostic in PC and no significant differences were observed between CDK4 expression and the patient's clinical status, though larger studies, especially concerning CDK4 protein expressions, are required for a more thorough understanding. The use of CDK4/6 inhibitors in PC is still in clinical trials. However, due to only modest improvements observed in the use of single-agent therapies, efforts have focused on combinatorial approaches.

MicroRNAs for Virus Pathogenicity and Host Responses, Identified in SARS-CoV-2 Genomes, May Play Roles in Viral-Host Co-Evolution in Putative Zoonotic Host Species.

Our recent study identified seven key microRNAs (miR-8066, 5197, 3611, 3934-3p, 1307-3p, 3691-3p, 1468-5p) similar between SARS-CoV-2 and the human genome, pointing at miR-related mechanisms in viral entry and the regulatory effects on host immunity. To identify the putative roles of these miRs in zoonosis, we assessed their conservation, compared with humans, in some key wild and domestic animal carriers of zoonotic viruses, including bat, pangolin, pig, cow, rat, and chicken. Out of the seven miRs under study, miR-3611 was the most strongly conserved across all species; miR-5197 was the most conserved in pangolin, pig, cow, bat, and rat; miR-1307 was most strongly conserved in pangolin, pig, cow, bat, and human; miR-3691-3p in pangolin, cow, and human; miR-3934-3p in pig and cow, followed by pangolin and bat; miR-1468 was most conserved in pangolin, pig, and bat; while miR-8066 was most conserved in pangolin and pig. In humans, miR-3611 and miR-1307 were most conserved, while miR-8066, miR-5197, miR-3334-3p and miR-1468 were least conserved, compared with pangolin, pig, cow, and bat. Furthermore, we identified that changes in the miR-5197 nucleotides between pangolin and human can generate three new miRs, with differing tissue distribution in the brain, lung, intestines, lymph nodes, and muscle, and with different downstream regulatory effects on KEGG pathways. This may be of considerable importance as miR-5197 is localized in the spike protein transcript area of the SARS-CoV-2 genome. Our findings may indicate roles for these miRs in viral-host co-evolution in zoonotic hosts, particularly highlighting pangolin, bat, cow, and pig as putative zoonotic carriers, while highlighting the miRs' roles in KEGG pathways linked to viral pathogenicity and host responses in humans. This in silico study paves the way for investigations into the roles of miRs in zoonotic disease.

Inhibition on JNK Mimics Silencing of Wnt-11 Mediated Cellular Response in Androgen-Independent Prostate Cancer Cells.

Prostate cancer (PCa) is one of the most common cancers among men, and one of the leading causes of cancer death for men. The c-Jun N-terminal kinase (JNK) pathway is required for several cellular functions, such as survival, proliferation, differentiation, and migration. Wnt-11, a member of the Wnt family, has been identified for its upregulation in PCa; however, downstream signalling of Wnt-11 remains to be fully characterized. In this study, we investigated the role of the JNK pathway as a potential downstream factor for Wnt-11 signalling. For this purpose, LNCaP, DU145, and PC-3 PCa cells and normal epithelial PNT1A cells were treated with a specific JNK kinase inhibitor: JNKVIII. Our results showed that JNK inhibition decreased mitochondrial membrane potential and promoted cell death in a cell type-dependent manner. We found that JNK inhibition led to an increase in autophagy and prevented epithelial-mesenchymal transition (EMT) in independently growing androgen cells. JNK inhibition and the silencing of Wnt-11 showed similar responses in DU145 and PC-3 cells and decreased metastasis-related biomarkers, cell migration, and invasion. Overall, our results suggest that JNK signalling plays a significant role in the pathophysiology of PCa by mediating Wnt-11 induced signals. Our data highlights that both the JNK pathway and Wnt-11 could be a useful therapeutic target for the combinatory application of current PCa.

The Prediction of miRNAs in SARS-CoV-2 Genomes: hsa-miR Databases Identify 7 Key miRs Linked to Host Responses and Virus Pathogenicity-Related KEGG Pathways Significant for Comorbidities.

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) is a member of the betacoronavirus family, which causes COVID-19 disease. SARS-CoV-2 pathogenicity in humans leads to increased mortality rates due to alterations of significant pathways, including some resulting in exacerbated inflammatory responses linked to the "cytokine storm" and extensive lung pathology, as well as being linked to a number of comorbidities. Our current study compared five SARS-CoV-2 sequences from different geographical regions to those from SARS, MERS and two cold viruses, OC43 and 229E, to identify the presence of miR-like sequences. We identified seven key miRs, which highlight considerable differences between the SARS-CoV-2 sequences, compared with the other viruses. The level of conservation between the five SARS-CoV-2 sequences was identical but poor compared with the other sequences, with SARS showing the highest degree of conservation. This decrease in similarity could result in reduced levels of transcriptional control, as well as a change in the physiological effect of the virus and associated host-pathogen responses. MERS and the milder symptom viruses showed greater differences and even significant sequence gaps. This divergence away from the SARS-CoV-2 sequences broadly mirrors the phylogenetic relationships obtained from the whole-genome alignments. Therefore, patterns of mutation, occurring during sequence divergence from the longer established human viruses to the more recent ones, may have led to the emergence of sequence motifs that can be related directly to the pathogenicity of SARS-CoV-2. Importantly, we identified 7 key-microRNAs (miRs 8066, 5197, 3611, 3934-3p, 1307-3p, 3691-3p, 1468-5p) with significant links to KEGG pathways linked to viral pathogenicity and host responses. According to Bioproject data (PRJNA615032), SARS-CoV-2 mediated transcriptomic alterations were similar to the target pathways of the selected 7 miRs identified in our study. This mechanism could have considerable significance in determining the symptom spectrum of future potential pandemics. KEGG pathway analysis revealed a number of critical pathways linked to the seven identified miRs that may provide insight into the interplay between the virus and comorbidities. Based on our reported findings, miRNAs may constitute potential and effective therapeutic approaches in COVID-19 and its pathological consequences.

Modelling, inference and big data in biophysics.

In recognition of the increasing importance of big data in biophysics, a new session called 'Modelling, inference, big data' is incorporated into the IUPAB/EBSA Congress on 18 July 2017 at Edinburgh, UK.

The interaction of Wnt-11 and signalling cascades in prostate cancer.

Prostate cancer (PCa) is the second most common cancer among the male population. Conventional therapies target androgen signalling, which drives tumour growth; however, they provide limited survival benefits for patients. It is essential, therefore, to develop a more specific biomarker than the current gold standard, PSA testing. The Wnt signalling pathway induces expression of target genes through cell surface receptors. A non-canonical member of this family, Wnt-11, is evolutionarily highly conserved and is normally expressed by various cells in the developing embryo, as well as in the heart, liver and skeletal muscle of adult humans. We comprehensively review several cell signalling pathways to explain how they interact with Wnt-11, demonstrating its use as a potential biomarker for PCa. Several studies have shown that the expression of Wnt-11 is associated with gastric, renal and colorectal adenocarcinomas and PCa. Moreover, Wnt-11 affects extracellular matrix composition and cytoskeletal rearrangement, and it is required for proliferation and/or survival during cell differentiation. It was found that PCa cell lines express high levels of Wnt-11, which allows differentiation of the epithelial prostate tumour cells to neuron-like (NE) cells. The NE cells produce additional factors that can cause regression after treatment. Accumulating evidence shows that Wnt-11 could be a potential biomarker in diagnosing PCa. Many studies have shown both non-canonical and canonical Wnts interact with several signalling cascades such as PKC, JNK, NF-κB, Rho, PKA and PI3K. In particular, evidence demonstrates Wnt-11 is involved in the progression of PCa, thus it could have the potential to become both a specific disease marker and an important therapeutic target.

Binding modes of diketo-acid inhibitors of HIV-1 integrase: a comparative molecular dynamics simulation study.

HIV-1 integrase (IN) has become an attractive target since drug resistance against HIV-1 reverse transcriptase (RT) and protease (PR) has appeared. Diketo acid (DKA) inhibitors are potent and selective inhibitors of HIV-1 IN: however the action mechanism is not well understood. Here, to study the inhibition mechanism of DKAs we performed 10 ns comparative molecular dynamics simulations on HIV-1 IN bound with three most representative DKA inhibitors: Shionogi inhibitor, S-1360 and two Merck inhibitors L-731,988 and L-708,906. Our simulations show that the acidic part of S-1360 formed salt bridge and cation-π interactions with Lys159. In addition, the catalytic Glu152 in S-1360 was pushed away from the active site to form an ion-pair interaction with Arg199. The Merck inhibitors can maintain either one or both of these ion-pair interaction features. The difference in potencies of the DKA inhibitors is thus attributed to the different binding modes at the catalytic site. Such structural information at atomic level, not only demonstrates the action modes of DKA inhibitors but also provides a novel starting point for structural-based design of HIV-1 IN inhibitors.

The distinct conformational dynamics of K-Ras and H-Ras A59G.

Ras proteins regulate signaling cascades crucial for cell proliferation and differentiation by switching between GTP- and GDP-bound conformations. Distinct Ras isoforms have unique physiological functions with individual isoforms associated with different cancers and developmental diseases. Given the small structural differences among isoforms and mutants, it is currently unclear how these functional differences and aberrant properties arise. Here we investigate whether the subtle differences among isoforms and mutants are associated with detectable dynamical differences. Extensive molecular dynamics simulations reveal that wild-type K-Ras and mutant H-Ras A59G are intrinsically more dynamic than wild-type H-Ras. The crucial switch 1 and switch 2 regions along with loop 3, helix 3, and loop 7 contribute to this enhanced flexibility. Removing the gamma-phosphate of the bound GTP from the structure of A59G led to a spontaneous GTP-to-GDP conformational transition in a 20-ns unbiased simulation. The switch 1 and 2 regions exhibit enhanced flexibility and correlated motion when compared to non-transitioning wild-type H-Ras over a similar timeframe. Correlated motions between loop 3 and helix 5 of wild-type H-Ras are absent in the mutant A59G reflecting the enhanced dynamics of the loop 3 region. Taken together with earlier findings, these results suggest the existence of a lower energetic barrier between GTP and GDP states of the mutant. Molecular dynamics simulations combined with principal component analysis of available Ras crystallographic structures can be used to discriminate ligand- and sequence-based dynamic perturbations with potential functional implications. Furthermore, the identification of specific conformations associated with distinct Ras isoforms and mutants provides useful information for efforts that attempt to selectively interfere with the aberrant functions of these species.

Computational analysis of phosphopeptide binding to the polo-box domain of the mitotic kinase PLK1 using molecular dynamics simulation.

The Polo-Like Kinase 1 (PLK1) acts as a central regulator of mitosis and is over-expressed in a wide range of human tumours where high levels of expression correlate with a poor prognosis. PLK1 comprises two structural elements, a kinase domain and a polo-box domain (PBD). The PBD binds phosphorylated substrates to control substrate phosphorylation by the kinase domain. Although the PBD preferentially binds to phosphopeptides, it has a relatively broad sequence specificity in comparison with other phosphopeptide binding domains. We analysed the molecular determinants of recognition by performing molecular dynamics simulations of the PBD with one of its natural substrates, CDC25c. Predicted binding free energies were calculated using a molecular mechanics, Poisson-Boltzmann surface area approach. We calculated the per-residue contributions to the binding free energy change, showing that the phosphothreonine residue and the mainchain account for the vast majority of the interaction energy. This explains the very broad sequence specificity with respect to other sidechain residues. Finally, we considered the key role of bridging water molecules at the binding interface. We employed inhomogeneous fluid solvation theory to consider the free energy of water molecules on the protein surface with respect to bulk water molecules. Such an analysis highlights binding hotspots created by elimination of water molecules from hydrophobic surfaces. It also predicts that a number of water molecules are stabilized by the presence of the charged phosphate group, and that this will have a significant effect on the binding affinity. Our findings suggest a molecular rationale for the promiscuous binding of the PBD and highlight a role for bridging water molecules at the interface. We expect that this method of analysis will be very useful for probing other protein surfaces to identify binding hotspots for natural binding partners and small molecule inhibitors.

Unraveling evolutionary constraints: a heterogeneous conservation in dynamics of the titin Ig domains.

The giant protein titin, which comprises immunoglobulin (Ig) domains, acts as a bidirectional spring in muscle. The unfolding of Ig domains has been extensively studied, but their dynamics under native states have not been well-characterized. We performed molecular dynamics simulation on a single titin Ig domain and multi-domains. Mobile regions displaying concerted motions were identified. The dynamics of Ig domains are constrained by evolutionary pressures, in such a way that global dominant motion is conserved, yet different flexibilities within Ig domains and in linkers connecting neighbouring domains were observed. We explain these heterogeneous conserved dynamics in relation to sequence conservation across species and the sequence diversity among neighbouring Ig domains.

A network of dynamically conserved residues deciphers the motions of maltose transporter.

The maltose transporter of Escherichia coli is a member of the ATP-binding cassette (ABC) transporter superfamily. The crystal structures of maltose transporter MalK have been determined for distinct conformations in the presence and absence of the ligand ATP, and other interacting proteins. Using the distinct MalK structures, normal mode analysis was performed to understand the dynamics behavior of the system. A network of dynamically important residues was obtained from the normal mode analysis and the analysis of point mutation on the normal modes. Our results suggest that the intradomain rotation occurs earlier than the interdomain rotation during the maltose-binding protein (MBP)-driven conformational changes of MalK. We inquire if protein motion and functional-driven evolutionary conservation are related. The sequence conservation of MalK was analyzed to derive a network of evolutionarily important residues. There are highly significant correlations between protein sequence and protein dynamics in many regions on the maltose transporter MalK, suggesting a link between the protein evolution and dynamics. The significant overlaps between the network of dynamically important residues and the network of evolutionarily important residues form a network of dynamically conserved residues.

Dynamics of conserved waters in human Hsp90: implications for drug design.

The flexibility of a promising protein target, human heat shock protein 90 (Hsp90), is investigated by molecular dynamics simulations. These simulations focus on: (i) the interactions between the protein and conserved water molecules; and (ii) the interactions of the ligand PU3, the conserved water molecules and the protein. This is followed by a virtual screening docking study of the PU3 family of compounds and Hsp90 incorporating several conserved water molecules.

Classification of proteins based on similarity of two-dimensional protein maps.

Data reduction techniques are now a vital part of numerical analysis and principal component analysis is often used to identify important molecular features from a set of descriptors. We now take a different approach and apply data reduction techniques directly to protein structure. With this we can reduce the three-dimensional structural data into two-dimensions while preserving the correct relationships. With two-dimensional representations, structural comparisons between proteins are accelerated significantly. This means that protein-protein similarity comparisons are now feasible on a large scale. We show how the approach can help to predict the function of kinase structures according to the Hanks' classification based on their structural similarity to different kinase classes.

The effect of solvation on biomolecular conformation: 2-amino-1-phenylethanol.

Small molecule neurotransmitters form one the most important classes of pharmaceutical molecules. While the behavior of these molecules in their neutral forms in the gas phase is well understood, their behavior in more biologically relevant scenarios (protonated and in aqueous solution) has received comparatively little attention. Here we address this problem by using molecular mechanics simulations to build up a detailed picture of the conformational behavior of 2-amino-1-phenylethanol, a noradrenaline analogue, in aqueous solution in both its neutral and protonated forms. For the sake of comparison, equivalent simulations are also performed on the gas-phase molecules and gas-phase hydrated clusters. These calculations reveal the important role that water has to play in determining the conformational preferences and dynamic behavior of the molecules. Water molecules are found to bridge between the various functional groups within the molecule, significantly affecting their relative stabilities in comparison to the gas-phase values. The reorganization of these solvation structures also provides a mechanism for conformational interconversion. The role of the solvent in mediating interactions between the various functional groups within the molecule suggests that in noradrenaline the catechol groups will be able to interact, albeit indirectly, with the other functional groups, thereby influencing the behavior of the molecule.

Role of aromatic amino acids in protein-nucleic acid recognition.

Statistical analysis of structures from the PBD has been used to examine the role that the aromatic amino acids play in protein-nucleic acid recognition. In protein-DNA complexes, the residues Phe and His are found to bind selectively to the DNA chain--Phe to A and T, and His to T and G. The preferred binding modes are identified, and the interactions involving Phe are shown to be important in the transcription process. In protein-RNA complexes, Phe is found to occur far less often and is instead replaced by Trp, which binds selectively to C and G, offering a possible mechanism for differentiation between the two nucleic acids. SASA analysis of the two sets of complexes suggests that all of the aromatic amino acids are more heavily involved in binding than would be expected on the balance of probability. Phe and Tyr occur approximately equal in both sets of data, whereas the proportions of His and Trp vary considerably, supporting the idea that these residues may be involved in differentiating between the two nucleic acids.

Modeling Aromatic Liquids: Toluene, Phenol, and Pyridine.

Aromatic groups are now acknowledged to play an important role in many systems of interest. However, existing molecular mechanics methods provide a poor representation of these groups. In a previous paper, we have shown that the molecular mechanics treatment of benzene can be improved by the incorporation of an explicit representation of the aromatic π electrons. Here, we develop this concept further, developing charge-separation models for toluene, phenol, and pyridine. Monte Carlo simulations are used to parametrize the models, via the reproduction of experimental thermodynamic data, and our models are shown to outperform an existing atom-centered model. The models are then used to make predictions about the structures of the liquids at the molecular level and are tested further through their application to the modeling of gas-phase dimers and cation-π interactions.

Characterization and tissue-specific expression of two lepidopteran farnesyl diphosphate synthase homologs: implications for the biosynthesis of ethyl-substituted juvenile hormones.

The sesquiterpenoid juvenile hormone (JH) regulates insect development and reproduction. Most insects produce only one chemical form of JH, but the Lepidoptera produce four derivatives featuring ethyl branches. The biogenesis of these JHs requires the synthesis of ethyl-substituted farnesyl diphosphate (FPP) by FPP synthase (FPPS). To determine if there exist more than one lepidopteran FPPS, and whether one FPPS homolog is better adapted for binding the bulkier ethyl-branched substrates/products, we cloned three lepidopteran FPPS cDNAs, two from Choristoneura fumiferana and one from Pseudaletia unipuncta. Amino acid sequence comparisons among these and other eukaryotic FPPSs led to the recognition of two lepidopteran FPPS types. Type-I FPPSs display unique active site substitutions, including several in and near the first aspartate-rich motif, whereas type-II proteins have a more "conventional" catalytic cavity. In a yeast assay, a Drosophila FPPS clone provided full complementation of an FPPS mutation, but lepidopteran FPPS clones of either type yielded only partial complementation, suggesting unusual catalytic features and/or requirements of these enzymes. Although a structural analysis of lepidopteran FPPS active sites suggested that type-I enzymes are better suited than type-II for generating ethyl-substituted products, a quantitative real-time PCR assessment of their relative abundance in insect tissues indicated that type-I expression is ubiquitous whereas that of type-II is essentially confined to the JH-producing glands, where its transcripts are approximately 20 times more abundant than those of type-I. These results suggest that type-II FPPS plays a leading role in lepidopteran JH biosynthesis in spite of its apparently more conventional catalytic cavity.

A solvent induced mechanism for conformational change.

Molecular dynamics simulations have been used to investigate the dynamic behaviour of two small molecule neurotransmitter analogues in aqueous solution, leading to the elucidation of a mechanism for conformational change which is driven by the presence of the solvent molecules.

The Structure of Liquid Benzene.

The interactions of aromatic groups have been identified as playing a crucial role in many systems of interest. Unfortunately, conventional atom-centered force fields provide only an approximate representation of these molecules owing to their failure to consider the quadrupole moment arising from the π electrons. In this paper the structure of liquid benzene, the prototypical aromatic system, is investigated using a novel approach to Monte Carlo simulation, parametrized against experimental thermodynamic data, which incorporates an explicit representation of the aromatic π electrons. In contrast to previous simulations of liquid benzene it is found that a perpendicular arrangement of benzene molecules is preferred to a parallel arrangement. This result is in good agreement with experimental data.