Coupled flexibility change in cytochrome P450cam substrate binding determined by neutron scattering, NMR, and molecular dynamics simulation

Joint neutron/molecular dynamics vibrational spectroscopy reveals softening of HIV-1 protease upon binding of a tight inhibitor

Biomacromolecules are inherently dynamic, and their dynamics are interwoven into function. The fast collective vibrational dynamics in proteins occurs in the low picosecond timescale corresponding to frequencies of ∼5-50 cm^-1. This sub-to-low THz frequency regime covers the low-amplitude collective breathing motions of a whole protein and vibrations of the constituent secondary structure elements, such as α-helices, β-sheets and loops. We have used inelastic neutron scattering experiments in combination with molecular dynamics simulations to demonstrate the vibrational dynamics softening of HIV-1 protease, a target of HIV/AIDS antivirals, upon binding of a tight clinical inhibitor darunavir. Changes in the vibrational density of states of matching structural elements in the two monomers of the homodimeric protein are not identical, indicating asymmetric effects of darunavir on the vibrational dynamics. Three of the 11 major secondary structure elements contribute over 40% to the overall changes in the vibrational density of states upon darunavir binding. Molecular dynamics simulations informed by experiments allowed us to estimate that the altered vibrational dynamics of the protease would contribute -3.6 kcal mol^-1 at 300 K, or 25%, to the free energy of darunavir binding. As HIV-1 protease drug resistance remains a concern, our results open a new avenue to help establish a direct quantitative link between protein vibrational dynamics and drug resistance.

Reconciling conformational heterogeneity and substrate recognition in cytochrome P450

Cytochrome P450, the ubiquitous metalloenzyme involved in detoxification of foreign components, has remained one of the most popular systems for substrate-recognition process. However, despite being known for its high substrate specificity, the mechanistic basis of substrate-binding by archetypal system cytochrome P450cam has remained at odds with the contrasting reports of multiple diverse crystallographic structures of its substrate-free form. Here, we address this issue by elucidating the probability of mutual dynamical transition to the other crystallographic pose of cytochrome P450cam and vice versa via unbiased all-atom computer simulation. A robust Markov state model, constructed using adaptively sampled 84-μs-long molecular dynamics simulation trajectories, maps the broad and heterogenous P450cam conformational landscape into five key substates. In particular, the Markov state model identifies an intermediate-assisted dynamic equilibrium between a pair of conformations of P450cam, in which the substrate-recognition sites remain "closed" and "open," respectively. However, the estimate of a significantly higher stationary population of closed conformation, coupled with faster rate of open → closed transition than its reverse process, dictates that the net conformational equilibrium would be swayed in favor of "closed" conformation. Together, the investigation quantitatively infers that although a potential substrate of cytochrome P450cam would, in principle, explore a diverse array of conformations of substrate-free protein, it would mostly encounter a "closed" or solvent-occluded conformation and hence would follow an induced-fit-based recognition process. Overall, the work reconciles multiple precedent crystallographic, spectroscopic investigations and establishes how a statistical elucidation of conformational heterogeneity in protein would provide crucial insights in the mechanism of potential substrate-recognition process.

Dynamics underlying hydroxylation selectivity of cytochrome P450cam

Structural heterogeneity and the dynamics of the complexes of enzymes with substrates can determine the selectivity of catalysis; however, fully characterizing how remains challenging as heterogeneity and dynamics can vary at the spatial level of an amino acid residue and involve rapid timescales. We demonstrate the nascent approach of site-specific two-dimensional infrared (IR) spectroscopy to investigate the archetypical cytochrome P450, P450cam, to better delineate the mechanism of the lower regioselectivity of hydroxylation of the substrate norcamphor in comparison to the native substrate camphor. Specific locations are targeted throughout the enzyme by selectively introducing cyano groups that have frequencies in a spectrally isolated region of the protein IR spectrum as local vibrational probes. Linear and two-dimensional IR spectroscopy were applied to measure the heterogeneity and dynamics at each probe and investigate how they differentiate camphor and norcamphor recognition. The IR data indicate that the norcamphor complex does not fully induce a large-scale conformational change to a closed state of the enzyme adopted in the camphor complex. Additionally, a probe directed at the bound substrate experiences rapidly interconverting states in the norcamphor complex that explain the hydroxylation product distribution. Altogether, the study reveals large- and small-scale structural heterogeneity and dynamics that could contribute to selectivity of a cytochrome P450 and illustrates the approach of site-selective IR spectroscopy to elucidate protein dynamics.

Substrate Binding Stiffens Aspartate Aminotransferase by Altering the Enzyme Picosecond Vibrational Dynamics

Protein dynamics on various time scales from femtoseconds to milliseconds impacts biological function by driving proteins to conformations conducive to ligand binding and creating functional states in enzyme catalysis. Neutron vibrational spectroscopy carried out by measuring inelastic neutron scattering from protein molecules in combination with molecular simulations has the unique ability of detecting and visualizing changes in the picosecond protein vibrational dynamics due to ligand binding. Here we present neutron vibrational spectra of a homodimeric pyridoxal 5'-phosphate-dependent enzyme, aspartate aminotransferase, obtained from the open internal aldimine and closed external aldimine conformational states. We observe that in the external aldimine state the protein structure stiffens relative to the internal aldimine state, indicating rigidified vibrational dynamics on the picosecond time scale in the low-frequency regime of 5-50 cm^-1. Our molecular dynamics simulations indicate substantial changes in the picosecond dynamics of the enzyme secondary structure elements upon substrate binding, with the largest contributions from just two helices and the β-sheet.

Neutron scattering in the biological sciences: progress and prospects

The scattering of neutrons can be used to provide information on the structure and dynamics of biological systems on multiple length and time scales. Pursuant to a National Science Foundation-funded workshop in February 2018, recent developments in this field are reviewed here, as well as future prospects that can be expected given recent advances in sources, instrumentation and computational power and methods. Crystallography, solution scattering, dynamics, membranes, labeling and imaging are examined. For the extraction of maximum information, the incorporation of judicious specific deuterium labeling, the integration of several types of experiment, and interpretation using high-performance computer simulation models are often found to be particularly powerful.

Entropic contribution to enhanced thermal stability in the thermostable P450 CYP119

The enhanced thermostability of thermophilic proteins with respect to their mesophilic counterparts is often attributed to the enthalpy effect, arising from strong interactions between protein residues. Intuitively, these strong interresidue interactions will rigidify the biomolecules. However, the present work utilizing neutron scattering and solution NMR spectroscopy measurements demonstrates a contrary example that the thermophilic cytochrome P450, CYP119, is much more flexible than its mesophilic counterpart, CYP101A1, something which is not apparent just from structural comparison of the two proteins. A mechanism to explain this apparent contradiction is that higher flexibility in the folded state of CYP119 increases its conformational entropy and thereby reduces the entropy gain during denaturation, which will increase the free energy needed for unfolding and thus stabilize the protein. This scenario is supported by thermodynamic data on the temperature dependence of unfolding free energy, which shows a significant entropic contribution to the thermostability of CYP119 and lends an added dimension to enhanced stability, previously attributed only to presence of aromatic stacking interactions and salt bridge networks. Our experimental data also support the notion that highly thermophilic P450s such as CYP119 may use a mechanism that partitions flexibility differently from mesophilic P450s between ligand binding and thermal stability.

Acceleration of biomolecular kinetics in Gaussian accelerated molecular dynamics

Recent studies demonstrated that Gaussian accelerated molecular dynamics (GaMD) is a robust computational technique, which provides simultaneous unconstrained enhanced sampling and free energy calculations of biomolecules. However, the exact acceleration of biomolecular dynamics or speedup of kinetic rates in GaMD simulations and, more broadly, in enhanced sampling methods, remains a challenging task to be determined. Here, the GaMD acceleration is examined using alanine dipeptide in explicit solvent as a biomolecular model system. Relative to long conventional molecular dynamics simulation, GaMD simulations exhibited ∼36-67 times speedup for sampling of the backbone dihedral transitions. The acceleration depended on level of the GaMD boost potential. Furthermore, Kramers' rate theory was applied to estimate GaMD acceleration using simulation-derived diffusion coefficients, curvatures and barriers of free energy profiles. In most cases, the calculations also showed significant speedup of dihedral transitions in GaMD, although the GaMD acceleration factors tended to be underestimated by ∼3-96 fold. Because greater boost potential can be applied in GaMD simulations of systems with increased sizes, which potentially leads to higher acceleration, it is subject to future studies on accelerating the dynamics and recovering kinetic rates of larger biomolecules such as proteins and protein-protein/nucleic acid complexes.

Dynamic Neutron Scattering by Biological Systems

Dynamic neutron scattering directly probes motions in biological systems on femtosecond to microsecond timescales. When combined with molecular dynamics simulation and normal mode analysis, detailed descriptions of the forms and frequencies of motions can be derived. We examine vibrations in proteins, the temperature dependence of protein motions, and concepts describing the rich variety of motions detectable using neutrons in biological systems at physiological temperatures. New techniques for deriving information on collective motions using coherent scattering are also reviewed.

Changes in dynamics of α-chymotrypsin due to covalent inhibitors investigated by elastic incoherent neutron scattering

The dynamics of chymotrypsin increases when bound to two different covalent inhibitors. These effects were analyzed by univariate and multivariate methods.

Determination of functional collective motions in a protein at atomic resolution using coherent neutron scattering

Protein function often depends on global, collective internal motions. However, the simultaneous quantitative experimental determination of the forms, amplitudes, and time scales of these motions has remained elusive. We demonstrate that a complete description of these large-scale dynamic modes can be obtained using coherent neutron-scattering experiments on perdeuterated samples. With this approach, a microscopic relationship between the structure, dynamics, and function in a protein, cytochrome P450cam, is established. The approach developed here should be of general applicability to protein systems.

Unconstrained Enhanced Sampling for Free Energy Calculations of Biomolecules: A Review

Free energy calculations are central to understanding the structure, dynamics and function of biomolecules. Yet insufficient sampling of biomolecular configurations is often regarded as one of the main sources of error. Many enhanced sampling techniques have been developed to address this issue. Notably, enhanced sampling methods based on biasing collective variables (CVs), including the widely used umbrella sampling, adaptive biasing force and metadynamics, have been discussed in a recent excellent review (Abrams and Bussi, Entropy, 2014). Here, we aim to review enhanced sampling methods that do not require predefined system-dependent CVs for biomolecular simulations and as such do not suffer from the hidden energy barrier problem as encountered in the CV-biasing methods. These methods include, but are not limited to, replica exchange/parallel tempering, self-guided molecular/Langevin dynamics, essential energy space random walk and accelerated molecular dynamics. While it is overwhelming to describe all details of each method, we provide a summary of the methods along with the applications and offer our perspectives. We conclude with challenges and prospects of the unconstrained enhanced sampling methods for accurate biomolecular free energy calculations.

Quasielastic neutron scattering in biology: Theory and applications

Neutrons scatter quasielastically from stochastic, diffusive processes, such as overdamped vibrations, localized diffusion and transitions between energy minima. In biological systems, such as proteins and membranes, these relaxation processes are of considerable physical interest. We review here recent methodological advances and applications of quasielastic neutron scattering (QENS) in biology, concentrating on the role of molecular dynamics simulation in generating data with which neutron profiles can be unambiguously interpreted. We examine the use of massively-parallel computers in calculating scattering functions, and the application of Markov state modeling. The decomposition of MD-derived neutron dynamic susceptibilities is described, and the use of this in combination with NMR spectroscopy. We discuss dynamics at very long times, including approximations to the infinite time mean-square displacement and nonequilibrium aspects of single-protein dynamics. Finally, we examine how neutron scattering and MD can be combined to provide information on lipid nanodomains. This article is part of a Special Issue entitled "Science for Life" Guest Editor: Dr. Austen Angell, Dr. Salvatore Magazù and Dr. Federica Migliardo.

Discovery of a regioselectivity switch in nitrating P450s guided by molecular dynamics simulations and Markov models

The dynamic motions of protein structural elements, particularly flexible loops, are intimately linked with diverse aspects of enzyme catalysis. Engineering of these loop regions can alter protein stability, substrate binding, and even dramatically impact enzyme function. When these flexible regions are structurally unresolvable, computational reconstruction in combination with large-scale molecular dynamics simulations can be used to guide the engineering strategy. Here, we present a collaborative approach consisting of both experiment and computation that led to the discovery of a single mutation in the F/G loop of the nitrating cytochrome P450 TxtE that simultaneously controls loop dynamics and completely shifts the enzyme's regioselectivity from the C4 to the C5 position of L-tryptophan. Furthermore, we find that this loop mutation is naturally present in a subset of homologous nitrating P450s and confirm that these uncharacterized enzymes exclusively produce 5-nitro-L-tryptophan, a previously unknown biosynthetic intermediate.

General trends of dihedral conformational transitions in a globular protein

Dihedral conformational transitions are analyzed systematically in a model globular protein, cytochrome P450cam, to examine their structural and chemical dependences through combined conventional molecular dynamics (cMD), accelerated molecular dynamics (aMD) and adaptive biasing force (ABF) simulations. The aMD simulations are performed at two acceleration levels, using dihedral and dual boost, respectively. In comparison with cMD, aMD samples protein dihedral transitions approximately two times faster on average using dihedral boost, and ∼ 3.5 times faster using dual boost. In the protein backbone, significantly higher dihedral transition rates are observed in the bend, coil, and turn flexible regions, followed by the β bridge and β sheet, and then the helices. Moreover, protein side chains of greater length exhibit higher transition rates on average in the aMD-enhanced sampling. Side chains of the same length (particularly Nχ = 2) exhibit decreasing transition rates with residues when going from hydrophobic to polar, then charged and aromatic chemical types. The reduction of dihedral transition rates is found to be correlated with increasing energy barriers as identified through ABF free energy calculations. These general trends of dihedral conformational transitions provide important insights into the hierarchical dynamics and complex free energy landscapes of functional proteins.

Comparison of intrinsic dynamics of cytochrome p450 proteins using normal mode analysis

Cytochrome P450 enzymes are hemeproteins that catalyze the monooxygenation of a wide-range of structurally diverse substrates of endogenous and exogenous origin. These heme monooxygenases receive electrons from NADH/NADPH via electron transfer proteins. The cytochrome P450 enzymes, which constitute a diverse superfamily of more than 8,700 proteins, share a common tertiary fold but < 25% sequence identity. Based on their electron transfer protein partner, cytochrome P450 proteins are classified into six broad classes. Traditional methods of pro are based on the canonical paradigm that attributes proteins' function to their three-dimensional structure, which is determined by their primary structure that is the amino acid sequence. It is increasingly recognized that protein dynamics play an important role in molecular recognition and catalytic activity. As the mobility of a protein is an intrinsic property that is encrypted in its primary structure, we examined if different classes of cytochrome P450 enzymes display any unique patterns of intrinsic mobility. Normal mode analysis was performed to characterize the intrinsic dynamics of five classes of cytochrome P450 proteins. The present study revealed that cytochrome P450 enzymes share a strong dynamic similarity (root mean squared inner product > 55% and Bhattacharyya coefficient > 80%), despite the low sequence identity (< 25%) and sequence similarity (< 50%) across the cytochrome P450 superfamily. Noticeable differences in Cα atom fluctuations of structural elements responsible for substrate binding were noticed. These differences in residue fluctuations might be crucial for substrate selectivity in these enzymes.

Delicate conformational balance of the redox enzyme cytochrome P450cam

The energy landscapes of proteins are highly complex and can be influenced by changes in physical and chemical conditions under which the protein is studied. The redox enzyme cytochrome P450cam undergoes a multistep catalytic cycle wherein two electrons are transferred to the heme group and the enzyme visits several conformational states. Using paramagnetic NMR spectroscopy with a lanthanoid tag, we show that the enzyme bound to its redox partner, putidaredoxin, is in a closed state at ambient temperature in solution. This result contrasts with recent crystal structures of the complex, which suggest that the enzyme opens up when bound to its partner. The closed state supports a model of catalysis in which the substrate is locked in the active site pocket and the enzyme acts as an insulator for the reactive intermediates of the reaction.

Insights into molecular interactions between CaM and its inhibitors from molecular dynamics simulations and experimental data

In order to contribute to the structural basis for rational design of calmodulin (CaM) inhibitors, we analyzed the interaction of CaM with 14 classic antagonists and two compounds that do not affect CaM, using docking and molecular dynamics (MD) simulations, and the data were compared to available experimental data. The Ca(2+)-CaM-Ligands complexes were simulated 20 ns, with CaM starting in the "open" and "closed" conformations. The analysis of the MD simulations provided insight into the conformational changes undergone by CaM during its interaction with these ligands. These simulations were used to predict the binding free energies (ΔG) from contributions ΔH and ΔS, giving useful information about CaM ligand binding thermodynamics. The ΔG predicted for the CaM's inhibitors correlated well with available experimental data as the r(2) obtained was 0.76 and 0.82 for the group of xanthones. Additionally, valuable information is presented here: I) CaM has two preferred ligand binding sites in the open conformation known as site 1 and 4, II) CaM can bind ligands of diverse structural nature, III) the flexibility of CaM is reduced by the union of its ligands, leading to a reduction in the Ca(2+)-CaM entropy, IV) enthalpy dominates the molecular recognition process in the system Ca(2+)-CaM-Ligand, and V) the ligands making more extensive contact with the protein have higher affinity for Ca(2+)-CaM. Despite their limitations, docking and MD simulations in combination with experimental data continue to be excellent tools for research in pharmacology, toward a rational design of new drugs.

Dynamics of CYP51: implications for function and inhibitor design

Sterol 14α-demethylase (cytochrome P450 family 51 (CYP51)) is an essential enzyme occurring in all biological kingdoms. In eukaryotes, it is located in the membrane of the endoplasmic reticulum. Selective inhibitors of trypanosomal CYP51s that do not affect the human CYP51 have been discovered in vitro and found to cure acute and chronic mouse Chagas disease without severe side effects in vivo. Crystal structures indicate that CYP51 may be more rigid than most CYPs, and it has been proposed that this property may facilitate antiparasitic drug design. Therefore, to investigate the dynamics of trypanosomal CYP51, we built a model of membrane-bound Trypanosoma brucei CYP51 and then performed molecular dynamics simulations of T. brucei CYP51 in membrane-bound and soluble forms. We compared the dynamics of T. brucei CYP51 with those of human CYP51, CYP2C9, and CYP2E1. In the simulations, the CYP51s display low mobility in the buried active site although overall mobility is similar in all the CYPs studied. The simulations suggest that in CYP51, pathway 2f serves as the major ligand access tunnel, and both pathways 2f (leading to membrane) and S (leading to solvent) can serve as ligand egress tunnels. Compared with the other CYPs, the residues at the entrance of the ligand access tunnels in CYP51 have higher mobility that may be necessary to facilitate the passage of its large sterol ligands. The water (W) tunnel is accessible to solvent during most of the simulations of CYP51, but its width is affected by the conformations of the heme's two propionate groups. These differ from those observed in the other CYPs studied because of differences in their hydrogen-bonding network. Our simulations give insights into the dynamics of CYP51 that complement the available experimental data and have implications for drug design against CYP51 enzymes.

P450cin active site water: implications for substrate binding and solvent accessibility

In P450cin, Tyr81, Asp241, Asn242, two water molecules, and the substrate participate in a complex H-bonded network. The role of this H-bonded network in substrate binding and catalysis has been probed by crystallography, spectroscopy, kinetics, isothermal titration calorimetry (ITC), and molecular dynamics. For the Y81F mutant, the substrate binds about 20-fold more weakly and Vmax decreases by about 30% in comparison to WT. The enhanced susceptibility of the heme to H₂O₂-mediated destruction in Y81F suggests that this mutant favors the open, low-spin conformational state. Asn242 H-bonds directly with the substrate, and replacing this residue with Ala results in water taking the place of the missing Asn side chain. This mutant exhibits a 70% decrease in activity. Crystal structures and molecular dynamics simulations of substrate-bound complexes show that the solvent has more ready access to the active site, especially for the N242A mutant. This accounts for about a 64% uncoupling of electron transfer from substrate hydroxylation. These data indicate the importance of the interconnected water network on substrate binding and on the open/closed conformational equilibrium, which are both critically important for maintaining high-coupling efficiency.

Zaccai neutron resilience and site-specific hydration dynamics in a globular protein

A discussion is presented of contributions of the Zaccai group to the understanding of flexibility in biological macromolecules using dynamic neutron scattering. The concept of resilience as introduced by Zaccai is discussed and investigated using molecular dynamics simulation on camphor-bound cytochrome P450. The resilience of hydrophilic residues is found to be more strongly affected by hydration than that of hydrophobic counterparts. The hydration-induced softening of protein propagates from the surface into the dry core. Moreover, buried hydrophilic residues behave more like those exposed on the protein surface, and are different from their hydrophobic counterparts.