Development and Validation of Fluorinated Amino Acid Parameters for use with the AMBER ff15ipq Protein Force Field

We developed force field parameters for fluorinated aromatic amino acids enabling molecular dynamics (MD) simulations of fluorinated proteins. These parameters are tailored to the AMBER ff15ipq protein force field and enable the modeling of 4, 5, 6, and 7F-tryptophan, 3F- and 3,5F-tyrosine, and 4F- or 4-CF3-phenylalanine. The parameters include 181 unique atomic charges derived using the Implicitly Polarized Charge (IPolQ) scheme in the presence of SPC/Eb explicit water molecules and 9 unique bond, angle, or torsion terms. Our simulations of benchmark peptides and proteins maintain expected conformational propensities on the μs-timescale. In addition, we have developed an open-source Python program to calculate fluorine relaxation rates from MD simulations. The extracted relaxation rates from protein simulations are in good agreement with experimental values determined by 19F NMR. Collectively, our results illustrate the power and robustness of the IPolQ lineage of force fields for modeling structure and dynamics of fluorine containing proteins at the atomic level.

The Omicron Variant Increases the Interactions of SARS-CoV-2 Spike Glycoprotein with ACE2

The Omicron variant (B.1.1.529) comprises 30 mutations on the spike glycoprotein (S), 15 of which are located on its receptor-binding domain (RBD_Omicron). RBD interacts with the peptidase domain (PD) of angiotensin-converting enzyme 2 (ACE2) receptors and plays a critical role in the host cell entry of the virus. We performed all-atom simulations of the RBD_Omicron-PD in the presence of explicit water and ions. Simulations showed a considerably more extensive interactions network between RBD_Omicron and PD compared to RBD_WT, comprising a 250%, 10% and -25% change in the total number of salt bridges, hydrophobic interactions, hydrogen bonds at the S-ACE2 interface, respectively. Using the conformations sampled in each our MD trajectories, binding energies of two sets of RBD_WT-PD and four sets of RBD_Omicron-PD simulations were calculated via the Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) method, estimating ~44% stronger binding energy for RBD_Omicron compared to RBD_WT. Our results suggest that an increase in the number of salt bridges in the S-ACE2 interface result in a higher binding strength of RBD to PD, which may result in a higher efficiency of the SARS-CoV-2 virus to infect host cells. Furthermore, RBD_Omicron exhibits a more dispersed interaction network on both sides of the RBD-PD interaction surface compared to WT.

Non-Markovian modeling of protein folding

We extract the folding free energy landscape and the time-dependent friction function, the two ingredients of the generalized Langevin equation (GLE), from explicit-water molecular dynamics (MD) simulations of the α-helix forming polypeptide alanine9 for a one-dimensional reaction coordinate based on the sum of the native H-bond distances. Folding and unfolding times from numerical integration of the GLE agree accurately with MD results, which demonstrate the robustness of our GLE-based non-Markovian model. In contrast, Markovian models do not accurately describe the peptide kinetics and in particular, cannot reproduce the folding and unfolding kinetics simultaneously, even if a spatially dependent friction profile is used. Analysis of the GLE demonstrates that memory effects in the friction significantly speed up peptide folding and unfolding kinetics, as predicted by the Grote–Hynes theory, and are the cause of anomalous diffusion in configuration space. Our methods are applicable to any reaction coordinate and in principle, also to experimental trajectories from single-molecule experiments. Our results demonstrate that a consistent description of protein-folding dynamics must account for memory friction effects.

Effect of Polyphosphorylation on Behavior of Protein Disordered Regions

Proteins interact with many charged biological macromolecules (polyelectrolytes), including inorganic polyphosphates. Recently a new protein post-translational modification, polyphosphorylation, or a covalent binding of polyphosphate chain to lysine, was demonstrated in human and yeast. Herein, we performed the first molecular modeling study of a possible effect of polyphosphorylation on behavior of the modified protein using replica exchange molecular dynamics simulations in atomistic force field with explicit water. Human endoplasmin (GRP-94), a member of heat shock protein 90 family, was selected as a model protein. Intrinsically disordered region in N-terminal domain serving as a charged linker between domains and containing a polyacidic serine and lysine-rich motif, was selected as a potent polyphosphorylation site according to literature data. Polyphosphorylation, depending on exact modification site, has been shown to influence on the disordered loop flexibility and induce its further expanding, as well as induce changes in interaction with ordered part of the molecule. As a result, polyphosphorylation in N-terminal domain might affect interaction of HSP90 with client proteins since these chaperones play a key role in protein folding.

Explanation of the Formation of Complexes between Representatives of Oxazolidinones and HDAS-β-CD Using Molecular Modeling as a Complementary Technique to cEKC and NMR

Molecular modeling (MM) results for tedizolid and radezolid with heptakis-(2,3-diacetyl-6-sulfo)-β-cyclodextrin (HDAS-β-CD) are presented and compared with the results previously obtained for linezolid and sutezolid. The mechanism of interaction of chiral oxazolidinone ligands belonging to a new class of antibacterial agents, such as linezolid, tedizolid, radezolid, and sutezolid, with HDAS-β-CD based on capillary electrokinetic chromatography (cEKC), nuclear magnetic resonance (NMR) spectroscopy, and MM methods was described. Principles of chiral separation of oxazolidinone analogues using charged single isomer derivatives of cyclodextrin by the cEKC method were presented, including the selection of the optimal chiral selector and separation conditions, complex stoichiometry, and binding constants, which provided a comprehensive basis for MM studies. In turn, NMR provided, where possible, direct information on the geometry of the inclusion complexes and also provided the necessary structural information to validate the MM calculations. Consequently, MM contributed to the understanding of the structure of diastereomeric complexes, the thermodynamics of complexation, and the visualization of their structures. The most probable mean geometries of the studied supramolecular complexes and their dynamics (geometry changes over time) were determined by molecular dynamics methods. Oxazolidinone ligands have been shown to complex mainly the inner part of cyclodextrin, while the external binding is less privileged, which is consistent with the conclusions of the NMR studies. Enthalpy values of binding of complexes were calculated using long-term molecular dynamics in explicit water as well as using molecular mechanics, the Poisson–Boltzmann or generalized Born, and surface area continuum solvation (MM/PBSA and MM/GBSA) methods. Computational methods predicted the effect of changes in pH and composition of the solution on the strength and complexation process, and it adapted the conditions selected as optimal during the cEKC study. By changing the dielectric constant in the MM/PBSA and MM/GBSA calculations, the effect of changing the solution to methanol/acetonitrile was investigated. A fairly successful attempt was made to predict the chiral separation of the oxazolidinones using the modified cyclodextrin by computational methods.