Biomedical subjects

James M Briggs

Publications and source records attributed to James M Briggs.

At least 19 recordsLinked to original sources

PubMed2006

Computed pore potentials of the nicotinic acetylcholine receptor.

Electrostatic surface potentials in the vestibule of the nicotinic acetylcholine receptor (nAChR) were computed from structural models using the University of Houston Brownian Dynamics program to determine their effect on ion conduction and ionic selectivity. To further determine whether computed potentials accurately reflect the electrostatic environment of the channel, the potentials were used to predict the rate constants for diffusion-enhanced fluorescence energy transfer; the calculated energy transfer rates are directly comparable with those determined experimentally (see companion article by Meltzer et al. in this issue). To include any effects on the local potentials by the bound acceptor fluorophore crystal violet, its binding site was first localized within the pore by fluorescence energy transfer measurements from dansyl-C6-choline bound to the agonist sites and also by simulations of binding using Autodock. To compare the computed potentials with those determined experimentally, we used the predicted energy transfer rates from Tb3+ chelates of varying charge to calculate an expected potential using the Boltzmann relationship. This expected potential (from -20 to -40 mV) overestimates the values determined experimentally (from -10 to -25 mV) by two- to fourfold at similar conditions of ionic strength. Although the results indicate a basic discrepancy between experimental and computed surface potentials, both methods demonstrate that the vestibular potential has a relatively small effect on conduction and selectivity.

Cell Membrane↗

PubMed2006

Electrostatic steering at acetylcholine binding sites.

The electrostatic environments near the acetylcholine binding sites on the nicotinic acetylcholine receptor (nAChR) and acetylcholinesterase were measured by diffusion-enhanced fluorescence energy transfer (DEFET) to determine the influence of long-range electrostatic interactions on ligand binding kinetics and net binding energy. Changes in DEFET from variously charged Tb3+ -chelates revealed net potentials of -20 mV at the nAChR agonist sites and -14 mV at the entrance to the AChE active site, in physiological ionic strength conditions. The potential at the alphadelta-binding site of the nAChR was determined independently in the presence of d-tubocurarine to be -14 mV; the calculated potential at the alphagamma-site was approximately threefold stronger than at the alphadelta-site. By determining the local potential in increasing ionic strength, Debye-Hückel theory predicted that the potentials near the nAChR agonist binding sites are constituted by one to three charges in close proximity to the binding site. Examination of the binding kinetics of the fluorescent acetylcholine analog dansyl-C6-choline at ionic strengths from 12.5 to 400 mM revealed a twofold decrease in association rate. Debye-Hückel analysis of the kinetics revealed a similar charge distribution as seen by changes in the potentials. To determine whether the experimentally determined potentials are reflected by continuum electrostatics calculations, solutions to the nonlinear Poisson-Boltzmann equation were used to compute the potentials expected from DEFET measurements from high-resolution models of the nAChR and AChE. These calculations are in good agreement with the DEFET measurements for AChE and for the alphagamma-site of the nAChR. We conclude that long-range electrostatic interactions contribute -0.3 and -1 kcal/mol to the binding energy at the nAChR alphadelta- and alphagamma-sites due to an increase in association rates.

Acetylcholine↗

PubMed2006

Dynamic pharmacophore model optimization: identification of novel HIV-1 integrase inhibitors.

We extended the previously described dynamic pharmacophore model studies of HIV-1 integrase (IN) by considering more key residues in the active site, including Mg2+. First, we applied a Monte Carlo sampling method to map the complementary features of the IN binding surface. Two types of dynamic pharmacophore models were generated. One considers Mg2+ as part of the IN and therefore as an excluded volume, and the other treats Mg2+ as a positively charged feature, representing a new type of pharmacophore model aimed to identify compounds potentially preventing Mg2+ binding. Second, we validated the models with 385 known active (IC50 < 20 microM) and 235 (IC50 > 100 microM) inactive IN inhibitors. Third, we used the derived models to screen our small molecule database. Twenty-two structurally novel compounds were tested in an in vitro assay specific for IN, and two of them showed IC50 < or = 10 microM for strand transfer reaction.

Algorithms↗

PubMed2005

Nucleotide modification at the gamma-phosphate leads to the improved fidelity of HIV-1 reverse transcriptase.

The mechanism by which HIV-1 reverse transcriptase (HIV-RT) discriminates between the correct and incorrect nucleotide is not clearly understood. Chemically modified nucleotides containing 1-aminonaphthalene-5-sulfonate (ANS) attached to their gamma-phosphate were synthesized and used to probe nucleotide selection by this error prone polymerase. Primer extension reactions provide direct evidence that the polymerase is able to incorporate the gamma-modified nucleotides. Forward mutation assays reveal a 6-fold reduction in the mutational frequency with the modified nucleotides, and specific base substitutions are dramatically reduced or eliminated. Molecular modeling illustrates potential interactions between critical residues within the polymerase active site and the modified nucleotides. Our data demonstrate that the fidelity of reverse transcriptase is improved using modified nucleotides, and we suggest that specific modifications to the gamma-phosphate may be useful in designing new antiviral therapeutics or, more generally, as a tool for defining the structural role that the polymerase active site has on nucleotide selectivity.

Binding Sites↗

PubMed2005

Comparative molecular dynamics simulations of HIV-1 integrase and the T66I/M154I mutant: binding modes and drug resistance to a diketo acid inhibitor.

HIV-1 IN is an essential enzyme for viral replication and an interesting target for the design of new pharmaceuticals for use in multidrug therapy of AIDS. L-731,988 is one of the most active molecules of the class of beta-diketo acids. Individual and combined mutations of HIV-1 IN at residues T66, S153, and M154 confer important degrees of resistance to one or more inhibitors belonging to this class. In an effort to understand the molecular mechanism of the resistance of T66I/M154I IN to the inhibitor L-731,988 and its specific binding modes, we have carried out docking studies, explicit solvent MD simulations, and binding free energy calculations. The inhibitor was docked against different protein conformations chosen from prior MD trajectories, resulting in 2 major orientations within the active site. MD simulations have been carried out for the T66I/M154I DM IN, DM IN in complex with L-731,988 in 2 different orientations, and 1QS4 IN in complex with L-731,988. The results of these simulations show a similar dynamical behavior between T66I/M154I IN alone and in complex with L-731,988, while significant differences are observed in the mobility of the IN catalytic loop (residues 138-149). Water molecules bridging the inhibitor to residues from the active site have been identified, and residue Gln62 has been found to play an important role in the interactions between the inhibitor and the protein. This work provides information about the binding modes of L-731,988, as well as insight into the mechanism of inhibitor-resistance in HIV-1 integrase.

Amino Acid Substitution↗

PubMed2005

Comparison of multiple molecular dynamics trajectories calculated for the drug-resistant HIV-1 integrase T66I/M154I catalytic domain.

HIV-1 integrase (IN) is an essential enzyme for the viral replication and an interesting target for the design of new pharmaceuticals for multidrug therapy of AIDS. Single and multiple mutations of IN at residues T66, S153, or M154 confer degrees of resistance to several inhibitors that prevent the enzyme from performing its normal strand transfer activity. Four different conformations of IN were chosen from a prior molecular dynamics (MD) simulation on the modeled IN T66I/M154I catalytic core domain as starting points for additional MD studies. The aim of this article is to understand the dynamic features that may play roles in the catalytic activity of the double mutant enzyme in the absence of any inhibitor. Moreover, we want to verify the influence of using different starting points on the MD trajectories and associated dynamical properties. By comparison of the trajectories obtained from these MD simulations we have demonstrated that the starting point does not affect the conformational space explored by this protein and that the time of the simulation is long enough to achieve convergence for this system.

Biophysics↗

PubMed2005

Dynamic receptor-based pharmacophore model development and its application in designing novel HIV-1 integrase inhibitors.

We present here a dynamic receptor-based pharmacophore model representing the complementary features of the active site region of HIV-1 integrase (IN), which was developed from a series of representative conformations of IN. Conformations of IN were sampled through a molecular dynamics study of the catalytic domain of an IN monomer, and an ensemble of representative IN structures were collected via a probability-based representative conformer sampling method that considers both the potential energy and the structural similarity of the protein conformations. The dynamic pharmacophore model was validated by a set of 128 known inhibitors, and the results showed that over 72% of the active inhibitors (IC(50) lower than 20 microM) could be successfully identified by the dynamic model. Therefore, we screened our in-house database of commercially available compounds against this model and successfully identified a set of structurally novel IN inhibitors. Compounds 7 and 18 with IC(50)s of 8 microM and 15 microM, respectively, against the strand transfer reaction were the most potent. Moreover, 7, 8 and 20 showed a 5-fold selectivity for the strand transfer reaction over 3'-processing.

Binding Sites↗

PubMed2005

Large-scale conformational dynamics of the HIV-1 integrase core domain and its catalytic loop mutants.

HIV-1 integrase is one of the three essential enzymes required for viral replication and has great potential as a novel target for anti-HIV drugs. Although tremendous efforts have been devoted to understanding this protein, the conformation of the catalytic core domain around the active site, particularly the catalytic loop overhanging the active site, is still not well characterized by experimental methods due to its high degree of flexibility. Recent studies have suggested that this conformational dynamics is directly correlated with enzymatic activity, but the details of this dynamics is not known. In this study, we conducted a series of extended-time molecular dynamics simulations and locally enhanced sampling simulations of the wild-type and three loop hinge mutants to investigate the conformational dynamics of the core domain. A combined total of >480 ns of simulation data was collected which allowed us to study the conformational changes that were not possible to observe in the previously reported short-time molecular dynamics simulations. Among the main findings are a major conformational change (>20 A) in the catalytic loop, which revealed a gatinglike dynamics, and a transient intraloop structure, which provided a rationale for the mutational effects of several residues on the loop including Q(148), P(145), and Y(143). Further, clustering analyses have identified seven major conformational states of the wild-type catalytic loop. Their implications for catalytic function and ligand interaction are discussed. The findings reported here provide a detailed view of the active site conformational dynamics and should be useful for structure-based inhibitor design for integrase.

Binding Sites↗

PubMed2005

The 1.9 A crystal structure of alanine racemase from Mycobacterium tuberculosis contains a conserved entryway into the active site.

We report the crystal structure of alanine racemase from Mycobacterium tuberculosis (Alr(Mtb)) at 1.9 A resolution. In our structure, Alr(Mtb) is found to be a dimer formed by two crystallographically different monomers, each comprising 384 residues. The domain makeup of each monomer is similar to that of Bacillus and Pseudomonas alanine racemases and includes both an alpha/beta-barrel at the N-terminus and a C-terminus primarily made of beta-strands. The hinge angle between these two domains is unique for Alr(Mtb), but the active site geometry is conserved. In Alr(Mtb), the PLP cofactor is covalently bound to the protein via an internal aldimine bond with Lys42. No guest substrate is noted in its active site, although some residual electron density is observed in the enzyme's active site pocket. Analysis of the active site pocket, in the context of other known alanine racemases, allows us to propose the inclusion of conserved residues found at the entrance to the binding pocket as additional targets in ongoing structure-aided drug design efforts. Also, as observed in other alanine racemase structures, PLP adopts a conformation that significantly distorts the planarity of the extended conjugated system between the PLP ring and the internal aldimine bond.

Alanine Racemase↗

PubMed2004

Cluster analysis of water molecules in alanine racemase and their putative structural role.

Conservation of water molecules was identified by a cluster analysis of seven crystal structures of alanine racemase from Bacillus stearothermophilus. A total of 47 clusters of consensus water sites were determined and found to be highly localized, as indicated by their low mobilities. These clusters are located in the region of the active sites as well as at the interface between the N-terminal domain (the alpha/beta-barrel) of the first monomer and the C-terminal domain of the second monomer. The clusters located at the dimer interface form extensive hydrogen-bonding networks linked to the protein backbone. These water-mediated hydrogen bonds, and also all hydrogen-bonding interactions at the dimer interface, were monitored during a 2 ns molecular dynamics simulation and showed that when the inhibitor propionate was bound to the enzyme, some of these interactions were disrupted. The data we present here indicate that the consensus water sites identified at the interface between the two monomers of alanine racemase may play a structural role, which is to maintain and stabilize the alanine racemase dimer. A second role might be to supply the active site continuously with water molecules in order to allow rapid equilibration of active site protons with the solvent.

Alanine↗

PubMed2004

HIV-1 integrase pharmacophore model derived from diverse classes of inhibitors.

A three-dimensional pharmacophore model has been generated for HIV-1 integrase (HIV-1 IN) from known inhibitors. A dataset consisting of 26 inhibitors was selected on the basis of the information content of the structures and activity data as required by the catalyst/HypoGen program. Our model was able to predict the activity of other known HIV-1 IN inhibitors not included in the model generation, and can be further used to identify structurally diverse compounds with desired biological activity by virtual screening.

HIV Integrase↗

PubMed2004

Molecular modeling study of the editing active site of Escherichia coli leucyl-tRNA synthetase: two amino acid binding sites in the editing domain.

Aminoacyl-tRNA synthetases (aaRSs) strictly discriminate their cognate amino acids. Some aaRSs accomplish this via proofreading and editing mechanisms. Mursinna and coworkers recently reported that substituting a highly conserved threonine (T252) with an alanine within the editing domain of Escherichia coli leucyl-tRNA synthetase (LeuRS) caused LeuRS to cleave its cognate aminoacylated leucine from tRNA(Leu) (Mursinna et al., Biochemistry 2001;40:5376-5381). To achieve atomic level insight into the role of T252 in LeuRS and the editing reaction of aaRSs, a series of molecular modeling studies including homology modeling and automated docking simulations were carried out. A 3D structure of E. coli LeuRS was constructed via homology modeling using the X-ray structure of Thermus thermophilus LeuRS as a template because the E. coli LeuRS structure is not available from X-ray or NMR studies. However, both the X-ray T. thermophilus and homology-modeled E. coli structures were used in our studies. Amino acid binding sites in the proposed editing domain, which is also called the connective polypeptide 1 (CP1) domain, were investigated by automated docking studies. The root mean square deviation (RMSD) for backbone atoms between the X-ray and homology-modeled structures was 1.18 A overall and 0.60 A for the editing (CP1) domain. Automated docking studies of a leucine ligand into the editing domain were performed for both structures: homology structure of E. coli LeuRS and X-ray structure of T. thermophilus LeuRS for comparison. The results of the docking studies suggested that there are two possible amino acid binding sites in the CP1 domain for both proteins. The first site lies near a threonine-rich region that includes the highly conserved T252 residue, which is important for amino acid discrimination. The second site is located in a flexible loop region surrounded by residues E292, A293, M295, A296, and M298. The important T252 residue is at the bottom of the first binding pocket.

Adenosine Monophosphate↗

PubMed2004

Prediction of HIV-1 integrase/viral DNA interactions in the catalytic domain by fast molecular docking.

This study details the separate analyses of binding specificity of HIV-1 integrase (IN) and viral B-DNA forms through ligand-receptor docking studies by means of a fast molecular docking method. The application of solvated electrostatics with the University of Houston Brownian Dynamics Program (UHBD) and configurational sampling by the Daughter of Turnip (DOT) docking program resulted in the computation of energies of more than 113 billion configurations for each ligand-receptor docking study, a procedure considered computationally intractable a few years ago. A specific binding pattern of viral DNA to the IN catalytic domain region has been predicted as a result of these calculations. In a representative docked configuration, we observe the 3'-hydroxyl of the conserved deoxyadenosine to be close to one of the two divalent metal ions that are necessary for catalysis. A superimposition of our energy-minimized docked complex on representative structures from a molecular dynamics (MD) simulation of a crystallographically resolved IN/inhibitor complex revealed an overlap of viral DNA with the inhibitor, indicating that the bound inhibitor might operate by blocking substrate binding. The DOT docking calculation also identified a second, adjacent DNA-binding site, which we believe is the nonspecific host DNA binding site. The binding pattern predicted by DOT complements previous electrostatics, MD simulation, photo-cross-linking, and mutagenesis studies and also provides a further refinement of the IN/viral DNA binding interaction as a basis for new structure-based design efforts.

Catalytic Domain↗

PubMed2004

Molecular dissection of a critical specificity determinant within the amino acid editing domain of leucyl-tRNA synthetase.

A highly conserved threonine residue marks the amino acid binding pocket within the editing active site of leucyl-tRNA synthetases (LeuRSs). It is essential to substrate specificity for the Escherichia coli enzyme in that it blocks the cognate leucine amino acid from binding in the hydrolytic editing active site. We combined mutagenesis and computational approaches to elucidate the molecular role of the critical side chain of this threonine residue. Removal of the terminal methyl group of the threonine side chain by replacement with serine yielded a mutant LeuRS that hydrolyzes Leu-tRNA(Leu). Substitution of valine for the conserved threonine conferred similar activities to the wild-type enzyme. However, an additional substitution within the editing active site suggested synergistic interactions with the conserved threonine site that significantly affected amino acid editing. On the basis of our combined biochemical and computational data, we propose that the threonine 252 side chain not only sterically hinders the cognate charged leucine from binding for hydrolysis but also plays a critical role in maintaining an active site geometry that is required for the fidelity of LeuRS.

Amino Acid Substitution↗

PubMed2003

Charge-charge interactions are key determinants of the pK values of ionizable groups in ribonuclease Sa (pI=3.5) and a basic variant (pI=10.2).

The pK values of the titratable groups in ribonuclease Sa (RNase Sa) (pI=3.5), and a charge-reversed variant with five carboxyl to lysine substitutions, 5K RNase Sa (pI=10.2), have been determined by NMR at 20 degrees C in 0.1M NaCl. In RNase Sa, 18 pK values and in 5K, 11 pK values were measured. The carboxyl group of Asp33, which is buried and forms three intramolecular hydrogen bonds in RNase Sa, has the lowest pK (2.4), whereas Asp79, which is also buried but does not form hydrogen bonds, has the most elevated pK (7.4). These results highlight the importance of desolvation and charge-dipole interactions in perturbing pK values of buried groups. Alkaline titration revealed that the terminal amine of RNase Sa and all eight tyrosine residues have significantly increased pK values relative to model compounds.A primary objective in this study was to investigate the influence of charge-charge interactions on the pK values by comparing results from RNase Sa with those from the 5K variant. The solution structures of the two proteins are very similar as revealed by NMR and other spectroscopic data, with only small changes at the N terminus and in the alpha-helix. Consequently, the ionizable groups will have similar environments in the two variants and desolvation and charge-dipole interactions will have comparable effects on the pK values of both. Their pK differences, therefore, are expected to be chiefly due to the different charge-charge interactions. As anticipated from its higher net charge, all measured pK values in 5K RNase are lowered relative to wild-type RNase Sa, with the largest decrease being 2.2 pH units for Glu14. The pK differences (pK(Sa)-pK(5K)) calculated using a simple model based on Coulomb's Law and a dielectric constant of 45 agree well with the experimental values. This demonstrates that the pK differences between wild-type and 5K RNase Sa are mainly due to changes in the electrostatic interactions between the ionizable groups. pK values calculated using Coulomb's Law also showed a good correlation (R=0.83) with experimental values. The more complex model based on a finite-difference solution to the Poisson-Boltzmann equation, which considers desolvation and charge-dipole interactions in addition to charge-charge interactions, was also used to calculate pK values. Surprisingly, these values are more poorly correlated (R=0.65) with the values from experiment. Taken together, the results are evidence that charge-charge interactions are the chief perturbant of the pK values of ionizable groups on the protein surface, which is where the majority of the ionizable groups are positioned in proteins.

Aspartic Acid↗

PubMed2003

pK values of histidine residues in ribonuclease Sa: effect of salt and net charge.

The primary goal of this study was to gain a better understanding of the effect of environment and ionic strength on the pK values of histidine residues in proteins. The salt-dependence of pK values for two histidine residues in ribonuclease Sa (RNase Sa) (pI=3.5) and a variant in which five acidic amino acids have been changed to lysine (5K) (pI=10.2) was measured and compared to pK values of model histidine-containing peptides. The pK of His53 is elevated by two pH units (pK=8.61) in RNase Sa and by nearly one pH unit (pK=7.39) in 5K at low salt relative to the pK of histidine in the model peptides (pK=6.6). The pK for His53 remains elevated in 1.5M NaCl (pK=7.89). The elevated pK for His53 is a result of screenable electrostatic interactions, particularly with Glu74, and a non-screenable hydrogen bond interaction with water. The pK of His85 in RNase Sa and 5K is slightly below the model pK at low salt and merges with this value at 1.5M NaCl. The pK of His85 reflects mainly effects of long-range Coulombic interactions that are screenable by salt. The tautomeric states of the neutral histidine residues are changed by charge reversal. The histidine pK values in RNase Sa are always higher than the pK values in the 5K variant. These results emphasize that the net charge of the protein influences the pK values of the histidine residues. Structure-based pK calculations capture the salt-dependence relatively well but are unable to predict absolute histidine pK values.

Electrophoresis, Gel, Two-Dimensional↗

PubMed2003

Molecular dynamics studies of alanine racemase: a structural model for drug design.

Alanine racemase (AlaR) is a bacterial enzyme that catalyzes the interconversion of L- and D-alanine, which is an essential constituent of the peptidoglycan layer of the bacterial cell wall and requires pyridoxal 5'-phosphate (PLP) as a cofactor. The enzyme is universal to bacteria, including mycobacteria, making it an attractive target for drug design. To investigate the effects of flexibility on the binding modes of the substrate and an inhibitor and to analyze how the active site is affected by the presence of the substrate versus inhibitor, a molecular dynamics simulation on the full AlaR dimer from Bacillus stearothermophilus (pdb code: 1SFT) with a D-alanine molecule in one active site and the noncovalent inhibitor, propionate, in the second site has been carried out. Within the time scale of the simulation, we show that the active site becomes more stabilized in the presence of substrate versus inhibitor. The results of this simulation are in agreement with the proposed mechanism of alanine racemase reaction in which the substrate carboxyl group directly participates in the catalysis by acting cooperatively with Tyr 265' and Lys 39. A structural water molecule in contact with both substrate and inhibitor (i.e., in both active sites) and bridging residues in both active sites was identified. It shows a remarkably low mobility and does not exchange with bulk water. This water molecule can be taken into account for the design of specific AlaR inhibitors by either utilizing it as a bridging group or displacing it with an inhibitor atom. The results presented here provide insights into the dynamics of the alanine racemase in the presence of substrate/inhibitor, which will be used for the rational design of novel inhibitors.

Alanine↗

PubMed2003

Molecular dynamics studies of the wild-type and double mutant HIV-1 integrase complexed with the 5CITEP inhibitor: mechanism for inhibition and drug resistance.

The human immunodeficiency virus type 1 (HIV-1) integrase (IN) is an essential enzyme in the life cycle of the virus and is an attractive target for the development of new drugs useful in acquired immunodeficiency syndrome multidrug therapy. Starting from the crystal structure of the 5CITEP inhibitor bound to the active site in the catalytic domain of the HIV-1 IN, two different molecular dynamics simulations in water have been carried out. In the first simulation the wild-type IN was used, whereas in the second one the double mutation T66I/M154I, described to lead to drug resistance, was introduced in the protein. Compelling differences have been observed in these two structures during analyses of the molecular dynamics trajectories, particularly in the inhibitor binding modes and in the conformational flexibility of the loop (residues 138-149) located near the three catalytic residues in the active site (Asp(64), Asp(116), Glu(152)). Because the conformational flexibility of this region is important for efficient biological activity and its behavior is quite different in the two models, we suggest a hypothetical mechanism for the inhibition and drug resistance of HIV-1 IN. These results can be useful for the rational design of more potent and selective integrase inhibitors and may allow for the design of inhibitors that will be more robust against known resistance mutations.

Binding Sites↗