101
|
|
102
|
Swanson JMJ, Adcock SA, McCammon JA. Optimized Radii for Poisson-Boltzmann Calculations with the AMBER Force Field. J Chem Theory Comput 2015; 1:484-93. [PMID: 26641515 DOI: 10.1021/ct049834o] [Citation(s) in RCA: 62] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Implicit solvent models are a standard tool for assessing the electrostatics of biomolecular systems. The accuracy of quantitative predictions, such as pKa values, transfer free energies, binding energies, and solvation forces, is strongly dependent on one's choice of continuum parameters: the solute charges, dielectric coefficient, and radii, which define the dielectric boundary. To ensure quantitative accuracy, these parameters can be benchmarked against explicit solvent simulations. Here we present two sets of optimized radii to define either abrupt or cubic-spline smoothed dielectric boundaries in Poisson-Boltzmann calculations of protein systems with AMBER (parm99) charges. Spline smoothing stabilizes the electrostatic potential at the molecular surface, allowing for continuum force calculations. Most implementations, however, require significantly different radii than the abrupt boundary surfaces. The optimal continuum radii are initially approximated from the solvent radial charge distribution surrounding each atom type. A genetic algorithm is then used to fine-tune the starting values to reproduce charging free energies measured from explicit solvent simulations. The optimized radii are tested on four protein-like polypeptides. The results show increased accuracy of molecular solvation energies and atomic forces relative to commonly used continuum parameter sets. These radii are suitable for Poisson-Boltzmann calculations with the AMBER force field and offer energetic congruence to any model that combines molecular mechanics and Poisson-Boltzmann solvation energies.
Collapse
|
103
|
Cerutti DS, Ten Eyck LF, McCammon JA. Rapid Estimation of Solvation Energy for Simulations of Protein-Protein Association. J Chem Theory Comput 2015; 1:143-52. [PMID: 26641125 DOI: 10.1021/ct049946f] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
We have formulated the Energy by Linear Superposition of Corrections Approximation (ELSCA) for estimating the electrostatic and apolar solvation energy of bringing two proteins into close proximity or into contact as defined by the linearized Poisson-Boltzmann model and a linear function of the solvent-accessible surface area. ELSCA utilizes potentials of mean force between atom types found in the AMBER ff99 force field, a uniform distance-dependent dielectric, and a potential that mimics the change in solvent accessible surface area for bringing two solvated spheres into contact. ELSCA was trained by a linear least-squares fit on more than 39 000 putative complexes, each formed from pairs of nonhomologous proteins with a range of shapes, sizes, and charges. The training set was also designed to capture various stages of complex formation. ELSCA was tested against over 8000 non-native complexes of 45 enzyme/inhibitor, antibody/antigen, and other systems that are known to form complexes and gives an overall correlation of 0.962 with PBSA-derived energies for these complexes. The predictions have a slope of 0.89 on the actual values with a bias of 11.1 kcal/mol. When applied to native complexes of these 45 protein systems, ELSCA reproduces PBSA results with a correlation of 0.787, a slope of 1.13, and a bias of 13.0 kcal/mol. We report parameters for ELSCA in the context of the AMBER ff99 parameter set. Our model is most useful in macromolecular docking and protein association simulations, where large portions of each molecule may be considered rigid.
Collapse
|
104
|
Swanson JMJ, Wagoner JA, Baker NA, McCammon JA. Optimizing the Poisson Dielectric Boundary with Explicit Solvent Forces and Energies: Lessons Learned with Atom-Centered Dielectric Functions. J Chem Theory Comput 2015; 3:170-83. [PMID: 26627162 DOI: 10.1021/ct600216k] [Citation(s) in RCA: 40] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Accurate implicit solvent models require parameters that have been optimized in a system- or atom-specific manner on the basis of experimental data or more rigorous explicit solvent simulations. Models based on the Poisson or Poisson-Boltzmann equation are particularly sensitive to the nature and location of the boundary which separates the low dielectric solute from the high dielectric solvent. Here, we present a novel method for optimizing the solute radii, which define the dielectric boundary, on the basis of forces and energies from explicit solvent simulations. We use this method to optimize radii for protein systems defined by AMBER ff99 partial charges and a spline-smoothed solute surface. The spline-smoothed surface is an atom-centered dielectric function that enables stable and efficient force calculations. We explore the relative performance of radii optimized with forces alone and those optimized with forces and energies. We show that our radii reproduce the explicit solvent forces and energies more accurately than four other parameter sets commonly used in conjunction with the AMBER force field, each of which has been appropriately scaled for spline-smoothed surfaces. Finally, we demonstrate that spline-smoothed surfaces show surprising accuracy for small, compact systems but may have limitations for highly solvated protein systems. The optimization method presented here is efficient and applicable to any system with explicit solvent parameters. It can be used to determine the optimal continuum parameters when experimental solvation energies are unavailable and the computational costs of explicit solvent charging free energies are prohibitive.
Collapse
|
105
|
Hamacher K, McCammon JA. Computing the Amino Acid Specificity of Fluctuations in Biomolecular Systems. J Chem Theory Comput 2015; 2:873-8. [PMID: 26626694 DOI: 10.1021/ct050247s] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We developed a new amino acid specific method for the computation of spatial fluctuations of proteins around their native structures. We show the consistency with experimental values and the increased performance in comparison to an established model, based on statistical estimates for a set of test proteins. We apply the new method to HIV-1 protease in its wild-type form and to a V82F-I84V mutant that shows resistance to protease inhibitors. We further show how the method can be successfully used to explain the molecular biophysics of drug resistance of the mutant.
Collapse
|
106
|
Lindert S, Cheng Y, Kekenes-Huskey P, Regnier M, McCammon JA. Effects of HCM cTnI mutation R145G on troponin structure and modulation by PKA phosphorylation elucidated by molecular dynamics simulations. Biophys J 2015; 108:395-407. [PMID: 25606687 DOI: 10.1016/j.bpj.2014.11.3461] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2014] [Revised: 11/21/2014] [Accepted: 11/21/2014] [Indexed: 10/24/2022] Open
Abstract
Cardiac troponin (cTn) is a key molecule in the regulation of human cardiac muscle contraction. The N-terminal cardiac-specific peptide of the inhibitory subunit of troponin, cTnI (cTnI(1-39)), is a target for phosphorylation by protein kinase A (PKA) during β-adrenergic stimulation. We recently presented evidence indicating that this peptide interacts with the inhibitory peptide (cTnl(137-147)) when S23 and S24 are phosphorylated. The inhibitory peptide is also the target of the point mutation cTnI-R145G, which is associated with hypertrophic cardiomyopathy (HCM), a disease associated with sudden death in apparently healthy young adults. It has been shown that both phosphorylation and this mutation alter the cTnC-cTnI (C-I) interaction, which plays a crucial role in modulating contractile activation. However, little is known about the molecular-level events underlying this modulation. Here, we computationally investigated the effects of the cTnI-R145G mutation on the dynamics of cTn, cTnC Ca(2+) handling, and the C-I interaction. Comparisons were made with the cTnI-R145G/S23D/S24D phosphomimic mutation, which has been used both experimentally and computationally to study the cTnI N-terminal specific effects of PKA phosphorylation. Additional comparisons between the phosphomimic mutations and the real phosphorylations were made. For this purpose, we ran triplicate 150 ns molecular dynamics simulations of cTnI-R145G Ca(2+)-bound cTnC(1-161)-cTnI(1-172)-cTnT(236-285), cTnI-R145G/S23D/S24D Ca(2+)-bound cTnC(1-161)-cTnI(1-172)-cTnT(236-285), and cTnI-R145G/PS23/PS24 Ca(2+)-bound cTnC(1-161)-cTnI(1-172)-cTnT(236-285), respectively. We found that the cTnI-R145G mutation did not impact the overall dynamics of cTn, but stabilized crucial Ca(2+)-coordinating interactions. However, the phosphomimic mutations increased overall cTn fluctuations and destabilized Ca(2+) coordination. Interestingly, cTnI-R145G blunted the intrasubunit interactions between the cTnI N-terminal extension and the cTnI inhibitory peptide, which have been suggested to play a crucial role in modulating troponin function during β-adrenergic stimulation. These findings offer a molecular-level explanation for how the HCM mutation cTnI-R145G reduces the modulation of cTn by phosphorylation of S23/S24 during β-adrenergic stimulation.
Collapse
|
107
|
Cheng Y, Rao V, Tu AY, Lindert S, Wang D, Oxenford L, McCulloch AD, McCammon JA, Regnier M. Troponin I Mutations R146G and R21C Alter Cardiac Troponin Function, Contractile Properties, and Modulation by Protein Kinase A (PKA)-mediated Phosphorylation. J Biol Chem 2015; 290:27749-66. [PMID: 26391394 DOI: 10.1074/jbc.m115.683045] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2015] [Indexed: 11/06/2022] Open
Abstract
Two hypertrophic cardiomyopathy-associated cardiac troponin I (cTnI) mutations, R146G and R21C, are located in different regions of cTnI, the inhibitory peptide and the cardiac-specific N terminus. We recently reported that these regions may interact when Ser-23/Ser-24 are phosphorylated, weakening the interaction of cTnI with cardiac TnC. Little is known about how these mutations influence the affinity of cardiac TnC for cTnI (KC-I) or contractile kinetics during β-adrenergic stimulation. Here, we tested how cTnI(R146G) or cTnI(R21C) influences contractile activation and relaxation and their response to protein kinase A (PKA). Both mutations significantly increased Ca(2+) binding affinity to cTn (KCa) and KC-I. PKA phosphorylation resulted in a similar reduction of KCa for all complexes, but KC-I was reduced only with cTnI(WT). cTnI(WT), cTnI(R146G), and cTnI(R21C) were complexed into cardiac troponin and exchanged into rat ventricular myofibrils, and contraction/relaxation kinetics were measured ± PKA phosphorylation. Maximal tension (Tmax) was maintained for cTnI(R146G)- and cTnI(R21C)-exchanged myofibrils, and Ca(2+) sensitivity of tension (pCa50) was increased. PKA phosphorylation decreased pCa50 for cTnI(WT)-exchanged myofibrils but not for either mutation. PKA phosphorylation accelerated the early slow phase relaxation for cTnI(WT) myofibrils, especially at Ca(2+) levels that the heart operates in vivo. Importantly, this effect was blunted for cTnI(R146G)- and cTnI(R21C)-exchanged myofibrils. Molecular dynamics simulations suggest both mutations inhibit formation of intra-subunit contacts between the N terminus and the inhibitory peptide of cTnI that is normally seen with WT-cTn upon PKA phosphorylation. Together, our results suggest that cTnI(R146G) and cTnI(R21C) blunt PKA modulation of activation and relaxation kinetics by prohibiting cardiac-specific N-terminal interaction with the cTnI inhibitory peptide.
Collapse
|
108
|
Huang YMM, Huber G, McCammon JA. Electrostatic steering enhances the rate of cAMP binding to phosphodiesterase: Brownian dynamics modeling. Protein Sci 2015; 24:1884-9. [PMID: 26346301 DOI: 10.1002/pro.2794] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2015] [Revised: 08/20/2015] [Accepted: 08/24/2015] [Indexed: 01/18/2023]
Abstract
Signaling in cells often involves co-localization of the signaling molecules. Most experimental evidence has shown that intracellular compartmentalization restricts the range of action of the second messenger, 3'-5'-cyclic adenosine monophosphate (cAMP), which is degraded by phosphodiesterases (PDEs). The objective of this study is to understand the details of molecular encounter that may play a role in efficient operation of the cAMP signaling apparatus. The results from electrostatic potential calculations and Brownian dynamics simulations suggest that positive potential of the active site from PDE enhances capture of diffusing cAMP molecules. This electrostatic steering between cAMP and the active site of a PDE plays a major role in the enzyme-substrate encounter, an effect that may be of significance in sequestering cAMP released from a nearby binding site or in attracting more freely diffusing cAMP molecules.
Collapse
|
109
|
Rao V, Cheng Y, Lindert S, Wang D, Oxenford L, McCulloch AD, McCammon JA, Regnier M. PKA phosphorylation of cardiac troponin I modulates activation and relaxation kinetics of ventricular myofibrils. Biophys J 2015; 107:1196-1204. [PMID: 25185555 DOI: 10.1016/j.bpj.2014.07.027] [Citation(s) in RCA: 42] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2014] [Revised: 07/10/2014] [Accepted: 07/15/2014] [Indexed: 10/24/2022] Open
Abstract
Protein kinase A (PKA) phosphorylation of myofibril proteins constitutes an important pathway for β-adrenergic modulation of cardiac contractility and relaxation. PKA targets the N-terminus (Ser-23/24) of cardiac troponin I (cTnI), cardiac myosin-binding protein C (cMyBP-C) and titin. The effect of PKA-mediated phosphorylation on the magnitude of contraction has been studied in some detail, but little is known about how it modulates the kinetics of thin filament activation and myofibril relaxation as Ca(2+) levels vary. Troponin C (cTnC) interaction with cTnI (C-I interaction) is a critical step in contractile activation that can be modulated by cTnI phosphorylation. We tested the hypothesis that altering C-I interactions by PKA, or by cTnI phosphomimetic mutations (S23D/S24D-cTnI), directly affects thin filament activation and myofilament relaxation kinetics. Rat ventricular myofibrils were isolated and endogenous cTn was exchanged with either wild-type cTnI, or S23D/S24D-cTnI recombinant cTn. Contractile mechanics were monitored at maximum and submaximal Ca(2+) concentrations. PKA treatment of wild-type cTn or exchange of cTn containing S23D/S24D-cTnI resulted in an increase in the rate of early, slow phase of relaxation (kREL,slow) and a decrease in its duration (tREL,slow). These effects were greater for submaximal Ca(2+) activated contractions. PKA treatment also reduced the rate of contractile activation (kACT) at maximal, but not submaximal Ca(2+), and reduced the Ca(2+) sensitivity of contraction. Using a fluorescent probe coupled to cTnC (C35S-IANBD), the Ca(2+)-cTn binding affinity and C-I interaction were monitored. Ca(2+) binding to cTn (pCa50) was significantly decreased when cTnI was phosphorylated by PKA (ΔpCa50 = 0.31). PKA phosphorylation of cTnI also weakened C-I interaction in the presence of Ca(2+). These data suggest that weakened C-I interaction, via PKA phosphorylation of cTnI, may slow thin filament activation and result in increased myofilament relaxation kinetics, the latter of which could enhance early phase diastolic relaxation during β-adrenergic stimulation.
Collapse
|
110
|
Miao Y, Feher VA, McCammon JA. Gaussian Accelerated Molecular Dynamics: Unconstrained Enhanced Sampling and Free Energy Calculation. J Chem Theory Comput 2015; 11:3584-3595. [PMID: 26300708 PMCID: PMC4535365 DOI: 10.1021/acs.jctc.5b00436] [Citation(s) in RCA: 439] [Impact Index Per Article: 48.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2015] [Indexed: 12/20/2022]
Abstract
A Gaussian accelerated molecular dynamics (GaMD) approach for simultaneous enhanced sampling and free energy calculation of biomolecules is presented. By constructing a boost potential that follows Gaussian distribution, accurate reweighting of the GaMD simulations is achieved using cumulant expansion to the second order. Here, GaMD is demonstrated on three biomolecular model systems: alanine dipeptide, chignolin folding, and ligand binding to the T4-lysozyme. Without the need to set predefined reaction coordinates, GaMD enables unconstrained enhanced sampling of these biomolecules. Furthermore, the free energy profiles obtained from reweighting of the GaMD simulations allow us to identify distinct low-energy states of the biomolecules and characterize the protein-folding and ligand-binding pathways quantitatively.
Collapse
|
111
|
Metzger VT, Eun C, Kekenes-Huskey PM, Huber G, McCammon JA. Electrostatic channeling in P. falciparum DHFR-TS: Brownian dynamics and Smoluchowski modeling. Biophys J 2015; 107:2394-402. [PMID: 25418308 DOI: 10.1016/j.bpj.2014.09.039] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2014] [Revised: 09/15/2014] [Accepted: 09/19/2014] [Indexed: 11/15/2022] Open
Abstract
We perform Brownian dynamics simulations and Smoluchowski continuum modeling of the bifunctional Plasmodium falciparum dihydrofolate reductase-thymidylate synthase (P. falciparum DHFR-TS) with the objective of understanding the electrostatic channeling of dihydrofolate generated at the TS active site to the DHFR active site. The results of Brownian dynamics simulations and Smoluchowski continuum modeling suggest that compared to Leishmania major DHFR-TS, P. falciparum DHFR-TS has a lower but significant electrostatic-mediated channeling efficiency (?15-25%) at physiological pH (7.0) and ionic strength (150 mM). We also find that removing the electric charges from key basic residues located between the DHFR and TS active sites significantly reduces the channeling efficiency of P. falciparum DHFR-TS. Although several protozoan DHFR-TS enzymes are known to have similar tertiary and quaternary structure, subtle differences in structure, active-site geometry, and charge distribution appear to influence both electrostatic-mediated and proximity-based substrate channeling.
Collapse
|
112
|
Wang N, McCammon JA. Substrate channeling between the human dihydrofolate reductase and thymidylate synthase. Protein Sci 2015; 25:79-86. [PMID: 26096018 DOI: 10.1002/pro.2720] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2015] [Revised: 06/02/2015] [Accepted: 06/05/2015] [Indexed: 12/17/2022]
Abstract
In vivo, as an advanced catalytic strategy, transient non-covalently bound multi-enzyme complexes can be formed to facilitate the relay of substrates, i. e. substrate channeling, between sequential enzymatic reactions and to enhance the throughput of multi-step enzymatic pathways. The human thymidylate synthase and dihydrofolate reductase catalyze two consecutive reactions in the folate metabolism pathway, and experiments have shown that they are very likely to bind in the same multi-enzyme complex in vivo. While reports on the protozoa thymidylate synthase-dihydrofolate reductase bifunctional enzyme give substantial evidences of substrate channeling along a surface "electrostatic highway," attention has not been paid to whether the human thymidylate synthase and dihydrofolate reductase, if they are in contact with each other in the multi-enzyme complex, are capable of substrate channeling employing surface electrostatics. This work utilizes protein-protein docking, electrostatics calculations, and Brownian dynamics to explore the existence and mechanism of the substrate channeling between the human thymidylate synthase and dihydrofolate reductase. The results show that the bound human thymidylate synthase and dihydrofolate reductase are capable of substrate channeling and the formation of the surface "electrostatic highway." The substrate channeling efficiency between the two can be reasonably high and comparable to that of the protozoa.
Collapse
|
113
|
Blachly PG, Sandala GM, Giammona DA, Bashford D, McCammon JA, Noodleman L. Broken-Symmetry DFT Computations for the Reaction Pathway of IspH, an Iron-Sulfur Enzyme in Pathogenic Bacteria. Inorg Chem 2015; 54:6439-61. [PMID: 26098647 PMCID: PMC4568833 DOI: 10.1021/acs.inorgchem.5b00751] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The recently discovered methylerythritol phosphate (MEP) pathway provides new targets for the development of antibacterial and antimalarial drugs. In the final step of the MEP pathway, the [4Fe-4S] IspH protein catalyzes the 2e(-)/2H(+) reductive dehydroxylation of (E)-4-hydroxy-3-methyl-but-2-enyl diphosphate (HMBPP) to afford the isoprenoid precursors isopentenyl pyrophosphate (IPP) and dimethylallyl pyrophosphate (DMAPP). Recent experiments have attempted to elucidate the IspH catalytic mechanism to drive inhibitor development. Two competing mechanisms have recently emerged, differentiated by their proposed HMBPP binding modes upon 1e(-) reduction of the [4Fe-4S] cluster: (1) a Birch reduction mechanism, in which HMBPP remains bound to the [4Fe-4S] cluster through its terminal C4-OH group (ROH-bound) until the -OH is cleaved as water; and (2) an organometallic mechanism, in which the C4-OH group rotates away from the [4Fe-4S] cluster, allowing the HMBPP olefin group to form a metallacycle complex with the apical iron (η(2)-bound). We perform broken-symmetry density functional theory computations to assess the energies and reduction potentials associated with the ROH- and η(2)-bound states implicated by these competing mechanisms. Reduction potentials obtained for ROH-bound states are more negative (-1.4 to -1.0 V) than what is typically expected of [4Fe-4S] ferredoxin proteins. Instead, we find that η(2)-bound states are lower in energy than ROH-bound states when the [4Fe-4S] cluster is 1e(-) reduced. Furthermore, η(2)-bound states can already be generated in the oxidized state, yielding reduction potentials of ca. -700 mV when electron addition occurs after rotation of the HMBPP C4-OH group. We demonstrate that such η(2)-bound states are kinetically accessible both when the IspH [4Fe-4S] cluster is oxidized and 1e(-) reduced. The energetically preferred pathway gives 1e(-) reduction of the cluster after substrate conformational change, generating the 1e(-) reduced intermediate proposed in the organometallic mechanism.
Collapse
|
114
|
Miao Y, Feixas F, Eun C, McCammon JA. Accelerated molecular dynamics simulations of protein folding. J Comput Chem 2015; 36:1536-49. [PMID: 26096263 DOI: 10.1002/jcc.23964] [Citation(s) in RCA: 110] [Impact Index Per Article: 12.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2015] [Revised: 05/11/2015] [Accepted: 05/19/2015] [Indexed: 02/02/2023]
Abstract
Folding of four fast-folding proteins, including chignolin, Trp-cage, villin headpiece and WW domain, was simulated via accelerated molecular dynamics (aMD). In comparison with hundred-of-microsecond timescale conventional molecular dynamics (cMD) simulations performed on the Anton supercomputer, aMD captured complete folding of the four proteins in significantly shorter simulation time. The folded protein conformations were found within 0.2-2.1 Å of the native NMR or X-ray crystal structures. Free energy profiles calculated through improved reweighting of the aMD simulations using cumulant expansion to the second-order are in good agreement with those obtained from cMD simulations. This allows us to identify distinct conformational states (e.g., unfolded and intermediate) other than the native structure and the protein folding energy barriers. Detailed analysis of protein secondary structures and local key residue interactions provided important insights into the protein folding pathways. Furthermore, the selections of force fields and aMD simulation parameters are discussed in detail. Our work shows usefulness and accuracy of aMD in studying protein folding, providing basic references in using aMD in future protein-folding studies.
Collapse
|
115
|
Caliman AD, Swift SE, Wang Y, Miao Y, McCammon JA. Investigation of the conformational dynamics of the apo A2A adenosine receptor. Protein Sci 2015; 24:1004-12. [PMID: 25761901 PMCID: PMC4456113 DOI: 10.1002/pro.2681] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2014] [Accepted: 03/08/2015] [Indexed: 12/19/2022]
Abstract
The activation/deactivation processes for G-protein coupled receptors (GPCRs) have been computationally studied for several different classes, including rhodopsin, the β2 adrenergic receptor, and the M2 muscarinic receptor. Despite determined cocrystal structures of the adenosine A2A receptor (A2A AR) in complex with antagonists, agonists and an antibody, the deactivation process of this GPCR is not completely understood. In this study, we investigate the convergence of two apo simulations, one starting with an agonist-bound conformation (PDB: 3QAK)(14) and the other starting with an antagonist-bound conformation (PDB: 3EML)(11) . Despite the two simulations not completely converging, we were able to identify distinct intermediate steps of the deactivation process characterized by the movement of Y288(7.53) in the NPxxY motif. We find that Y288(7.53) contributes to the process by forming hydrogen bonds to residues in transmembrane helices 2 and 7 and losing these interactions upon full deactivation. Y197(5.58) also plays a role in the process by forming a hydrogen bond only once the side chain moves from the lipid interface to the middle of the helical bundle.
Collapse
|
116
|
Zhou S, Cheng LT, Sun H, Che J, Dzubiella J, Li B, McCammon JA. LS-VISM: A software package for analysis of biomolecular solvation. J Comput Chem 2015; 36:1047-59. [PMID: 25766844 PMCID: PMC4412808 DOI: 10.1002/jcc.23890] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2014] [Accepted: 02/08/2015] [Indexed: 12/30/2022]
Abstract
We introduce a software package for the analysis of biomolecular solvation. The package collects computer codes that implement numerical methods for a variational implicit-solvent model (VISM). The input of the package includes the atomic data of biomolecules under consideration and the macroscopic parameters such as solute-solvent surface tension, bulk solvent density and ionic concentrations, and the dielectric coefficients. The output includes estimated solvation free energies and optimal macroscopic solute-solvent interfaces that are obtained by minimizing the VISM solvation free-energy functional among all possible solute-solvent interfaces enclosing the solute atoms. We review the VISM with various descriptions of electrostatics. We also review our numerical methods that consist mainly of the level-set method for relaxing the VISM free-energy functional and a compact coupling interface method for the dielectric Poisson-Boltzmann equation. Such numerical methods and algorithms constitute the central modules of the software package. We detail the structure of the package, format of input and output files, workflow of the codes, and the postprocessing of output data. Our demo application to a host-guest system illustrates how to use the package to perform solvation analysis for biomolecules, including ligand-receptor binding systems. The package is simple and flexible with respect to minimum adjustable parameters and a wide range of applications. Future extensions of the package use can include the efficient identification of ligand binding pockets on protein surfaces.
Collapse
|
117
|
Kaus JW, McCammon JA. Enhanced ligand sampling for relative protein-ligand binding free energy calculations. J Phys Chem B 2015; 119:6190-7. [PMID: 25906170 PMCID: PMC4442669 DOI: 10.1021/acs.jpcb.5b02348] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
Free
energy calculations are used to study how strongly potential
drug molecules interact with their target receptors. The accuracy
of these calculations depends on the accuracy of the molecular dynamics
(MD) force field as well as proper sampling of the major conformations
of each molecule. However, proper sampling of ligand conformations
can be difficult when there are large barriers separating the major
ligand conformations. An example of this is for ligands with an asymmetrically
substituted phenyl ring, where the presence of protein loops hinders
the proper sampling of the different ring conformations. These ring
conformations become more difficult to sample when the size of the
functional groups attached to the ring increases. The Adaptive Integration
Method (AIM) has been developed, which adaptively changes the alchemical
coupling parameter λ during the MD simulation so that conformations
sampled at one λ can aid sampling at the other λ values.
The Accelerated Adaptive Integration Method (AcclAIM) builds on AIM
by lowering potential barriers for specific degrees of freedom at
intermediate λ values. However, these methods may not work when
there are very large barriers separating the major ligand conformations.
In this work, we describe a modification to AIM that improves sampling
of the different ring conformations, even when there is a very large
barrier between them. This method combines AIM with conformational
Monte Carlo sampling, giving improved convergence of ring populations
and the resulting free energy. This method, called AIM/MC, is applied
to study the relative binding free energy for a pair of ligands that
bind to thrombin and a different pair of ligands that bind to aspartyl
protease β-APP cleaving enzyme 1 (BACE1). These protein–ligand
binding free energy calculations illustrate the improvements in conformational
sampling and the convergence of the free energy compared to both AIM
and AcclAIM.
Collapse
|
118
|
Kekenes-Huskey PM, Gillette AK, McCammon JA. Predicting the influence of long-range molecular interactions on macroscopic-scale diffusion by homogenization of the Smoluchowski equation. J Chem Phys 2015; 140:174106. [PMID: 24811624 DOI: 10.1063/1.4873382] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The macroscopic diffusion constant for a charged diffuser is in part dependent on (1) the volume excluded by solute "obstacles" and (2) long-range interactions between those obstacles and the diffuser. Increasing excluded volume reduces transport of the diffuser, while long-range interactions can either increase or decrease diffusivity, depending on the nature of the potential. We previously demonstrated [P. M. Kekenes-Huskey et al., Biophys. J. 105, 2130 (2013)] using homogenization theory that the configuration of molecular-scale obstacles can both hinder diffusion and induce diffusional anisotropy for small ions. As the density of molecular obstacles increases, van der Waals (vdW) and electrostatic interactions between obstacle and a diffuser become significant and can strongly influence the latter's diffusivity, which was neglected in our original model. Here, we extend this methodology to include a fixed (time-independent) potential of mean force, through homogenization of the Smoluchowski equation. We consider the diffusion of ions in crowded, hydrophilic environments at physiological ionic strengths and find that electrostatic and vdW interactions can enhance or depress effective diffusion rates for attractive or repulsive forces, respectively. Additionally, we show that the observed diffusion rate may be reduced independent of non-specific electrostatic and vdW interactions by treating obstacles that exhibit specific binding interactions as "buffers" that absorb free diffusers. Finally, we demonstrate that effective diffusion rates are sensitive to distribution of surface charge on a globular protein, Troponin C, suggesting that the use of molecular structures with atomistic-scale resolution can account for electrostatic influences on substrate transport. This approach offers new insight into the influence of molecular-scale, long-range interactions on transport of charged species, particularly for diffusion-influenced signaling events occurring in crowded cellular environments.
Collapse
|
119
|
Bertolini R, Goepfert C, Andrieu T, Nichols S, Walter MA, Frey FJ, McCammon JA, Frey BM. 18F-RB390: innovative ligand for imaging the T877A androgen receptor mutant in prostate cancer via positron emission tomography (PET). Prostate 2015; 75:348-59. [PMID: 25358634 DOI: 10.1002/pros.22919] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 03/13/2014] [Accepted: 09/17/2014] [Indexed: 11/09/2022]
Abstract
BACKGROUND Detecting prostate cancer before spreading or predicting a favorable therapy are challenging issues for impacting patient's survival. Presently, 2-[(18) F]-fluoro-2-deoxy-D-glucose ((18) F-FDG) and/or (18) F-fluorocholine ((18) F-FCH) are the generally used PET-tracers in oncology yet do not emphasize the T877A androgen receptor (AR) mutation being exclusively present in cancerous tissue and escaping androgen deprivation treatment. METHODS We designed and synthesized fluorinated 5α-dihydrotestosterone (DHT) derivatives to target T877A-AR. We performed binding assays to select suitable candidates using COS-7 cells transfected with wild-type or T877A AR (WT-AR, T877A-AR) expressing plasmids and investigated cellular uptake of candidate (18) F-RB390. Stability, biodistribution analyses and PET-Imaging were assessed by injecting (18) F-RB390 (10MBq), with and without co-injection of an excess of unlabeled DHT in C4-2 and PC-3 tumor bearing male SCID mice (n = 12). RESULTS RB390 presented a higher relative binding affinity (RBA) (28.1%, IC50 = 32 nM) for T877A-AR than for WT-AR (1.7%, IC50 = 357 nM) related to DHT (RBA = 100%). A small fraction of (18) F-RB390 was metabolized when incubated with murine liver homogenate or human blood for 3 hr. The metabolite of RB390, 3-hydroxysteroid RB448, presented similar binding characteristics as RB390. (18) F-RB390 but not (18) F-FDG or (18) F-FCH accumulated 2.5× more in COS-7 cells transfected with pSG5AR-T877A than with control plasmid. Accumulation was reduced with an excess of DHT. PET/CT imaging and biodistribution studies revealed a significantly higher uptake of (18) F-RB390 in T877A mutation positive xenografts compared to PC-3 control tumors. This effect was blunted with DHT. CONCLUSION Given the differential binding capacity and the favorable radioactivity pattern, (18) F-RB390 represents the portrayal of the first imaging ligand with predictive potential for mutant T877A-AR in prostate cancer for guiding therapy. Prostate 75:348-359, 2015. © 2014 Wiley Periodicals, Inc.
Collapse
|
120
|
Mouchlis VD, Bucher D, McCammon JA, Dennis EA. Membranes serve as allosteric activators of phospholipase A2, enabling it to extract, bind, and hydrolyze phospholipid substrates. Proc Natl Acad Sci U S A 2015; 112:E516-25. [PMID: 25624474 PMCID: PMC4330758 DOI: 10.1073/pnas.1424651112] [Citation(s) in RCA: 68] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023] Open
Abstract
Defining the molecular details and consequences of the association of water-soluble proteins with membranes is fundamental to understanding protein-lipid interactions and membrane functioning. Phospholipase A2 (PLA2) enzymes, which catalyze the hydrolysis of phospholipid substrates that compose the membrane bilayers, provide the ideal system for studying protein-lipid interactions. Our study focuses on understanding the catalytic cycle of two different human PLA2s: the cytosolic Group IVA cPLA2 and calcium-independent Group VIA iPLA2. Computer-aided techniques guided by deuterium exchange mass spectrometry data, were used to create structural complexes of each enzyme with a single phospholipid substrate molecule, whereas the substrate extraction process was studied using steered molecular dynamics simulations. Molecular dynamic simulations of the enzyme-substrate-membrane systems revealed important information about the mechanisms by which these enzymes associate with the membrane and then extract and bind their phospholipid substrate. Our data support the hypothesis that the membrane acts as an allosteric ligand that binds at the allosteric site of the enzyme's interfacial surface, shifting its conformation from a closed (inactive) state in water to an open (active) state at the membrane interface.
Collapse
|
121
|
Lindert S, Li MX, Sykes B, McCammon JA. Computer-Aided Drug Discovery Approach Finds Calcium Sensitizer of Cardiac Troponin. Biophys J 2015. [DOI: 10.1016/j.bpj.2014.11.2074] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
|
122
|
Wang N, Zhou S, Kekenes-Huskey PM, Li B, McCammon JA. Poisson-Boltzmann versus Size-Modified Poisson-Boltzmann Electrostatics Applied to Lipid Bilayers. J Phys Chem B 2014; 118:14827-32. [PMID: 25426875 PMCID: PMC4280115 DOI: 10.1021/jp511702w] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
![]()
Mean-field
methods, such as the Poisson–Boltzmann equation
(PBE), are often used to calculate the electrostatic properties of
molecular systems. In the past two decades, an enhancement of the
PBE, the size-modified Poisson–Boltzmann equation (SMPBE),
has been reported. Here, the PBE and the SMPBE are reevaluated for
realistic molecular systems, namely, lipid bilayers, under eight different
sets of input parameters. The SMPBE appears to reproduce the molecular
dynamics simulation results better than the PBE only under specific
parameter sets, but in general, it performs no better than the Stern
layer correction of the PBE. These results emphasize the need for
careful discussions of the accuracy of mean-field calculations on
realistic systems with respect to the choice of parameters and call
for reconsideration of the cost-efficiency and the significance of
the current SMPBE formulation.
Collapse
|
123
|
Eun C, Kekenes-Huskey PM, Metzger VT, McCammon JA. A model study of sequential enzyme reactions and electrostatic channeling. J Chem Phys 2014; 140:105101. [PMID: 24628210 DOI: 10.1063/1.4867286] [Citation(s) in RCA: 42] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We study models of two sequential enzyme-catalyzed reactions as a basic functional building block for coupled biochemical networks. We investigate the influence of enzyme distributions and long-range molecular interactions on reaction kinetics, which have been exploited in biological systems to maximize metabolic efficiency and signaling effects. Specifically, we examine how the maximal rate of product generation in a series of sequential reactions is dependent on the enzyme distribution and the electrostatic composition of its participant enzymes and substrates. We find that close proximity between enzymes does not guarantee optimal reaction rates, as the benefit of decreasing enzyme separation is countered by the volume excluded by adjacent enzymes. We further quantify the extent to which the electrostatic potential increases the efficiency of transferring substrate between enzymes, which supports the existence of electrostatic channeling in nature. Here, a major finding is that the role of attractive electrostatic interactions in confining intermediate substrates in the vicinity of the enzymes can contribute more to net reactive throughput than the directional properties of the electrostatic fields. These findings shed light on the interplay of long-range interactions and enzyme distributions in coupled enzyme-catalyzed reactions, and their influence on signaling in biological systems.
Collapse
|
124
|
Tallorin LC, Durrant JD, Nguyen QG, McCammon JA, Burkart MD. Celastrol inhibits Plasmodium falciparum enoyl-acyl carrier protein reductase. Bioorg Med Chem 2014; 22:6053-6061. [PMID: 25284249 PMCID: PMC4807855 DOI: 10.1016/j.bmc.2014.09.002] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2014] [Revised: 08/25/2014] [Accepted: 09/02/2014] [Indexed: 01/29/2023]
Abstract
Enoyl-acyl carrier protein reductase (ENR), a critical enzyme in type II fatty acid biosynthesis, is a promising target for drug discovery against hepatocyte-stage Plasmodium falciparum. In order to identify PfENR-specific inhibitors, we docked 70 FDA-approved, bioactive, and/or natural product small molecules known to inhibit the growth of whole-cell blood-stage P. falciparum into several PfENR crystallographic structures. Subsequent in vitro activity assays identified a noncompetitive low-micromolar PfENR inhibitor, celastrol, from this set of compounds.
Collapse
|
125
|
Cheng Y, Lindert S, Kekenes-Huskey P, Rao VS, Solaro RJ, Rosevear PR, Amaro R, McCulloch AD, McCammon JA, Regnier M. Computational studies of the effect of the S23D/S24D troponin I mutation on cardiac troponin structural dynamics. Biophys J 2014; 107:1675-85. [PMID: 25296321 PMCID: PMC4190606 DOI: 10.1016/j.bpj.2014.08.008] [Citation(s) in RCA: 43] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2014] [Revised: 07/25/2014] [Accepted: 08/11/2014] [Indexed: 11/27/2022] Open
Abstract
During β-adrenergic stimulation, cardiac troponin I (cTnI) is phosphorylated by protein kinase A (PKA) at sites S23/S24, located at the N-terminus of cTnI. This phosphorylation has been shown to decrease KCa and pCa50, and weaken the cTnC-cTnI (C-I) interaction. We recently reported that phosphorylation results in an increase in the rate of early, slow phase of relaxation (kREL,slow) and a decrease in its duration (tREL,slow), which speeds up the overall relaxation. However, as the N-terminus of cTnI (residues 1-40) has not been resolved in the whole cardiac troponin (cTn) structure, little is known about the molecular-level behavior within the whole cTn complex upon phosphorylation of the S23/S24 residues of cTnI that results in these changes in function. In this study, we built up the cTn complex structure (including residues cTnC 1-161, cTnI 1-172, and cTnT 236-285) with the N-terminus of cTnI. We performed molecular-dynamics (MD) simulations to elucidate the structural basis of PKA phosphorylation-induced changes in cTn structure and Ca(2+) binding. We found that introducing two phosphomimic mutations into sites S23/S24 had no significant effect on the coordinating residues of Ca(2+) binding site II. However, the overall fluctuation of cTn was increased and the C-I interaction was altered relative to the wild-type model. The most significant changes involved interactions with the N-terminus of cTnI. Interestingly, the phosphomimic mutations led to the formation of intrasubunit interactions between the N-terminus and the inhibitory peptide of cTnI. This may result in altered interactions with cTnC and could explain the increased rate and decreased duration of slow-phase relaxation seen in myofibrils.
Collapse
|