Tetrameric mouse acetylcholinesterase: continuum diffusion rate calculations by solving the steady-state Smoluchowski equation using finite element methods

The tetramer is the most important form for acetylcholinesterase in physiological conditions, i.e., in the neuromuscular junction and the nervous system. It is important to study the diffusion of acetylcholine to the active sites of the tetrameric enzyme to understand the overall signal transduction process in these cellular components. Crystallographic studies revealed two different forms of tetramers, suggesting a flexible tetramer model for acetylcholinesterase. Using a recently developed finite element solver for the steady-state Smoluchowski equation, we have calculated the reaction rate for three mouse acetylcholinesterase tetramers using these two crystal structures and an intermediate structure as templates. Our results show that the reaction rates differ for different individual active sites in the compact tetramer crystal structure, and the rates are similar for different individual active sites in the other crystal structure and the intermediate structure. In the limit of zero salt, the reaction rates per active site for the tetramers are the same as that for the monomer, whereas at higher ionic strength, the rates per active site for the tetramers are approximately 67%-75% of the rate for the monomer. By analyzing the effect of electrostatic forces on ACh diffusion, we find that electrostatic forces play an even more important role for the tetramers than for the monomer. This study also shows that the finite element solver is well suited for solving the diffusion problem within complicated geometries.     

Ligand-induced conformational change in the alpha7 nicotinic receptor ligand binding domain

Molecular dynamics simulations of a homology model of the ligand binding domain of the alpha7 nicotinic receptor are conducted with a range of bound ligands to induce different conformational states. Four simulations of 15 ns each are run with no ligand, antagonist d-tubocurarine (dTC), agonist acetylcholine (ACh), and agonist ACh with potentiator Ca(2+), to give insight into the conformations of the active and inactive states of the receptor and suggest the mechanism for conformational change. The main structural factor distinguishing the active and inactive states is that a more open, symmetric arrangement of the five subunits arises for the two agonist simulations, whereas a more closed and asymmetric arrangement results for the apo and dTC cases. Most of the difference arises in the lower portion of the ligand binding domain near its connection to the adjacent transmembrane domain. The transfer of the more open state to the transmembrane domain could then promote ion flow through the channel. Variation in how subunits pack together with no ligand bound appears to give rise to asymmetry in the apo case. The presence of dTC expands the receptor but induces rotations in alternate directions in adjacent subunits that lead to an asymmetric arrangement as in the apo case. Ca(2+) appears to promote a slightly greater expansion in the subunits than ACh alone by stabilizing the C-loop and ACh positions. Although the simulations are unlikely to be long enough to view the full conformational changes between open and closed states, a collection of different motions at a range of length scales are observed that are likely to participate in the conformational change.     

Finite element solution of the steady-state Smoluchowski equation for rate constant calculations

This article describes the development and implementation of algorithms to study diffusion in biomolecular systems using continuum mechanics equations. Specifically, finite element methods have been developed to solve the steady-state Smoluchowski equation to calculate ligand binding rate constants for large biomolecules. The resulting software has been validated and applied to mouse acetylcholinesterase. Rates for inhibitor binding to mAChE were calculated at various ionic strengths with several different reaction criteria. The calculated rates were compared with experimental data and show very good agreement when the correct reaction criterion is used. Additionally, these finite element methods require significantly less computational resources than existing particle-based Brownian dynamics methods.     

Proliferating cell nuclear antigen loaded onto double-stranded DNA: dynamics, minor groove interactions and functional implications

Proliferating cell nuclear antigen (PCNA) acts as a biologically essential processivity factor that encircles DNA and provides binding sites for polymerase, flap endonuclease-1 (FEN-1) and ligase during DNA replication and repair. We have computationally characterized the interactions of human and Archaeoglobus fulgidus PCNA trimer with double-stranded DNA (ds DNA) using multi-nanosecond classical molecular dynamics simulations. The results reveal the interactions of DNA passing through the PCNA trimeric ring including the contacts formed, overall orientation and motion with respect to the sliding clamp. Notably, we observe pronounced tilting of the axis of dsDNA with respect to the PCNA ring plane reflecting interactions between the DNA phosphodiester backbone and positively charged arginine and lysine residues lining the PCNA inner surface. Covariance matrix analysis revealed a pattern of correlated motions within and between the three equivalent subunits involving the PCNA C-terminal region and linker strand associated with partner protein binding sites. Additionally, principal component analysis identified low frequency global PCNA subunit motions suitable for translocation along duplex DNA. The PCNA motions and interactions with the DNA minor groove, identified here computationally, provide an unexpected basis for PCNA to act in the coordinated handoff of intermediates from polymerase to FEN-1 to ligase during DNA replication and repair.     

E230Q mutation of the catalytic subunit of cAMP-dependent protein kinase affects local structure and the binding of peptide inhibitor

The active site of the mammalian cAMP-dependent protein kinase catalytic subunit (C-subunit) has a cluster of nonconserved acidic residues-Glu127, Glu170, Glu203, Glu230, and Asp241-that are crucial for substrate recognition and binding. Studies have shown that the Glu230 to Gln mutant (E230Q) of the enzyme has physical properties similar to the wild-type enzyme and has decreased affinity for a short peptide substrate, Kemptide. However, recent experiments intended to crystallize ternary complex of the E230Q mutant with MgATP and protein kinase inhibitor (PKI) could only obtain crystals of the apo-enzyme of E230Q mutant. To deduce the possible mechanism that prevented ternary complex formation, we used the relaxed-complex method (Lin, J.-H., et al. J Am Chem Soc 2002, 24, 5632-5633) to study PKI binding to the E230Q mutant C-subunit. In the E230Q mutant, we observed local structural changes of the peptide binding site that correlated closely to the reduced PKI affinity. The structural changes occurred in the F-to-G helix loop and appeared to hinder PKI binding. Reduced electrostatic potential repulsion among Asp241 from the helix loop section and the other acidic residues in the peptide binding site appear to be responsible for the structural change.     

Sensitive fluorescence in situ hybridization signal detection in maize using directly labeled probes produced by high concentration DNA polymerase nick translation

The signal produced by fluorescence in situ hybridization (FISH) often is inconsistent among cells and sensitivity is low. Small DNA targets on the chromatin are difficult to detect. We report here an improved nick translation procedure for Texas red and Alexa Fluor 488 direct labeling of FISH probes. Brighter probes can be obtained by adding excess DNA polymerase I. Using such probes, a 30 kb yeast transgene, and the rp1, rp3 and zein multigene clusters were clearly detected.     

Simulations of the p97 complex suggest novel conformational states of hydrolysis intermediates

The vitally important AAA (ATPases associated with various cellular activities) protein p97 is involved in cellular functions ranging from replication to degradation of misfolded proteins and has recently been proposed as a novel chemotherapeutic target. p97 is a large molecular machine that has been shown to hexamerize in vitro, with each monomer consisting of an N domain responsible for binding to effector proteins and two AAA repeats (D1 and D2). However, structural studies are inconclusive or in disagreement with one another on several important features such as the locations of the N domains, the relative orientations of the D1 and D2 rings, and the dimensions of the central pore. Here, we present atomic-scale simulations of the p97 hexamer in the prehydrolysis, transition, and post-hydrolysis states. To improve the agreement between low- and high-resolution experimental studies, we first use a biased simulation technique, molecular dynamics flexible fitting (MDFF), to improve the correlation between the structures described in these experiments. We follow this with extended, classical molecular dynamics simulations, which not only show that structures generated in the MDFF phase are stable, but reveal insights into the dynamics important to each state. Simulation results suggest a hybrid model for hydrolysis, in which the N and D2 domains are dynamic while the D1 domains are relatively static, salt bridges stabilize the position of the N domains in the pre-hydrolysis state, and the rings formed by D1 and D2 rotate relative to one another.     

Identification of potential small molecule binding pockets on Rho family GTPases

Rho GTPases are conformational switches that control a wide variety of signaling pathways critical for eukaryotic cell development and proliferation. They represent attractive targets for drug design as their aberrant function and deregulated activity is associated with many human diseases including cancer. Extensive high-resolution structures (>100) and recent mutagenesis studies have laid the foundation for the design of new structure-based chemotherapeutic strategies. Although the inhibition of Rho signaling with drug-like compounds is an active area of current research, very little attention has been devoted to directly inhibiting Rho by targeting potential allosteric non-nucleotide binding sites. By avoiding the nucleotide binding site, compounds may minimize the potential for undesirable off-target interactions with other ubiquitous GTP and ATP binding proteins. Here we describe the application of molecular dynamics simulations, principal component analysis, sequence conservation analysis, and ensemble small-molecule fragment mapping to provide an extensive mapping of potential small-molecule binding pockets on Rho family members. Characterized sites include novel pockets in the vicinity of the conformationaly responsive switch regions as well as distal sites that appear to be related to the conformations of the nucleotide binding region. Furthermore the use of accelerated molecular dynamics simulation, an advanced sampling method that extends the accessible time-scale of conventional simulations, is found to enhance the characterization of novel binding sites when conformational changes are important for the protein mechanism.     

Inhibition of cathepsin B by Au(I) complexes: a kinetic and computational study

Gold(I) compounds have been used in the treatment of rheumatoid arthritis for over 80 years, but the biological targets and the structure–activity relationships of these drugs are not well understood. Of particular interest is the molecular mechanism behind the antiarthritic activity of the orally available drug triethylphosphine(2,3,4,6-tetra-O-acetyl-β-1-d-thiopyranosato-S) gold(I) (auranofin, Ridaura). The cathepsin family of lysosomal, cysteine-dependent enzymes is an attractive biological target of Au(I) and is inhibited by auranofin and auranofin analogs with reasonable potency. Here we employ a combination of experimental and computational investigations into the effect of changes in the phosphine ligand of auranofin on its in vitro inhibition of cathepsin B. Sequential replacement of the ethyl substituents of triethylphosphine by phenyl groups leads to increasing potency in the resultant Au(I) complexes, due in large part to favorable interactions of the more sterically bulky Au(I)–PR₃ fragments with the enzyme active site.     

Hot-spot residues at the E9/Im9 interface help binding via different mechanisms

Protein-protein association involves many interface interactions, but they do not contribute equally. Ala scanning experiments reveal that only a few mutations significantly lower binding affinity. These key residues, which appear to drive protein-protein association, are called hot-spot residues. Molecular dynamics simulations of the Colicin E9/Im9 complex show Im9 Glu41 and Im9 Ser50, both hot-spots, bind via different mechanisms. The results suggest that Im9 Ser50 restricts Glu41 in a conformation auspicious for salt-bridge formation across the interface. This type of model may be helpful in engineering hot-spot clusters at protein-protein interfaces and, consequently, the design of specificity. (c) 2008 Wiley Periodicals, Inc. Biopolymers 89: 916-920, 2008.This article was originally published online as an accepted preprint. The "Published Online" date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com.