Structural Dynamics of the Methyl-Coenzyme M Reductase Active Site Are Influenced by Coenzyme F430 Modifications.

Methyl-coenzyme M reductase (MCR) is a central player in methane biogeochemistry, governing methanogenesis and the anaerobic oxidation of methane (AOM) in methanogens and anaerobic methanotrophs (ANME), respectively. The prosthetic group of MCR is coenzyme F430, a nickel-containing tetrahydrocorphin. Several modified versions of F430 have been discovered, including the 172-methylthio-F430 (mtF430) used by ANME-1 MCR. Here, we employ molecular dynamics (MD) simulations to investigate the active site dynamics of MCR from Methanosarcina acetivorans and ANME-1 when bound to the canonical F430 compared to 172-thioether coenzyme F430 variants and substrates (methyl-coenzyme M and coenzyme B) for methane formation. Our simulations highlight the importance of the Gln to Val substitution in accommodating the 172 methylthio modification in ANME-1 MCR. Modifications at the 172 position disrupt the canonical substrate positioning in M. acetivorans MCR. However, in some replicates, active site reorganization to maintain substrate positioning suggests that the modified F430 variants could be accommodated in a methanogenic MCR. We additionally report the first quantitative estimate of MCR intrinsic electric fields that are pivotal in driving methane formation. Our results suggest that the electric field aligned along the CH3-S-CoM thioether bond facilitates homolytic bond cleavage, coinciding with the proposed catalytic mechanism. Structural perturbations, however, weaken and misalign these electric fields, emphasizing the importance of the active site structure in maintaining their integrity. In conclusion, our results deepen the understanding of MCR active site dynamics, the enzyme's organizational role in intrinsic electric fields for catalysis, and the interplay between active site structure and electrostatics.

Differences in Conformational Sampling and Intrinsic Electric Fields Drive Ion Binding in Telomeric and TERRA G-Quadruplexes.

The formation of G-quadruplexes (GQs) occurs in guanine-rich sequences of DNA and RNA, producing highly stable and structurally diverse noncanonical nucleic acid structures. GQs play crucial roles in regulating transcription, translation, and replication and maintaining the genome, among others; thus, changes to their structures can lead to diseases such as cancer. Previous studies using polarizable molecular dynamics simulations have shown differences in ion binding properties between telomeric and telomeric repeat-containing RNA GQs despite architectural similarities. Here, we used volume-based metadynamics and repulsive potential simulations in conjunction with polarizable force fields to quantify the impact of ion binding on the GQ dynamics and ion binding free energies. Furthermore, we describe how GQs exert electric fields on their surroundings to link dynamics with variations in the electronic structure. Our findings provide new insights into the energetic, physical, and conformational properties of GQs and expose subtle but important differences between DNA and RNA GQs with the same fold.

Differences in Conformational Sampling and Intrinsic Electric Fields Drive Ion Binding in Telomeric and TERRA G-Quadruplexes.

The formation of G-quadruplexes (GQs) occurs in guanine-rich sequences of DNA and RNA, producing highly stable and structurally diverse noncanonical nucleic acid structures. GQs play crucial roles in regulating transcription, translation, and replication; and maintaining the genome, among others, thus changes to their structures can lead to diseases such as cancer. Previous studies using polarizable molecular dynamics simulations have shown differences in ion binding properties between telomeric and TERRA GQs despite architectural similarities. Here, we used volume-based metadynamics and repulsive potential simulations in conjunction with polarizable force fields to quantify the impact of ion binding on GQ dynamics and ion binding free energies. Furthermore, we describe how GQs exert electric fields on their surroundings to link dynamics with variations in electronic structure. Our findings provide new insights into the energetic, physical, and conformational properties of GQs and expose subtle, but important, differences between DNA and RNA GQs with the same fold.

Monomeric and dimeric states of human ZO1-PDZ2 are functional partners of the SARS-CoV-2 E protein.

The Envelope (E) protein of SARS-CoV-2 plays a key role in virus maturation, assembly, and virulence mechanisms. The E protein is characterized by the presence of a PDZ-binding motif (PBM) at its C-terminus that allows it to interact with several PDZ-containing proteins in the intracellular environment. One of the main binding partners of the SARS-CoV-2 E protein is the PDZ2 domain of ZO1, a protein with a crucial role in the formation of epithelial and endothelial tight junctions (TJs). In this work, through a combination of analytical ultracentrifugation analysis and equilibrium and kinetic folding experiments, we show that ZO1-PDZ2 domain is able to fold in a monomeric state, an alternative form to the dimeric conformation that is reported to be functional in the cell for TJs assembly. Importantly, surface plasmon resonance (SPR) data indicate that the PDZ2 monomer is fully functional and capable of binding the C-terminal portion of the E protein of SARS-CoV-2, with a measured affinity in the micromolar range. Moreover, we present a detailed computational analysis of the complex between the C-terminal portion of E protein with ZO1-PDZ2, both in its monomeric conformation (computed as a high confidence AlphaFold2 model) and dimeric conformation (obtained from the Protein Data Bank), by using both polarizable and nonpolarizable simulations. Together, our results indicate both the monomeric and dimeric states of PDZ2 to be functional partners of the E protein, with similar binding mechanisms, and provide mechanistic and structural information about a fundamental interaction required for the replication of SARS-CoV-2.

TUPÃ: Electric field analyses for molecular simulations.

We introduce TUPÃ, a Python-based algorithm to calculate and analyze electric fields in molecular simulations. To demonstrate the features in TUPÃ, we present three test cases in which the orientation and magnitude of the electric field exerted by biomolecules help explain biological phenomena or observed kinetics. As part of TUPÃ, we also provide a PyMOL plugin to help researchers visualize how electric fields are organized within the simulation system. The code is freely available and can be obtained at https://mdpoleto.github.io/tupa/.

Integration of Experimental Data and Use of Automated Fitting Methods in Developing Protein Force Fields.

The development of accurate protein force fields has been the cornerstone of molecular simulations for the past 50 years. During this period, many lessons have been learned regarding the use of experimental target data and parameter fitting procedures. Here, we review recent advances in protein force field development. We discuss the recent emergence of polarizable force fields and the role of electronic polarization and areas in which additive force fields fall short. The use of automated fitting methods and the inclusion of additional experimental solution data during parametrization is discussed as a means to highlight possible routes to improve the accuracy of force fields even further.

Making it Rain: Cloud-Based Molecular Simulations for Everyone.

We present a user-friendly front-end for running molecular dynamics (MD) simulations using the OpenMM toolkit on the Google Colab framework. Our goals are (1) to highlight the usage of a cloud-computing scheme for educational purposes for a hands-on approach when learning MD simulations and (2) to exemplify how low-income research groups can perform MD simulations in the microsecond time scale. We hope this work facilitates teaching and learning of molecular simulation throughout the community.

Rotational Profiler: A Fast, Automated, and Interactive Server to Derive Torsional Dihedral Potentials for Classical Molecular Simulations.

Rotational Profiler provides an analytical algorithm to compute sets of classical torsional dihedral parameters by fitting an empirical energy profile to a reference one that can be obtained experimentally or by quantum-mechanical methods. The resulting profiles are compatible with the functional forms in the most widely used biomolecular force fields (e.g., GROMOS, AMBER, OPLS, and CHARMM). The linear least-squares regression method is used to generate sets of parameters that best satisfy the fitting. Rotational Profiler is free to use, analytical, and force field/package independent. The formalism is herein described, and its usage, in an interactive and automated manner, is made available as a Web server at http://rotprof.lncc.br.

ConfID: an analytical method for conformational characterization of small molecules using molecular dynamics trajectories.

MOTIVATION: The conformational space of small molecules can be vast and difficult to assess. Molecular dynamics (MD) simulations of free ligands in solution have been applied to predict conformational populations, but their characterization is often based on clustering algorithms or manual efforts. RESULTS: Here, we introduce ConfID, an analytical tool for conformational characterization of small molecules using MD trajectories. The evolution of conformational sampling and population frequencies throughout trajectories is calculated to check for sampling convergence while allowing to map relevant conformational transitions. The tool is designed to track conformational transition events and calculate time-dependent properties for each conformational population detected. AVAILABILITY AND IMPLEMENTATION: Toolkit and documentation are freely available at http://sbcb.inf.ufrgs.br/confid. CONTACT: marcelo.poleto@ufv.br or bigrisci@inf.ufrgs.br. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The Lazy Life of Lipid-Linked Oligosaccharides in All Life Domains.

Lipid-linked oligosaccharides (LLOs) play an important role in the N-glycosylation pathway as the donor substrate of oligosaccharyltransferases (OSTs), which are responsible for the en bloc transfer of glycan chains onto a nascent polypeptide. The lipid component of LLO in both eukarya and archaea consists of a dolichol, and an undecaprenol in prokarya, whereas the number of isoprene units may change between species. Given the potential relevance of LLOs and their related enzymes to diverse biotechnological applications, obtaining reliable LLO models from distinct domains of life could support further studies on complex formation and their processing by OSTs, as well as protein engineering on such systems. In this work, molecular modeling techniques, such as quantum mechanics calculations, molecular dynamics simulations, and metadynamics were employed to study eukaryotic (Glc3-Man9-GlcNAc2-PP-Dolichol), bacterial (Glc1-GalNAc5-Bac1-PP-Undecaprenol), and archaeal (Glc1-Man1-Gal1-Man1-Glc1-Gal1-Glc1-P-Dolichol) LLOs in membrane bilayers. Microsecond molecular dynamics simulations and metadynamics calculations of LLOs revealed that glycan chains are more prone to interact with the membrane lipid head groups, while the PP linkages are positioned at the lipid phosphate head groups level. The dynamics of isoprenoid chains embedded within the bilayer are described, and membrane dynamics and related properties are also investigated. Overall, there are similarities regarding the structure and dynamics of the eukaryotic, the bacterial, and the archaeal LLOs in bilayers, which can support the comprehension of their association with OSTs. These data may support future studies on the transferring mechanism of the oligosaccharide chain to an acceptor protein.

Discovery of novel West Nile Virus protease inhibitor based on isobenzonafuranone and triazolic derivatives of eugenol and indan-1,3-dione scaffolds.

The West Nile Virus (WNV) NS2B-NS3 protease is an attractive target for the development of therapeutics against this arboviral pathogen. In the present investigation, the screening of a small library of fifty-eight synthetic compounds against the NS2-NB3 protease of WNV is described. The following groups of compounds were evaluated: 3-(2-aryl-2-oxoethyl)isobenzofuran-1(3H)-ones; eugenol derivatives bearing 1,2,3-triazolic functionalities; and indan-1,3-diones with 1,2,3-triazolic functionalities. The most promising of these was a eugenol derivative, namely 4-(3-(4-allyl-2-methoxyphenoxy)-propyl)-1-(2-bromobenzyl)-1H-1,2,3-triazole (35), which inhibited the protease with IC50 of 6.86 µmol L-1. Enzyme kinetic assays showed that this derivative of eugenol presents competitive inhibition behaviour. Molecular docking calculations predicted a recognition pattern involving the residues His51 and Ser135, which are members of the catalytic triad of the WNV NS2B-NS3 protease.

Role of structural ions on the dynamics of the Pseudomonas fluorescens 07A metalloprotease.

The molecular dynamics of the Pseudomonas fluorescens 07A metalloprotease in the presence of structural Ca2+ and Mn2+ ions was evaluated. Seven Ca2+ ions are primarily bound to the C-terminus, while a divalent cation is located at the catalytic site, acting as a cofactor. The observed enzyme's experimental activity suggests that Mn2+ could compete for the active site of the enzyme with Ca2+, Zn2+ or other divalent cations, thus providing greater catalytic power to the enzyme. Our molecular dynamics simulations suggest that these ions partially protect the enzyme's structure from thermal denaturation. Moreover, our simulations have shown a collective movement of opening-closing of the active-site in simulations with structural Ca2+ and Mn2+ ions bound, leading to a proposal of a dynamical model of P. fluorescens 07A metalloprotease active and inactive conformations. These findings can support the development of measures to control the activity of P. fluorescens and other spoilage microorganism proteases.

Development of GROMOS-Compatible Parameter Set for Simulations of Chalcones and Flavonoids.

Chalcones and flavonoids constitute a large family of plant secondary metabolites that have been explored as a potential source of novel pharmaceutical products. While the simulation of these compounds by molecular dynamics (MD) can be a valuable strategy to assess their conformational properties and so further develop their role in drug discovery, there are no set of force field parameters specifically designed and experimentally validated for their conformational description in condensed phase. So the current work developed a new parameter set for MD simulations of these compounds' main scaffolds under GROMOS force field. We employed a protocol adjusting the atomic charges and torsional parameters to the respective quantum mechanical derived dipole moments and dihedrals rotational profiles, respectively. Experimental properties of organic liquids were used as references to the calculated values to validate the parameters. Additionally, metadynamics simulations were performed to evaluate the conformational space of complex chalcones and flavonoids, while NOE contacts during simulations were measured and compared to experimental data. Accordingly, the employed protocol allowed us to obtain force field parameters that reproduce well the target data and may be expected to contribute in more accurate computational studies on the biological/therapeutical role of such molecules.

An Unusual Intramolecular Halogen Bond Guides Conformational Selection.

PIK-75 is a phosphoinositide-3-kinase (PI3K) α-isoform-selective inhibitor with high potency. Although published structure-activity relationship data show the importance of the NO2 and the Br substituents in PIK-75, none of the published studies could correctly determine the underlying reason for their importance. In this publication, we report the first X-ray crystal structure of PIK-75 in complex with the kinase GSK-3ß. The structure shows an unusual U-shaped conformation of PIK-75 within the active site of GSK-3ß that is likely stabilized by an atypical intramolecular Br⋅⋅⋅NO2 halogen bond. NMR and MD simulations show that this conformation presumably also exists in solution and leads to a binding-competent preorganization of the PIK-75 molecule, thus explaining its high potency. We therefore suggest that the site-specific incorporation of halogen bonds could be generally used to design conformationally restricted bioactive substances with increased potencies.

Aromatic Rings Commonly Used in Medicinal Chemistry: Force Fields Comparison and Interactions With Water Toward the Design of New Chemical Entities.

The identification of lead compounds usually includes a step of chemical diversity generation. Its rationale may be supported by both qualitative (SAR) and quantitative (QSAR) approaches, offering models of the putative ligand-receptor interactions. In both scenarios, our understanding of which interactions functional groups can perform is mostly based on their chemical nature (such as electronegativity, volume, melting point, lipophilicity etc.) instead of their dynamics in aqueous, biological solutions (solvent accessibility, lifetime of hydrogen bonds, solvent structure etc.). As a consequence, it is challenging to predict from 2D structures which functional groups will be able to perform interactions with the target receptor, at which intensity and relative abundance in the biological environment, all of which will contribute to ligand potency and intrinsic activity. With this in mind, the aim of this work is to assess properties of aromatic rings, commonly used for drug design, in aqueous solution through molecular dynamics simulations in order to characterize their chemical features and infer their impact in complexation dynamics. For this, common aromatic and heteroaromatic rings were selected and received new atomic charge set based on the direction and module of the dipole moment from MP2/6-31G* calculations, while other topological terms were taken from GROMOS53A6 force field. Afterwards, liquid physicochemical properties were simulated for a calibration set composed by nearly 40 molecules and compared to their respective experimental data, in order to validate each topology. Based on the reliance of the employed strategy, we expanded the dataset to more than 100 aromatic rings. Properties in aqueous solution such as solvent accessible surface area, H-bonds availability, H-bonds residence time, and water structure around heteroatoms were calculated for each ring, creating a database of potential interactions, shedding light on features of drugs in biological solutions, on the structural basis for bioisosterism and on the enthalpic/entropic costs for ligand-receptor complexation dynamics.

Influence of Na+ and Mg2+ ions on RNA structures studied with molecular dynamics simulations.

The structure of ribonucleic acid (RNA) polymers is strongly dependent on the presence of, in particular Mg2+ cations to stabilize structural features. Only in high-resolution X-ray crystallography structures can ions be identified reliably. Here, we perform molecular dynamics simulations of 24 RNA structures with varying ion concentrations. Twelve of the structures were helical and the others complex folded. The aim of the study is to predict ion positions but also to evaluate the impact of different types of ions (Na+ or Mg2+) and the ionic strength on structural stability and variations of RNA. As a general conclusion Mg2+ is found to conserve the experimental structure better than Na+ and, where experimental ion positions are available, they can be reproduced with reasonable accuracy. If a large surplus of ions is present the added electrostatic screening makes prediction of binding-sites less reproducible. Distinct differences in ion-binding between helical and complex folded structures are found. The strength of binding (ΔG‡ for breaking RNA atom-ion interactions) is found to differ between roughly 10 and 26 kJ/mol for the different RNA atoms. Differences in stability between helical and complex folded structures and of the influence of metal ions on either are discussed.

ESI-MS/MS of expanded porphyrins: a look into their structure and aromaticity.

Electrospray mass spectrometry/mass spectrometry was used to investigate the gas-phase properties of protonated expanded porphyrins, in order to correlate those with their structure and conformation. We have selected five expanded meso-pentafluorophenyl porphyrins, respectively, a pair of oxidized/reduced fused pentaphyrins (22 and 24 π electrons), a pair of oxidized/reduced regular hexaphyrins (26 and 28 π electrons) and a regular doubly N-fused hexaphyrin (28 π electrons). The gas-phase behavior of the protonated species of oxidized and reduced expanded porphyrins is different. The oxidized species (aromatic Hückel systems) fragment more extensively, mainly by the loss of two HF molecules. The reduced species (Möbius aromatic or Möbius-like aromatic systems) fragment less than their oxidized counterparts because of their increased flexibility. The protonated regular doubly fused hexaphyrin (non-aromatic Hückel system) shows the least fragmentation even at higher collision energies. In general, cyclization through losses of HF molecules decreases from the aromatic Hückel systems to Möbius aromatic or Möbius-like aromatic systems to non-aromatic Hückel systems and is related to an increase in conformational distortion. Copyright © 2016 John Wiley & Sons, Ltd.

[28]Hexaphyrin derivatives for anion recognition in organic and aqueous media.

Hexaphyrin-based anion chemosensors are reported for the first time. The meso-hexakis(ethylenediamine)-substituted [28]hexaphyrins 2 and 3 revealed strong affinity for F(-), AcO(-) and H2PO4(-). Adsorption constants in aqueous media were determined on a gold piezoelectric crystal coated with 2 and 3. (1)H NMR titrations and molecular dynamics simulations showed the main interactions between hosts and guests.

A porcine circovirus-2 mutant isolated in Brazil contains low-frequency substitutions in regions of immunoprotective epitopes in the capsid protein.

Porcine circovirus-2 (PCV2) is the etiologic agent of several diseases in pigs, including multi-systemic wasting syndrome (PMWS). In this work, a new mutant PCV2b was isolated from PMWS-affected pigs on a Brazilian farm. Its genome showed high sequence similarity (>99% identity) to those from a group of emerging mutants isolated from cases of PMWS outbreaks in vaccinated pigs in China, the USA and South Korea. Here, we show that these isolates share a combination of low-frequency substitutions (single amino acid polymorphisms with a frequency of ≤25%) in the viral capsid protein, mainly in regions of immunoprotective epitopes, and an additional lysine residue at position 234. These isolates were phylogenetically grouped in the PCV2b clade, reinforcing the idea of the emergence of a new group of mutants PCV2b associated with outbreaks worldwide. The identification of these polymorphisms in the viral capsid highlights the importance of considering these isolates for the development of more-effective vaccines.