The gating mechanism of the bacterial mechanosensitive channel MscL revealed by molecular dynamics simulations

One of the ultimate goals of the study on mechanosensitive (MS) channels is to understand the biophysical mechanisms of how the MS channel protein senses forces and how the sensed force induces channel gating. The bacterial MS channel MscL is an ideal subject to reach this goal owing to its resolved 3D protein structure in the closed state on the atomic scale and large amounts of electrophysiological data on its gating kinetics. However, the structural basis of the dynamic process from the closed to open states in MscL is not fully understood. In this study, we performed molecular dynamics (MD) simulations on the initial process of MscL opening in response to a tension increase in the lipid bilayer. To identify the tension-sensing site(s) in the channel protein, we calculated interaction energy between membrane lipids and candidate amino acids (AAs) facing the lipids. We found that Phe78 has a conspicuous interaction with the lipids, suggesting that Phe78 is the primary tension sensor of MscL. Increased membrane tension by membrane stretch dragged radially the inner (TM1) and outer (TM2) helices of MscL at Phe78, and the force was transmitted to the pentagon-shaped gate that is formed by the crossing of the neighboring TM1 helices in the inner leaflet of the bilayer. The radial dragging force induced radial sliding of the crossing portions, leading to a gate expansion. Calculated energy for this expansion is comparable to an experimentally estimated energy difference between the closed and the first subconductance state, suggesting that our model simulates the initial step toward the full opening of MscL. The model also successfully mimicked the behaviors of a gain of function mutant (G22N) and a loss of function mutant (F78N), strongly supporting that our MD model did simulate some essential biophysical aspects of the mechano-gating in MscL.     

The dynamics of protein-protein interactions between domains of MscL at the cytoplasmic-lipid interface

The bacterial mechanosensitive channel of large conductance, MscL, is one of the best characterized mechanosensitive channels serving as a paradigm for how proteins can sense and transduce mechanical forces. The physiological role of MscL is that of an emergency release valve that opens a large pore upon a sudden drop in the osmolarity of the environment. A crystal structure of a closed state of MscL shows it as a homopentamer, with each subunit consisting of two transmembrane domains (TM). There is consensus that the TM helices move in an iris like manner tilting in the plane of the membrane while gating. An N-terminal amphipathic helix that lies along the cytoplasmic membrane (S1), and the portion of TM2 near the cytoplasmic interface (TM2_ci), are relatively close in the crystal structure, yet predicted to be dynamic upon gating. Here we determine how these two regions interact in the channel complex, and study how these interactions change as the channel opens. We have screened 143 double-cysteine mutants of E. coli MscL for their efficiency in disulfide bridging and generated a map of protein-protein interactions between these two regions. Interesting candidates have been further studied by patch clamp and show differences in channel activity under different redox potentials; the results suggest a model for the dynamics of these two domains during MscL gating.     

The dynamics of protein-protein interactions between domains of MscL at the cytoplasmic-lipid interface

The bacterial mechanosensitive channel of large conductance, MscL, is one of the best characterized mechanosensitive channels serving as a paradigm for how proteins can sense and transduce mechanical forces. The physiological role of MscL is that of an emergency release valve that opens a large pore upon a sudden drop in the osmolarity of the environment. A crystal structure of a closed state of MscL shows it as a homopentamer, with each subunit consisting of two transmembrane domains (TM). There is consensus that the TM helices move in an iris like manner tilting in the plane of the membrane while gating. An N-terminal amphipathic helix that lies along the cytoplasmic membrane (S1), and the portion of TM2 near the cytoplasmic interface (TM2_ci), are relatively close in the crystal structure, yet predicted to be dynamic upon gating. Here we determine how these two regions interact in the channel complex, and study how these interactions change as the channel opens. We have screened 143 double-cysteine mutants of E. coli MscL for their efficiency in disulfide bridging and generated a map of protein-protein interactions between these two regions. Interesting candidates have been further studied by patch clamp and show differences in channel activity under different redox potentials; the results suggest a model for the dynamics of these two domains during MscL gating.     

Energetic and spatial parameters for gating of the bacterial large conductance mechanosensitive channel, MscL

MscL is multimeric protein that forms a large conductance mechanosensitive channel in the inner membrane of Escherichia coli. Since MscL is gated by tension transmitted through the lipid bilayer, we have been able to measure its gating parameters as a function of absolute tension. Using purified MscL reconstituted in liposomes, we recorded single channel currents and varied the pressure gradient (P) to vary the tension (T). The tension was calculated from P and the radius of curvature was obtained using video microscopy of the patch. The probability of being open (Po) has a steep sigmoidal dependence on T, with a midpoint (T1/2) of 11.8 dyn/cm. The maximal slope sensitivity of Po/Pc was 0.63 dyn/cm per e-fold. Assuming a Boltzmann distribution, the energy difference between the closed and fully open states in the unstressed membrane was DeltaE = 18.6 kBT. If the mechanosensitivity arises from tension acting on a change of in-plane area (DeltaA), the free energy, TDeltaA, would correspond to DeltaA = 6.5 nm2. MscL is not a binary channel, but has four conducting states and a closed state. Most transition rates are independent of tension, but the rate-limiting step to opening is the transition between the closed state and the lowest conductance substate. This transition thus involves the greatest DeltaA. When summed over all transitions, the in-plane area change from closed to fully open was 6 nm2, agreeing with the value obtained in the two-state analysis. Assuming a cylindrical channel, the dimensions of the (fully open) pore were comparable to DeltaA. Thus, the tension dependence of channel gating is primarily one of increasing the external channel area to accommodate the pore of the smallest conducting state. The higher conducting states appear to involve conformational changes internal to the channel that don't involve changes in area.     

Toward elucidating the heat activation mechanism of the TRPV1 channel gating by molecular dynamics simulation

As a key cellular sensor, the TRPV1 cation channel undergoes a gating transition from a closed state to an open state in response to various physical and chemical stimuli including noxious heat. Despite years of study, the heat activation mechanism of TRPV1 gating remains enigmatic at the molecular level. Toward elucidating the structural and energetic basis of TRPV1 gating, we have performed molecular dynamics (MD) simulations (with cumulative simulation time of 3 μs), starting from the high‐resolution closed and open structures of TRPV1 solved by cryo‐electron microscopy. In the closed‐state simulations at 30°C, we observed a stably closed channel constricted at the lower gate (near residue I679), while the upper gate (near residues G643 and M644) is dynamic and undergoes flickery opening/closing. In the open‐state simulations at 60°C, we found higher conformational variation consistent with a large entropy increase required for the heat activation, and both the lower and upper gates are dynamic with transient opening/closing. Through ensemble‐based structural analyses of the closed state versus the open state, we revealed pronounced closed‐to‐open conformational changes involving the membrane proximal domain (MPD) linker, the outer pore, and the TRP helix, which are accompanied by breaking/forming of a network of closed/open‐state specific hydrogen bonds. By comparing the closed‐state simulations at 30°C and 60°C, we observed heat‐activated conformational changes in the MPD linker, the outer pore, and the TRP helix that resemble the closed‐to‐open conformational changes, along with partial formation of the open‐state specific hydrogen bonds. Some of the residues involved in the above key hydrogen bonds were validated by previous mutational studies. Taken together, our MD simulations have offered rich structural and dynamic details beyond the static structures of TRPV1, and promising targets for future mutagenesis and functional studies of the TRPV1 channel. Proteins 2016; 84:1938–1949. © 2016 Wiley Periodicals, Inc.     

On the conformation of the COOH-terminal domain of the large mechanosensitive channel MscL

COOH-terminal (S3) domains are conserved within the MscL family of bacterial mechanosensitive channels, but their function remains unclear. The X-ray structure of MscL from Mycobacterium tuberculosis (TbMscL) revealed cytoplasmic domains forming a pentameric bundle (Chang, G., R.H. Spencer, A.T. Lee, M.T. Barclay, and D.C. Rees. 1998. SCIENCE: 282:2220-2226). The helices, however, have an unusual orientation in which hydrophobic sidechains face outside while charged residues face inside, possibly due to specific crystallization conditions. Based on the structure of pentameric cartilage protein, we modeled the COOH-terminal region of E. coli MscL to better satisfy the hydrophobicity criteria, with sidechains of conserved aliphatic residues all inside the bundle. Molecular dynamic simulations predicted higher stability for this conformation compared with one modeled after the crystal structure of TbMscL, and suggested distances for disulfide trapping experiments. The single cysteine mutants L121C and I125C formed dimers under ambient conditions and more so in the presence of an oxidant. The double-cysteine mutants, L121C/L122C and L128C/L129C, often cross-link into tetrameric and pentameric structures, consistent with the new model. Patch-clamp examination of these double mutants under moderately oxidizing or reducing conditions indicated that the bundle cross-linking neither prevents the channel from opening nor changes thermodynamic parameters of gating. Destabilization of the bundle by replacing conservative leucines with small polar residues, or complete removal of COOH-terminal domain (Delta110-136 mutation), increased the occupancy of subconducting states but did not change gating parameters substantially. The Delta110-136 truncation mutant was functional in in vivo osmotic shock assays; however, the amount of ATP released into the shock medium was considerably larger than in controls. The data strongly suggest that in contrast to previous gating models (Sukharev, S., M. Betanzos, C.S. Chiang, and H.R. Guy. 2001a. NATURE: 409:720-724.), S3 domains are stably associated in both closed and open conformations. The bundle-like assembly of cytoplasmic helices provides stability to the open conformation, and may function as a size-exclusion filter at the cytoplasmic entrance to the MscL pore, preventing loss of essential metabolites.     

Conversion of a mechanosensitive channel protein from a membrane-embedded to a water-soluble form by covalent modification with amphiphiles

Covalent modification of integral membrane proteins with amphiphiles may provide a general approach to the conversion of membrane proteins into water-soluble forms for biophysical and high-resolution structural studies. To test this approach, we mutated four surface residues of the pentameric Mycobacterium tuberculosis mechanosensitive channel of large conductance (MscL) to cysteine residues as anchors for amphiphile attachment. A series of modified ion channels with four amphiphile groups attached per channel subunit was prepared. One construct showed the highest water solubility to a concentration of up to 4mg/ml in the absence of detergent. This analog also formed native-like, alpha-helical homo-pentamers in the absence of detergent as judged by circular dichroism spectroscopy, size-exclusion chromatography and various light-scattering techniques. Proteins with longer, or shorter polymers attached, or proteins modified exclusively with polar cysteine-reactive small molecules, exhibited reduced to no solubility and higher-order aggregation. Electron microscopy revealed a homogeneous population of particles consistent with a pentameric channel. Solubilization of membrane proteins by covalent attachment of amphiphiles results in homogeneous particles that may prove useful for crystallization, solution NMR spectroscopy, and electron microscopy.     

Fundamental gating mechanism of nicotinic receptor channel revealed by mutation causing a congenital myasthenic syndrome

We describe the genetic and kinetic defects in a congenital myasthenic syndrome due to the mutation epsilonA411P in the amphipathic helix of the acetylcholine receptor (AChR) epsilon subunit. Myasthenic patients from three unrelated families are either homozygous for epsilonA411P or are heterozygous and harbor a null mutation in the second epsilon allele, indicating that epsilonA411P is recessive. We expressed human AChRs containing wild-type or A411P epsilon subunits in 293HEK cells, recorded single channel currents at high bandwidth, and determined microscopic rate constants for individual channels using hidden Markov modeling. For individual wild-type and mutant channels, each rate constant distributes as a Gaussian function, but the spread in the distributions for channel opening and closing rate constants is greatly expanded by epsilonA411P. Prolines engineered into positions flanking residue 411 of the epsilon subunit greatly increase the range of activation kinetics similar to epsilonA411P, whereas prolines engineered into positions equivalent to epsilonA411 in beta and delta subunits are without effect. Thus, the amphipathic helix of the epsilon subunit stabilizes the channel, minimizing the number and range of kinetic modes accessible to individual AChRs. The findings suggest that analogous stabilizing structures are present in other ion channels, and possibly allosteric proteins in general, and that they evolved to maintain uniformity of activation episodes. The findings further suggest that the fundamental gating mechanism of the AChR channel can be explained by a corrugated energy landscape superimposed on a steeply sloped energy well.     

Computer simulation of the KvAP voltage-gated potassium channel: steered molecular dynamics of the voltage sensor

The recent crystal structures of the voltage-gated potassium channel KvAP and its isolated voltage-sensing 'paddle' (composed of segments S1-S4) challenge existing models of voltage gating and raise a number of questions about the structure of the physiologically relevant state. We investigate a possible gating mechanism based on the crystal structures in a 10 ns steered molecular dynamics simulation of KvAP in a membrane-mimetic octane layer. The structure of the full KvAP protein has been modified by restraining the S2-S4 domain to the conformation of the isolated high-resolution paddle structure. After an initial relaxation, the paddle tips are pulled through the membrane from the intracellular to the extracellular side, corresponding to a putative change from closed to open. We describe the effect of this large-scale motion on the central pore domain, which remains largely unchanged, on the protein hydrogen-bonding network and on solvent. We analyze the motion of the S3b-S4 portion of the protein and propose a possible coupling mechanism between the paddle motion and the opening of the channel. Interactions between the arginine residues in S4, solvent and chloride ions are likely to play a role in the gating charge.     

Insight into the key features for ligand binding in Y1230 mutated c-Met kinase domain by molecular dynamics simulations

c-Met kinase has been considered as an attractive target for developing antitumor agents. The strong interactions between Tyr1230 and the inhibitors emphasized its importance for ligand binding. The clinically related Tyr1230 mutations have made negative impacts on current c-Met kinase inhibitors, especially the exquisitely selective ones, like PF-04217903, while the multi-targeted inhibitors, like Crizotinib, were not affected so much. In this study, the protein–ligand interactions between c-Met kinase domain (wild, Y1230C and Y1230H) and these inhibitors were compared. The binding site was expanded and the post-mutated regions became solvent accessible. The heavy dependency of PF-04217903 on the interactions with Tyr1230 resulted in the steep decrease of its potency against the Y1230 mutants. It was found that the ligand entrance region contributed consistently to the binding of Crizotinib, but not PF-04217903. Additional groups substituted in the ligand entrance region with stable interactions should be beneficial for improving the inhibitory activity of PF-04217903 against the Y1230 mutants. These findings will facilitate the discovery of potent inhibitors against Y1230 mutated c-Met kinase.     

The mechanism of high affinity pentasaccharide binding to antithrombin,insights from Gaussian accelerated molecular dynamics simulations

Abstract

The activity of antithrombin (AT), a serpin protease inhibitor, is enhanced by heparin and heparin analogs against its target proteases, mainly thrombin, factors Xa and IXa. Considerable amount of information is available on the multistep mechanism of the heparin pentasaccharide binding and conformational activation. However, much of the details were inferred from ‘static’ structures obtained by X-ray diffraction. Moreover, limited information is available for the early steps of binding mechanism other than kinetic studies with various ligands. To gain insights into these processes, we performed enhanced sampling molecular dynamics (MD) simulations using the Gaussian Accelerated Molecular Dynamics (GAMD) method, applied previously in drug binding studies. We were able to observe the binding of the pentasaccharide idraparinux to a ‘non-activated’ AT conformation in two separate trajectories with low root mean square deviation (RMSD) values compared to X-ray structures of the bound state. These trajectories along with further simulations of the AT-pentasaccharide complex provided insights into the mechanisms of multiple conformational transitions, including the expulsion of the hinge region, the extension of helix D and the conformational behavior of the reactive center loop (RCL). We could also confirm the high stability of helix P in non-activated AT conformations, such states might play an important role in heparin binding. ‘Generalized correlation’ matrices revealed possible paths of allosteric signal propagation to the binding sites for the target proteases, factors Xa and IXa. Enhanced MD simulations of ligand binding to AT may assist the design of new anticoagulant drugs.

Communicated by Ramaswamy H. Sarma     

Investigation of molecular dissolution mechanism of ketoprofen binary and ternary solid dispersions by molecular dynamics simulations

Abstract

The research aimed to investigate the molecular dissolution mechanism of both binary and ternary solid dispersions by molecular dynamics simulations. The simulation results indicated that the drug molecules were much easier to be released from surfactant-containing ternary systems than from binary ones. Moreover, sodium dodecyl sulfate as an additive in ternary systems had better effects than Tween 60. The simulation results were in well agreement with the experimental results. This research presented a reasonable explanation of molecular dissolution mechanism for both binary and ternary solid dispersions, which may benefit the future development of solid dispersion formulations.     

Investigating the size and dynamics of voltage-gated sodium channel fenestrations: A molecular dynamics study

Eukaryotic voltage-gated sodium channels (VGSCs) are essential for the initiation and propagation of action potentials in electrically excitable cells, and are important pharmaceutical targets for the treatment of neurological disorders such as epilepsy, cardiac arrhythmias, and chronic pain. Evidence suggests that small, hydrophobic, VGSC-blocking drugs can gain access to binding residues within the central cavity of these channels by passing through lateral, lipid-filled “fenestrations” which run between the exterior of the protein and its central pore. Here, we use molecular dynamics simulations to investigate how the size and shape of fenestrations change over time in several bacterial VGSC models and a homology model of Nav1.4. We show that over the course of the simulations, the size of the fenestrations is primarily influenced by rapid protein motions, such as amino acid side-chain rotation, and highlight that differences between fenestration bottleneck-contributing residues are the primary cause of variations in fenestration size between the 6 bacterial models. In the eukaryotic channel model, 2 fenestrations are wide, but 2 are narrow due to differences in the amino acid sequence in the 4 domains. Lipid molecules are found to influence the size of the fenestrations by protruding acyl chains into the fenestrations and displacing amino acid side-chains. Together, the results suggest that fenestrations provide viable pathways for small, flexible, hydrophobic drugs.     

Gating of the bacterial sodium channel, NaChBac: voltage-dependent charge movement and gating currents

The bacterial sodium channel, NaChBac, from Bacillus halodurans provides an excellent model to study structure-function relationships of voltage-gated ion channels. It can be expressed in mammalian cells for functional studies as well as in bacterial cultures as starting material for protein purification for fine biochemical and biophysical studies. Macroscopic functional properties of NaChBac have been described previously (Ren, D., B. Navarro, H. Xu, L. Yue, Q. Shi, and D.E. Clapham. 2001. Science. 294:2372-2375). In this study, we report gating current properties of NaChBac expressed in COS-1 cells. Upon depolarization of the membrane, gating currents appeared as upward inflections preceding the ionic currents. Gating currents were detectable at -90 mV while holding at -150 mV. Charge-voltage (Q-V) curves showed sigmoidal dependence on voltage with gating charge saturating at -10 mV. Charge movement was shifted by -22 mV relative to the conductance-voltage curve, indicating the presence of more than one closed state. Consistent with this was the Cole-Moore shift of 533 micros observed for a change in preconditioning voltage from -160 to -80 mV. The total gating charge was estimated to be 16 elementary charges per channel. Charge immobilization caused by prolonged depolarization was also observed; Q-V curves were shifted by approximately -60 mV to hyperpolarized potentials when cells were held at 0 mV. The kinetic properties of NaChBac were simulated by simultaneous fit of sodium currents at various voltages to a sequential kinetic model. Gating current kinetics predicted from ionic current experiments resembled the experimental data, indicating that gating currents are coupled to activation of NaChBac and confirming the assertion that this channel undergoes several transitions between closed states before channel opening. The results indicate that NaChBac has several closed states with voltage-dependent transitions between them realized by translocation of gating charge that causes activation of the channel.     

Structure and interactions of RecA: plasticity revealed by molecular dynamics simulations

Eleven independent simulations, each involving three consecutive molecules in the RecA filament, carried out on the protein from Mycobacterium tuberculosis, Mycobacterium smegmatis and Escherichia coli and their Adenosine triphosphate (ATP) complexes, provide valuable information which is complementary to that obtained from crystal structures, in addition to confirming the robust common structural framework within which RecA molecules from different eubacteria function. Functionally important loops, which are largely disordered in crystal structures, appear to adopt in each simulation subsets of conformations from larger ensembles. The simulations indicate the possibility of additional interactions involving the P-loop which remains largely invariant. The phosphate tail of the ATP is firmly anchored on the loop while the nucleoside moiety exhibits substantial structural variability. The most important consequence of ATP binding is the movement of the ‘switch’ residue. The relevant simulations indicate the feasibility of a second nucleotide binding site, but the pathway between adjacent molecules in the filament involving the two nucleotide binding sites appears to be possible only in the mycobacterial proteins.     

Molecular dissection of the large mechanosensitive ion channel (MscL) of E. coli: Mutants with altered channel gating and pressure sensitivity

In the search for the essential functional domains of the large mechanosensitive ion channel (MscL) of E. coli, we have cloned several mutants of the mscL gene into a glutathione S-transferase fusion protein expression system. The resulting mutated MscL proteins had either amino acid additions, substitutions or deletions in the amphipathic N-terminal region, and/or deletions in the amphipathic central or hydrophilic C-terminal regions. Proteolytic digestion of the isolated fusion proteins by thrombin yielded virtually pure recombinant MscL proteins that were reconstituted into artificial liposomes and examined for function by the patch-clamp technique. The addition of amino acid residues to the N-terminus of the MscL did not affect channel activity, whereas N-terminal deletions or changes to the N-terminal amino acid sequence were poorly tolerated and resulted in channels exhibiting altered pressure sensitivity and gating. Deletion of 27 amino acids from the C-terminus resulted in MscL protein that formed channels similar to the wild-type, while deletion of 33 C-terminal amino acids extinguished channel activity. Similarly, deletion of the internal amphipathic region of the MscL abolished activity. In accordance with a recently proposed spatial model of the MscL, our results suggest that (i) the N-terminal portion participates in the channel activation by pressure, and (ii) the essential channel functions are associated with both, the putative central amphipathic α-helical portion of the protein and the six C-terminal residues RKKEEP forming a charge cluster following the putative M2 membrane spanning α-helix. Received: 25 September 1996/Revised: 21 November 1996     

分子动力学模拟尿素对水溶液中蛋白质构象转变的影响

采用分子动力学方法和全原子模型研究尿素和水分子对模型蛋白S-肽链结构转化的影响。模拟结果显示S-肽链的变性速率常数k值随着尿素浓度的增加而先降低后升高,在尿素浓度为2.9 mol/L时达到最低值。模拟了不同尿素浓度下尿素-肽链、水-肽链以及肽链分子氢键的形成状况。结果表明:尿素浓度较低时,尿素分子与S-肽链的极性氨基酸侧链形成氢键,但不破坏其分子内的骨架氢键,尿素在S-肽链水化层外形成限制性空间,增强了S-肽链的稳定性。随着尿素的升高,尿素分子进入S-肽链内部并与其内部氨基酸残基形成氢键,导致S-肽链的骨架氢键丧失,S-肽链发生去折叠。上述模拟结果与文献报道的实验结果一致,从分子水平上揭示了尿素对蛋白质分子结构变化的影响机制,对于研究和发展蛋白质折叠及稳定化技术具有指导意义。     

Conformational dynamics of the ligand-binding domain of inward rectifier K channels as revealed by molecular dynamics simulations: toward an understanding of Kir channel gating

Inward rectifier (Kir) potassium channels are characterized by two transmembrane helices per subunit, plus an intracellular C-terminal domain that controls channel gating in response to changes in concentration of various ligands. Based on the crystal structure of the tetrameric C-terminal domain of Kir3.1, it is possible to build a homology model of the ATP-binding C-terminal domain of Kir6.2. Molecular dynamics simulations have been used to probe the dynamics of Kir C-terminal domains and to explore the relationship between their dynamics and possible mechanisms of channel gating. Multiple simulations, each of 10 ns duration, have been performed for Kir3.1 (crystal structure) and Kir6.2 (homology model), in both their monomeric and tetrameric forms. The Kir6.2 simulations were performed with and without bound ATP. The results of the simulations reveal comparable conformational stability for the crystal structure and the homology model. There is some decrease in conformational flexibility when comparing the monomers with the tetramers, corresponding mainly to the subunit interfaces in the tetramer. The beta-phosphate of ATP interacts with the side chain of K185 in the Kir6.2 model and simulations. The flexibility of the Kir6.2 tetramer is not changed greatly by the presence of bound ATP, other than in two loop regions. Principal components analysis of the simulated dynamics suggests loss of symmetry in both the Kir3.1 and Kir6.2 tetramers, consistent with "dimer-of-dimers" motion of subunits in C-terminal domains of the corresponding Kir channels. This is suggestive of a gating model in which a transition between exact tetrameric symmetry and dimer-of-dimers symmetry is associated with a change in transmembrane helix packing coupled to gating of the channel. Dimer-of-dimers motion of the C-terminal domain tetramer is also supported by coarse-grained (anisotropic network model) calculations. It is of interest that loss of exact rotational symmetry has also been suggested to play a role in gating in the bacterial Kir homolog, KirBac1.1, and in the nicotinic acetylcholine receptor channel.     

Dimerization of the p53 oligomerization domain: identification of a folding nucleus by molecular dynamics simulations

Dimerization of the p53 oligomerization domain involves coupled folding and binding of monomers. To examine the dimerization, we have performed molecular dynamics (MD) simulations of dimer folding from the rate-limiting transition state ensemble (TSE). Among 799 putative transition state structures that were selected from a large ensemble of high-temperature unfolding trajectories, 129 were identified as members of the TSE via calculation of a 50% transmission coefficient from at least 20 room-temperature simulations. This study is the first to examine the refolding of a protein dimer using MD simulations in explicit water, revealing a folding nucleus for dimerization. Our atomistic simulations are consistent with experiment and offer insight that was previously unobtainable.     

Investigation of ligand selectivity in CYP3A7 by molecular dynamics simulations

Cytochrome P450 (CYP) 3A7 plays a crucial role in the biotransformation of the metabolized endogenous and exogenous steroids. To compare the metabolic capabilities of CYP3A7–ligands complexes, three endogenous ligands were selected, namely dehydroepiandrosterone (DHEA), estrone, and estradiol. In this study, a three-dimensional model of CYP3A7 was constructed by homology modeling using the crystal structure of CYP3A4 as the template and refined by molecular dynamics simulation (MD). The docking method was adopted, combined with MD simulation and the molecular mechanics generalized born surface area method, to probe the ligand selectivity of CYP3A7. These results demonstrate that DHEA has the highest binding affinity, and the results of the binding free energy were in accordance with the experimental conclusion that estrone is better than estradiol. Moreover, several key residues responsible for substrate specificity were identified on the enzyme. Arg372 may be the most important residue due to the low interaction energies and the existence of hydrogen bond with DHEA throughout simulation. In addition, a cluster of Phe residues provides a hydrophobic environment to stabilize ligands. This study provides insights into the structural features of CYP3A7, which could contribute to further understanding of related protein structures and dynamics.