Dynamic multiple-scattering treatment of X-ray absorption: Parameterization of a new molecular dynamics force field for myoglobin

Understanding the role of R266K mutation in cystathionine β-synthase (CBS) enzyme: an in silico study

Human cystathionine β-synthase (hCBS) is a Heme-containing, unique pyridoxal 5'-phosphate (PLP) dependent enzyme. CBS catalyzes the bio-chemical condensation reactions in the transsulfuration pathway. The role of Heme in the catalytic activities of the hCBS enzyme is still unknown, even though various experimental studies indicated its participation in the bi-directional electronic communication with the PLP center. The hypothesis is, Heme acts as an electron density reservoir for the catalytic reaction center rather than a redox electron source. In this work, we have investigated In Silico dynamical aspects of the bi-directional communications by performing classical molecular dynamics (MD) simulations upon developing the necessary force field parameters for the cysteine and histidine bound hexa-coordinated Heme. The comparative aspects, of electron density overlap across the communicating pathways, were investigated adopting the Density Functional Theory (DFT) in conjunction with the hybrid exchange-correlation functional for the CBS^WT (wild-type) and CBS^R266K (mutated) enzymes. The molecular dynamics simulations and subsequent explorations of the electronic structures confirm the reported observations. It also provides an in-depth mechanistic understanding of how the non-covalent hydrogen bonding interactions with Cys52 control such long-distance communication. Our study also provides a convincing answer to the reduced enzymatic activities in the R266K mutated hCBS compared to the wild-type enzymes. The difference in hydrogen-bonding patterns and salt-bridge interactions play the pivotal roles in such long distant bi-directional communications.Communicated by Ramaswamy H. Sarma.

Stationary and Time-Dependent Carbon Monoxide Stretching Mode Features in Carboxy Myoglobin: A Theoretical-Computational Reappraisal

The stationary and time-dependent infrared spectrum (IR) of the CO stretching mode (ν_CO) in carboxymyoglobin (MbCO), a longstanding problem of biophysical chemistry, has been modeled through a theoretical-computational method specifically designed for simulating quantum observables in complex atomic-molecular systems and based on a combined application of long time scale molecular dynamics simulations and quantum-chemical calculations. This study is basically focused on two aspects: (i) the origin of the stationary IR substates (termed as A₀, A₁, and A₃) and (ii) the modeling and the interpretation of the ν_CO energy relaxation. The results, strengthened by a more than satisfactory agreement with the experimental data, concisely indicate that (i) the conformational His64-FeCO relevant substates, i.e., characterized by the formation-disruption of the H-bond between the above moieties, are the main responsible of the presence of two distinct and well separated (A₀ and A₁/A₃) spectroscopic regions; (ii) the characteristic bimodal shape of the A₁/A₃ spectral region, according to our model, is the result of the fluctuation of the electric field pattern as provided by the protein-solvent framework perturbing the bound His64-CO-Heme complex; and (iii) the electric field pattern, in conjunction with the relatively high density of MbCO vibrational states, is also the main determinant of the ν_CO energy relaxation, characterizing its kinetic efficiency.

Understanding conformational diversity of heat shock protein 90 (HSP90) and binding features of inhibitors to HSP90 via molecular dynamics simulations

Heat shock protein 90 (HSP90) is a promising target for treatment of cancer, and inhibitor bindings can generate efficient suppression on tumor in multiple ways. In this work, 140-ns molecular dynamics simulations were performed on six systems. Principal component analysis was subsequently carried out to explore the conformational diversity of HSP90. The results suggest that inhibitor bindings induce large conformational changes of HSP90, which tends to enlarge the volume of the binding pocket to facilitate the entrance of inhibitors. Hierarchical clustering analyses, the calculation of the energy contribution of each atom, and the analyses of hydrogen-bonding interactions were performed. The results indicate that 20 residues in group A of the hierarchical tree are responsible for major contributions, and van der Waals interactions as well as hydrogen-bonding interactions between important residues in HSP90 and key regions of inhibitors are the main force for promoting inhibitor bindings. We expect that this work can provide useful theoretical information for development of efficient inhibitors targeting HSP90.

Ligand pathways in neuroglobin revealed by low-temperature photodissociation and docking experiments

A combined biophysical approach was applied to map gas-docking sites within murine neuroglobin (Ngb), revealing snapshots of events that might govern activity and dynamics in this unique hexacoordinate globin, which is most likely to be involved in gas-sensing in the central nervous system and for which a precise mechanism of action remains to be elucidated. The application of UV-visible microspectroscopy in crystallo, solution X-ray absorption near-edge spectroscopy and X-ray diffraction experiments at 15-40 K provided the structural characterization of an Ngb photolytic intermediate by cryo-trapping and allowed direct observation of the relocation of carbon monoxide within the distal heme pocket after photodissociation. Moreover, X-ray diffraction at 100 K under a high pressure of dioxygen, a physiological ligand of Ngb, unravelled the existence of a storage site for O2 in Ngb which coincides with Xe-III, a previously described docking site for xenon or krypton. Notably, no other secondary sites were observed under our experimental conditions.