Role of methyl groups in dynamics and evolution of biomolecules

De novo amyloid peptides with subtle sequence variations differ in their self-assembly and nanomechanical properties

Proteinaceous amyloids are well known for their widespread pathological roles but lately have emerged also as key components in several biological functions. The remarkable ability of amyloid fibers to form tightly packed conformations in a cross β-sheet arrangement manifests in their robust enzymatic and structural stabilities. These characteristics of amyloids make them attractive for designing proteinaceous biomaterials for various biomedical and pharmaceutical applications. In order to design customizable and tunable amyloid nanomaterials, it is imperative to understand the sensitivity of the peptide sequence for subtle changes based on amino acid position and chemistry. Here we report our results from four rationally-designed amyloidogenic decapeptides that subtly differ in hydrophobicity and polarity at positions 5 and 6. We show that making the two positions hydrophobic renders the peptide with enhanced aggregation and material properties while introducing polar residues in position 5 dramatically changes the structure and nanomechanical properties of the fibrils formed. A charged residue at position 6, however, abrogates amyloid formation. In sum, we show that subtle changes in the sequence do not make the peptide innocuous but rather sensitive to aggregation, reflected in the biophysical and nanomechanical properties of the fibrils. We conclude that tolerance of peptide amyloid for changes in the sequence, however small they may be, should not be neglected for the effective design of customizable amyloid nanomaterials.

The effect of a methyl group on structure and function: Serine vs. threonine glycosylation and phosphorylation

A variety of glycan structures cover the surface of all cells and are involved in myriad biological processes, including but not limited to, cell adhesion and communication, protein quality control, signal transduction and metabolism, while also being intimately involved in innate and adaptive immune functions. Immune surveillance and responses to foreign carbohydrate antigens, such as capsular polysaccharides on bacteria and surface protein glycosylation of viruses, are the basis of microbial clearance, and most antimicrobial vaccines target these structures. In addition, aberrant glycans on tumors called Tumor-Associated Carbohydrate Antigens (TACAs) elicit immune responses to cancer, and TACAs have been used in the design of many antitumor vaccine constructs. A majority of mammalian TACAs are derived from what are referred to as mucin-type O-linked glycans on cell-surface proteins and are linked to the protein backbone through the hydroxyl group of either serine or threonine residues. A small group of structural studies that have compared mono- and oligosaccharides attached to each of these residues have shown that there are distinct differences in conformational preferences assumed by glycans attached to either "unmethylated" serine or ß-methylated threonine. This suggests that the linkage point of antigenic glycans will affect their presentation to the immune system as well as to various carbohydrate binding molecules (e.g., lectins). This short review, followed by our hypothesis, will examine this possibility and extend the concept to the presentation of glycans on surfaces and in assay systems where recognition of glycans by proteins and other binding partners can be defined by different attachment points that allow for a range of conformational presentations.

A concept of a broadband inverted geometry spectrometer for the Second Target Station at the Spallation Neutron Source

BWAVES is an acronym for Broadband Wide-Angle VElocity Selector spectrometer, indicating that a novel WAVES (Wide-Angle VElocity Selector) device will be used to select the velocity/wavelength of the detected neutrons after they are scattered by the sample. We describe a conceptual design of BWAVES, a time-of-flight broadband inverted-geometry neutron spectrometer for the Second Target Station at the Spallation Neutron Source operated by Oak Ridge National Laboratory. Being the first inverted geometry spectrometer where the energy of the detected neutrons can be chosen by a WAVES device mechanically, irrespective of the limitations imposed by the crystal analyzers or filters, BWAVES will feature a uniquely broad, continuous dynamic range of measurable energy transfers, spanning 4.5 decades. This will enable measurements of both vibrational and relaxational excitations within the same, continuous scattering spectra. Novel approaches that are necessary for the implementation of a WAVES device at the BWAVES spectrometer will result in a spectrometer with the design and characteristics much different from those displayed by the neutron spectrometers in existence today.

Low rotational barriers for the most dynamically active methyl groups in the proposed antiviral drugs for treatment of SARS-CoV-2, apilimod and tetrandrine

A recent screening study highlighted a molecular compound, apilimod, for its efficacy against the SARS-CoV-2 virus, while another compound, tetrandrine, demonstrated a remarkable synergy with the benchmark antiviral drug, remdesivir. Here, we find that because of significantly reduced potential energy barriers, which also give rise to pronounced quantum effects, the rotational dynamics of the most dynamically active methyl groups in apilimod and tetrandrine are much faster than those in remdesivir. Because dynamics of methyl groups are essential for biochemical activity, screening studies based on the computed potential energy profiles may help identify promising candidates within a given class of drugs.

Hydration-Induced Disorder Lowers the Energy Barriers for Methyl Rotation in Drug Molecules

The thermally activated dynamics of methyl groups are important for biochemical activity as they allow for a more efficient sampling of the energy landscape. Here, we compare methyl rotations in the dry and variously hydrated states of three primary drugs under consideration to treat the recent coronavirus disease (COVID-19), namely, hydroxychloroquine and its sulfate, dexamethasone and its sodium diphosphate, and remdesivir. We find that the main driving force behind the considerable reduction in the activation energy for methyl rotations in the hydrated state is the hydration-induced disorder in the methyl group local environments. Furthermore, the activation energy for methyl rotations in the hydration-induced disordered state is much lower than that in an isolated drug molecule, indicating that neither isolated molecules nor periodic crystalline structures can be used to analyze the potential landscape governing the side group dynamics in drug molecules. Instead, only the explicitly considered disordered structures can provide insight.

Effect of Hydration on the Molecular Dynamics of Hydroxychloroquine Sulfate

Chloroquine and its derivative hydroxychloroquine are primarily known as antimalaria drugs. Here, we investigate the influence of hydration water on the molecular dynamics in hydroxychloroquine sulfate, a commonly used solubilized drug form. When hydration, even at a low level, results in a disordered structure, as opposed to the highly ordered structure of dry hydroxychloroquine sulfate, the activation barriers for the rotation of methyl groups in the drug molecules become randomized and, on average, significantly reduced. The facilitated stochastic motions of the methyl groups may benefit the biomolecular activity due to the more efficient sampling of the energy landscape in the disordered hydration environment experienced by the drug molecules in vivo.

Water in protein hydration and ligand recognition

This review describes selected basics of water in biomolecular recognition. We focus on a qualitative understanding of the most important physical aspects, how these change in magnitude between bulk water and protein environment, and how the roles that water plays for proteins arise from them. These roles include mechanical support, thermal coupling, dielectric screening, mass and charge transport, and the competition with a ligand for the occupation of a binding site. The presence or absence of water has ramifications that range from the thermodynamic binding signature of a single ligand up to cellular survival. The large inhomogeneity in water density, polarity and mobility around a solute is hard to assess in experiment. This is a source of many difficulties in the solvation of protein models and computational studies that attempt to elucidate or predict ligand recognition. The influence of water in a protein binding site on the experimental enthalpic and entropic signature of ligand binding is still a point of much debate. The strong water‐water interaction in enthalpic terms is counteracted by a water molecule's high mobility in entropic terms. The complete arrest of a water molecule's mobility sets a limit on the entropic contribution of a water displacement process, while the solvent environment sets limits on ligand reactivity.

Anticancer drug impact on DNA – a study by neutron spectroscopy coupled with synchrotron-based FTIR and EXAFS

Complementary information on drug–DNA interplay has been achieved for Pt/Pd anticancer agents, by a combined QENS, SR-FTIR-ATR and EXAFS approach.

Structural effects of modified ribonucleotides and magnesium in transfer RNAs

Modified nucleotides are ubiquitous and important to tRNA structure and function. To understand their effect on tRNA conformation, we performed a series of molecular dynamics simulations on yeast tRNA^Phe and tRNA_init, Escherichia coli tRNA_init and HIV tRNA^Lys. Simulations were performed with the wild type modified nucleotides, using the recently developed CHARMM compatible force field parameter set for modified nucleotides (J. Comput. Chem.2016, 37, 896), or with the corresponding unmodified nucleotides, and in the presence or absence of Mg²⁺. Results showed a stabilizing effect associated with the presence of the modifications and Mg²⁺ for some important positions, such as modified guanosine in position 37 and dihydrouridines in 16/17 including both structural properties and base interactions. Some other modifications were also found to make subtle contributions to the structural properties of local domains. While we were not able to investigate the effect of adenosine 37 in tRNA_init and limitations were observed in the conformation of E. coli tRNA_init, the presence of the modified nucleotides and of Mg²⁺ better maintained the structural features and base interactions of the tRNA systems than in their absence indicating the utility of incorporating the modified nucleotides in simulations of tRNA and other RNAs.

Atomistic details of protein dynamics and the role of hydration water

BACKGROUND

The importance of protein dynamics for their biological activity is now well recognized. Different experimental and computational techniques have been employed to study protein dynamics, hierarchy of different processes and the coupling between protein and hydration water dynamics. Yet, understanding the atomistic details of protein dynamics and the role of hydration water remains rather limited.

SCOOP OF REVIEW

Based on overview of neutron scattering, molecular dynamic simulations, NMR and dielectric spectroscopy results we present a general picture of protein dynamics covering time scales from faster than ps to microseconds and the influence of hydration water on different relaxation processes.

MAJOR CONCLUSIONS

Internal protein dynamics spread over a wide time range from faster than picosecond to longer than microseconds. We suggest that the structural relaxation in hydrated proteins appears on the microsecond time scale, while faster processes present mostly motion of side groups and some domains. Hydration water plays a crucial role in protein dynamics on all time scales. It controls the coupled protein-hydration water relaxation on 10-100ps time scale. This process defines the friction for slower protein dynamics. Analysis suggests that changes in amount of hydration water affect not only general friction, but also influence significantly the protein's energy landscape.

GENERAL SIGNIFICANCE

The proposed atomistic picture of protein dynamics provides deeper understanding of various relaxation processes and their hierarchy, similarity and differences between various biological macromolecules, including proteins, DNA and RNA. This article is part of a Special Issue entitled "Science for Life" Guest Editor: Dr. Austen Angell, Dr. Salvatore Magazù and Dr. Federica Migliardo".

Protein dynamics: from rattling in a cage to structural relaxation

We present an overview of protein dynamics based mostly on results of neutron scattering, dielectric relaxation spectroscopy and molecular dynamics simulations. We identify several major classes of protein motions on the time scale from faster than picoseconds to several microseconds, and discuss the coupling of these processes to solvent dynamics. Our analysis suggests that the microsecond backbone relaxation process might be the main structural relaxation of the protein that defines its glass transition temperature, while faster processes present some localized secondary relaxations. Based on the overview, we formulate a general picture of protein dynamics and discuss the challenges in this field.

Painting biological low-frequency vibrational modes from small peptides to proteins

We use experiments and simulation to investigate the validity of different model systems used to study the low-frequency vibrations of proteins.

Dynamics and rigidity in an intrinsically disordered protein, β-casein

The emergence of intrinsically disordered proteins (IDPs) as a recognized structural class has forced the community to confront a new paradigm of structure, dynamics, and mechanical properties for proteins. We present novel data on the similarities and differences in the dynamics and nanomechanical properties of IDPs and other biomacromolecules on the picosecond time scale. An IDP, β-casein (CAS), has been studied in a calcium bound and unbound state using neutron and light scattering techniques. We show that CAS partially folds and stiffens upon calcium binding, but in the unfolded state, it is softer than folded proteins such as green fluorescence protein (GFP). We also see that some localized diffusive motions in CAS have a larger amplitude than in GFP at this time scale but are still smaller than those observed in tRNA. In spite of these differences, CAS dynamics are consistent with the classes of motions seen in folded protein on this time scale.

Coherent neutron scattering and collective dynamics in the protein, GFP

Collective dynamics are considered to be one of the major properties of soft materials, including biological macromolecules. We present coherent neutron scattering studies of the low-frequency vibrations, the so-called boson peak, in fully deuterated green fluorescent protein (GFP). Our analysis revealed unexpectedly low coherence of the atomic motions in GFP. This result implies a low amount of in-phase collective motion of the secondary structural units contributing to the boson peak vibrations and fast conformational fluctuations on the picosecond timescale. These observations are in contrast to earlier studies of polymers and glass-forming systems, and suggest that random or out-of-phase motions of the β-strands contribute greater than two-thirds of the intensity to the low-frequency vibrational spectra of GFP.

Dynamics of protein and its hydration water: neutron scattering studies on fully deuterated GFP

We present a detailed analysis of the picosecond-to-nanosecond motions of green fluorescent protein (GFP) and its hydration water using neutron scattering spectroscopy and hydrogen/deuterium contrast. The analysis reveals that hydration water suppresses protein motions at lower temperatures (<~ 200 K), and facilitates protein dynamics at high temperatures. Experimental data demonstrate that the hydration water is harmonic at temperatures <~ 180-190 K and is not affected by the proteins' methyl group rotations. The dynamics of the hydration water exhibits changes at ~ 180-190 K that we ascribe to the glass transition in the hydrated protein. Our results confirm significant differences in the dynamics of protein and its hydration water at high temperatures: on the picosecond-to-nanosecond timescale, the hydration water exhibits diffusive dynamics, while the protein motions are localized to <~3 Å. The diffusion of the GFP hydration water is similar to the behavior of hydration water previously observed for other proteins. Comparison with other globular proteins (e.g., lysozyme) reveals that on the timescale of 1 ns and at equivalent hydration level, GFP dynamics (mean-square displacements and quasielastic intensity) are of much smaller amplitude. Moreover, the suppression of the protein dynamics by the hydration water at low temperatures appears to be stronger in GFP than in other globular proteins. We ascribe this observation to the barrellike structure of GFP.