| 1. |
- Barrozo, Alexandre, et al.
(författare)
-
Computer simulations of the catalytic mechanism of wild-type and mutant beta-phosphoglucomutase
- 2018
-
Ingår i: Organic and biomolecular chemistry. - : Royal Society of Chemistry. - 1477-0520 .- 1477-0539. ; 16:12, s. 2060-2073
-
Tidskriftsartikel (refereegranskat)abstract
- beta-Phosphoglucomutase (beta-PGM) has served as an important model system for understanding biological phosphoryl transfer. This enzyme catalyzes the isomerization of beta-glucose-1-phosphate to -glucose-6-phosphate in a two-step process proceeding via a bisphosphate intermediate. The conventionally accepted mechanism is that both steps are concerted processes involving acid-base catalysis from a nearby aspartate (D10) side chain. This argument is supported by the observation that mutation of D10 leaves the enzyme with no detectable activity. However, computational studies have suggested that a substrate-assisted mechanism is viable for many phosphotransferases. Therefore, we carried out empirical valence bond (EVB) simulations to address the plausibility of this mechanistic alternative, including its role in the abolished catalytic activity of the D10S, D10C and D10N point mutants of beta-PGM. In addition, we considered both of these mechanisms when performing EVB calculations of the catalysis of the wild type (WT), H20A, H20Q, T16P, K76A, D170A and E169A/D170A protein variants. Our calculated activation free energies confirm that D10 is likely to serve as the general base/acid for the reaction catalyzed by the WT enzyme and all its variants, in which D10 is not chemically altered. Our calculations also suggest that D10 plays a dual role in structural organization and maintaining electrostatic balance in the active site. The correct positioning of this residue in a catalytically competent conformation is provided by a functionally important conformational change in this enzyme and by the extensive network of H-bonding interactions that appear to be exquisitely preorganized for the transition state stabilization.
|
|
| 2. |
|
|
| 3. |
|
|
| 4. |
- Bauer, Paul, et al.
(författare)
-
Q6 : A comprehensive toolkit for empirical valence bond and related free energy calculations
- 2018
-
Ingår i: SoftwareX. - : Elsevier. - 2352-7110. ; 7, s. 388-395
-
Tidskriftsartikel (refereegranskat)abstract
- Atomistic simulations have become one of the main approaches to study the chemistry and dynamicsof biomolecular systems in solution. Chemical modelling is a powerful way to understand biochemistry,with a number of different programs available to perform specialized calculations. We present here Q6, anew version of the Q software package, which is a generalized package for empirical valence bond, linearinteraction energy, and other free energy calculations. In addition to general technical improvements, Q6extends the reach of the EVB implementation to fast approximations of quantum effects, extended solventdescriptions and quick estimation of the contributions of individual residues to changes in the activationfree energy of reactions.
|
|
| 5. |
- Esguerra, Mauricio, et al.
(författare)
-
GPCR-ModSim : A comprehensive web based solution for modeling G-protein coupled receptors
- 2016
-
Ingår i: Nucleic Acids Research. - : Oxford University Press (OUP). - 0305-1048 .- 1362-4962. ; 44:W1, s. W455-W462
-
Tidskriftsartikel (refereegranskat)abstract
- GPCR-ModSim (http://open.gpcr-modsim.org) is a centralized and easy to use service dedicated to the structural modeling of G-protein Coupled Receptors (GPCRs). 3D molecular models can be generated from amino acid sequence by homology-modeling techniques, considering different receptor conformations. GPCR-ModSim includes a membrane insertion and molecular dynamics (MD) equilibration protocol, which can be used to refine the generated model or any GPCR structure uploaded to the server, including if desired non-protein elements such as orthosteric or allosteric ligands, structural waters or ions. We herein revise the main characteristics of GPCR-ModSim and present new functionalities. The templates used for homology modeling have been updated considering the latest structural data, with separate profile structural alignments built for inactive, partially-active and active groups of templates. We have also added the possibility to perform multiple-template homology modeling in a unique and flexible way. Finally, our new MD protocol considers a series of distance restraints derived from a recently identified conserved network of helical contacts, allowing for a smoother refinement of the generated models which is particularly advised when there is low homology to the available templates. GPCR- ModSim has been tested on the GPCR Dock 2013 competition with satisfactory results.
|
|
| 6. |
- Jespers, Willem, et al.
(författare)
-
QligFEP : an automated workflow for small molecule free energy calculations in Q
- 2019
-
Ingår i: Journal of Cheminformatics. - : BMC. - 1758-2946. ; 11
-
Tidskriftsartikel (refereegranskat)abstract
- The process of ligand binding to a biological target can be represented as the equilibrium between the relevant solvated and bound states of the ligand. This which is the basis of structure-based, rigorous methods such as the estimation of relative binding affinities by free energy perturbation (FEP). Despite the growing capacity of computing power and the development of more accurate force fields, a high throughput application of FEP is currently hampered due to the need, in the current schemes, of an expert user definition of the alchemical transformations between molecules in the series explored. Here, we present QligFEP, a solution to this problem using an automated workflow for FEP calculations based on a dual topology approach. In this scheme, the starting poses of each of the two ligands, for which the relative affinity is to be calculated, are explicitly present in the MD simulations associated with the (dual topology) FEP transformation, making the perturbation pathway between the two ligands univocal. We show that this generalized method can be applied to accurately estimate solvation free energies for amino acid sidechain mimics, as well as the binding affinity shifts due to the chemical changes typical of lead optimization processes. This is illustrated in a number of protein systems extracted from other FEP studies in the literature: inhibitors of CDK2 kinase and a series of A(2A) adenosine G protein-coupled receptor antagonists, where the results obtained with QligFEP are in excellent agreement with experimental data. In addition, our protocol allows for scaffold hopping perturbations to identify the binding affinities between different core scaffolds, which we illustrate with a series of Chk1 kinase inhibitors. QligFEP is implemented in the open-source MD package Q, and works with the most common family of force fields: OPLS, CHARMM and AMBER.
|
|
| 7. |
- Lind, Christoffer, et al.
(författare)
-
A close-up view of codon selection in eukaryotic initiation
- 2017
-
Ingår i: RNA Biology. - : TAYLOR & FRANCIS INC. - 1547-6286 .- 1555-8584. ; 14:7, s. 815-819
-
Tidskriftsartikel (refereegranskat)abstract
- When given an option to choose among a set of alternatives and only one selection is right, one might stop and reflect over which one is best. However, the ribosome has no time to stop and make such reflections, proteins need to be produced and very fast. Eukaryotic translation initiation is an example of such a conundrum. Here, scanning for the correct codon match must be fast, efficient and accurate. We highlight our recent computational findings, which show how the initiation machinery manages to recognize one specific codon among many possible challengers, by fine-tuning the energetic landscape of base-pairing with the aid of the initiation factors eIF1 and eIF1A. Using a recent 3-dimensional structure of the eukaryotic initiation complex we have performed simulations of codon recognition in atomic detail. These calculations provide an in-depth energetic and structural view of how discrimination against near-cognate codons is achieved by the initiation complex.
|
|
| 8. |
- Lind, Christoffer, et al.
(författare)
-
Free energy calculations of RNA interactions
- 2019
-
Ingår i: Methods. - : Elsevier. - 1046-2023 .- 1095-9130. ; 162-163, s. 85-95
-
Tidskriftsartikel (refereegranskat)abstract
- This review discusses the use of molecular dynamics free energy calculations for characterizing RNA interactions, with particular emphasis on molecular recognition events involved in mRNA translation on the ribosome. The general methodology for efficient free energy calculations is outlined and our specific implementation for binding free energy changes due to base mutations in mRNA and tRNA is described, We show that there are a number of key problems related to the accuracy of protein synthesis that can be addressed with this type of computational approach and several such examples are discussed in detail. These include the decoding of mRNA during peptide chain elongation, initiation and termination of translation, as well as the energetic effects of base tautomerization and tRNA modifications. It is shown that free energy calculations can be made sufficiently reliable to allow quantitative conclusions to be drawn regarding the energetics of cognate versus non-cognate interactions and its structural origins.
|
|
| 9. |
- Vasile, Silvana, et al.
(författare)
-
Characterization of Ligand Binding to GPCRs Through Computational Methods
- 2018
-
Ingår i: Computational Methods for GPCR Drug Discovery. - New York, NY : Humana Press. - 9781493974641 - 9781493984947 - 9781493974658 ; , s. 23-44
-
Bokkapitel (refereegranskat)abstract
- The recent increase in available G protein-coupled receptor structures now contributes decisively to the structure-based ligand design. In this context, computational approaches in combination with medicinal chemistry and pharmacology are extremely helpful. Here, we provide an update on our structure-based computational protocols, used to answer key questions related to GPCR-ligand binding. All combined, these techniques can shed light on ligand binding modes, determine the molecular basis of conformational selection, for agonists and antagonists, as well as of subtype selectivity. To illustrate each of these questions, we will consider examples from existing projects on three families of class A (rhodopsin-like) GPCRs: one small-molecule (nucleotide-like) family, i.e., the adenosine receptors, and two peptide-binding receptors: neuropeptide-Y and angiotensin II receptors. The successful application of the same computational protocols to investigate this diverse group of receptor families gives an idea of the general applicability of our methodology in the characterization of GPCR-ligand binding.
|
|
