Biomedical subjects

J Andrew McCammon

Publications and source records attributed to J Andrew McCammon.

At least 19 recordsLinked to original sources

PubMed2007

SR protein kinase 1 is resilient to inactivation.

SR protein kinase 1 (SRPK1) is a constitutively active kinase, which processively phosphorylates multiple serines within its substrates, ASF/SF2. We describe crystallographic, molecular dynamics, and biochemical results that shed light on how SRPK1 preserves its constitutive active conformation. Our structure reveals that unlike other known active kinase structures, the activation loop remains in an active state without any specific intraprotein interactions. Moreover, SRPK1 remains active despite extensive mutation to the activation segment. Molecular dynamics simulations reveal that SRPK1 partially absorbs the effect of mutations by forming compensatory interactions that maintain a catalytically competent chemical environment. Furthermore, SRPK1 is similarly resistant to deletion of its spacer loop region. Based upon a model of SRPK1 bound to a segment encompassing the docking motif and active-site peptide of ASF/SF2, we suggest a mechanism for processive phosphorylation and propose that the atypical resiliency we observed is critical for SRPK1's processive activity.

Crystallography, X-Ray↗

PubMed2006

Order N algorithm for computation of electrostatic interactions in biomolecular systems.

Poisson-Boltzmann electrostatics is a well established model in biophysics; however, its application to large-scale biomolecular processes such as protein-protein encounter is still limited by the efficiency and memory constraints of existing numerical techniques. In this article, we present an efficient and accurate scheme that incorporates recently developed numerical techniques to enhance our computational ability. In particular, a boundary integral equation approach is applied to discretize the linearized Poisson-Boltzmann equation; the resulting integral formulas are well conditioned and are extended to systems with arbitrary numbers of biomolecules. The solution process is accelerated by Krylov subspace methods and a new version of the fast multipole method. In addition to the electrostatic energy, fast calculations of the forces and torques are made possible by using an interpolation procedure. Numerical experiments show that the implemented algorithm is asymptotically optimal O(N) in both CPU time and required memory, and application to the acetylcholinesterase-fasciculin complex is illustrated.

Acetylcholinesterase↗

PubMed2006

On the application of accelerated molecular dynamics to liquid water simulations.

Our group recently proposed a robust bias potential function that can be used in an efficient all-atom accelerated molecular dynamics (MD) approach to simulate the transition of high energy barriers without any advance knowledge of the potential-energy landscape. The main idea is to modify the potential-energy surface by adding a bias, or boost, potential in regions close to the local minima, such that all transitions rates are increased. By applying the accelerated MD simulation method to liquid water, we observed that this new simulation technique accelerates the molecular motion without losing its microscopic structure and equilibrium properties. Our results showed that the application of a small boost energy on the potential-energy surface significantly reduces the statistical inefficiency of the simulation while keeping all the other calculated properties unchanged. On the other hand, although aggressive acceleration of the dynamics simulation increases the self-diffusion coefficient of water molecules greatly and dramatically reduces the correlation time of the simulation, configurations representative of the true structure of liquid water are poorly sampled. Our results also showed the strength and robustness of this simulation technique, which confirm this approach as a very useful and promising tool to extend the time scale of the all-atom simulations of biological system with explicit solvent models. However, we should keep in mind that there is a compromise between the strength of the boost applied in the simulation and the reproduction of the ensemble average properties.

Algorithms↗

PubMed2006

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.

Computational Biology↗

PubMed2006

In-situ synthesis of a tacrine-triazole-based inhibitor of acetylcholinesterase: configurational selection imposed by steric interactions.

Recently, researchers have used acetylcholinesterase (AChE) as a reaction vessel to synthesize its own inhibitors. Thus, 1 (syn-TZ2PA6), a femtomolar AChE inhibitor, which is formed in a 1:1 mixture with its anti-isomer by solution phase reaction from 3 (TZ2) and 4 (PA6), can be synthesized exclusively inside the AChE gorge. Our computational approach based on quantum mechanical/molecular mechanical (QM/MM) calculations, molecular dynamics (MD), and targeted molecular dynamics (TMD) studies answers why 1 is the sole product in the AChE environment. Ab initio QM/MM results show that the reaction in the AChE gorge occurs when 3/azide and 4/acetylene are extended in a parallel orientation. An MD simulation started from the final structure of QM/MM calculations keeps the azide's and acetylene's parallel orientations intact for 10 ns of simulation time. A TMD simulation applied on an antiparallel azide-acetylene conformation flips the acetylene easily to bring it to a position that is parallel to azide. A second set of QM/MM calculations performed on this flipped structure generates a similar minimum-energy path as obtained previously. Even a TMD simulation carried out on a parallel azide-acetylene conformation could not deform their parallel arrangement. All of these results, thus, imply that inside the AChE gorge, the azide group of 3 and the acetylene group of 4 always remain parallel, with the consequence that 1 is the only product. The architecture of the gorge plays an important role in this selective formation of 1.

Acetylcholinesterase↗

PubMed2006

Protein complex formation by acetylcholinesterase and the neurotoxin fasciculin-2 appears to involve an induced-fit mechanism.

Specific, rapid association of protein complexes is essential for all forms of cellular existence. The initial association of two molecules in diffusion-controlled reactions is often influenced by the electrostatic potential. Yet, the detailed binding mechanisms of proteins highly depend on the particular system. A complete protein complex formation pathway has been delineated by using structural information sampled over the course of the transformation reaction. The pathway begins at an encounter complex that is formed by one of the apo forms of neurotoxin fasciculin-2 (FAS2) and its high-affinity binding protein, acetylcholinesterase (AChE), followed by rapid conformational rearrangements into an intermediate complex that subsequently converts to the final complex as observed in crystal structures. Formation of the intermediate complex has also been independently captured in a separate 20-ns molecular dynamics simulation of the encounter complex. Conformational transitions between the apo and liganded states of FAS2 in the presence and absence of AChE are described in terms of their relative free energy profiles that link these two states. The transitions of FAS2 after binding to AChE are significantly faster than in the absence of AChE; the energy barrier between the two conformational states is reduced by half. Conformational rearrangements of FAS2 to the final liganded form not only bring the FAS2/AChE complex to lower energy states, but by controlling transient motions that lead to opening or closing one of the alternative passages to the active site of the enzyme also maximize the ligand's inhibition of the enzyme.

Acetylcholinesterase↗

PubMed2006

Binding of aminoglycosidic antibiotics to the oligonucleotide A-site model and 30S ribosomal subunit: Poisson-Boltzmann model, thermal denaturation, and fluorescence studies.

The binding of paromomycin and similar antibiotics to the oligonucleotide A-site model and the small (30S) ribosomal subunit has been studied using continuum electrostatics methods. Crystallographic information from complexes of paromomycin, tobramycin, and Geneticin bound to an A-site oligonucleotide, and paromomycin and streptomycin complexed to the 30S subunit was used as a foundation to develop structures of similar antibiotics in the same ribosomal binding site. Relative binding free energies were calculated by combining the electrostatic contribution, which was obtained by solving the Poisson-Boltzmann equation, with a surface-area-dependent apolar term and contributions from conformational changes. These computed results showed good correlation with the experimental data resulting from fluorescence binding assays and thermal denaturation studies, demonstrating the ability of the Poisson-Boltzmann model to provide insight into the electrostatic mechanisms for aminoglycoside binding and direction for designing more effective antibiotics.

Aminoglycosides↗

PubMed2006

Insight into the role of hydration on protein dynamics.

The potential energy surface of a protein is rough. This intrinsic energetic roughness affects diffusion, and hence the kinetics. The dynamics of a system undergoing Brownian motion on this surface in an implicit continuum solvent simulation can be tuned via the frictional drag or collision frequency to be comparable to that of experiments or explicit solvent simulations. We show that the kinetic rate constant for a local rotational isomerization in stochastic simulations with continuum solvent and a collision frequency of 2 ps(-1) is about 10(4) times faster than that in explicit water and experiments. A further increase in the collision frequency to 60 ps(-1) slows down the dynamics, but does not fully compensate for the lack of explicit water. We also show that the addition of explicit water does not only slow down the dynamics by increasing the frictional drag, but also increases the local energetic roughness of the energy landscape by as much as 1.0 kcal/mol.

Computer Simulation↗

PubMed2006

Bio3d: an R package for the comparative analysis of protein structures.

UNLABELLED: An automated procedure for the analysis of homologous protein structures has been developed. The method facilitates the characterization of internal conformational differences and inter-conformer relationships and provides a framework for the analysis of protein structural evolution. The method is implemented in bio3d, an R package for the exploratory analysis of structure and sequence data. AVAILABILITY: The bio3d package is distributed with full source code as a platform-independent R package under a GPL2 license from: http://mccammon.ucsd.edu/~bgrant/bio3d/

Amino Acid Sequence↗

PubMed2006

Flap opening dynamics in HIV-1 protease explored with a coarse-grained model.

We present a one-bead coarse-grained model that enables dynamical simulations of proteins on the time scale of tens of microseconds. The parameterization of the force field includes accurate conformational terms that allow for fast and reliable exploration of the configurational space. The model is applied to the dynamics of flap opening in HIV-1 protease. The experimental structure of the recently crystallized semi-open conformation of HIV-1 protease is well reproduced in the simulation, which supports the accuracy of our model. Thanks to very long simulations and extensive sampling of opening and closing events, we also investigate the thermodynamics and kinetics of the opening process. We have shown that the effect of the solvent slows down the dynamics to the experimentally observed time scales. The model is found to be reliable for application to substrate docking simulations, which are currently in progress.

HIV Protease↗

PubMed2006

Targeted molecular dynamics study of C-loop closure and channel gating in nicotinic receptors.

The initial coupling between ligand binding and channel gating in the human alpha7 nicotinic acetylcholine receptor (nAChR) has been investigated with targeted molecular dynamics (TMD) simulation. During the simulation, eight residues at the tip of the C-loop in two alternating subunits were forced to move toward a ligand-bound conformation as captured in the crystallographic structure of acetylcholine binding protein (AChBP) in complex with carbamoylcholine. Comparison of apo- and ligand-bound AChBP structures shows only minor rearrangements distal from the ligand-binding site. In contrast, comparison of apo and TMD simulation structures of the nAChR reveals significant changes toward the bottom of the ligand-binding domain. These structural rearrangements are subsequently translated to the pore domain, leading to a partly open channel within 4 ns of TMD simulation. Furthermore, we confirmed that two highly conserved residue pairs, one located near the ligand-binding pocket (Lys145 and Tyr188), and the other located toward the bottom of the ligand-binding domain (Arg206 and Glu45), are likely to play important roles in coupling agonist binding to channel gating. Overall, our simulations suggest that gating movements of the alpha7 receptor may involve relatively small structural changes within the ligand-binding domain, implying that the gating transition is energy-efficient and can be easily modulated by agonist binding/unbinding.

Amino Acid Sequence↗

PubMed2006

Restrained molecular dynamics simulations of HIV-1 protease: the first step in validating a new target for drug design.

To test the anticorrelated relationship that was recently displayed in conventional molecular dynamics (MD) simulations, several different restrained MD simulations on a wild type and on the V82F/I84V drug-resistant mutant of HIV-1 protease were performed. This anticorrelated relationship refers to the observation that compression of the peripheral ear-to-cheek region of HIV protease (i.e., the elbow of the flap to the fulcrum and the cantilever) occurred as the active site flaps were opening, and, conversely, expansion of that ear-to-cheek region occurred as both flaps were closing. An additional examination of this anticorrelated relationship was necessary to determine whether it can be harnessed in a useful manner. Consequently, six different MD experiments were performed that incorporated pairwise distance restraints in that ear-to-cheek region (i.e., the distance between the alpha-carbons of Gly40 and Gln61 was restrained to either 7.7 or 10.5 A, in both monomers). Pushing the backbones of the ear and the cheek regions away from each other slightly did force the flaps that guard the active site to remain closed in both the wild type and the mutant systems-even though there were no ligands in the active sites. Thus, these restrained MD simulations provided evidence that the anticorrelated relationship can be exploited to affect the dynamic behavior of the flaps that guard the active site of HIV-1 protease. These simulations supported our hypothesis of the mechanism governing flap motion, and they are the first step towards validating that peripheral surface as a new target for drug design.

Binding Sites↗

PubMed2006

The influence of macromolecular crowding on HIV-1 protease internal dynamics.

High macromolecular concentrations, or crowded conditions, have been shown to affect a wide variety of molecular processes, including diffusion, association and dissociation, and protein folding and stability. Here, we model the effect of macromolecular crowding on the internal dynamics of a protein, HIV-1 protease, using Brownian dynamics simulations. HIV-1 protease possesses a pair of flaps which are postulated to open in the early stages of its catalytic mechanism. Compared to low concentrations, close-packed concentrations of repulsive crowding agents are found to significantly reduce the fraction of time that the protease flaps are open. Macromolecular crowding is likely to have a major effect on in vivo enzyme activity, and may play an important regulatory role in the viral life cycle.

Amino Acids↗

PubMed2006

Elasticity of peptide omega bonds.

We calculated the changes of the free energy profile of the peptidyl-prolyl torsional angle of the dipeptide valine-proline under pulling forces by simulations. Using a dynamic model built on the equilibrium properties of this system and previously studied dynamic properties of cis-trans isomerization of other dipeptides, we calculated the dynamic viscoelasticity of this degree of freedom. The results show significant differences between how thermal and mechanical forces alter the equilibrium and the dynamics of the isomerization transition. The former does not change the barrier heights but changes the prefactor of the kinetics owing to temperature effects, while the latter changes minima and thus the population. The force that is required to "excite" this degree of freedom is small. Compared to other systems, we found that this degree of freedom becomes already quite rigid at several hertz, which is a much lower value due to the high barrier of the cis-trans isomerization. We also found that the tensile elastic modulus of densely packed omega bonds is at the order of GPa, which is comparable to that of polymer materials. These results give mechanical properties of polyproline elasticity of a local nature and provide guidance for future experimental designs.

Binding Sites↗

PubMed2006

Gated binding of ligands to HIV-1 protease: Brownian dynamics simulations in a coarse-grained model.

The internal motions of proteins may serve as a "gate" in some systems, which controls ligand-protein association. This study applies Brownian dynamics simulations in a coarse-grained model to study the gated association rate constants of HIV-1 proteases and drugs. The computed gated association rate constants of three protease mutants, G48V/V82A/I84V/L90M, G48V, and L90M with three drugs, amprenavir, indinavir, and saquinavir, yield good agreements with experiments. The work shows that the flap dynamics leads to "slow gating". The simulations suggest that the flap flexibility and the opening frequency of the wild-type, the G48V and L90M mutants are similar, but the flaps of the variant G48V/V82A/I84V/L90M open less frequently, resulting in a lower gated rate constant. The developed methodology is fast and provides an efficient way to predict the gated association rate constants for various protease mutants and ligands.

Carbamates↗

PubMed2006

Potentials of mean force for acetylcholine unbinding from the alpha7 nicotinic acetylcholine receptor ligand-binding domain.

The nicotinic acetylcholine receptor is a prototype ligand-gated ion channel that mediates signal transduction in the neuromuscular junctions and other cholinergic synapses. The molecular basis for the energetics of ligand binding and unbinding is critical to our understanding of the pharmacology of this class of receptors. Here, we used steered molecular dynamics to investigate the unbinding of acetylcholine from the ligand-binding domain of human alpha7 nicotinic acetylcholine receptor along four different predetermined pathways. Pulling forces were found to correlate well with interactions between acetylcholine and residues in the binding site during the unbinding process. From multiple trajectories along these unbinding pathways, we calculated the potentials of mean force for acetylcholine unbinding. Four available methods based on Jarzynski's equality were used and compared for their efficiencies. The most probable pathway was identified to be along a direction approximately parallel to the membrane. The derived binding energy for acetylcholine was in good agreement with that derived from the experimental binding constant for acetylcholine binding protein, but significantly higher than that for the complete human alpha7 nicotinic acetylcholine receptor. In addition, it is likely that several intermediate states exist along the unbinding pathways.

Acetylcholine↗

PubMed2006

How does activation loop phosphorylation modulate catalytic activity in the cAMP-dependent protein kinase: a theoretical study.

Phosphorylation mediates the function of many proteins and enzymes. In the catalytic subunit of cAMP-dependent protein kinase, phosphorylation of Thr 197 in the activation loop strongly influences its catalytic activity. In order to provide theoretical understanding about this important regulatory process, classical molecular dynamics simulations and ab initio QM/MM calculations have been carried out on the wild-type PKA-Mg(2) ATP-substrate complex and its dephosphorylated mutant, T197A. It was found that pThr 197 not only facilitates the phosphoryl transfer reaction by stabilizing the transition state through electrostatic interactions but also strongly affects its essential protein dynamics as well as the active site conformation.

Catalysis↗

PubMed2006

Dependency map of proteins in the small ribosomal subunit.

The assembly of the ribosome has recently become an interesting target for antibiotics in several bacteria. In this work, we extended an analytical procedure to determine native state fluctuations and contact breaking to investigate the protein stability dependence in the 30S small ribosomal subunit of Thermus thermophilus. We determined the causal influence of the presence and absence of proteins in the 30S complex on the binding free energies of other proteins. The predicted dependencies are in overall agreement with the experimentally determined assembly map for another organism, Escherichia coli. We found that the causal influences result from two distinct mechanisms: one is pure internal energy change, the other originates from the entropy change. We discuss the implications on how to target the ribosomal assembly most effectively by suggesting six proteins as targets for mutations or other hindering of their binding. Our results show that by blocking one out of this set of proteins, the association of other proteins is eventually reduced, thus reducing the translation efficiency even more. We could additionally determine the binding dependency of THX--a peptide not present in the ribosome of E. coli--and suggest its assembly path.

Computational Biology↗