Homology modeling in the time of collective and artificial intelligence

Conformational Stability and Denaturation Processes of Proteins Investigated by Electrophoresis under Extreme Conditions

The functional structure of proteins results from marginally stable folded conformations. Reversible unfolding, irreversible denaturation, and deterioration can be caused by chemical and physical agents due to changes in the physicochemical conditions of pH, ionic strength, temperature, pressure, and electric field or due to the presence of a cosolvent that perturbs the delicate balance between stabilizing and destabilizing interactions and eventually induces chemical modifications. For most proteins, denaturation is a complex process involving transient intermediates in several reversible and eventually irreversible steps. Knowledge of protein stability and denaturation processes is mandatory for the development of enzymes as industrial catalysts, biopharmaceuticals, analytical and medical bioreagents, and safe industrial food. Electrophoresis techniques operating under extreme conditions are convenient tools for analyzing unfolding transitions, trapping transient intermediates, and gaining insight into the mechanisms of denaturation processes. Moreover, quantitative analysis of electrophoretic mobility transition curves allows the estimation of the conformational stability of proteins. These approaches include polyacrylamide gel electrophoresis and capillary zone electrophoresis under cold, heat, and hydrostatic pressure and in the presence of non-ionic denaturing agents or stabilizers such as polyols and heavy water. Lastly, after exposure to extremes of physical conditions, electrophoresis under standard conditions provides information on irreversible processes, slow conformational drifts, and slow renaturation processes. The impressive developments of enzyme technology with multiple applications in fine chemistry, biopharmaceutics, and nanomedicine prompted us to revisit the potentialities of these electrophoretic approaches. This feature review is illustrated with published and unpublished results obtained by the authors on cholinesterases and paraoxonase, two physiologically and toxicologically important enzymes.

Molecular Modeling of ABHD5 Structure and Ligand Recognition

Alpha/beta hydrolase domain-containing 5 (ABHD5), also termed CGI-58, is the key upstream activator of adipose triglyceride lipase (ATGL), which plays an essential role in lipid metabolism and energy storage. Mutations in ABHD5 disrupt lipolysis and are known to cause the Chanarin-Dorfman syndrome. Despite its importance, the structure of ABHD5 remains unknown. In this work, we combine computational and experimental methods to build a 3D structure of ABHD5. Multiple comparative and machine learning-based homology modeling methods are used to obtain possible models of ABHD5. The results from Gaussian accelerated molecular dynamics and experimental data of the apo models and their mutants are used to select the most likely model. Moreover, ensemble docking is performed on representative conformations of ABHD5 to reveal the binding mechanism of ABHD5 and a series of synthetic ligands. Our study suggests that the ABHD5 models created by deep learning-based methods are the best candidate structures for the ABHD5 protein. The mutations of E41, R116, and G328 disturb the hydrogen bonding network with nearby residues and suppress membrane targeting or ATGL activation. The simulations also reveal that the hydrophobic interactions are responsible for binding sulfonyl piperazine ligands to ABHD5. Our work provides fundamental insight into the structure of ABHD5 and its ligand-binding mode, which can be further applied to develop ABHD5 as a therapeutic target for metabolic disease and cancer.

Molecular Docking Study of IPBCC.08.610 Glucose Oxidase Mutant for Increasing Gluconic Acid Production

Glucose oxidase (GOD) is an oxidoreductase enzyme that catalyzes the oxidation of glucose to gluconolactone and hydrogen peroxide. Then, gluconolactone will be hydrolyzed to gluconic acid. The wide application of gluconic acid in various industries has increased production demand. However, glucose concentrations higher than 40% (w/w) inhibited the conversion of glucose to gluconic acid due to a decrease in the oxygen solubility concentration at pH 6, 30℃, and 1 bar pressure. Therefore, decreasing the value of Km is predicted to reduce saturation and enhance gluconic acid production. This study aimed to analyze the interaction between the IPBCC.08.610 GOD mutant with β-D-Glucose in improving gluconic acid production by decreasing the Km value. Mutations were performed in silico using Chimera and then docked using AutoDock Vina. The mutations resulted in distinct ligand poses in the binding pocket, different -OH conformations of the ligands, and changes in the T554M/D578P mutant’s hydrophobicity index (554 mutated from threonine to methionine, and 578 mutated from aspartate to proline), and decreased ΔG and Km values in the H559D mutant (559 mutated from histidine to aspartate), D578P and T554M/D578P. This decrease might strengthen the ligand-receptor interaction, increasing gluconic acid production. The H559D was the best mutant to increase production based on the ΔG, Km value, and stability due to the addition of hydrogen bonds.

How much metagenome data is needed for protein structure prediction: The advantages of targeted approach from the ecological and evolutionary perspectives

It has been proven that three-dimensional protein structures could be modeled by supplementing homologous sequences with metagenome sequences. Even though a large volume of metagenome data is utilized for such purposes, a significant proportion of proteins remain unsolved. In this review, we focus on identifying ecological and evolutionary patterns in metagenome data, decoding the complicated relationships of these patterns with protein structures, and investigating how these patterns can be effectively used to improve protein structure prediction. First, we proposed the metagenome utilization efficiency and marginal effect model to quantify the divergent distribution of homologous sequences for the protein family. Second, we proposed that the targeted approach effectively identifies homologous sequences from specified biomes compared with the untargeted approach's blind search. Finally, we determined the lower bound for metagenome data required for predicting all the protein structures in the Pfam database and showed that the present metagenome data is insufficient for this purpose. In summary, we discovered ecological and evolutionary patterns in the metagenome data that may be used to predict protein structures effectively. The targeted approach is promising in terms of effectively extracting homologous sequences and predicting protein structures using these patterns.

Peroxide-Mediated Oxygenation of Organic Compounds by Fungal Peroxygenases

Unspecific peroxygenases (UPOs), whose sequences can be found in the genomes of thousands of filamentous fungi, many yeasts and certain fungus-like protists, are fascinating biocatalysts that transfer peroxide-borne oxygen (from H2O2 or R-OOH) with high efficiency to a wide range of organic substrates, including less or unactivated carbons and heteroatoms. A twice-proline-flanked cysteine (PCP motif) typically ligates the heme that forms the heart of the active site of UPOs and enables various types of relevant oxygenation reactions (hydroxylation, epoxidation, subsequent dealkylations, deacylation, or aromatization) together with less specific one-electron oxidations (e.g., phenoxy radical formation). In consequence, the substrate portfolio of a UPO enzyme always combines prototypical monooxygenase and peroxidase activities. Here, we briefly review nearly 20 years of peroxygenase research, considering basic mechanistic, molecular, phylogenetic, and biotechnological aspects.

AI-based protein structure databases have the potential to accelerate rare diseases research: AlphaFoldDB and the case of IAHSP/Alsin

Artificial intelligence (AI)-based protein structure databases are expected to have an impact on drug discovery. Here, we show how AlphaFold could support rare diseases research programs. We focus on Alsin, a protein responsible for rare motor neuron diseases, such as infantile-onset ascending hereditary spastic paralysis (IAHSP) and juvenile primary lateral sclerosis (JPLS), and involved in some cases of amyotrophic lateral sclerosis (ALS). First, we compared the AlphaFoldDB human Alsin model with homology models of alsin domains. We then evaluated the flexibility profile of Alsin and of experimentally characterized mutants present in patients with IAHSP. Next, we compared preliminary models of dimeric/tetrameric Alsin responsible for its physiological action with hypothetical models reported in the literature. Finally, we suggest the best animal model for drug candidates testing. Overall, we computationally show that drug discovery efforts toward Alsin-involving diseases should be pursued.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

Characterization of a Stable Form of Carboxypeptidase G2 (Glucarpidase), a Potential Biobetter Variant, From Acinetobacter sp. 263903-1

Carboxypeptidase G2 (CPG2) is a bacterial enzyme widely used to detoxify methotrexate (MTX) and in enzyme/prodrug therapy for cancer treatment. However, several drawbacks, such as instability, have limited its efficiency. Herein, we have evaluated the properties of a putative CPG2 from Acinetobacter sp. 263903-1 (AcCPG2). AcCPG2 is compared with a CPG2 derived from Pseudomonas sp. strain RS-16 (PsCPG2), available as an FDA-approved medication called glucarpidase. After modeling AcCPG2 using the I-TASSER program, the refined model was validated by PROCHECK, VERIFY 3D and according to the Z score of the model. Using computational analyses, AcCPG2 displayed higher thermodynamic stability and a lower aggregation propensity than PsCPG2. AcCPG2 showed an optimum pH of 7.5 against MTX and was stable over a pH range of 5-10. AcCPG2 exhibited optimum activity at 50 °C and higher thermal stability at a temperature range of 20-70 °C compared to PsCPG2. The K_m value of the purified AcCPG2 toward folate and MTX was 31.36 µM and 44.99 µM, respectively. The V_max value of AcCPG2 for folate and MTX was 125.80 µmol/min/mg and 48.90  µmol/min/mg, respectively. Accordingly, thermostability and pH versatility makes AcCPG2 a potential biobetter variant for therapeutic applications.

Incorporating Machine Learning into Established Bioinformatics Frameworks

The exponential growth of biomedical data in recent years has urged the application of numerous machine learning techniques to address emerging problems in biology and clinical research. By enabling the automatic feature extraction, selection, and generation of predictive models, these methods can be used to efficiently study complex biological systems. Machine learning techniques are frequently integrated with bioinformatic methods, as well as curated databases and biological networks, to enhance training and validation, identify the best interpretable features, and enable feature and model investigation. Here, we review recently developed methods that incorporate machine learning within the same framework with techniques from molecular evolution, protein structure analysis, systems biology, and disease genomics. We outline the challenges posed for machine learning, and, in particular, deep learning in biomedicine, and suggest unique opportunities for machine learning techniques integrated with established bioinformatics approaches to overcome some of these challenges.