ABSTRACT
Cholinesterases are well-known and widely studied enzymes crucial to human health and involved in neurology, Alzheimer's, and lipid metabolism. The protonation pattern of active sites of cholinesterases influences all the chemical processes within, including reaction, covalent inhibition by nerve agents, and reactivation. Despite its significance, our comprehension of the fine structure of cholinesterases remains limited. In this study, we employed enhanced-sampling quantum-mechanical/molecular-mechanical calculations to show that cholinesterases predominantly operate as dynamic mixtures of two protonation states. The proton transfer between two non-catalytic glutamate residues follows the Grotthuss mechanism facilitated by a mediator water molecule. We show that this uncovered complexity of active sites presents a challenge for classical molecular dynamics simulations and calls for special treatment. The calculated proton transfer barrier of 1.65 kcal/mol initiates a discussion on the potential existence of two coupled low-barrier hydrogen bonds in the inhibited form of butyrylcholinesterase. These findings expand our understanding of structural features expressed by highly evolved enzymes and guide future advances in cholinesterase-related protein and drug design studies.
Subject(s)
Butyrylcholinesterase , Catalytic Domain , Molecular Dynamics Simulation , Protons , Butyrylcholinesterase/chemistry , Butyrylcholinesterase/metabolism , Humans , Hydrogen Bonding , Acetylcholinesterase/chemistry , Acetylcholinesterase/metabolism , Cholinesterases/chemistry , Cholinesterases/metabolismABSTRACT
We employed a structural bioinformatics approach to develop novel peptides with predicted affinity to the binding site for negative allosteric modulators (NAMs) of metabotropic glutamate receptor 5 (mGluR5). Primary screening in zebrafish (Danio rerio) revealed a stimulatory effect of two peptides, LCGM-10 and LCGM-15. Target validation studies using calcium ion flux imaging and a luciferase reporter assay confirmed mGluR5 as the target. LCGM-10 showed greater potency than LCGM-15; it was comparable to that of the mGluR5 NAM 2-methyl-6-(phenylethynyl) pyridine (MPEP). Rodent behavioral screening in the open field and elevated plus maze revealed increased locomotor activity in both tests after acute LCGM-10 treatment, supported by further analysis of home cage spontaneous locomotor activity (SLA). The stimulating effect of a single LCGM-10 administration on SLA was evident up to 60 min after administration and was not accompanied by hypokinetic rebound observed for caffeine. According to our results, LCGM-10 has therapeutic potential to treat hypo- and dyskinesias of various etiologies. Further investigation of LCGM-10 effects in the delay discounting model of impulsive choice in rats revealed reduced trait impulsivity after single and chronic administrations, suggesting potential implication for attention deficit hyperactivity disorder, obsessive compulsive disorder, and addictions.
ABSTRACT
Soybean is a leguminous plant cultivated in many countries and is considered important in the food industry due to the high levels of oil and protein content in the beans. The high demand for soybeans and its products in the industry requires the expansion of cultivation areas. Despite climatic restrictions, West Siberia is gradually expanding its area of soybean cultivation. In this study, we present the first analysis of the population structure and genetic diversity of the 175 soybean Glycine max breeding lines and varieties cultivated in West Siberia (103 accessions) and other regions of Russia (72 accessions), and we compare them with the cultivated soybean varieties from other geographical locations. Principal component analysis revealed several genetic clusters with different levels of genetic heterogeneity. Studied accessions are genetically similar to varieties from China, Japan, and the USA and are genetically distant to varieties from South Korea. Admixture analysis revealed four ancestry groups based on genetic ancestry and geographical origin, which are consistent with the regions of cultivation and origin of accessions and correspond to the principal component analysis result. Population statistics, including nucleotide diversity, Tajima's D, and linkage disequilibrium, are comparatively similar to those observed for studied accessions of a different origin. This study provides essential population and genetic information about the unique collection of breeding lines and varieties cultivated in West Siberia and other Russian regions to foster further evolutionary, genome-wide associations and functional breeding studies.
ABSTRACT
Classical molecular dynamics simulation is a powerful and established method of modern computational chemistry. Being able to obtain accurate information on molecular behavior is crucial to get valuable insights into structure-function relationships that translate into fundamental findings and practical applications. Active sites of enzymes are known to be particularly intricate, therefore, simpler non-polarizable force fields may provide an inaccurate description. In this work, we addressed this hypothesis in a case of a canonical serine triad protease trypsin in its complex with a substrate-mimicking inhibitor. We tested six modern and popular force fields to find that significantly diverging results may be obtained. Amber FB-15 and OPLS-AA/M turned out to model the active site incorrectly. Amber ff19sb and ff15ipq demonstrated mixed performance. The best performing force fields were CHARMM36m and Amber ff99sb-ildn, therefore, they are recommended for use with this and related systems. We speculate that a similar lack of cross-force field convergence may be characteristic of other enzymatic systems. Therefore, we advocate for careful consideration of different force fields in any study within the field of computational enzymology.
Subject(s)
Amber , Serine Proteases , Molecular Dynamics SimulationABSTRACT
Negative heterosis can occur on different economically important traits, but the exact biological mechanisms of this phenomenon are still unknown. The present study focuses on determining the genetic factors associated with negative heterosis in interspecific hybrids between domestic sheep (Ovis aries) and argali (Ovis ammon). One locus (rs417431015) associated with viability and two loci (rs413302370, rs402808951) associated with meat productivity were identified. One gene (ARAP2) was prioritized for viability and three for meat productivity (PDE2A, ARAP1, and PCDH15). The loci associated with meat productivity were demonstrated to fit the overdominant inheritance model and could potentially be involved int negative heterosis mechanisms.
ABSTRACT
Hybrid quantum mechanical/molecular mechanical (QM/MM) simulations fuel discoveries in many fields of science including computational biochemistry and enzymology. Development of more convenient tools leads to an increase in the number of works in which mechanical insights into enzymes' mode of operation are obtained. Most commonly, these tools feature hydrogen-capping (link atom) approach to provide coupling between QM and MM subsystems across a covalent bond. Extensive studies were conducted to provide a solid foundation for the correctness of such an approach when a bond to a nonpolar MM atom is considered. However, not every task may be accomplished this way. Certain scenarios of using QM/MM in computational enzymology encourage or even necessitate the incorporation of backbone atoms into the QM region. Two out of three backbone atoms are polar, and in QM/MM with electrostatic embedding, a neighboring link atom will be hyperpolarized. Several schemes to mitigate this effect were previously proposed alongside a rigorous assessment of quantitative effects on model systems. However, it was not clear whether they may translate into qualitatively different results and how link atom hyperpolarization may manifest itself in a real-life enzymological scenario. Here, we show that the consequences of such an artifact may be severe and may completely overturn the conclusions drawn from the simulations. Our case advocates for the use of charge redistribution schemes whenever intra-backbone QM/MM boundaries are considered. Moreover, we addressed how different boundary types and charge redistribution schemes influence backbone dynamics. We showed that the results are heavily dependent on which boundary MM terms are retained, with charge alteration being of secondary importance. In the worst case, only three intra-backbone boundaries may be used with relative confidence in the adequacy of resulting simulations, irrespective of the hyperpolarization mitigation scheme. Thus, advances in the field are certainly needed to fuel new discoveries. As of now, we believe that issues raised in this work might encourage authors in the field to report what boundaries, boundary MM terms, and charge redistribution schemes they are using, so their results may be correctly interpreted.
Subject(s)
Proteins , Quantum Theory , Proteins/chemistry , HydrogenABSTRACT
The cysteine protease from the tobacco etch virus (TEVp) is a well-known and widely utilized enzyme. TEVp's chymotrypsin-like fold is generally associated with serine catalytic triads that differ in terms of a reaction mechanism from the most well-studied papain-like cysteine proteases. The question of what dominates the TEVp mechanism, nucleophile identity, or structural composition has never been previously addressed. Here, we use enhanced sampling multiscale modeling to uncover that TEVp combines the features of two worlds in such a way that potentially hampers its activity. We show that TEVp cysteine is strictly in the anionic form in a free enzyme similar to papain. Peptide binding shifts the equilibrium toward the nucleophile's protonated form, characteristic of chymotrypsin-like proteases, although the cysteinyl anion form is still present and interconversion is rapid. This way cysteine protonation generates enzyme states that are a diversion from the most effective course of action, with only 13.2% of Michaelis complex sub-states able to initiate the reaction. As a result, we propose an updated view on the reaction mechanism catalyzed by TEVp. We also demonstrate that AlphaFold is able to construct protease-substrate complexes with high accuracy. We propose that our findings open a way for its industrious use in enzymological tasks. Unique features of TEVp discovered in this work open a discussion on the evolutionary history and trade-offs of optimizing serine triad-associated folds to cysteine as a nucleophile.
ABSTRACT
Organophosphate hydrolases are promising as potential biotherapeutic agents to treat poisoning with pesticides or nerve gases. However, these enzymes often need to be further engineered in order to become useful in practice. One example of such enhancement is the alteration of enantioselectivity of diisopropyl fluorophosphatase (DFPase). Molecular modeling techniques offer a unique opportunity to address this task rationally by providing a physical description of the substrate-binding process. However, DFPase is a metalloenzyme, and correct modeling of metal cations is a challenging task generally coming with a tradeoff between simulation speed and accuracy. Here, we probe several molecular mechanical parameter combinations for their ability to empower long simulations needed to achieve a quantitative description of substrate binding. We demonstrate that a combination of the Amber19sb force field with the recently developed 12-6 Ca2+ models allows us to both correctly model DFPase and obtain new insights into the DFP binding process.
Subject(s)
Calcium/chemistry , Calcium/metabolism , Molecular Dynamics Simulation , Phosphoric Triester Hydrolases/chemistry , Phosphoric Triester Hydrolases/metabolism , Catalytic Domain , Models, Molecular , Protein ConformationABSTRACT
The aim of the study was to develop better anxiolytics and antidepressants. We focused on GABA A receptors and the α2δ auxiliary subunit of V-gated Ca2+ channels as putative targets because they are established as mediators of efficacious anxiolytics, antidepressants, and anticonvulsants. We further focused on short peptides as candidate ligands because of their high safety and tolerability profiles. We employed a structural bioinformatics approach to develop novel tetrapeptides with predicted affinity to GABA A receptors and α2δ. In silico docking studies of one of these peptides, LCGA-17, showed a high binding score for both GABA A receptors and α2δ, combined with anxiolytic-like properties in a Danio rerio behavioral screen. LCGA-17 showed anxiolytic-like effects in the novel tank test, the light-dark box, and the social preference test, with efficacy comparable to fluvoxamine and diazepam. In binding assays using rat brain membranes, [3H]-LCGA-17 was competed more effectively by gabapentinoid ligands of α2δ than ligands of GABA A receptors, suggesting that α2δ represents a likely target for LCGA-17. [3H]-LCGA-17 binding to brain lysates was unaffected by competition with ligands for GABAB, glutamate, dopamine, serotonin, and other receptors, suggesting specific interaction with α2δ. Dose-finding studies in mice using acute administration of LCGA-17 (i.p.) demonstrated anxiolytic-like effects in the open field test, elevated plus maze, and marble burying tests, as well as antidepressant-like properties in the forced swim test. The anxiolytic effects were effectively blocked by bicuculline. Therefore, LCGA-17 is a novel candidate anxiolytic and antidepressant that may act through α2δ, with possible synergism by GABA A receptors.
ABSTRACT
Despite their economic value, sheep remain relatively poorly studied animals in terms of the number of known loci and genes associated with commercially important traits. This gap in our knowledge can be filled in by performing new genome-wide association studies (GWAS) or by re-analyzing previously documented data using novel powerful statistical methods. This study is focused on the search for new loci associated with meat productivity and carcass traits in sheep. With a multivariate approach applied to publicly available GWAS results, we identified eight novel loci associated with the meat productivity and carcass traits in sheep. Using an in silico follow-up approach, we prioritized 13 genes in these loci. One of eight novel loci near the FAM3C and WNT16 genes has been replicated in an independent sample of Russian sheep populations (N = 108). The novel loci were added to our regularly updated database increasing the number of known loci to more than 140.
Subject(s)
Genome-Wide Association Study/veterinary , Quantitative Trait Loci , Sheep/genetics , Animals , Computer Simulation , Cytokines/genetics , Meat Products , Multivariate Analysis , Phenotype , Wnt Proteins/geneticsABSTRACT
The N-glycosylation of immunoglobulin G (IgG) affects its structure and function. It has been demonstrated that IgG N-glycosylation patterns are inherited as complex quantitative traits. Genome-wide association studies identified loci harboring genes encoding enzymes directly involved in protein glycosylation as well as loci likely to be involved in regulation of glycosylation biochemical pathways. Many of these loci could be linked to immune functions and risk of inflammatory and autoimmune diseases. The aim of the present study was to discover and replicate new loci associated with IgG N-glycosylation and to investigate possible pleiotropic effects of these loci onto immune function and the risk of inflammatory and autoimmune diseases. We conducted a multivariate genome-wide association analysis of 23 IgG N-glycosylation traits measured in 8090 individuals of European ancestry. The discovery stage was followed up by replication in 3147 people and in silico functional analysis. Our study increased the total number of replicated loci from 22 to 29. For the discovered loci, we suggest a number of genes potentially involved in the control of IgG N-glycosylation. Among the new loci, two (near RNF168 and TNFRSF13B) were previously implicated in rare immune deficiencies and were associated with levels of circulating immunoglobulins. For one new locus (near AP5B1/OVOL1), we demonstrated a potential pleiotropic effect on the risk of asthma. Our findings underline an important link between IgG N-glycosylation and immune function and provide new clues to understanding their interplay.
Subject(s)
Genetic Loci/genetics , Genetic Pleiotropy/genetics , Genome-Wide Association Study/methods , Immunity/genetics , Immunoglobulin G/genetics , Alleles , Autoimmune Diseases/genetics , Cohort Studies , Computer Simulation , Gene Frequency , Genome-Wide Association Study/statistics & numerical data , Genotype , Glycosylation , Humans , Immunoglobulin G/metabolism , Inflammation/genetics , Multivariate Analysis , Phenotype , Polymorphism, Single Nucleotide , Quantitative Trait Loci/geneticsABSTRACT
One of the hallmarks of Alzheimer's disease (AD) is the deposition of amyloid plaques in the brain parenchyma, which occurs 7-15 years before the onset of cognitive symptoms of the pathology. Timely diagnostics of amyloid formations allows identifying AD at an early stage and initiating inhibitor therapy, delaying the progression of the disease. However, clinically used radiopharmaceuticals based on 11C and 18F are synchrotron-dependent and short-lived. The design of new metal-containing radiopharmaceuticals for AD visualization is of interest. The development of coordination compounds capable of effectively crossing the blood-brain barrier (BBB) requires careful selection of a ligand moiety, a metal chelating scaffold, and a metal cation, defining the method of supposed Aß visualization. In this review, we have summarized metal-containing drugs for positron emission tomography (PET), magnetic resonance imaging (MRI), and single-photon emission computed tomography (SPECT) imaging of Alzheimer's disease. The obtained data allow assessing the structure-ability to cross the BBB ratio.
Subject(s)
Alzheimer Disease/diagnosis , Brain/diagnostic imaging , Brain/metabolism , Diagnostic Imaging , Metals , Diagnostic Imaging/methods , Female , Humans , Metals/chemistry , Metals/metabolism , Multimodal Imaging/methods , Neuroimaging , Structure-Activity RelationshipABSTRACT
Quantum mechanics/molecular mechanics (QM/MM) maturation of an immunoglobulin (Ig) powered by supercomputation delivers novel functionality to this catalytic template and facilitates artificial evolution of biocatalysts. We here employ density functional theory-based (DFT-b) tight binding and funnel metadynamics to advance our earlier QM/MM maturation of A17 Ig-paraoxonase (WTIgP) as a reactibody for organophosphorus toxins. It enables regulation of biocatalytic activity for tyrosine nucleophilic attack on phosphorus. The single amino acid substitution l-Leu47Lys results in 340-fold enhanced reactivity for paraoxon. The computed ground-state complex shows substrate-induced ionization of the nucleophilic l-Tyr37, now H-bonded to l-Lys47, resulting from repositioning of l-Lys47. Multiple antibody structural homologs, selected by phenylphosphonate covalent capture, show contrasting enantioselectivities for a P-chiral phenylphosphonate toxin. That is defined by crystallographic analysis of phenylphosphonylated reaction products for antibodies A5 and WTIgP. DFT-b analysis using QM regions based on these structures identifies transition states for the favored and disfavored reactions with surprising results. This stereoselection analysis is extended by funnel metadynamics to a range of WTIgP variants whose predicted stereoselectivity is endorsed by experimental analysis. The algorithms used here offer prospects for tailored design of highly evolved, genetically encoded organophosphorus scavengers and for broader functionalities of members of the Ig superfamily, including cell surface-exposed receptors.
ABSTRACT
Microbial communities are self-controlled by repertoires of lethal agents, the antibiotics. In their turn, these antibiotics are regulated by bioscavengers that are selected in the course of evolution. Kinase-mediated phosphorylation represents one of the general strategies for the emergence of antibiotic resistance. A new subfamily of AmiN-like kinases, isolated from the Siberian bear microbiome, inactivates antibiotic amicoumacin by phosphorylation. The nanomolar substrate affinity defines AmiN as a phosphotransferase with a unique catalytic efficiency proximal to the diffusion limit. Crystallographic analysis and multiscale simulations revealed a catalytically perfect mechanism providing phosphorylation exclusively in the case of a closed active site that counteracts substrate promiscuity. AmiN kinase is a member of the previously unknown subfamily representing the first evidence of a specialized phosphotransferase bioscavenger.
ABSTRACT
Peptides are promising drug candidates due to high specificity and standout safety. Identification of bioactive peptides de novo using molecular docking is a widely used approach. However, current scoring functions are poorly optimized for peptide ligands. In this work, we present a novel algorithm PeptoGrid that rescores poses predicted by AutoDock Vina according to frequency information of ligand atoms with particular properties appearing at different positions in the target protein's ligand binding site. We explored the relevance of PeptoGrid ranking with a virtual screening of peptide libraries using angiotensin-converting enzyme and GABAB receptor as targets. A reasonable agreement between the computational and experimental data suggests that PeptoGrid is suitable for discovering functional leads.
Subject(s)
Drug Discovery , Molecular Docking Simulation , Peptide Library , Peptides/chemistry , Peptides/pharmacology , Algorithms , Animals , Computer Simulation , Molecular Dynamics Simulation , Reproducibility of Results , Structure-Activity Relationship , ZebrafishABSTRACT
Identification of quantitative trait loci (QTLs) and candidate genes that affect growth intensity is a prerequisite for the marker-assisted selection of economically important traits. The number of QTL studies on sheep is relatively small in comparison to those on cattle and pigs. The current QTL sheep database - Sheep QTLdb - contains information on 1658 QTLs for 225 different traits. A few genes and markers associated with growth, carcass and meat productivity traits have been reported. The information about QTLs from the Sheep QTLdb cannot be directly used in marker-assisted selection due to the lack of essential information such as effective and reference alleles, the effect direction etc., and it requires manual curation and validation. In this study we performed a comprehensive search for QTLs focusing on single nucleotide polymorphisms (SNPs) associated with growth and meat traits in sheep. The database contains information about 156 SNP-trait associations (123 unique SNPs) and a list of 165 associated genes. The updated information is freely available at https://github.com/Defrag1236/Ovines_2018 (last access: 18 September 2019). This information can be useful for further association studies and preliminary estimation of genetic variability for economically important traits in different breeds.
ABSTRACT
Butyrylcholinesterase (BChE) is considered as an efficient stoichiometric antidote against organophosphorus (OP) poisons. Recently we utilized combination of calculations and ultrahigh-throughput screening (uHTS) to select BChE variants capable of catalytic destruction of OP pesticide paraoxon. The purpose of this study was to elucidate the molecular mechanism underlying enzymatic hydrolysis of paraoxon by BChE variants using hybrid quantum mechanical/molecular mechanical (QM/MM) calculations. Detailed analysis of accomplished QM/MM runs revealed that histidine residues introduced into the acyl-binding loop are always located in close proximity with aspartate residue at position 70. Histidine residue acts as general base thus leading to attacking water molecule activation and subsequent SN2 inline hydrolysis resulting in BChE reactivation. This combination resembles canonical catalytic triad found in active centers of various proteases. Carboxyl group activates histidine residue by altering its pKa, which in turn promotes the activation of water molecule in terms of its nucleophilicity. Observed re-protonation of catalytic serine residue at position 198 from histidine residue at position 438 recovers initial configuration of the enzyme's active center, facilitating next catalytic cycle. We therefore suggest that utilization of uHTS platform in combination with deciphering of molecular mechanisms by QM/MM calculations may significantly improve our knowledge of enzyme function, propose new strategies for enzyme design and open new horizons in generation of catalytic bioscavengers against OP poisons.