Validated ligand geometries for macromolecular refinement restraints and molecular mechanics force fields.

In macromolecular structure refinement the low observation-to-parameter ratio and the lack of high-resolution data is countered by using information in the form of restraints. Having accurate geometries of the chemical entities in the sample is paramount for generating accurate chemical restraints and, therefore, accurate macromolecular structures. In particular, it is desirable to have accurate restraints for known and novel ligand entities.

Recent Developments in Amber Biomolecular Simulations.

Amber is a molecular dynamics (MD) software package first conceived by Peter Kollman, his lab and collaborators to simulate biomolecular systems. The module is available as a serial version for central processing units (CPUs), NVIDIA and Advanced Micro Devices (AMD) graphics processing unit (GPU) versions as well as Message Passing Interface (MPI) parallel versions. Advanced capabilities include thermodynamic integration, replica exchange MD and accelerated MD methods.

Wavelength dependent photochemistry of an iron dinitrogen hydride complex multiple spectroscopies - competing ejection of axial ligands.

Nitrogenase (Nase) is a critical enzyme which catalyzes the reaction of N → NH in nature. Studies on the spectroscopy and photochemistry of -[Fe(DMeOPrPE)(N)H][BPh] (1) and its isotopologues (2-6) provide a possible first step to evaluate the geometries and properties of the real Nase-N structure(s). In this article, we have used FT-IR, FT-Raman, synchrotron-based nuclear resonant vibrational spectroscopy (NRVS) and DFT calculations to examine and assign the normal modes of these complexes.

Selective [9-N] Guanosine for Nuclear Magnetic Resonance Studies of Large Ribonucleic Acids.

RNAs regulate various cellular processes using malleable 3D structures, and understanding the factors that control RNA structure and dynamics is critical for understanding their mechanisms of action. To mitigate factors that have limited studies of large, functionally relevant RNAs by solution nuclear magnetic resonance (NMR) spectroscopy, we have extended a recently described H-enhanced, H-N correlation approach that used uniformly 15N-labeled guanosine triphosphate (GTP) by developing a chemoenzymatic labeling technology that grafts selectively labeled [9-N]-Guanine on to any labeled ribose to make [9-N]-GTP. The approach exploits advantageous NMR properties of the N9 nucleus which, when combined with extensive ribose deuteration and optimized NMR pulse sequences, affords sharp signals without complications that can arise using uniform [N]-guanine labeling.

Structure-Based Experimental Datasets for Benchmarking Protein Simulation Force Fields [Article v0.1].

This review article provides an overview of structurally oriented experimental datasets that can be used to benchmark protein force fields, focusing on data generated by nuclear magnetic resonance (NMR) spectroscopy and room temperature (RT) protein crystallography. We discuss what the observables are, what they tell us about structure and dynamics, what makes them useful for assessing force field accuracy, and how they can be connected to molecular dynamics simulations carried out using the force field one wishes to benchmark. We also touch on statistical issues that arise when comparing simulations with experiment.

Complex water networks visualized by cryogenic electron microscopy of RNA.

The stability and function of biomolecules are directly influenced by their myriad interactions with water. Here we investigated water through cryogenic electron microscopy (cryo-EM) on a highly solvated molecule: the Tetrahymena ribozyme. By using segmentation-guided water and ion modelling (SWIM), an approach combining resolvability and chemical parameters, we automatically modelled and cross-validated water molecules and Mg ions in the ribozyme core, revealing the extensive involvement of water in mediating RNA non-canonical interactions.

Complex Water Networks Visualized through 2.2-2.3 Å Cryogenic Electron Microscopy of RNA.

The stability and function of biomolecules are directly influenced by their myriad interactions with water. In this study, we investigated water through cryogenic electron microscopy (cryo-EM) on a highly solvated molecule, the ribozyme, determined at 2.2 and 2.

Simulating the Motion Underlying the Mechanism of Thioredoxin Reductase.

Thioredoxin reductase (TrxR) is an essential antioxidant in most cells; it reduces thioredoxin (Trx) and several more substrates, utilizing NADPH. However, the enzyme's internal active site is too small to accommodate the Trx substrate. Thus, TrxR evolved a disulfide shuttle that can carry reducing equivalents from the active site to the docking site of thioredoxin on the enzyme surface.

Enhanced TROSY Effect in [2- F, 2- C] Adenosine and ATP Analogs Facilitates NMR Spectroscopy of Very Large Biological RNAs in Solution.

Large RNAs are central to cellular functions, but characterizing such RNAs remains challenging by solution NMR. We present two labeling technologies based on [2- F, 2- C]-adenosine, which allow the incorporation of aromatic F- C spin pairs. The labels when coupled with the transverse relaxation optimized spectroscopy (TROSY) enable us to probe RNAs comprising up to 124 nucleotides.

X-ray Crystallography Module in MD Simulation Program Amber 2023. Refining the Models of Protein Crystals.

The MD simulation package Amber offers an attractive platform to refine crystallographic structures of proteins: (i) state-of-the-art force fields help to regularize protein coordinates and reconstruct the poorly diffracting elements of the structure, such as flexible loops; (ii) MD simulations restrained by the experimental diffraction data provide an effective strategy to optimize structural models of protein crystals, including explicitly modeled interstitial solvent as well as crystal contacts. Here, we present the new crystallography module , released as a part of the Amber 2023 package. This module contains functions to calculate and scale structure factors (including the contributions from bulk solvent), evaluate the maximum-likelihood-type crystallographic potential, and compute its derivative forces.

An RNA excited conformational state at atomic resolution.

Sparse and short-lived excited RNA conformational states are essential players in cell physiology, disease, and therapeutic development, yet determining their 3D structures remains challenging. Combining mutagenesis, NMR spectroscopy, and computational modeling, we determined the 3D structural ensemble formed by a short-lived (lifetime ~2.1 ms) lowly-populated (~0.

AmberTools.

AmberTools is a free and open-source collection of programs used to set up, run, and analyze molecular simulations. The newer features contained within AmberTools23 are briefly described in this Application note.

MD simulations of macromolecular crystals: Implications for the analysis of Bragg and diffuse scattering.

Some of our most detailed information about structure and dynamics of macromolecules comes from X-ray-diffraction studies in crystalline environments. More than 170,000 atomic models have been deposited in the Protein Data Bank, and the number of observations (typically of intensities of Bragg diffraction peaks) is generally quite large, when compared to other experimental methods. Nevertheless, the general agreement between calculated and observed intensities is far outside the experimental precision, and the majority of scattered photons fall between the sharp Bragg peaks, and are rarely taken into account.

New experimental evidence for pervasive dynamics in proteins.

There is ample computational, but only sparse experimental data suggesting that pico-ns motions with 1 Å amplitude are pervasive in proteins in solution. Such motions, if present in reality, must deeply affect protein function and protein entropy. Several NMR relaxation experiments have provided insights into motions of proteins in solution, but they primarily report on azimuthal angle variations of vectors of covalently-linked atoms.

Robust total X-ray scattering workflow to study correlated motion of proteins in crystals.

The breathing motions of proteins are thought to play a critical role in function. However, current techniques to study key collective motions are limited to spectroscopy and computation. We present a high-resolution experimental approach based on the total scattering from protein crystals at room temperature (TS/RT-MX) that captures both structure and collective motions.

RNA Electrostatics: How Ribozymes Engineer Active Sites to Enable Catalysis.

Electrostatic interactions are fundamental to RNA structure and function, and intimately influenced by solvation and the ion atmosphere. RNA enzymes, or ribozymes, are catalytic RNAs that are able to enhance reaction rates over a million-fold, despite having only a limited repertoire of building blocks and available set of chemical functional groups. Ribozyme active sites usually occur at junctions where negatively charged helices come together, and in many cases leverage this strained electrostatic environment to recruit metal ions in solution that can assist in catalysis.

Colicin E1 opens its hinge to plug TolC.

The double membrane architecture of Gram-negative bacteria forms a barrier that is impermeable to most extracellular threats. Bacteriocin proteins evolved to exploit the accessible, surface-exposed proteins embedded in the outer membrane to deliver cytotoxic cargo. Colicin E1 is a bacteriocin produced by, and lethal to, that hijacks the outer membrane proteins (OMPs) TolC and BtuB to enter the cell.

Integral equation models for solvent in macromolecular crystals.

The solvent can occupy up to ∼70% of macromolecular crystals, and hence, having models that predict solvent distributions in periodic systems could improve the interpretation of crystallographic data. Yet, there are few implicit solvent models applicable to periodic solutes, and crystallographic structures are commonly solved assuming a flat solvent model. Here, we present a newly developed periodic version of the 3D-reference interaction site model (RISM) integral equation method that is able to solve efficiently and describe accurately water and ion distributions in periodic systems; the code can compute accurate gradients that can be used in minimizations or molecular dynamics simulations.

Simulations of Kindlin-2 PIP binding domains reveal protonation-dependent membrane binding modes.

Kindlin-2, a member of the Kindlin family of peripheral membrane proteins, is important for integrin activation and stabilization of epidermal growth factor receptor. It associates with the cytoplasmic face of the plasma membrane via dedicated phosphatidylinositol phosphate binding domains located in the N-terminal F0 and Pleckstrin Homology domains. These domains have binding affinity for phosphatidylinositol 4,5-bisphosphate and, to a greater degree, phosphatidylinositol 3,4,5-trisphosphate.

Refinement of RNA Structures Using Amber Force Fields.

Atomic models for nucleic acids derived from X-ray diffraction data at low resolution provide much useful information, but the observed scattering intensities can be fit with models that can differ in structural detail. Tradtional geometric restraints favor models that have bond length and angle terms derived from small molecule crystal structures. Here we explore replacing these restraints with energy gradients derived from force fields, including recently developed integral equation models to account for the effects of water molecules and ions that are not part of the explicit model.

Molecular dynamics analysis of a flexible loop at the binding interface of the SARS-CoV-2 spike protein receptor-binding domain.

Since the identification of the SARS-CoV-2 virus as the causative agent of the current COVID-19 pandemic, considerable effort has been spent characterizing the interaction between the Spike protein receptor-binding domain (RBD) and the human angiotensin converting enzyme 2 (ACE2) receptor. This has provided a detailed picture of the end point structure of the RBD-ACE2 binding event, but what remains to be elucidated is the conformation and dynamics of the RBD prior to its interaction with ACE2. In this work, we utilize molecular dynamics simulations to probe the flexibility and conformational ensemble of the unbound state of the receptor-binding domain from SARS-CoV-2 and SARS-CoV.

Rapid and accurate determination of atomistic RNA dynamic ensemble models using NMR and structure prediction.

Biomolecules form dynamic ensembles of many inter-converting conformations which are key for understanding how they fold and function. However, determining ensembles is challenging because the information required to specify atomic structures for thousands of conformations far exceeds that of experimental measurements. We addressed this data gap and dramatically simplified and accelerated RNA ensemble determination by using structure prediction tools that leverage the growing database of RNA structures to generate a conformation library.