A Computational Study

The dilemma of reconciling the contradictory evidence regarding the conformation of long solvated peptide chains is the so-called “reconciliation problem”. Clues regarding the stability of certain conformations likely lie in the electronic structure at the peptide–solvent interface, but the peptide–solvent interaction is not fully understood. Here, we study the influence of aqueous solvent on peptide conformations by using classical molecular dynamics (MD) and quantum mechanical/molecular mechanical (QM/MM) energy calculations. The model systems include an 11-residue peptide, X 2 A 7 O 2 (XAO), where X, A, and O denote diaminobutyric acid, alanine, and ornithine, respectively, and a 9-mer (Arg-Pro-Pro-Gly-Phe-Ser-Ala-Phe-Lys). Spectroscopic and MD data present conflicting evidence regarding the structure of XAO in water; some results indicate that XAO adopts a polyproline II (P II ) conformation, whereas other findings suggest that XAO explores a range of conformations. To investigate this contradiction, we present here the results of MD simulations of XAO and the 9-mer in aqueous solution, combined with QM/MM energy calculations

Probing the Structure of [NiFeSe] Hydrogenase with QM/MM Computations

The geometry and vibrational behavior of selenocysteine [NiFeSe] hydrogenase isolated from Desulfovibrio vulgaris Hildenborough have been investigated using a hybrid quantum mechanical (QM)/ molecular mechanical (MM) approach. Structural models have been built based on the three conformers identified in the recent crystal structure resolved at 1.3 Å from X-ray crystallography. In the models, a diamagnetic Ni2+ atom was modeled in combination with both Fe2+ and Fe3+ to investigate the effect of iron oxidation on geometry and vibrational frequency of the nonproteic ligands, CO and CN-, coordinated to the Fe atom. Overall, the QM/MM optimized geometries are in good agreement with the experimentally resolved geometries. Analysis of computed vibrational frequencies, in comparison with experimental Fourier-transform infrared (FTIR) frequencies, suggests that a mixture of conformers as well as Fe2+ and Fe3+ oxidation states may be responsible for the acquired vibrational spectra.DFG, 390540038, EXC 2008: UniSysCatDFG, 414044773, Open Access Publizieren 2019 - 2020 / Technische Universität BerlinEC/H2020/810856/EU/Twin to Illuminate Metals in Biology and Biocatalysis through Biospectroscopy/TIMB

Structural and electronic properties of the active site of [ZnFe] SulE

The function of the recently isolated sulerythrin (SulE) has been investigated using a combination of structural and electronic analyses based on quantum mechanical calculations. In the SulE structure of Fushinobu et al. (2003), isolated from a strictly aerobic archaeon, Sulfolobus tokadaii, a dioxygen-containing species was tentatively included at the active site during crystallographic refinement although the substrate specificity of SulE remains unclear. Studies have suggested that a structurally related enzyme, rubrerythrin, functions as a hydrogen peroxide reductase. Since SulE is a truncated version of rubrerythrin, the enzymes are hypothesized to function similarly. Hence, using available X-ray crystallography data (1.7 Å), we constructed various models of SulE containing a ZnII–Fe active site, differing in the nature of the substrate specificity (O2, H2O2), the oxidation level and the spin state of the iron ion, and the protonation states of the coordinating glutamate residues. Also, the substrate H2O2 is modeled in two possible configurations, differing in the orientation of the hydrogen atoms. Overall, the optimized geometries with an O2 substrate do not show good agreement with the experimentally resolved geometry. In contrast, excellent agreement between crystal structure arrangement and optimized geometries is achieved considering a H2O2 substrate and FeII in both spin states, when Glu92 is protonated. These results suggest that the dioxo species detected at the [ZnFe] active site of sulerythrin is H2O2, rather than an O2 molecule in agreement with experimental data indicating that only the diferrous oxidation state of the dimetal site in rubrerythrin reacts rapidly with H2O2. Based on our computations, we proposed a possible reaction pathway for substrate binding at the ZnFeII site of SulE with a H2O2 substrate. In this reaction pathway, Fe or another electron donor, such as NAD(P)H, catalyzes the reduction of H2O2 to water at the zinc–iron site

Cholesterol Transport in Wild-Type NPC1 and P691S: Molecular Dynamics Simulations Reveal Changes in Dynamical Behavior

The Niemann–Pick C1 (NPC1) protein is the main protein involved in NPC disease, a fatal lysosomal lipid storage disease. NPC1, containing 1278 amino acids, is comprised of three lumenal domains (N-terminal, middle lumenal, C-terminal) and a transmembrane (TM) domain that contains a five helix bundle referred to as the sterol-sensing domain (SSD). The exact purpose of the SSD is not known, but it is believed that the SSD may bind cholesterol, either as a part of the lipid trafficking pathway or as part of a signaling mechanism. A recent cryo-EM structure has revealed an itraconazole binding site (IBS) in the SSD of human NPC1. Using this structural data, we constructed a model of cholesterol-bound wild-type (WT) and mutant P691S and performed molecular dynamics (MD) simulations of each cholesterol-bound protein. For WT NPC1, cholesterol migrates laterally, in the direction of the lipid bilayer. In the case of P691S, cholesterol is observed for the first time to migrate away from the SSD toward the N-terminal domain via a putative tunnel that connects the IBS with the lumenal domains. Structural features of the IBS are analyzed to identify the causes for different dynamical behavior between cholesterol-bound WT and cholesterol-bound P691S. The side chain of Ser691 in the P691S mutant introduces a hydrogen bond network that is not present in the WT protein. This change is likely responsible for the altered dynamical behavior observed in the P691S mutant and helps explain the disrupted cholesterol trafficking behavior observed in experiments

Cholesterol transport in wild-type NPC1 and P691S: Molecular dynamics simulations reveal changes in dynamical behavior

The Niemann–Pick C1 (NPC1) protein is the main protein involved in NPC disease, a fatal lysosomal lipid storage disease. NPC1, containing 1278 amino acids, is comprised of three lumenal domains (N-terminal, middle lumenal, C-terminal) and a transmembrane (TM) domain that contains a five helix bundle referred to as the sterol-sensing domain (SSD). The exact purpose of the SSD is not known, but it is believed that the SSD may bind cholesterol, either as a part of the lipid trafficking pathway or as part of a signaling mechanism. A recent cryo-EM structure has revealed an itraconazole binding site (IBS) in the SSD of human NPC1. Using this structural data, we constructed a model of cholesterol-bound wild-type (WT) and mutant P691S and performed molecular dynamics (MD) simulations of each cholesterol-bound protein. For WT NPC1, cholesterol migrates laterally, in the direction of the lipid bilayer. In the case of P691S, cholesterol is observed for the first time to migrate away from the SSD toward the N-terminal domain via a putative tunnel that connects the IBS with the lumenal domains. Structural features of the IBS are analyzed to identify the causes for different dynamical behavior between cholesterol-bound WT and cholesterol-bound P691S. The side chain of Ser691 in the P691S mutant introduces a hydrogen bond network that is not present in the WT protein. This change is likely responsible for the altered dynamical behavior observed in the P691S mutant and helps explain the disrupted cholesterol trafficking behavior observed in experiments.DFG, 414044773, Open Access Publizieren 2019 - 2020 / Technische Universität Berli

Insight into Inhibitor Binding in the Eukaryotic Proteasome: Computations of the 20S CP

A combination of molecular dynamics (MD) simulations and computational analyses uncovers structural features that may influence substrate passage and exposure to the active sites within the proteolytic chamber of the 20S proteasome core particle (CP). MD simulations of the CP reveal relaxation dynamics in which the CP slowly contracts over the 54 ns sampling period. MD simulations of the SyringolinA (SylA) inhibitor within the proteolytic B 1 ring chamber of the CP indicate that favorable van der Waals and electrostatic interactions account for the predominant association of the inhibitor with the walls of the proteolytic chamber. The time scale required for the inhibitor to travel from the center of the proteolytic chamber to the chamber wall is on the order of 4 ns, accompanied by an average energetic stabilization of approximately −20 kcal/mol

Simulations of NPC1(NTD):NPC2 Protein Complex Reveal Cholesterol Transfer Pathways

The Niemann Pick type C (NPC) proteins, NPC1 and NPC2, are involved in the lysosomal storage disease, NPC disease. The formation of a NPC1–NPC2 protein–protein complex is believed to be necessary for the transfer of cholesterol and lipids out of the late endosomal (LE)/lysosomal (Lys) compartments. Mutations in either NPC1 or NPC2 can lead to an accumulation of cholesterol and lipids in the LE/Lys, the primary phenotype of the NPC disease. We investigated the NPC1(NTD)–NPC2 protein–protein complex computationally using two putative binding interfaces. A combination of molecular modeling and molecular dynamics simulations reveals atomic details that are responsible for interface stability. Cholesterol binding energies associated with each of the binding pockets for the two models are calculated. Analyses of the cholesterol binding in the two models support bidirectional ligand transfer when a particular interface is established. Based on the results, we propose that, depending on the location of the cholesterol ligand, a dynamical interface between the NPC2 and NPC1(NTD) proteins exists. Structural features of a particular interface can lower the energy barrier and stabilize the passage of the cholesterol substrate from NPC2 to NPC1(NTD)

Electronic and Structural Properties of the Double Cubane Iron-Sulfur Cluster

The double-cubane cluster (DCC) refers to an [Fe8S9] iron-sulfur complex that is otherwise only known to exist in nitrogenases. Containing a bridging µ2-S ligand, the DCC in the DCC-containing protein (DCCP) is covalently linked to the protein scaffold via six coordinating cysteine residues. In this study, the nature of spin coupling and the effect of spin states on the cluster’s geometry are investigated computationally. Using density functional theory (DFT) and a broken symmetry (BS) approach to study the electronic ground state of the system, we computed the exchange interaction between the spin-coupled spins of the four FeFe dimers contained in the DCC. This treatment yields results that are in excellent agreement with both computed and experimentally determined exchange parameters for analogously coupled di-iron complexes. Hybrid quantum mechanical (QM)/molecular mechanical (MM) geometry optimizations show that cubane cluster A closest to charged amino acid side chains (Arg312, Glu140, Lys146) is less compact than cluster B, indicating that electrons of the same spin in a charged environment seek maximum separation. Overall, this study provides the community with a fundamental reference for subsequent studies of DCCP, as well as for investigations of other [Fe8S9]-containing enzymes.DFG, 390540038, EXC 2008: Unifying Systems in Catalysis "UniSysCat"DFG, 414044773, Open Access Publizieren 2021 - 2022 / Technische Universität Berli

Modelling the active SARS-CoV-2 helicase complex as a basis for structure-based inhibitor design

The RNA helicase (non-structural protein 13, NSP13) of SARS-CoV-2 is essential for viral replication, and it is highly conserved among the coronaviridae family, thus a prominent drug target to treat COVID-19. We present here structural models and dynamics of the helicase in complex with its native substrates based on thorough analysis of homologous sequences and existing experimental structures. We performed and analysed microseconds of molecular dynamics (MD) simulations, and our model provides valuable insights to the binding of the ATP and ssRNA at the atomic level. We identify the principal motions characterising the enzyme and highlight the effect of the natural substrates on this dynamics. Furthermore, allosteric binding sites are suggested by our pocket analysis. Our obtained structural and dynamical insights are important for subsequent studies of the catalytic function and for the development of specific inhibitors at our characterised binding pockets for this promising COVID-19 drug target