Dynamics of Hydrogen Bonds between Water and Intrinsically Disordered and Structured Regions of Proteins.

Recent studies indicate more restricted dynamics of water around intrinsically disordered proteins (IDPs) than structured proteins. We examine here the dynamics of hydrogen bonds between water molecules and two proteins, small ubiquitin-related modifier-1 (SUMO-1) and ubiquitin-conjugating enzyme E2I (UBC9), which we compare around intrinsically disordered regions (IDRs) and structured regions of these proteins. It has been recognized since some time that excluded volume effects, which influence access of water molecules to hydrogen-bonding sites, and the strength of hydrogen bonds between water and protein affect hydrogen bond lifetimes. While we find those two properties to mediate lifetimes of hydrogen bonds between water and protein residues in this study, we also find that the lifetimes are affected by the concentration of charged groups on other nearby residues. These factors are more important in determining the hydrogen bond lifetimes than whether a residue hydrogen bonding with water belongs to an IDR or to a structured region.

Locating dynamic contributions to allostery via determining rates of vibrational energy transfer.

Determining rates of energy transfer across non-covalent contacts for different states of a protein can provide information about dynamic and associated entropy changes during transitions between states. We investigate the relationship between rates of energy transfer across polar and nonpolar contacts and contact dynamics for the ß2-adrenergic receptor, a rhodopsin-like G-protein coupled receptor, in an antagonist-bound inactive state and agonist-bound active state. From structures sampled during molecular dynamics (MD) simulations, we find the active state to have, on average, a lower packing density, corresponding to generally more flexibility and greater entropy than the inactive state. Energy exchange networks (EENs) are computed for the inactive and active states from the results of the MD simulations. From the EENs, changes in the rates of energy transfer across polar and nonpolar contacts are found for contacts that remain largely intact during activation. Change in dynamics of the contact, and entropy associated with the dynamics, can be estimated from the change in rates of energy transfer across the contacts. Measurement of change in the rates of energy transfer before and after the transition between states thereby provides information about dynamic contributions to activation and allostery.

Energy Transport in Class B GPCRs: Role of Protein-Water Dynamics and Activation.

We compute energy exchange networks (EENs) through glucagon-like peptide-1 receptor (GLP-1R), a class B G-protein-coupled receptor (GPCR), in inactive and two active states, one activated by a peptide ligand and the other by a small molecule agonist, from results of molecular dynamics simulations. The reorganized network upon activation contains contributions from structural as well as from dynamic changes and corresponding entropic contributions to the free energy of activation, which are estimated in terms of the change in rates of energy transfer across non-covalent contacts. The role of water in the EENs and in activation of GLP-1R is also investigated. The dynamics of water in contact with the central polar network of the transmembrane region is found to be significantly slower for both activated states compared to the inactive state. This result is consistent with the contribution of water molecules to activation of GLP-1R previously suggested and resembles water dynamics in parts of the transmembrane region found in earlier studies of rhodopsin-like GPCRs.

A Statistical Journey through the Topological Determinants of the ß2 Adrenergic Receptor Dynamics.

Activation of G-protein-coupled receptors (GPCRs) is mediated by molecular switches throughout the transmembrane region of the receptor. In this work, we continued along the path of a previous computational study wherein energy transport in the ß2 Adrenergic Receptor (ß2-AR) was examined and allosteric switches were identified in the molecular structure through the reorganization of energy transport networks during activation. In this work, we further investigated the allosteric properties of ß2-AR, using Protein Contact Networks (PCNs). In this paper, we report an extensive statistical analysis of the topological and structural properties of ß2-AR along its molecular dynamics trajectory to identify the activation pattern of this molecular system. The results show a distinct character to the activation that both helps to understand the allosteric switching previously identified and confirms the relevance of the network formalism to uncover relevant functional features of protein molecules.

Biophysical Insight into the SARS-CoV2 Spike-ACE2 Interaction and Its Modulation by Hepcidin through a Multifaceted Computational Approach.

At the center of the SARS-CoV2 infection, the spike protein and its interaction with the human receptor ACE2 play a central role in the molecular machinery of SARS-CoV2 infection of human cells. Vaccine therapies are a valuable barrier to the worst effects of the virus and to its diffusion, but the need of purposed drugs is emerging as a core target of the fight against COVID19. In this respect, the repurposing of drugs has already led to discovery of drugs thought to reduce the effects of the cytokine storm, but still a drug targeting the spike protein, in the infection stage, is missing. In this work, we present a multifaceted computational approach strongly grounded on a biophysical modeling of biological systems, so to disclose the interaction of the SARS-CoV2 spike protein with ACE2 with a special focus to an allosteric regulation of the spike-ACE2 interaction. Our approach includes the following methodologies: Protein Contact Networks and Network Clustering, Targeted Molecular Dynamics, Elastic Network Modeling, Perturbation Response Scanning, and a computational analysis of energy flow and SEPAS as a protein-softness and monomer-based affinity predictor. We applied this approach to free (closed and open) states of spike protein and spike-ACE2 complexes. Eventually, we analyzed the interactions of free and bound forms of spike with hepcidin (HPC), the major hormone in iron regulation, recently addressed as a central player in the COVID19 pathogenesis, with a special emphasis to the most severe outcomes. Our results demonstrate that, compared with closed and open states, the spike protein in the ACE2-bound state shows higher allosteric potential. The correspondence between hinge sites and the Allosteric Modulation Region (AMR) in the S-ACE complex suggests a molecular basis for hepcidin involvement in COVID19 pathogenesis. We verify the importance of AMR in different states of spike and then study its interactions with HPC and the consequence of the HPC-AMR interaction on spike dynamics and its affinity for ACE2. We propose two complementary mechanisms for HPC effects on spike of SARS-CoV-2; (a) HPC acts as a competitive inhibitor when spike is in a preinfection state (open and with no ACE2), (b) the HPC-AMR interaction pushes the spike structure into the safer closed state. These findings need clear molecular in vivo verification beside clinical observations.

Enhanced Mobility during Diels-Alder Reaction: Results of Molecular Simulations.

Recent measurements indicate enhanced mobility of solvent molecules during Diels-Alder (DA) and other common chemical reactions. We present results of molecular dynamics simulations of the last stages of the DA cycloaddition reaction, from the transition state configuration to product, of furfurylamine and maleimide in acetonitrile at reactant concentrations studied experimentally. We find enhanced mobility of solvent and reactant molecules up to at least a nanometer from the DA product over hundreds of picoseconds. Local heating is ruled out as a factor in the enhanced mobility observed in the simulations, which is instead found to be due to solvent relaxation following the formation of the DA product.

The origin and impact of bound water around intrinsically disordered proteins.

Proteins and water couple dynamically over a wide range of time scales. Motivated by their central role in protein function, protein-water dynamics and thermodynamics have been extensively studied for structured proteins, where correspondence to structural features has been made. However, properties controlling intrinsically disordered protein (IDP)-water dynamics are not yet known. We report results of megahertz-to-terahertz dielectric spectroscopy and molecular dynamics simulations of a group of IDPs with varying charge content along with structured proteins of similar size. Hydration water around IDPs is found to exhibit more heterogeneous rotational and translational dynamics compared with water around structured proteins of similar size, yielding on average more restricted dynamics around individual residues of IDPs, charged or neutral, compared with structured proteins. The on-average slower water dynamics is found to arise from excess tightly bound water in the first hydration layer, which is related to greater exposure to charged groups. The more tightly bound water to IDPs correlates with the smaller hydration shell found experimentally, and affects entropy associated with protein-water interactions, the contribution of which we estimate based on the dielectric measurements and simulations. Water-IDP dynamic coupling at terahertz frequencies is characterized by the dielectric measurements and simulations.

Electric Fields Influence Intramolecular Vibrational Energy Relaxation and Line Widths.

Intramolecular vibrational energy relaxation (IVR) is fundamentally important to chemical dynamics. We show that externally applied electric fields affect IVR and vibrational line widths by changing the anharmonic couplings and frequency detunings between modes. We demonstrate this effect in benzonitrile for which prior experimental results show a decrease in vibrational line width as a function of applied electric field. We identify three major channels for IVR that depend on electric field. In the dominant channel, the electric field affects the frequency detuning, while in the other two channels, variation of anharmonic couplings as a function of field is the underlying mechanism. Consistent with experimental results, we show that the combination of all channels gives rise to reduced line widths with increasing electric field in benzonitrile. Our results are relevant for controlling IVR with external or internal fields and for gaining a more complete interpretation of line widths of vibrational Stark probes.

Activation-Induced Reorganization of Energy Transport Networks in the ß2 Adrenergic Receptor.

We compute energy exchange networks (EENs) through the ß2 adrenergic receptor (ß2AR), a G-protein coupled receptor (GPCR), in inactive and active states, based on the results of molecular dynamics simulations of this membrane bound protein. We introduce a new definition for the reorganization of EENs upon activation that depends on the relative change in rates of energy transfer across noncovalent contacts throughout the protein. On the basis of the reorganized network that we obtain for ß2AR upon activation, we identify a branched pathway between the agonist binding site and the cytoplasmic region, where a G-protein binds to the receptor when activated. The pathway includes all of the motifs containing molecular switches previously identified as contributing to the allosteric transition of ß2AR upon agonist binding. EENs and their reorganization upon activation are compared with structure-based contact networks computed for the inactive and active states of ß2AR.

Change in vibrational entropy with change in protein volume estimated with mode Grüneisen parameters.

For a small adjustment in average volume, due to a change in state of a protein or other macromolecule at constant temperature, the change in vibrational entropy is related to the mode Grüneisen parameters, which relate shifts in frequency to a small volume change. We report here values of mode Grüneisen parameters computed for two hydrated proteins, cytochrome c and myoglobin, which exhibit trends with mode frequency resembling those of glassy systems. We use the mode Grüneisen parameters to relate volumetric thermal expansion to previously computed values of the isothermal compressibility for several proteins. We also estimate changes in vibrational entropy resulting from the change in volume upon ligand bonding of myoglobin and the homodimeric hemoglobin from Scapharca inaequivalvis (HbI). We compare estimates of the change in entropy upon ligation obtained in terms of mode Grüneisen parameters with the results of normal mode analysis for myoglobin and earlier molecular dynamics simulations of HbI. The results illustrate how small changes in average volume can yield changes in entropy that contribute to ligand binding and allostery.

Locating and Navigating Energy Transport Networks in Proteins.

We review computational methods to locate energy transport networks in proteins that are based on the calculation of local energy diffusion in nanoscale systems. As an illustrative example, we discuss energy transport networks computed for the homodimeric hemoglobin from Scapharca inaequivalvis, where channels for facile energy transport, which include the cluster of water molecules at the interface of the globules, have been found to lie along pathways that experiments reveal are important in allosteric processes. We also review recent work on master equation simulations to model energy transport dynamics, including efforts to relate rate constants in the master equation to protein structural dynamics. Results for apomyoglobin involving relations between fluctuations in the length of hydrogen bonds and the energy flux between them are presented.

Energy Transfer across Nonpolar and Polar Contacts in Proteins: Role of Contact Fluctuations.

Molecular dynamics simulations of the villin headpiece subdomain HP36 have been carried out to examine relations between rates of vibrational energy transfer across non-covalently bonded contacts and equilibrium structural fluctuations, with focus on van der Waals contacts. Rates of energy transfer across van der Waals contacts vary inversely with the variance of the contact length, with the same constant of proportionality for all nonpolar contacts of HP36. A similar relation is observed for hydrogen bonds, but the proportionality depends on contact pairs, with hydrogen bonds stabilizing the α-helices all exhibiting the same constant of proportionality, one that is distinct from those computed for other polar contacts. Rates of energy transfer across van der Waals contacts are found to be up to 2 orders of magnitude smaller than rates of energy transfer across polar contacts.

Water-mediated biomolecular dynamics and allostery.

Dynamic coupling with water contributes to regulating the functional dynamics of a biomolecule. We discuss protein-water dynamics, with emphasis on water that is partially confined, and the role of protein-confined water dynamics in allosteric regulation. These properties are illustrated with two systems, a homodimeric hemoglobin from Scapharca inaequivalvis (HbI) and an A2A adenosine receptor (A2AAR). For HbI, water-protein interactions, long known to contribute to the thermodynamics of cooperativity, are seen to influence the dynamics of the protein not only around the protein-water interface but also into the core of each globule, where dynamic and entropic changes upon ligand binding are coupled to protein-water contact dynamics. Similarly, hydration waters trapped deep inside the core region of A2AAR enable the formation of an allosteric network made of water-mediated inter-residue contacts. Extending from the ligand binding pocket to the G-protein binding site, this allosteric network plays key roles in regulating the activity of the receptor.

Recent developments in the computational study of protein structural and vibrational energy dynamics.

Recent developments in the computational study of energy transport in proteins are reviewed, including advances in both methodology and applications. The concept of energy exchange network (EEN) is discussed, and a recent calculation of EENs for the allosteric protein FixL is reviewed, which illustrates how residues and protein regions involved in the allosteric transition can be identified. Recent work has examined relations between EENs and protein dynamics as well as structure. We review some of the computational studies carried out on several proteins that explore connections between energy conductivity across polar contacts in proteins and between proteins and water and equilibrium dynamics of the contacts, and we discuss some of the recent experimental work that addresses this topic.

Variation of Energy Transfer Rates across Protein-Water Contacts with Equilibrium Structural Fluctuations of a Homodimeric Hemoglobin.

Molecular dynamics simulations of the homodimeric hemoglobin from Scapharca inaequivalvis (HbI) have been carried out to examine relations between rates of vibrational energy transfer across nonbonded contacts and equilibrium structural fluctuations, with emphasis on protein-water contacts. The scaling of rates of energy transfer with equilibrium fluctuations of the contact length is found to hold up well for contacts between residues and hemes at the interface and the cluster of 17 interface water molecules in the unliganded state of HbI, as well as for the liganded state, for which the cluster contains on average 11 water molecules. In both states, the rate of energy transfer is also found to satisfy a diffusion relation. Within each globule, the scaling for polar contacts is similar to that found in an earlier analysis of myoglobin. Entropy associated with dynamics of polar contacts within each globule and with contacts between the hemes and water cluster is found to increase upon ligation. Energy exchange networks (EENs) for liganded and unliganded states obtained from the simulations are also presented and discussed. Energy transport networks through which nonbonded contacts transport energy in HbI, referred to as nonbonded networks (NBNs), are determined from the EENs and compared for the two states.

Tuning Molecular Vibrational Energy Flow within an Aromatic Scaffold via Anharmonic Coupling.

From guiding chemical reactivity in synthesis or protein folding to the design of energy diodes, intramolecular vibrational energy redistribution harnesses the power to influence the underlying fundamental principles of chemistry. To evaluate the ability to steer these processes, the mechanism and time scales of intramolecular vibrational energy redistribution through aromatic molecular scaffolds have been assessed by utilizing two-dimensional infrared (2D IR) spectroscopy. 2D IR cross peaks reveal energy relaxation through an aromatic scaffold from the azido- to the cyano-vibrational reporters in para-azidobenzonitrile (PAB) and para-(azidomethyl)benzonitrile (PAMB) prior to energy relaxation into the solvent. The rates of energy transfer are modulated by Fermi resonances, which are apparent by the coupling cross peaks identified within the 2D IR spectrum. Theoretical vibrational mode analysis allowed the determination of the origins of the energy flow, the transfer pathway, and a direct comparison of the associated transfer rates, which were in good agreement with the experimental results. Large variations in energy-transfer rates, approximately 1.9 ps for PAB and 23 ps for PAMB, illustrate the importance of strong anharmonic coupling, i.e., Fermi resonance, on the transfer pathways. In particular, vibrational energy rectification is altered by Fermi resonances of the cyano- and azido-modes allowing control of the propensity for energy flow.

Energy Transport across Interfaces in Biomolecular Systems.

Energy transport during chemical reactions or following photoexcitation in systems of biological molecules is mediated by numerous interfaces that separate chemical groups and molecules. Describing and predicting energy transport has been complicated by the inhomogeneous environment through which it occurs, and general rules are still lacking. We discuss recent work on identification of networks for vibrational energy transport in biomolecules and their environment, with focus on the nature of energy transfer across interfaces. Energy transport is influenced both by structure of the biomolecular system as well as by equilibrium fluctuations of nonbonded contacts between chemical groups, biomolecules, and water along the network. We also discuss recent theoretical and computational work on the related topic of thermal transport through molecular interfaces, with focus on systems important in biology as well as relevant experimental studies.

Scaling of Rates of Vibrational Energy Transfer in Proteins with Equilibrium Dynamics and Entropy.

Theoretical arguments and results of molecular dynamics (MD) simulations of myoglobin at 300 K are presented to relate rates of vibrational energy transfer across nonbonded contacts interacting via short-range potentials to dynamics of the contact. Both theory and the results of the simulations support a scaling relation between the energy transfer rate and the inverse of the variance in the distance between hydrogen-bonded contacts. The results of the MD simulations do not support such a relation for longer-range charged contacts. Instead, the energy transfer rate is found to scale as a power law in the distance between charged groups. The scaling between rates of vibrational energy transfer across nonbonded contacts interacting via short-range potentials and conformational dynamics suggests a relation between vibrational energy transfer rates and entropy associated with the dynamics of interacting residues. The use of time-resolved vibrational spectroscopy to determine change in conformational entropy with change in protein functional state is discussed, and an expression quantifying the connection is provided.