Delineating the mechanistic relevance of the TP53 gene and its mutational impact on gene expression and patients' survival in bladder cancer

Bladder carcinoma (BLCA) is a widespread urological malignancy causing significant global mortality, often hindered by delayed diagnosis and limited treatments. BLCA frequently exhibits TP53 mutations, playing a pivotal role in its pathogenesis and underscoring the potential of targeting TP53 as a therapeutic approach for this prevalent urological malignancy. Tumor tissues from 50 bladder cancer patients were used for mutational analysis in TP53's mutation-rich exons (5, 7, & 8). The gene expression of the TP53 gene, along with a TP53-target gene B-cell translocation gene 2 (BTG2) was also assessed in the cDNA samples from the same BLCA tissues and 15 urine controls of healthy people. The analysis revealed 22 % of patients with somatic hotspot mutations, 18 % with pathogenic missense mutations, and 12 % with intronic variants. Patients with somatic mutations exhibited the worst prognosis, supported by survival analysis from The Cancer Genome Atlas (TCGA) BLCA data. Interestingly, H296Y missense mutation correlated with higher TP53 expression and improved survival, while intronic SNPs were linked to worse outcomes. Additionally, upregulated BTG2 expression in mutated patients was observed which was correlated with poor prognosis, emphasizing the role of TP53 mutations in bladder cancer progression. The multivariate analysis highlighted the predictive power of TP53 mutations, with a high frequency of high-grade tumors (78.57 %) in mutated patients, underscoring their role in cancer progression. In conclusion, this study emphasizes the crucial role of TP53 mutations in bladder cancer patients from Bangladesh.

Designer high-density lipoprotein particles enhance endothelial barrier function and suppress inflammation

High-density lipoprotein (HDL) nanoparticles promote endothelial cell (EC) function and suppress inflammation, but their utility in treating EC dysfunction has not been fully explored. Here, we describe a fusion protein named ApoA1-ApoM (A1M) consisting of apolipoprotein A1 (ApoA1), the principal structural protein of HDL that forms lipid nanoparticles, and ApoM, a chaperone for the bioactive lipid sphingosine 1-phosphate (S1P). A1M forms HDL-like particles, binds to S1P, and is signaling competent. Molecular dynamics simulations showed that the S1P-bound ApoM moiety in A1M efficiently activated EC surface receptors. Treatment of human umbilical vein ECs with A1M-S1P stimulated barrier function either alone or cooperatively with other barrier-enhancing molecules, including the stable prostacyclin analog iloprost, and suppressed cytokine-induced inflammation. A1M-S1P injection into mice during sterile inflammation suppressed neutrophil influx and inflammatory mediator secretion. Moreover, systemic A1M administration led to a sustained increase in circulating HDL-bound S1P and suppressed inflammation in a murine model of LPS-induced endotoxemia. We propose that A1M administration may enhance vascular endothelial barrier function, suppress cytokine storm, and promote resilience of the vascular endothelium.

Uncovering water effects in protein-ligand recognition: importance in the second hydration shell and binding kinetics

Developing a ligand with high affinity for a specific protein target is essential for drug design, and water molecules are well known to play a key role in protein-drug recognition. However, predicting the role of particularly ordered water molecules in drug binding remains challenging. Furthermore, hydration free energy contributed from the water network, including the second shell of water molecules, is far from being well studied. In this research we focused on these aspects to accurately and efficiently evaluate water effects in protein-ligand binding affinity. We developed a new strategy using a free-energy calculation method, VM2. We successfully predicted the stable ordered water molecules in a number of protein systems: PDE 10a, HSP90, tryptophan synthase (TRPS), CDK2 and Factor Xa. In some of these, the second shell of water molecules appeared to be critical in protein-ligand binding. We also applied the strategy to largely improve binding free energy calculation using the MM/PBSA method. When applying MM/PBSA alone for two systems, CDK2 and Factor Xa, the computed binding free energy resulted in poor to moderate R² values with experimental data. However, including water free energy correction greatly improved the free energy calculation. Furthermore, our work helped to explain how xk263 is a 1000 times faster binder to HIVp than ritonavir, a potentially useful tool for investigating binding kinetics. Our studies reveal the importance of fully considering water effects in therapeutic developments in pharmaceutical and biotechnology industries and for fundamental research in protein-ligand recognition.

Water molecules at protein-drug interfaces: computational prediction and analysis methods

The fundamental importance of water molecules at drug-protein interfaces is now widely recognised and a significant feature in structure-based drug design. Experimental methods for analysing the role of water in drug binding have many challenges, including the accurate location of bound water molecules in crystal structures, and problems in resolving specific water contributions to binding thermodynamics. Computational analyses of binding site water molecules provide an alternative, and in principle complete, structural and thermodynamic picture, and their use is now commonplace in the pharmaceutical industry. In this review, we describe the computational methodologies that are available and discuss their strengths and weaknesses. Additionally, we provide a critical analysis of the experimental data used to validate the methods, regarding the type and quality of experimental structural data. We also discuss some of the fundamental difficulties of each method and suggest directions for future study.

Pyridazinones containing dithiocarbamoyl moieties as a new class of selective MAO-B inhibitors

A novel class of potential MAO-B inhibitors was designed and synthesized in good yield by combining the pyridazinone moiety with the dithiocarbamate framework, two relevant pharmacophores for drug discovery. The biological results obtained for the different pyridazinone/dithiocarbamate hybrids (compounds 8-14) indicated that most of them reversibly and selectively inhibit the hMAO-B in vitro with IC₅₀ values in the µM range and exhibit not significant cellular toxicity. The analogues 9a₁, 11a₁, 12a₂, 12b₁ and 12b₂, which present the dithiocarbamate fragment derivatized with a piperidin-1-yl or pyrrolidin-1-yl group and placed at C3 or C4 of the diazine ring, were the most attractive compounds of these series. Molecular modeling studies were performed to analyze the binding mode to the enzyme and the structure activity relationships of the titled compounds, as well as to predict their drug-like properties.

Assessment of CYP2C9 Structural Models for Site of Metabolism Prediction

Structure-based prediction of a compound's potential sites of metabolism (SOMs) mediated by cytochromes P450 (CYPs) is highly advantageous in the early stage of drug discovery. However, the accuracy of the SOMs prediction can be influenced by several factors. CYP2C9 is one of the major drug-metabolizing enzymes in humans and is responsible for the metabolism of ∼13 % of clinically used drugs. In this study, we systematically evaluated the effects of protein crystal structure models, scoring functions, heme forms, conserved active-site water molecules, and protein flexibility on SOMs prediction of CYP2C9 substrates. Our results demonstrated that, on average, ChemScore and GlideScore outperformed four other scoring functions: Vina, GoldScore, ChemPLP, and ASP. The performance of the crystal structure models with pentacoordinated heme was generally superior to that of the hexacoordinated iron-oxo heme (referred to as Compound I) models. Inclusion of the conserved active-site water molecule improved the prediction accuracy of GlideScore, but reduced the accuracy of ChemScore. In addition, the effect of the conserved water on SOMs prediction was found to be dependent on the receptor model and the substrate. We further found that one of snapshots from molecular dynamics simulations on the apo form can improve the prediction accuracy when compared to the crystal structural model.

Comparing Fragment Binding Poses Prediction Using HSP90 as a Key Study: When Bound Water Makes the Difference

Fragment-Based Drug Discovery (FBDD) approaches have gained popularitynot only in industry but also in academic research institutes. However, the computational prediction of the binding mode adopted by fragment-like molecules within a protein binding site is still a very challenging task. One of the most crucial aspects of fragment binding is related to the large amounts of bound waters in the targeted binding pocket. The binding affinity of fragmentsmay not be sufficientto displace the bound water molecules. In the present work, we confirmed the importance of the bound water molecules in the correct prediction of the fragment binding mode. Moreover, we investigate whether the use of methods based on explicit solvent molecular dynamics simulations can improve the accuracy of fragment posing. The protein chosen for this study is HSP-90.

PyVibMS: a PyMOL plugin for visualizing vibrations in molecules and solids

Visualizing vibrational motions calculated with different ab initio packages requires dedicated post-processing tools. Here, we present a PyMOL plugin called PyVibMS for visualizing the vibrational motions for both molecular and solid systems calculated by mainstream quantum chemical computer programs including Gaussian, Q-Chem, VASP, and CRYSTAL. Benefiting from the continuing development of the PyMOL platform, PyVibMS provides powerful functionalities and user-friendly interface. PyVibMS was written in Python and its open-source nature makes it flexible and sustainable. As an example, the motions of the Konkoli-Cremer local vibrational modes are shown in this work for the first time. PyVibMS is freely available at https://github.com/smutao/PyVibMS . Graphical abstract In this work, a PyMOL plugin named PyVibMS is developed to visualize molecular and lattice vibrations.

HomolWat: a web server tool to incorporate 'homologous' water molecules into GPCR structures

Internal water molecules play an essential role in the structure and function of membrane proteins including G protein-coupled receptors (GPCRs). However, technical limitations severely influence the number and certainty of observed water molecules in 3D structures. This may compromise the accuracy of further structural studies such as docking calculations or molecular dynamics simulations. Here we present HomolWat, a web application for incorporating water molecules into GPCR structures by using template-based modelling of homologous water molecules obtained from high-resolution structures. While there are various tools available to predict the positions of internal waters using energy-based methods, the approach of borrowing lacking water molecules from homologous GPCR structures makes HomolWat unique. The tool can incorporate water molecules into a protein structure in about a minute with around 85% of water recovery. The web server is freely available at http://lmc.uab.es/homolwat.

Applications of water molecules for analysis of macromolecule properties

Water molecules maintain proteins' structures, functions, stabilities and dynamics. They can occupy certain positions or pass quickly via a protein's interior. Regardless of their behaviour, water molecules can be used for the analysis of proteins' structural features and biochemical properties. Here, we present a list of several software programs that use the information provided by water molecules to: i) analyse protein structures and provide the optimal positions of water molecules for protein hydration, ii) identify high-occupancy water sites in order to analyse ligand binding modes, and iii) detect and describe tunnels and cavities. The analysis of water molecules' distribution and trajectories sheds a light on proteins' interactions with small molecules, on the dynamics of tunnels and cavities, on protein composition and also on the functionality, transportation network and location of functionally relevant residues. Finally, the correct placement of water molecules in protein crystal structures can significantly improve the reliability of molecular dynamics simulations.

Solvation Free Energy as a Measure of Hydrophobicity: Application to Serine Protease Binding Interfaces

Solvation and hydrophobicity play a key role in a variety of biological mechanisms. In substrate binding, but also in structure-based drug design, the thermodynamic properties of water molecules surrounding a given protein are of high interest. One of the main algorithms devised in recent years to quantify thermodynamic properties of water is the grid inhomogeneous solvation theory (GIST), which calculates these features on a grid surrounding the protein. Despite the inherent advantages of GIST, the computational demand is a major drawback, as calculations for larger systems can take days or even weeks. Here, we present a GPU accelerated version of the GIST algorithm, which facilitates efficient estimates of solvation free energy even of large biomolecular interfaces. Furthermore, we show that GIST can be used as a reliable tool to evaluate protein surface hydrophobicity. We apply the approach on a set of nine different proteases calculating localized solvation free energies on the surface of the binding interfaces as a measure of their hydrophobicity. We find a compelling agreement with the hydrophobicity of their substrates, i.e., peptides, binding into the binding cleft, and thus our approach provides a reliable description of hydrophobicity characteristics of these biological interfaces.

Prediction of Deleterious Non-synonymous SNPs of Human STK11 Gene by Combining Algorithms, Molecular Docking, and Molecular Dynamics Simulation

Serine-threonine kinase11 (STK11) is a tumor suppressor gene which plays a key role in regulating cell growth and apoptosis. It is widely known as a multitasking kinase and engaged in cell polarity, cell cycle arrest, chromatin remodeling, energy metabolism, and Wnt signaling. The substitutions of single amino acids in highly conserved regions of the STK11 protein are associated with Peutz-Jeghers syndrome (PJS), which is an autosomal dominant inherited disorder. The abnormal function of the STK11 protein is still not well understood. In this study, we classified disease susceptible single nucleotide polymorphisms (SNPs) in STK11 by using different computational algorithms. We identified the deleterious nsSNPs, constructed mutant protein structures, and evaluated the impact of mutation by employing molecular docking and molecular dynamics analysis. Our results show that W239R and W308C variants are likely to be highly deleterious mutations found in the catalytic kinase domain, which may destabilize structure and disrupt the activation of the STK11 protein as well as reduce its catalytic efficiency. The W239R mutant is likely to have a greater impact on destabilizing the protein structure compared to the W308C mutant. In conclusion, these mutants can help to further realize the large pool of disease susceptibilities linked with catalytic kinase domain activation of STK11 and assist to develop an effective drug for associated diseases.

The role of hydration effects in 5-fluorouridine binding to SOD1: insight from a new 3D-RISM-KH based protocol for including structural water in docking simulations

Misfolded Cu/Zn superoxide dismutase enzyme (SOD1) shows prion-like propagation in neuronal cells leading to neurotoxic aggregates that are implicated in amyotrophic lateral sclerosis (ALS). Tryptophan-32 (W32) in SOD1 is part of a potential site for templated conversion of wild type SOD1. This W32 binding site is located on a convex, solvent exposed surface of the SOD1 suggesting that hydration effects can play an important role in ligand recognition and binding. A recent X-ray crystal structure has revealed that 5-Fluorouridine (5-FUrd) binds at the W32 binding site and can act as a pharmacophore scaffold for the development of anti-ALS drugs. In this study, a new protocol is developed to account for structural (non-displaceable) water molecules in docking simulations and successfully applied to predict the correct docked conformation binding modes of 5-FUrd at the W32 binding site. The docked configuration is within 0.58 Å (RMSD) of the observed configuration. The docking protocol involved calculating a hydration structure around SOD1 using molecular theory of solvation (3D-RISM-KH, 3D-Reference Interaction Site Model-Kovalenko-Hirata) whereby, non-displaceable water molecules are identified for docking simulations. This protocol was also used to analyze the hydrated structure of the W32 binding site and to explain the role of solvation in ligand recognition and binding to SOD1. Structural water molecules mediate hydrogen bonds between 5-FUrd and the receptor, and create an environment favoring optimal placement of 5-FUrd in the W32 binding site.

Metabolic engineering of Methylobacterium extorquens AM1 for the production of butadiene precursor

BACKGROUND

Butadiene is a platform chemical used as an industrial feedstock for the manufacture of automobile tires, synthetic resins, latex and engineering plastics. Currently, butadiene is predominantly synthesized as a byproduct of ethylene production from non-renewable petroleum resources. Although the idea of biological synthesis of butadiene from sugars has been discussed in the literature, success for that goal has so far not been reported. As a model system for methanol assimilation, Methylobacterium extorquens AM1 can produce several unique metabolic intermediates for the production of value-added chemicals, including crotonyl-CoA as a potential precursor for butadiene synthesis.

RESULTS

In this work, we focused on constructing a metabolic pathway to convert crotonyl-CoA into crotyl diphosphate, a direct precursor of butadiene. The engineered pathway consists of three identified enzymes, a hydroxyethylthiazole kinase (THK) from Escherichia coli, an isopentenyl phosphate kinase (IPK) from Methanothermobacter thermautotrophicus and an aldehyde/alcohol dehydrogenase (ADHE2) from Clostridium acetobutylicum. The Km and kcat of THK, IPK and ADHE2 were determined as 8.35 mM and 1.24 s-1, 1.28 mM and 153.14 s-1, and 2.34 mM and 1.15 s-1 towards crotonol, crotyl monophosphate and crotonyl-CoA, respectively. Then, the activity of one of rate-limiting enzymes, THK, was optimized by random mutagenesis coupled with a developed high-throughput screening colorimetric assay. The resulting variant (THKM82V) isolated from over 3000 colonies showed 8.6-fold higher activity than wild-type, which helped increase the titer of crotyl diphosphate to 0.76 mM, corresponding to a 7.6% conversion from crotonol in the one-pot in vitro reaction. Overexpression of native ADHE2, IPK with THKM82V under a strong promoter mxaF in M. extorquens AM1 did not produce crotyl diphosphate from crotonyl-CoA, but the engineered strain did generate 0.60 μg/mL of intracellular crotyl diphosphate from exogenously supplied crotonol at mid-exponential phase.

CONCLUSIONS

These results represent the first step in producing a butadiene precursor in recombinant M. extorquens AM1. It not only demonstrates the feasibility of converting crotonol to key intermediates for butadiene biosynthesis, it also suggests future directions for improving catalytic efficiency of aldehyde/alcohol dehydrogenase to produce butadiene precursor from methanol.

Ligand Binding Free Energies with Adaptive Water Networks: Two-Dimensional Grand Canonical Alchemical Perturbations

Computational methods to calculate ligand binding affinities are increasing in popularity, due to improvements in simulation algorithms, computational resources, and easy-to-use software. However, issues can arise in relative ligand binding free energy simulations if the ligands considered have different active site water networks, as simulations are typically performed with a predetermined number of water molecules (fixed N ensembles) in preassigned locations. If an alchemical perturbation is attempted where the change should result in a different active site water network, the water molecules may not be able to adapt appropriately within the time scales of the simulations—particularly if the active site is occluded. By combining the grand canonical ensemble (μVT) to sample active site water molecules, with conventional alchemical free energy methods, the water network is able to dynamically adapt to the changing ligand. We refer to this approach as grand canonical alchemical perturbation (GCAP). In this work we demonstrate GCAP for two systems; Scytalone Dehydratase (SD) and Adenosine A_2A receptor. For both systems, GCAP is shown to perform well at reproducing experimental binding affinities. Calculating the relative binding affinities with a naïve, conventional attempt to solvate the active site illustrates how poor results can be if proper consideration of water molecules in occluded pockets is neglected. GCAP results are shown to be consistent with time-consuming double decoupling simulations. In addition, by obtaining the free energy surface for ligand perturbations, as a function of both the free energy coupling parameter and water chemical potential, it is possible to directly deconvolute the binding energetics in terms of protein–ligand direct interactions and protein binding site hydration.

Slow Dynamics of Tryptophan-Water Networks in Proteins

Water has a profound effect on the dynamics of biomolecules and governs many biological processes, leading to the concept that function is slaved to solvent dynamics within and surrounding the biomolecule. Protein conformational changes on μs-ms time scales are frequently associated with protein function, but little is known about the behavior of protein-bound water on these time scales. Here we have used NMR relaxation dispersion measurements to probe the tryptophan indoles in the enzyme dihydrofolate reductase (DHFR). We find that during structural changes on the μs-ms time scale, large chemical shift changes are often observed for the NH proton on the indole ring, while relatively smaller chemical shift changes are observed for the ring nitrogen atom. Comparison with experimental chemical shifts and density functional theory-based chemical shift predictions show that during the structural change the tryptophan indole NHs remain bound to water, but the geometry of the protein-bound water networks changes. These results establish that relaxation dispersion measurements can indirectly probe water dynamics and indicate that water can influence, or be influenced by, protein conformational changes on the μs-ms time scale. Our data show that structurally conserved bound water molecules can play a critical role in transmitting information between functionally important regions of the protein and provide evidence that internal protein motions can be coupled through the mediation of hydrogen-bonded water bound in the protein structure.

In Silico Studies of Small Molecule Interactions with Enzymes Reveal Aspects of Catalytic Function

Small molecules, such as solvent, substrate, and cofactor molecules, are key players in enzyme catalysis. Computational methods are powerful tools for exploring the dynamics and thermodynamics of these small molecules as they participate in or contribute to enzymatic processes. In-depth knowledge of how small molecule interactions and dynamics influence protein conformational dynamics and function is critical for progress in the field of enzyme catalysis. Although numerous computational studies have focused on enzyme-substrate complexes to gain insight into catalytic mechanisms, transition states and reaction rates, the dynamics of solvents, substrates, and cofactors are generally less well studied. Also, solvent dynamics within the biomolecular solvation layer play an important part in enzyme catalysis, but a full understanding of its role is hampered by its complexity. Moreover, passive substrate transport has been identified in certain enzymes, and the underlying principles of molecular recognition are an area of active investigation. Enzymes are highly dynamic entities that undergo different conformational changes, which range from side chain rearrangement of a residue to larger-scale conformational dynamics involving domains. These events may happen nearby or far away from the catalytic site, and may occur on different time scales, yet many are related to biological and catalytic function. Computational studies, primarily molecular dynamics (MD) simulations, provide atomistic-level insight and site-specific information on small molecule interactions, and their role in conformational pre-reorganization and dynamics in enzyme catalysis. The review is focused on MD simulation studies of small molecule interactions and dynamics to characterize and comprehend protein dynamics and function in catalyzed reactions. Experimental and theoretical methods available to complement and expand insight from MD simulations are discussed briefly.

The Roles of Water in the Protein Matrix: A Largely Untapped Resource for Drug Discovery

The value of thoroughly understanding the thermodynamics specific to a drug discovery/design study is well known. Over the past decade, the crucial roles of water molecules in protein structure, function, and dynamics have also become increasingly appreciated. This Perspective explores water in the biological environment by adopting its point of view in such phenomena. The prevailing thermodynamic models of the past, where water was seen largely in terms of an entropic gain after its displacement by a ligand, are now known to be much too simplistic. We adopt a set of terminology that describes water molecules as being "hot" and "cold", which we have defined as being easy and difficult to displace, respectively. The basis of these designations, which involve both enthalpic and entropic water contributions, are explored in several classes of biomolecules and structural motifs. The hallmarks for characterizing water molecules are examined, and computational tools for evaluating water-centric thermodynamics are reviewed. This Perspective's summary features guidelines for exploiting water molecules in drug discovery.

WONKA: objective novel complex analysis for ensembles of protein-ligand structures

WONKA is a tool for the systematic analysis of an ensemble of protein-ligand structures. It makes the identification of conserved and unusual features within such an ensemble straightforward. WONKA uses an intuitive workflow to process structural co-ordinates. Ligand and protein features are summarised and then presented within an interactive web application. WONKA's power in consolidating and summarising large amounts of data is described through the analysis of three bromodomain datasets. Furthermore, and in contrast to many current methods, WONKA relates analysis to individual ligands, from which we find unusual and erroneous binding modes. Finally the use of WONKA as an annotation tool to share observations about structures is demonstrated. WONKA is freely available to download and install locally or can be used online at http://wonka.sgc.ox.ac.uk.

Molecular Characterization of Legionellosis Drug Target Candidate Enzyme Phosphoglucosamine Mutase from Legionella pneumophila (strain Paris): An In Silico Approach

The harshness of legionellosis differs from mild Pontiac fever to potentially fatal Legionnaire's disease. The increasing development of drug resistance against legionellosis has led to explore new novel drug targets. It has been found that phosphoglucosamine mutase, phosphomannomutase, and phosphoglyceromutase enzymes can be used as the most probable therapeutic drug targets through extensive data mining. Phosphoglucosamine mutase is involved in amino sugar and nucleotide sugar metabolism. The purpose of this study was to predict the potential target of that specific drug. For this, the 3D structure of phosphoglucosamine mutase of Legionella pneumophila (strain Paris) was determined by means of homology modeling through Phyre2 and refined by ModRefiner. Then, the designed model was evaluated with a structure validation program, for instance, PROCHECK, ERRAT, Verify3D, and QMEAN, for further structural analysis. Secondary structural features were determined through self-optimized prediction method with alignment (SOPMA) and interacting networks by STRING. Consequently, we performed molecular docking studies. The analytical result of PROCHECK showed that 95.0% of the residues are in the most favored region, 4.50% are in the additional allowed region and 0.50% are in the generously allowed region of the Ramachandran plot. Verify3D graph value indicates a score of 0.71 and 89.791, 1.11 for ERRAT and QMEAN respectively. Arg419, Thr414, Ser412, and Thr9 were found to dock the substrate for the most favorable binding of S-mercaptocysteine. However, these findings from this current study will pave the way for further extensive investigation of this enzyme in wet lab experiments and in that way assist drug design against legionellosis.