Log-periodic oscillations as real-time signatures of hierarchical dynamics in proteins.

The time-dependent relaxation of a dynamical system may exhibit a power-law behavior that is superimposed by log-periodic oscillations. D. Sornette [Phys. Rep. 297, 239 (1998)] showed that this behavior can be explained by a discrete scale invariance of the system, which is associated with discrete and equidistant timescales on a logarithmic scale. Examples include such diverse fields as financial crashes, random diffusion, and quantum topological materials. Recent time-resolved experiments and molecular dynamics simulations suggest that discrete scale invariance may also apply to hierarchical dynamics in proteins, where several fast local conformational changes are a prerequisite for a slow global transition to occur. Employing entropy-based timescale analysis and Markov state modeling to a simple one-dimensional hierarchical model and biomolecular simulation data, it is found that hierarchical systems quite generally give rise to logarithmically spaced discrete timescales. By introducing a one-dimensional reaction coordinate that collectively accounts for the hierarchically coupled degrees of freedom, the free energy landscape exhibits a characteristic staircase shape with two metastable end states, which causes the log-periodic time evolution of the system. The period of the log-oscillations reflects the effective roughness of the energy landscape and can, in simple cases, be interpreted in terms of the barriers of the staircase landscape.

Lifestyle and Dietary Habits Change Before and During Quarantine and Subsequent Weight Gain.

Background: Pandemics greatly affected human health due to changes in dietary habits and lifestyle. Methods: A multi-centric comparative cross-sectional study was conducted online using a structured questionnaire with 323 respondents from two medical colleges of Lahore, Pakistan. The collected data was analyzed including various variables of dietary habits, lifestyle changes and were compared by using paired sample t-test. Chi-square test was applied to measure the relation of weight gain with lifestyle and dietary habits changes. Result: Out of 323 participants, 64.5% of them observed an increase in their weights and 64.1% of respondents noticed an increase in caloric intake. An increase in chicken, milk, oil, sugar and legumes intake was increased significantly during the quarantine. The lifestyle changes were noted in 91.6% of participants and decreased physical activity, working duration and outdoor games were found significant. Out of 7.4% of students who were smokers, 6.8% of them quit smoking during the quarantine. Conclusion: A significant increase in weight is indicated during quarantine as compared to before quarantine probably due to the changes in dietary habits and lifestyle. If the pandemic prolongs, the problem of weight gain may worsen and lead to obesity and other health problems.

Metallofullerenes as Robust Single-Atom Catalysts for Adsorption and Dissociation of Hydrogen Molecules: A Density Functional Study.

Hydrogen is currently considered as the best alternative for traditional fuels due to its sustainable and ecofriendly nature. Additionally, hydrogen dissociation is a critical step in almost all hydrogenation reactions, which is crucial in industrial chemical production. A cost-effective and efficient catalyst with favorable activity for this step is highly desirable. Herein, transition-metal-doped fullerene (TM@C60) complexes are designed and investigated as single-atom catalysts for the hydrogen splitting process. Interaction energy analysis (Eint) is also carried out to demonstrate the stability of designed TM@C60 metallofullerenes, which reveals that all the designed complexes have higher thermodynamic stability. Furthermore, among all the studied metallofullerenes, the best catalytic efficiency for hydrogen dissociation is seen for the Sc@C60 catalyst Ea = 0.13 eV followed by the V@C60 catalyst Ea = 0.19 eV. The hydrogen activation and dissociation processes over TM@C60 metallofullerenes is further elaborated by analyzing charge transfer via the natural bond orbital and electron density difference analyses. Additionally, quantum theory of atoms in molecule analysis is carried out to investigate the nature of interatomic interactions between hydrogen molecules and TMs@C60 metallofullerenes. Overall, results of the current study declare that the Sc@C60 catalyst can act as a low cost, highly efficient, and noble metal-free single-atom catalyst to efficiently catalyze hydrogen dissociation reaction.

Elucidating the Self-cleavage Dynamics of Hairpin Ribozyme by Mode-decomposed Infrared Spectroscopy.

While catalytic reactions of biomolecular processes play an indispensable role in life, extracting the underlying molecular picture often remains challenging. Based on ab initio simulations of the self-cleavage reaction of hairpin ribozyme, mode-decomposed infrared spectra, and cosine similarity analysis to correlate the product with reactant IR spectra, we demonstrate a strategy to extract molecular details from characteristic spectral changes. Our results are in almost quantitative agreement with the experimental IR band library of nucleic acids and suggest that the spectral range of 800-1200 cm-1 is particularly valuable to monitor self-cleavage. Importantly, the cosine similarities also disclose that IR peaks subject to slight shifts due to self-cleavage might be unrelated, while strongly shifting resonances can correspond to the same structural dynamics. This framework of correlating complex IR spectra at the molecular level along biocatalytic reaction pathways is broadly applicable.

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 α3-helix and ligand binding. In an intriguing experiment, Hamm and co-workers employed photoswitching of the α3-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 α3-helix and the core of the protein, which move in a concerted manner. Representing a one-step process and involving direct α3-ligand contacts, this cooperative transition is considered as the elementary step in the propagation of conformational change.

Through bonds or contacts? Mapping protein vibrational energy transfer using non-canonical amino acids.

Vibrational energy transfer (VET) is essential for protein function. It is responsible for efficient energy dissipation in reaction sites, and has been linked to pathways of allosteric communication. While it is understood that VET occurs via backbone as well as via non-covalent contacts, little is known about the competition of these two transport channels, which determines the VET pathways. To tackle this problem, we equipped the ß-hairpin fold of a tryptophan zipper with pairs of non-canonical amino acids, one serving as a VET injector and one as a VET sensor in a femtosecond pump probe experiment. Accompanying extensive non-equilibrium molecular dynamics simulations combined with a master equation analysis unravel the VET pathways. Our joint experimental/computational endeavor reveals the efficiency of backbone vs. contact transport, showing that even if cutting short backbone stretches of only 3 to 4 amino acids in a protein, hydrogen bonds are the dominant VET pathway.

Real-time observation of ligand-induced allosteric transitions in a PDZ domain.

While allostery is of paramount importance for protein regulation, the underlying dynamical process of ligand (un)binding at one site, resulting time evolution of the protein structure, and change of the binding affinity at a remote site are not well understood. Here the ligand-induced conformational transition in a widely studied model system of allostery, the PDZ2 domain, is investigated by transient infrared spectroscopy accompanied by molecular dynamics simulations. To this end, an azobenzene-derived photoswitch is linked to a peptide ligand in a way that its binding affinity to the PDZ2 domain changes upon switching, thus initiating an allosteric transition in the PDZ2 domain protein. The subsequent response of the protein, covering four decades of time, ranging from ∼1 ns to ∼µs, can be rationalized by a remodeling of its rugged free-energy landscape, with very subtle shifts in the populations of a small number of structurally well-defined states. It is proposed that structurally and dynamically driven allostery, often discussed as limiting scenarios of allosteric communication, actually go hand-in-hand, allowing the protein to adapt its free-energy landscape to incoming signals.

Master equation model to predict energy transport pathways in proteins.

Recent time-resolved experiments and accompanying molecular dynamics simulations allow us to monitor the flow of vibrational energy in biomolecules. As a simple means to describe these experimental and simulated data, Buchenberg et al. [J. Phys. Chem. Lett. 7, 25 (2016)] suggested a master equation model that accounts for the energy transport from an initially excited residue to some target residue. The transfer rates of the model were obtained from two scaling rules, which account for the energy transport through the backbone and via tertiary contacts, respectively, and were parameterized using simulation data of a small α-helical protein at low temperatures. To extend the applicability of the model to general proteins at room temperature, here a new parameterization is presented, which is based on extensive nonequilibrium molecular dynamics simulations of a number of model systems. With typical transfer times of 0.5-1 ps between adjacent residues, backbone transport represents the fastest channel of energy flow. It is well described by a diffusive-type scaling rule, which requires only an overall backbone diffusion coefficient and interatom distances as input. Contact transport, e.g., via hydrogen bonds, is considerably slower (6-30 ps) at room temperature. A new scaling rule depending on the inverse square contact distance is suggested, which is shown to successfully describe the energy transport in the allosteric protein PDZ3. Since both scaling rules require only the structure of the considered system, the model provides a simple and general means to predict energy transport in proteins. To identify the pathways of energy transport, Monte Carlo Markov chain simulations are performed, which highlight the competition between backbone and contact transport channels.

COVID-19 and Its Psychological Impacts on Healthcare Staff - A Multi-Centric Comparative Cross-Sectional Study.

Background Since the first case of coronavirus disease-19 (COVID-19) in Pakistan was reported in February 2020, the medical and paramedical staff has been working on the frontlines to deal with this disease. They have been facing significant strain and stress due to the pandemic, affecting their social, mental, and personal life. The purpose of this study is to investigate the psychological effects of the COVID-19 pandemic, etiology, personal coping mechanisms, and the strategies that are being adopted to reduce stress by the healthcare workers (HCWs) working in COVID-19 dedicated wards (group 2) and compare it with staff working in other departments but not in COVID-19 wards amid this pandemic (group 1) in various hospitals of Lahore, Pakistan. Methods The comparative cross-sectional study was designed which included doctors, nurses, and allied health professionals from various hospitals of Lahore, Pakistan. A questionnaire was designed which consisted of five sections, and 51 questions. A Chi-square test was used to compare the responses between these two groups. Results The study questionnaire was submitted by 200 participants, 100 responses for each group (see the Appendix). In group 1, HCWs not working in COVID-19 dedicated floors were afraid of getting infected, transmitting the infection to their families and concerned about using personal protective equipment (PPE) improperly. They reported a lack of concentration and tense muscles. The coping mechanisms of this group were exercise, strict precautions at work, and social distancing measures. While HCWs serving in COVID-19 dedicated wards were concerned and afraid of putting their families at risk by working in the high-risk environment; the major stresses in this group were: lack of knowledge about proper strategies for treatment, they faced insecurity due to physical and verbal violence by caretakers of COVID-19 patients, and lack of concentration. The coping mechanism was the support of their families and taking strict precautions, with self-isolation if required, to avoid any disease transmission to their families. The proposed strategies to be implemented included teaching skills for self-rescue as well as the implementation of policies at the administrative level to reduce working hours and frequent shift rotation. Conclusion The COVID-19 outbreak posed a great deal of mental stress among HCWs working on the COVID-19 floor as well as those serving in other departments of the hospital. The HCWs from group 1 were most afraid of getting infected and putting family members at risk, experienced tense muscles and lack of concentration, coped their stress by exercise and being more vigilant, and suggested the strategies of teaching skills for self-rescue and better community awareness. While the staff from the second group were most afraid of being the source of infection and violence from the caretakers of patients, experienced tense muscles, used family support, and strict isolation measures as coping mechanisms and suggested the strategies of self-rescue and increase in wages of directly exposed healthcare workers to deal with such pandemics in future in a better way.

Energy Transport Pathways in Proteins: A Non-equilibrium Molecular Dynamics Simulation Study.

To facilitate the observation of biomolecular energy transport in real time and with single-residue resolution, recent experiments by Baumann et al. ( Angew. Chem. Int. Ed. 2019 , 58 , 2899 , DOI: 10.1002/anie.201812995 ) have used unnatural amino acids ß-(1-azulenyl)alanine (Azu) and azidohomoalanine (Aha) to site-specifically inject and probe vibrational energy in proteins. To aid the interpretation of such experiments, non-equilibrium molecular dynamics simulations of the anisotropic energy flow in proteins TrpZip2 and PDZ3 domains are presented. On this account, an efficient simulation protocol is established that accurately mimics the excitation and probing steps of Azu and Aha. The simulations quantitatively reproduce the experimentally found cooling times of the solvated proteins at room temperature and predict that the cooling slows by a factor 2 below the glass temperature of water. In PDZ3, vibrational energy is shown to travel from the initially excited peptide ligand via a complex network of inter-residue contacts and backbone transport to distal regions of the protein. The supposed connection of these energy transport pathways with pathways of allosteric communication is discussed.

Photocontrolling Protein-Peptide Interactions: From Minimal Perturbation to Complete Unbinding.

An azobenzene-derived photoswitch has been covalently cross-linked to two sites of the S-peptide in the RNase S complex in a manner that the α-helical content of the S-peptide reduces upon cis-to-trans isomerization of the photoswitch. Three complementary experimental techniques have been employed, isothermal titration calorimetry, circular dichroism spectroscopy and intrinsic tyrosine fluorescence quenching, to determine the binding affinity of the S-peptide to the S-protein in the two states of the photoswitch. Five mutants with the photoswitch attached to different sites of the S-peptide have been explored, with the goal to maximize the change in binding affinity upon photoswitching, and to identify the mechanisms that determine the binding affinity. With regard to the first goal, one mutant has been identified, which binds with reasonable affinity in the one state of the photoswitch, while specific binding is completely switched off in the other state. With regard to the second goal, accompanying molecular dynamics simulations combined with a quantitative structure activity relationship revealed that the α-helicity of the S-peptide in the binding pocket correlates surprisingly well with measured dissociation constants. Moreover, the simulations show that both configurations of all S-peptides exhibit quite well-defined structures, even in apparently disordered states.

Azidohomoalanine: A Minimally Invasive, Versatile, and Sensitive Infrared Label in Proteins To Study Ligand Binding.

The noncanonical amino acid azidohomoalanine (Aha) is known to be an environment-sensitive infrared probe for the site-specific investigation of protein structure and dynamics. Here, the capability of that label is explored to detect protein-ligand interactions by incorporating it in the vicinity of the binding groove of a PDZ2 domain. Circular dichroism and isothermal titration calorimetry measurements reveal that the perturbation of the protein system by mutation is negligible, with minimal influence on protein stability and binding affinity. Two-dimensional infrared spectra exhibit small (1-3 cm-1) but clearly measurable red shifts of the Aha vibrational frequency upon binding of two different peptide ligands, while accompanying molecular dynamics simulations suggest that these red shifts are induced by polar contacts with side chains of the peptide ligands. Hence, Aha is a versatile and minimally invasive vibrational label that is not only able to report on large structural changes during, e.g., protein folding, but also on very subtle changes of the electrostatic environment upon ligand binding.

Molecular dynamics simulation unveils the conformational flexibility of the interdomain linker in the bacterial transcriptional regulator GabR from Bacillus subtilis bound to pyridoxal 5'-phosphate.

GabR from Bacillus subtilis is a transcriptional regulator belonging to the MocR subfamily of the GntR regulators. The structure of the MocR regulators is characterized by the presence of two domains: i) a N-terminal domain, about 60 residue long, possessing the winged-Helix-Turn-Helix (wHTH) architecture with DNA recognition and binding capability; ii) a C-terminal domain (about 350 residue) folded as the pyridoxal 5'-phosphate (PLP) dependent aspartate aminotransferase (AAT) with dimerization and effector binding functions. The two domains are linked to each other by a peptide bridge. Although structural and functional characterization of MocRs is proceeding at a fast pace, virtually nothing is know about the molecular changes induced by the effector binding and on how these modifications influence the properties of the regulator. An extensive molecular dynamics simulation on the crystallographic structure of the homodimeric B. subtilis GabR has been undertaken with the aim to envisage the role and the importance of conformational flexibility in the action of GabR. Molecular dynamics has been calculated for the apo (without PLP) and holo (with PLP bound) forms of the GabR. A comparison between the molecular dynamics trajectories calculated for the two GabR forms suggested that one of the wHTH domain detaches from the AAT-like domain in the GabR PLP-bound form. The most evident conformational change in the holo PLP-bound form is represented by the rotation and the subsequent detachment from the subunit surface of one of the wHTH domains. The movement is mediated by a rearrangement of the linker connecting the AAT domain possibly triggered by the presence of the negative charge of the PLP cofactor. This is the second most significant conformational modification. The C-terminal section of the linker docks into the "active site" pocket and establish stabilizing contacts consisting of hydrogen-bonds, salt-bridges and hydrophobic interactions.

2D-IR Spectroscopy of an AHA Labeled Photoswitchable PDZ2 Domain.

We explore the capability of the non-natural amino acid azidohomoalanine (AHA) as an IR label to sense relatively small structural changes in proteins with the help of 2D IR difference spectroscopy. To that end, we AHA-labeled an allosteric protein (the PDZ2 domain from human tyrosine-phosphatase 1E) and furthermore covalently linked it to an azobenzene-derived photoswitch as to mimic its conformational transition upon ligand binding. To determine the strengths and limitations of the AHA label, in total six mutants have been investigated with the label at sites with varying properties. Only one mutant revealed a measurable 2D IR difference signal. In contrast to the commonly observed frequency shifts that report on the degree of solvation, in this case we observe an intensity change. To understand this spectral response, we performed classical MD simulations, evaluating local contacts of the AHA labels to water molecules and protein side chains and calculating the vibrational frequency on the basis of an electrostatic model. Although these simulations revealed in part significant and complex changes of the number of intraprotein and water contacts upon trans-cis photoisomerization, they could not provide a clear explanation of why this one label would stick out. Subsequent quantum-chemistry calculations suggest that the response is the result of an electronic interaction involving charge transfer of the azido group with sulfonate groups from the photoswitch. To the best of our knowledge, such an effect has not been described before.