Collective dynamics underlying allosteric transitions in hemoglobin

Identification of histidine residues that affect the T/R-state conformations of human hemoglobin using constant pH molecular dynamics simulations

Human hemoglobin (Hb) is a tetrameric protein consisting of two α and two β subunits that can adopt a low-affinity T- and high-affinity R-state conformations. Under physiological pH conditions, histidine (His) residues are the main sites for proton binding or release, and their protonation states can affect the T/R-state conformation of Hb. However, it remains unclear which His residues can effectively affect the Hb conformation. Herein, the impact of the 38 His residues of Hb on its T/R-state conformations was evaluated using constant-pH molecular dynamics (CpHMD) simulations at physiological pH while focusing on the His protonation states. Overall, the protonation states of some His residues were found to be correlated with the Hb conformation state. These residues were mainly located in the proximity of the heme (α87 and β92), and at the α1β2 and α2β1 interfaces (α89 and β97). This correlation may be partly explained by how easily hydrogen bonds can be formed, which depends on the protonation states of the His residues. Taken together, these CpHMD-based findings provide new insights into the identification of titratable His residues α87, α89, β92, and β97 that can affect Hb conformational switching under physiological pH conditions.

Wordom update 2: A user-friendly program for the analysis of molecular structures and conformational ensembles

We present the second update of Wordom, a user-friendly and efficient program for manipulation and analysis of conformational ensembles from molecular simulations. The actual update expands some of the existing modules and adds 21 new modules to the update 1 published in 2011. The new adds can be divided into three sets that: 1) analyze atomic fluctuations and structural communication; 2) explore ion-channel conformational dynamics and ionic translocation; and 3) compute geometrical indices of structural deformation. Set 1 serves to compute correlations of motions, find geometrically stable domains, identify a dynamically invariant core, find changes in domain-domain separation and mutual orientation, perform wavelet analysis of large-scale simulations, process the output of principal component analysis of atomic fluctuations, perform functional mode analysis, infer regions of mechanical rigidity, analyze overall fluctuations, and perform the perturbation response scanning. Set 2 includes modules specific for ion channels, which serve to monitor the pore radius as well as water or ion fluxes, and measure functional collective motions like receptor twisting or tilting angles. Finally, set 3 includes tools to monitor structural deformations by computing angles, perimeter, area, volume, β-sheet curvature, radial distribution function, and center of mass. The ring perception module is also included, helpful to monitor supramolecular self-assemblies. This update places Wordom among the most suitable, complete, user-friendly, and efficient software for the analysis of biomolecular simulations. The source code of Wordom and the relative documentation are available under the GNU general public license at http://wordom.sf.net.

On the Uses of PCA to Characterise Molecular Dynamics Simulations of Biological Macromolecules: Basics and Tips for an Effective Use

Principal Component Analysis (PCA) is a procedure widely used to examine data collected from molecular dynamics simulations of biological macromolecules. It allows for greatly reducing the dimensionality of their configurational space, facilitating further qualitative and quantitative analysis. Its simplicity and relatively low computational cost explain its extended use. However, a judicious implementation of PCA requires the knowledge of its theoretical grounds as well as its weaknesses and capabilities. In this article, we review these issues and discuss several strategies developed over the last years to mitigate the main PCA flaws and enhance the reproducibility of its results.

Nonequilibrium Modeling of the Elementary Step in PDZ3 Allosteric Communication

While allostery is of paramount importance for protein signaling and regulation, the underlying dynamical process of allosteric communication is not well understood. The PDZ3 domain represents a prime example of an allosteric single-domain protein, as it features a well-established long-range coupling between the C-terminal α₃-helix and ligand binding. In an intriguing experiment, Hamm and co-workers employed photoswitching of the α₃-helix to initiate a conformational change of PDZ3 that propagates from the C-terminus to the bound ligand within 200 ns. Performing extensive nonequilibrium molecular dynamics simulations, the modeling of the experiment reproduces the measured time scales and reveals a detailed picture of the allosteric communication in PDZ3. In particular, a correlation analysis identifies a network of contacts connecting the α₃-helix and the core of the protein, which move in a concerted manner. Representing a one-step process and involving direct α₃-ligand contacts, this cooperative transition is considered as the elementary step in the propagation of conformational change.

Atomic-Level View of the Functional Transition in Vertebrate Hemoglobins: The Case of Antarctic Fish Hbs

Tetrameric hemoglobins (Hbs) are prototypal systems for studies aimed at unveiling basic structure-function relationships as well as investigating the molecular/structural basis of adaptation of living organisms to extreme conditions. However, a chronological analysis of decade-long studies conducted on Hbs is illuminating on the difficulties associated with the attempts of gaining functional insights from static structures. Here, we applied molecular dynamics (MD) simulations to explore the functional transition from the T to the R state of the hemoglobin of the Antarctic fish Trematomus bernacchii (HbTb). Our study clearly demonstrates the ability of the MD technique to accurately describe the transition of HbTb from the T to R-like states, as shown by a number of global and local structural indicators. A comparative analysis of the structural states that HbTb assumes in the simulations with those detected in previous MD analyses conducted on HbA (human Hb) highlights interesting analogies (similarity of the transition pathway) and differences (distinct population of intermediate states). In particular, the ability of HbTb to significantly populate intermediate states along the functional pathway explains the observed propensity of this protein to assume these structures in the crystalline state. It also explains some functional data reported on the protein that indicate the occurrence of other functional states in addition to the canonical R and T ones. These findings are in line with the emerging idea that the classical two-state view underlying tetrameric Hb functionality is probably an oversimplification and that other structural states play important roles in these proteins. The ability of MD simulations to accurately describe the functional pathway in tetrameric Hbs suggests that this approach may be effectively applied to unravel the molecular and structural basis of Hbs exhibiting peculiar functional properties as a consequence of the environmental adaptation of the host organism.

Cooperative Protein Allosteric Transition Mediated by a Fluctuating Transmission Network

Allosteric communication between distant protein sites represents a key mechanism of biomolecular regulation and signal transduction. Compared to other processes such as protein folding, however, the dynamical evolution of allosteric transitions is still not well understood. As an example of allosteric coupling between distant protein regions, we consider the global open-closed motion of the two domains of T4 lysozyme, which is triggered by local motion in the hinge region. Combining extensive molecular dynamics simulations with a correlation analysis of interresidue contacts, we identify a network of interresidue distances that move in a concerted manner. The cooperative process originates from a cogwheel-like motion of the hydrophobic core in the hinge region, which constitutes an evolutionary conserved and flexible transmission network. Through rigid contacts and the protein backbone, the small local changes of the hydrophobic core are passed on to the distant terminal domains and lead to the emergence of a rare global conformational transition. As in an Ising-type model, the cooperativity of the allosteric transition can be explained via the interaction of local fluctuations.

Quaternary Structure Transitions of Human Hemoglobin: An Atomic-Level View of the Functional Intermediate States

Human hemoglobin (HbA) is one of the prototypal systems used to investigate structure-function relationships in proteins. Indeed, HbA has been used to develop the basic concepts of protein allostery, although the atomic-level mechanism underlying the HbA functionality is still highly debated. This is due to the fact that most of the three-dimensional structural information collected over the decades refers to the endpoints of HbA functional transition with little data available for the intermediate states. Here, we report molecular dynamics (MD) simulations by focusing on the relevance of the intermediate states of the protein functional transition unraveled by the crystallographic studies carried out on vertebrate Hbs. Fully atomistic simulations of the HbA T-state indicate that the protein undergoes a spontaneous transition toward the R-state. The inspection of the trajectory structures indicates that the protein significantly populates the intermediate HL-(C) state previously unraveled by crystallography. In the structural transition, it also assumes the intermediate states crystallographically detected in Antarctic fish Hbs. This finding suggests that HbA and Antarctic fish Hbs, in addition to the endpoints of the transitions, also share a similar deoxygenation pathway despite a distace of hundreds of millions of years in the evolution scale. Finally, using the essential dynamic sampling methodology, we gained some insights into the reverse R to T transition that is not spontaneously observed in classic MD simulations.

Exposing the distinctive modular behavior of β-strands and α-helices in folded proteins

Although folded proteins are commonly depicted as simplistic combinations of β-strands and α-helices, the actual properties and functions of these secondary-structure elements in their native contexts are just partly understood. The principal reason is that the behavior of individual β- and α-elements is obscured by the global folding cooperativity. In this study, we have circumvented this problem by designing frustrated variants of the mixed α/β-protein S6, which allow the structural behavior of individual β-strands and α-helices to be targeted selectively by stopped-flow kinetics, X-ray crystallography, and solution-state NMR. Essentially, our approach is based on provoking intramolecular "domain swap." The results show that the α- and β-elements have quite different characteristics: The swaps of β-strands proceed via global unfolding, whereas the α-helices are free to swap locally in the native basin. Moreover, the α-helices tend to hybridize and to promote protein association by gliding over to neighboring molecules. This difference in structural behavior follows directly from hydrogen-bonding restrictions and suggests that the protein secondary structure defines not only tertiary geometry, but also maintains control in function and structural evolution. Finally, our alternative approach to protein folding and native-state dynamics presents a generally applicable strategy for in silico design of protein models that are computationally testable in the microsecond-millisecond regime.

Use of Molecular Dynamics Simulations in Structure-Based Drug Discovery

BACKGROUND

Traditional drug discovery is a lengthy process which involves a huge amount of resources. Modern-day drug discovers various multidisciplinary approaches amongst which, computational ligand and structure-based drug designing methods contribute significantly. Structure-based drug designing techniques require the knowledge of structural information of drug target and drug-target complexes. Proper understanding of drug-target binding requires the flexibility of both ligand and receptor to be incorporated. Molecular docking refers to the static picture of the drug-target complex(es). Molecular dynamics, on the other hand, introduces flexibility to understand the drug binding process.

OBJECTIVE

The aim of the present study is to provide a systematic review on the usage of molecular dynamics simulations to aid the process of structure-based drug design.

METHOD

This review discussed findings from various research articles and review papers on the use of molecular dynamics in drug discovery. All efforts highlight the practical grounds for which molecular dynamics simulations are used in drug designing program. In summary, various aspects of the use of molecular dynamics simulations that underline the basis of studying drug-target complexes were thoroughly explained.

RESULTS

This review is the result of reviewing more than a hundred papers. It summarizes various problems that use molecular dynamics simulations.

CONCLUSION

The findings of this review highlight how molecular dynamics simulations have been successfully implemented to study the structure-function details of specific drug-target complexes. It also identifies the key areas such as stability of drug-target complexes, ligand binding kinetics and identification of allosteric sites which have been elucidated using molecular dynamics simulations.

Graph-theoretical analysis for energy landscape reveals the organization of state transitions in the resting-state human cerebral cortex

The resting-state brain is often considered a nonlinear dynamic system transitioning among multiple coexisting stable states. Despite the increasing number of studies on the multistability of the brain system, the processes of state transitions have rarely been systematically explored. Thus, we investigated the state transition processes of the human cerebral cortex system at rest by introducing a graph-theoretical analysis of the state transition network. The energy landscape analysis of brain state occurrences, estimated using the pairwise maximum entropy model for resting-state fMRI data, identified multiple local minima, some of which mediate multi-step transitions toward the global minimum. The state transition among local minima is clustered into two groups according to state transition rates and most inter-group state transitions were mediated by a hub transition state. The distance to the hub transition state determined the path length of the inter-group transition. The cortical system appeared to have redundancy in inter-group transitions when the hub transition state was removed. Such a hub-like organization of transition processes disappeared when the connectivity of the cortical system was altered from the resting-state configuration. In the state transition, the default mode network acts as a transition hub, while coactivation of the prefrontal cortex and default mode network is captured as the global minimum. In summary, the resting-state cerebral cortex has a well-organized architecture of state transitions among stable states, when evaluated by a graph-theoretical analysis of the nonlinear state transition network of the brain.

Response to comment on 'Valid molecular dynamics simulations of human hemoglobin require a surprisingly large box size'

We recently reported that molecular dynamics simulations for hemoglobin require a surprisingly large box size to stabilize the T(0) state relative to R(0), as observed in experiments (El Hage et al., 2018). Gapsys and de Groot have commented on this work but do not provide convincing evidence that the conclusions of El Hage et al., 2018 are incorrect. Here we respond to these concerns, argue that our original conclusions remain valid, and raise our own concerns about some of the results reported in the comment by Gapsys and de Groot that require clarification.

Comment on 'Valid molecular dynamics simulations of human hemoglobin require a surprisingly large box size'

A recent molecular dynamics investigation into the stability of hemoglobin concluded that the unliganded protein is only stable in the T state when a solvent box is used in the simulations that is ten times larger than what is usually employed (El Hage et al., 2018). Here, we express three main concerns about that study. In addition, we find that with an order of magnitude more statistics, the reported box size dependence is not reproducible. Overall, no significant effects on the kinetics or thermodynamics of conformational transitions were observed.

Identifying coupled clusters of allostery participants through chemical shift perturbations

Allosteric couplings underlie many cellular signaling processes and provide an exciting avenue for development of new diagnostics and therapeutics. A general method for identifying important residues in allosteric mechanisms would be very useful, but remains elusive due to the complexity of long-range phenomena. Here, we introduce an NMR method to identify residues involved in allosteric coupling between two ligand-binding sites in a protein, which we call chemical shift detection of allostery participants (CAP). Networks of functional groups responding to each ligand are defined through correlated NMR perturbations. In this process, we also identify allostery participants, groups that respond to both binding events and likely play a role in the coupling between the binding sites. Such residues exhibit multiple functional states with distinct NMR chemical shifts, depending on binding status at both binding sites. Such a strategy was applied to the prototypical ion channel KcsA. We had previously shown that the potassium affinity at the extracellular selectivity filter is strongly dependent on proton binding at the intracellular pH sensor. Here, we analyzed proton and potassium binding networks and identified groups that depend on both proton and potassium binding (allostery participants). These groups are viewed as candidates for transmitting information between functional units. The vital role of one such identified amino acid was validated through site-specific mutagenesis, electrophysiology functional studies, and NMR-detected thermodynamic analysis of allosteric coupling. This strategy for identifying allostery participants is likely to have applications for many other systems.

Dissimilar flexibility of α and β subunits of human adult hemoglobin influences the protein dynamics and its alteration induced by allosteric effectors

The general question by what mechanism an "effector" molecule and the hemes of hemoglobin interact over widely separated intramolecular distances to change the oxygen affinity has been extensively investigated, and still has remained of central interest. In the present work we were interested in clarifying the general role of the protein matrix and its dynamics in the regulation of human adult hemoglobin (HbA). We used a spectroscopy approach that yields the compressibility (κ) of the protein matrix around the hemes of the subunits in HbA and studied how the binding of heterotropic allosteric effectors modify this parameter. κ is directly related to the variance of volume fluctuation, therefore it characterizes the molecular dynamics of the protein structure. For the experiments the heme groups either in the α or in the β subunits of HbA were replaced by fluorescent Zn-protoporphyrinIX, and series of fluorescence line narrowed spectra were measured at varied pressures. The evaluation of the spectra yielded the compressibility that showed significant dynamic asymmetry between the subunits: κ of the α subunit was 0.17±0.05/GPa, while for the β subunit it was much higher, 0.36±0.07/GPa. The heterotropic effectors, chloride ions, inositol hexaphosphate and bezafibrate did not cause significant changes in κ of the α subunits, while in the β subunits the effectors lead to a significant reduction down to 0.15±0.04/GPa. We relate our results to structural data, to results of recent functional studies and to those of molecular dynamics simulations, and find good agreements. The observed asymmetry in the flexibility suggests a distinct role of the subunits in the regulation of Hb that results in the observed changes of the oxygen binding capability.

The Dynamic Behavior of the P2X4 Ion Channel in the Closed Conformation

We present the results of a detailed molecular dynamics study of the closed form of the P2X₄ receptor. The fluctuations observed in the simulations were compared with the changes that occur in the transition from the closed to the open structure. To get further insight on the opening mechanism, the actual displacements were decomposed into interchain motions and intrachain deformations. This analysis revealed that the iris-like expansion of the transmembrane helices mainly results from interchain motions that already take place in the closed conformation. However, these movements cannot reach the amplitude required for the opening of the channel because they are impeded by interactions occurring around the ATP binding pocket. This suggests that the union of ATP produces distortions in the chains that eliminate the restrictions on the interchain displacements, leading to the opening of the pore.

Consistent Principal Component Modes from Molecular Dynamics Simulations of Proteins

Principal component analysis is a technique widely used for studying the movements of proteins using data collected from molecular dynamics simulations. In spite of its extensive use, the technique has a serious drawback: equivalent simulations do not afford the same PC-modes. In this article, we show that concatenating equivalent trajectories and calculating the PC-modes from the concatenated one significantly enhances the reproducibility of the results. Moreover, the consistency of the modes can be systematically improved by adding more individual trajectories to the concatenated one.

Comparative Study of Elastic Network Model and Protein Contact Network for Protein Complexes: The Hemoglobin Case

The overall topology and interfacial interactions play key roles in understanding structural and functional principles of protein complexes. Elastic Network Model (ENM) and Protein Contact Network (PCN) are two widely used methods for high throughput investigation of structures and interactions within protein complexes. In this work, the comparative analysis of ENM and PCN relative to hemoglobin (Hb) was taken as case study. We examine four types of structural and dynamical paradigms, namely, conformational change between different states of Hbs, modular analysis, allosteric mechanisms studies, and interface characterization of an Hb. The comparative study shows that ENM has an advantage in studying dynamical properties and protein-protein interfaces, while PCN is better for describing protein structures quantitatively both from local and from global levels. We suggest that the integration of ENM and PCN would give a potential but powerful tool in structural systems biology.

Induced Fit in Protein Multimerization: The HFBI Case

Hydrophobins, produced by filamentous fungi, are small amphipathic proteins whose biological functions rely on their unique surface-activity properties. Understanding the mechanistic details of the multimerization process is of primary importance to clarify the interfacial activity of hydrophobins. We used free energy calculations to study the role of a flexible β-hairpin in the multimerization process in hydrophobin II from Trichoderma reesei (HFBI). We characterized how the displacement of this β-hairpin controls the stability of the monomers/dimers/tetramers in solution. The regulation of the oligomerization equilibrium of HFBI will necessarily affect its interfacial properties, fundamental for its biological function and for technological applications. Moreover, we propose possible routes for the multimerization process of HFBI in solution. This is the first case where a mechanism by which a flexible loop flanking a rigid patch controls the protein-protein binding equilibrium, already known for proteins with charged binding hot-spots, is described within a hydrophobic patch.

Fungi proliferate by creating a complex hyphal network growing within a wet environment. However, for most fungi to colonize new territories, they must produce spores carried by aerial hyphae and spread them into the air. Aerial structures need to overcome the surface tension of the surrounding water in order to grow into the air. This process requires hydrophobins, a remarkable class of self-associating fungal proteins which lower the surface tension at the air/water interface by creating a thin amphipathic layer. In solution they form multimers in equilibrium with the interfacial layer. Due to their unique surface-activity properties, hydrophobins have been used for a variety of biotechnical applications. We used enhanced sampling molecular dynamics simulations methods to study the multimerization process in solution of a hydrophobin from Trichoderma reesei (HFBI). We clarified the fundamental role of a small flexible region within the HFBI monomer involved in the formation of multimers. A flexible loop flanking a rigid interaction patch is able to fine-tune the interaction energy. This mechanism, already known for charged binding patches, is described here for hydrophobic hot-spots. This result is remarkably important in order to clarify the mechanism of arranging at the interface and enhancing hydrophobin-based technological applications.

Mechanisms for allosteric activation of protease DegS by ligand binding and oligomerization as revealed from molecular dynamics simulations

A local perturbation of a protein may lead to functional changes at some distal site, a phenomenon denoted as allostery. Here, we study the allosteric control of a protease using molecular dynamics simulations. The system considered is the bacterial protein DegS which includes a protease domain activated on ligand binding to an adjacent PDZ domain. Starting from crystallographic structures of DegS homo-trimers, we perform simulations of the ligand-free and -bound state of DegS at equilibrium. Considering a single protomer only, the trimeric state was mimicked by applying restraints on the residues in contact with other protomers in the DegS trimer. In addition, the bound state was also simulated without any restraints to mimic the monomer. Our results suggest that not only ligand release but also disassembly of a DegS trimer inhibits proteolytic activity. Considering various observables for structural changes, we infer allosteric pathways from the interface with other protomers to the active site. Moreover, we study how ligand release leads to (i) catalytically relevant changes involving residues 199-201 and (ii) a transition from a stretched to a bent conformation for residues 217-219 (which prohibits proper substrate binding). Finally, based on ligand-induced C_α shifts we identify residues in contact with other protomers in the DegS trimer that likely transduce the perturbation from ligand release from a given protomer to adjacent protomers. These residues likely play a key role in the experimentally known effect of ligand release from a protomer on the proteolytic activity of the other protomers. Proteins 2016; 84:1690-1705. © 2016 Wiley Periodicals, Inc.

Near infrared light induces post-translational modifications of human red blood cell proteins

There is a growing body of evidence that near infrared (NIR) light exerts beneficial effects on cells. Its usefulness in the treatment of cancer, acute brain injuries, strokes and neurodegenerative disorders has been proposed. The mechanism of the NIR action is probably of photochemical nature, however it is not fully understood. Here, using a relatively simple biological model, human red blood cells (RBCs), and a polychromatic non-polarized light source, we investigate the impact of NIR radiation on the oxygen carrier, hemoglobin (Hb), and anion exchanger (AE1, Band 3). The exposure of intact RBCs to NIR light causes quaternary transitions in Hb, dehydration of proteins and decreases the amount of physiologically inactive methemoglobin, as detected by Raman spectroscopy. These effects are accompanied by a lowering of the intracellular pH (pHi) and changes in the cell membrane topography, as documented by atomic force microscopy (AFM). All those changes are in line with our previous studies where alterations of the membrane fluidity and membrane potential were attributed to NIR action on RBCs. The rate of the above listed changes depends strictly on the dose of NIR light that the cells receive, nonetheless it should not be considered as a thermal effect.

Long-Range Conformational Response of a PDZ Domain to Ligand Binding and Release: A Molecular Dynamics Study

The binding of a ligand to a protein may induce long-range structural or dynamical changes in the biomacromolecule even at sites physically well separated from the binding pocket. A system for which such behavior has been widely discussed is the PDZ2 domain of human tyrosine phosphatase 1E. Here, we present results from equilibrium trajectories of the PDZ2 domain in the free and ligand-bound state, as well as nonequilibrium simulations of the relaxation of PDZ2 after removal of its peptide ligand. The study reveals changes in inter-residue contacts, backbone dihedral angles, and C(α) positions upon ligand release. Our findings show a long-range conformational response of the PDZ2 domain to ligand release in the form of a collective shift of the secondary structure elements α2, β2, β3, α1-β4, and the C terminal loop relative to the rest of the protein away from the N-terminus, and a shift of the loops β2-β3 and β1-β2 in the opposite direction. The shifts lead to conformational changes in the backbone, especially in the β2-β3 loop but also in the β5-α2 and the α2-β6 loop, and are accompanied by changes of inter-residue contacts mainly within the β2-β3 loop as well as between the α2 helix and other segments. The residues showing substantial changes of inter-residue contacts, backbone conformations, or C(α) positions are considered "key residues" for the long-range conformational response of PDZ2. By comparing these residues with various sets of residues highlighted by previous studies of PDZ2, we investigate the statistical correlation of the various approaches. Interestingly, we find a considerable correlation of our findings with several works considering structural changes but no significant correlations with approaches considering energy flow or networks based on inter-residue energies.

Long-range conformational transition of a photoswitchable allosteric protein: molecular dynamics simulation study

A local perturbation of a protein may lead to functional changes at some distal site. An example is the PDZ2 domain of human tyrosine phosphatase 1E, which shows an allosteric transition upon binding to a peptide ligand. Recently Buchli et al. presented a time-resolved study of this transition by covalently linking an azobenzene photoswitch across the binding groove and using a femtosecond laser pulse that triggers the cis-trans photoisomerization of azobenzene. To aid the interpretation of these experiments, in this work seven microsecond runs of all-atom molecular dynamics simulations each for the wild-type PDZ2 in the ligand-bound and -free state, as well as the photoswitchable protein (PDZ2S) in the cis and trans states of the photoswitch, in explicit water were conducted. First the theoretical model is validated by recalculating the available NMR data from the simulations. By comparing the results for PDZ2 and PDZ2S, it is analyzed to what extent the photoswitch indeed mimics the free-bound transition. A detailed description of the conformational rearrangement following the cis-trans photoisomerization of PDZ2S reveals a series of photoinduced structural changes that propagate from the anchor residues of the photoswitch via intermediate secondary structure segments to the C-terminus of PDZ2S. The changes of the conformational distribution of the C-terminal region is considered as the distal response of the isolated allosteric protein.