A simple clustering algorithm can be accurate enough for use in calculations of pKs in macromolecules

Effect of Calcium Ion Removal, Ionic Strength, and Temperature on the Conformation Change in Calmodulin Protein at Physiological pH

The response of the calmodulin (CaM) protein as a function of calcium ion removal, ionic strength, and temperature at physiological pH condition was investigated using classical molecular dynamics simulations. Changing the ionic strength and temperature came out to be two of the possible routes for observing a conformation change in the protein. This behavior is similar to the conformation change observed in our previous study where a change in the pH was observed to trigger a conformation change in this protein. In the present study, as the calcium ions are removed from the protein, the protein is observed to acquire more flexibility. This flexibility is observed to be more prominent at a higher ionic strength. At a lower ionic strength of 150 mM with all the four calcium ions intact, the N- and C-lobes are observed to come close to a distance of 30 Å starting from an initial separation distance of 48 Å. This conformation change is observed to take place around 50 ns in a simulation of 100 ns. As a second parameter, temperature is observed to play a key role in the conformation change of the protein. With an increase in the temperature, the protein is observed to acquire a more compact form with the formation of different salt bridges between the residues of the N- and the C-lobes. The salt bridge formation leads to an overall lowering of the energy of the protein thus favoring the bending of the two lobes towards each other. The improper and dihedral terms show a significant shift thus leading to a more compact form on increasing the temperature. Another set of simulations is also performed at an increased temperature of 500 K to verify the reproducibility of the results. Thus a set of three possible alterations in the environmental conditions of the protein CaM are studied, with two of them giving rise to a conformation change and one adding flexibility to the protein.

Unstable, metastable, or stable halogen bonding interaction involving negatively charged donors? A statistical and computational chemistry study

The noncovalent halogen bonding could be attributed to the attraction between the positively charged σ-hole and a nucleophile. Quantum mechanics (QM) calculation indicated that the negatively charged organohalogens have no positively charged σ-hole on their molecular surface, leading to a postulation of repulsion between negatively charged organohalogens and nucleophiles in vacuum. However, PDB survey revealed that 24% of the ligands with halogen bonding geometry could be negatively charged. Moreover, 36% of ionizable drugs in CMC (Comprehensive Medicinal Chemistry) are possibly negatively charged at pH 7.0. QM energy scan showed that the negatively charged halogen bonding is probably metastable in vacuum. However, the QM calculated bonding energy turned negative in various solvents, suggesting that halogen bonding with negatively charged donors should be stable in reality. Indeed, QM/MM calculation on three crystal structures with negatively charged ligands revealed that the negatively charged halogen bonding was stable. Hence, we concluded that halogen bonding with negatively charged donors is unstable or metastable in vacuum but stable in protein environment, and possesses similar geometric and energetic characteristics as conventional halogen bonding. Therefore, negatively charged organohalogens are still effective halogen bonding donors for medicinal chemistry and other applications.

Visualizing phosphodiester-bond hydrolysis by an endonuclease

The enzymatic hydrolysis of DNA phosphodiester bonds has been widely studied, but the chemical reaction has not yet been observed. Here we follow the generation of a DNA double-strand break (DSB) by the Desulfurococcus mobilis homing endonuclease I-DmoI, trapping sequential stages of a two-metal-ion cleavage mechanism. We captured intermediates of the different catalytic steps, and this allowed us to watch the reaction by 'freezing' multiple states. We observed the successive entry of two metals involved in the reaction and the arrival of a third cation in a central position of the active site. This third metal ion has a crucial role, triggering the consecutive hydrolysis of the targeted phosphodiester bonds in the DNA strands and leaving its position once the DSB is generated. The multiple structures show the orchestrated conformational changes in the protein residues, nucleotides and metals during catalysis.

A Molecular Dynamics Investigation of Mycobacterium tuberculosis Prenyl Synthases: Conformational Flexibility and Implications for Computer-aided Drug Discovery

With the rise in antibiotic resistance, there is interest in discovering new drugs active against new targets. Here, we investigate the dynamic structures of three isoprenoid synthases from Mycobacterium tuberculosis using molecular dynamics (MD) methods with a view to discovering new drug leads. Two of the enzymes, cis-farnesyl diphosphate synthase (cis-FPPS) and cis-decaprenyl diphosphate synthase (cis-DPPS), are involved in bacterial cell wall biosynthesis, while the third, tuberculosinyl adenosine synthase (Rv3378c), is involved in virulence factor formation. The MD results for these three enzymes were then compared with previous results on undecaprenyl diphosphate synthase (UPPS) by means of active site volume fluctuation and principal component analyses. In addition, an analysis of the binding of prenyl diphosphates to cis-FPPS, cis-DPPS, and UPPS utilizing the new MD results is reported. We also screened libraries of inhibitors against cis-DPPS, finding ~1 μm inhibitors, and used the receiver operating characteristic-area under the curve (ROC-AUC) method to test the predictive power of X-ray and MD-derived cis-DPPS receptors. We found that one compound with potent M. tuberculosis cell growth inhibition activity was an IC(50) ~0.5- to 20-μm inhibitor (depending on substrate) of cis-DPPS, a ~660-nm inhibitor of Rv3378c as well as a 4.8-μm inhibitor of cis-FPPS, opening up the possibility of multitarget inhibition involving both cell wall biosynthesis and virulence factor formation.

Design of sweet protein based sweeteners: hints from structure-function relationships

Sweet proteins represent a class of natural molecules, which are extremely interesting regarding their potential use as safe low-calories sweeteners for individuals who need to control sugar intake, such as obese or diabetic subjects. Punctual mutations of amino acid residues of MNEI, a single chain derivative of the natural sweet protein monellin, allow the modulation of its taste. In this study we present a structural and functional comparison between MNEI and a sweeter mutant Y65R, containing an extra positive charge on the protein surface, in conditions mimicking those of typical beverages. Y65R exhibits superior sweetness in all the experimental conditions tested, has a better solubility at mild acidic pH and preserves a significant thermal stability in a wide range of pH conditions, although slightly lower than MNEI. Our findings confirm the advantages of structure-guided protein engineering to design improved low-calorie sweeteners and excipients for food and pharmaceutical preparations.

In silico insights into protein-protein interactions and folding dynamics of the saposin-like domain of Solanum tuberosum aspartic protease

The plant-specific insert is an approximately 100-residue domain found exclusively within the C-terminal lobe of some plant aspartic proteases. Structurally, this domain is a member of the saposin-like protein family, and is involved in plant pathogen defense as well as vacuolar targeting of the parent protease molecule. Similar to other members of the saposin-like protein family, most notably saposins A and C, the recently resolved crystal structure of potato (Solanum tuberosum) plant-specific insert has been shown to exist in a substrate-bound open conformation in which the plant-specific insert oligomerizes to form homodimers. In addition to the open structure, a closed conformation also exists having the classic saposin fold of the saposin-like protein family as observed in the crystal structure of barley (Hordeum vulgare L.) plant-specific insert. In the present study, the mechanisms of tertiary and quaternary conformation changes of potato plant-specific insert were investigated in silico as a function of pH. Umbrella sampling and determination of the free energy change of dissociation of the plant-specific insert homodimer revealed that increasing the pH of the system to near physiological levels reduced the free energy barrier to dissociation. Furthermore, principal component analysis was used to characterize conformational changes at both acidic and neutral pH. The results indicated that the plant-specific insert may adopt a tertiary structure similar to the characteristic saposin fold and suggest a potential new structural motif among saposin-like proteins. To our knowledge, this acidified PSI structure presents the first example of an alternative saposin-fold motif for any member of the large and diverse SAPLIP family.

Structural basis for the recruitment and activation of the Legionella phospholipase VipD by the host GTPase Rab5

A challenge for microbial pathogens is to assure that their translocated effector proteins target only the correct host cell compartment during infection. The Legionella pneumophila effector vacuolar protein sorting inhibitor protein D (VipD) localizes to early endosomal membranes and alters their lipid and protein composition, thereby protecting the pathogen from endosomal fusion. This process requires the phospholipase A1 (PLA1) activity of VipD that is triggered specifically on VipD binding to the host cell GTPase Rab5, a key regulator of endosomes. Here, we present the crystal structure of VipD in complex with constitutively active Rab5 and reveal the molecular mechanism underlying PLA1 activation. An active site-obstructing loop that originates from the C-terminal domain of VipD is repositioned on Rab5 binding, thereby exposing the catalytic pocket within the N-terminal PLA1 domain. Substitution of amino acid residues located within the VipD-Rab5 interface prevented Rab5 binding and PLA1 activation and caused a failure of VipD mutant proteins to target to Rab5-enriched endosomal structures within cells. Experimental and computational analyses confirmed an extended VipD-binding interface on Rab5, explaining why this L. pneumophila effector can compete with cellular ligands for Rab5 binding. Together, our data explain how the catalytic activity of a microbial effector can be precisely linked to its subcellular localization.

Protocols utilizing constant pH molecular dynamics to compute pH-dependent binding free energies

In protein–ligand binding, the electrostatic environments of the two binding partners may vary significantly in bound and unbound states, which may lead to protonation changes upon binding. In cases where ligand binding results in a net uptake or release of protons, the free energy of binding is pH-dependent. Nevertheless, conventional free energy calculations and molecular docking protocols typically do not rigorously account for changes in protonation that may occur upon ligand binding. To address these shortcomings, we present a simple methodology based on Wyman’s binding polynomial formalism to account for the pH dependence of binding free energies and demonstrate its use on cucurbit[7]uril (CB[7]) host–guest systems. Using constant pH molecular dynamics and a reference binding free energy that is taken either from experiment or from thermodynamic integration computations, the pH-dependent binding free energy is determined. This computational protocol accurately captures the large pK_a shifts observed experimentally upon CB[7]:guest association and reproduces experimental binding free energies at different levels of pH. We show that incorrect assignment of fixed protonation states in free energy computations can give errors of >2 kcal/mol in these host–guest systems. Use of the methods presented here avoids such errors, thus suggesting their utility in computing proton-linked binding free energies for protein–ligand complexes.

Structure-based engineering of a minimal porin reveals loop-independent channel closure

Porins, like outer membrane protein G (OmpG) of Escherichia coli, are ideal templates among ion channels for protein and chemical engineering because of their robustness and simple architecture. OmpG shows fast transitions between open and closed states, which were attributed to loop 6 (L6). As flickering limits single-channel-based applications, we pruned L6 by either 8 or 12 amino acids. While the open probabilities of both L6 variants resemble that of native OmpG, their gating frequencies were reduced by 63 and 81%, respectively. Using the 3.2 Å structure of the shorter L6 variant in the open state, we engineered a minimal porin (220 amino acids), where all remaining extramembranous loops were truncated. Unexpectedly, this minimized porin still exhibited gating, but it was 5-fold less frequent than in OmpG. The residual gating of the minimal pore is hence independent of L6 rearrangements and involves narrowing of the ion conductance pathway most probably driven by global stretching-flexing deformations of the membrane-embedded β-barrel.

Cysteine-mediated dynamic hydrogen-bonding network in the active site of Pin1

Enzymes catalyze a plethora of chemical reactions that are tightly regulated and intricately coupled in biology. Catalysis of phosphorylation-dependent cis-trans isomerization of peptidyl-prolyl bonds, which act as conformational switches in regulating many post-phosphorylation processes, is considered to be one of the most critical. Pin1 is a cis-trans isomerase of peptidyl-prolyl(ω-) bonds of phosphorylated-Ser/Thr-Pro motifs and has been implicated in many diseases. Structural and experimental studies are still unable to resolve the mechanistic role and protonation states of two adjacent histidines (His59 and His157) and a cysteine (Cys113) in the active site of Pin1. Here, we show that the protonation state of Cys113 mediates a dynamic hydrogen-bonding network in the active site of Pin1, involving the two adjacent histidines and several other residues that are highly conserved and necessary for catalysis. We have used detailed free energy calculations and molecular dynamics simulations, complementing previous experiments, to resolve the ambiguities in the orientations of the histidines and protonation states of these key active site residues, details that are critical for fully understanding the mechanism of Pin1 and necessary for developing potent inhibitors. Importantly, Cys113 is shown to alternate between the unprotonated and neutral states, unprotonated in free Pin1 and neutral in substrate-bound Pin1. Our results are consistent with experiments and provide an explanation for the chemical reactivity of free Pin1 that is suggested to be necessary for the regulation of the enzyme.

Molecular Dynamics Study of Helicobacter pylori Urease

Helicobacter pylori have been implicated in an array of gastrointestinal disorders including, but not limited to, gastric and duodenal ulcers and adenocarcinoma. This bacterium utilizes an enzyme, urease, to produce copious amounts of ammonia through urea hydrolysis in order to survive the harsh acidic conditions of the stomach. Molecular dynamics (MD) studies on the H. pylori urease enzyme have been employed in order to study structural features of this enzyme that may shed light on the hydrolysis mechanism. A total of 400 ns of MD simulation time were collected and analyzed in this study. A wide-open flap state previously observed in MD simulations on Klebsiella aerogenes [Roberts et al. J. Am. Chem. Soc.2012, 134, 9934] urease has been identified in the H. pylori enzyme that has yet to be experimentally observed. Critical distances between residues on the flap, contact points in the closed state, and the separation between the active site Ni²⁺ ions and the critical histidine α322 residue were used to characterize flap motion. An additional flap in the active site was elaborated upon that we postulate may serve as an exit conduit for hydrolysis products. Finally we discuss the internal hollow cavity and present analysis of the distribution of sodium ions over the course of the simulation.

Utility of computational structural biology in mass spectrometry

Recent developments of mass spectrometry (MS) allow us to identify, estimate, and characterize proteins and protein complexes. At the same time, structural biology helps to determine the protein structure and its structure-function relationship. Together, they aid to understand the protein structure, property, function, protein-complex assembly, protein-protein interaction and dynamics. The present chapter is organized with illustrative results to demonstrate how experimental mass spectrometry can be combined with computational structural biology for detailed studies of protein's structures. We have used tumor differentiation factor protein/peptide as ligand and Hsp70/Hsp90 as receptor protein as examples to study ligand-protein interaction. To investigate possible protein conformation we will describe two proteins, lysozyme and myoglobin.

New insights into the binding and catalytic mechanisms of Bacillus thuringiensis lactonase: insights into B. thuringiensis AiiA mechanism

The lactonase enzyme (AiiA) produced by Bacillus thuringiensis serves to degrade autoinducer-1 (AI-1) signaling molecules in what is an evolved mechanism by which to compete with other bacteria. Bioassays have been previously performed to determine whether the AI-1 aliphatic tail lengths have any effect on AiiA's bioactivity, however, data to date are conflicting. Additionally, specific residue contributions to the catalytic activity of AiiA provide for some interesting questions. For example, it has been proposed that Y194 serves to provide an oxyanion hole to AI-1 which is curious given the fact the substrate spans two Zn(2+) ions. These ions might conceivably provide enough charge to promote both ligand stability and the carbonyl activation necessary to drive a nucleophilic attack. To investigate these questions, multiple molecular dynamics simulations were performed across a family of seven acylated homoserine lactones (AHL) along with their associated intermediate and product states. Distance analyses and interaction energy analyses were performed to investigate current bioassay data. Our simulations are consistent with experimental studies showing that AiiA degrades AHLs in a tail length independent manner. However, the presence of the tail is required for activity. Also, the putative oxyanion hole function of Y194 toward the substrate is not observed in any of the reactant or product state simulation trajectories, but does seem to show efficacy in stabilizing the intermediate state. Last, we argue through ionization state analyses, that the proton shuttling necessary for catalytic activity might be mediated by both water and substrate-based intra-molecular proton transfer. Based on this argument, an alternate catalytic mechanism is proposed.

Phenotypic expression of a novel C282Y/R226G compound heterozygous state in HFE hemochromatosis: molecular dynamics and biochemical studies

Most adults affected with hereditary hemochromatosis are homozygous for a single point mutation of HFE (p.Cys282Tyr). Apart from the compound heterozygous state for the p.Cys282Tyr mutant and the widespread p.His63Asp variant allele, other rare HFE mutations can be found in trans and may have clinical impact. In the present report we describe the structural and functional consequences of a new mutation, namely the p.Arg226Gly which was inherited in trans with the p.Cys282Tyr allele in a patient affected with a mild iron overload. Because the R226G substitution is located in the vicinity of the normal Cys225S-S282Cys disulfide bond we initially investigated the structure of the variant by molecular dynamics techniques in order to estimate the effect of the mutation on the global structure of HFE domain α3. We found that the solvation free energy, hydrophobicity and formation of salt bridges are slightly modified with the global secondary structure of the α3 domain being conserved. In a previous paper, we demonstrated that the Q283P substitution leads to the loss of the normal Cys225S-S282Cys disulfide bridge. Similar to the Q283P substitution, the R226G substitution does not substitute a residue directly involved in the formation of the disulfide bridge. However, unlike the p.Gln283Pro variant which destroyed the normal disulfide bridge, the R226G mutation does not affect the normal Cys225S-S282Cys bridge. Furthermore based on cell line studies we clearly show that the mutation does not prevent cell surface localization, β2-microglobulin association and binding to transferrin receptor 1. This new compound heterozygous phenotype is very close to those of the C282Y/H63D compound heterozygous patients who display the biochemical hemochromatosis phenotype but with lower body iron stores than C282Y homozygotes. Our results do not exclude unknown genetic and/or metabolic factors that may act synergistically to increase the ferritin level.

Improving the Efficiency of Free Energy Calculations in the Amber Molecular Dynamics Package

Alchemical transformations are widely used methods to calculate free energies. Amber has traditionally included support for alchemical transformations as part of the sander molecular dynamics (MD) engine. Here we describe the implementation of a more efficient approach to alchemical transformations in the Amber MD package. Specifically we have implemented this new approach within the more computational efficient and scalable pmemd MD engine that is included with the Amber MD package. The majority of the gain in efficiency comes from the improved design of the calculation, which includes better parallel scaling and reduction in the calculation of redundant terms. This new implementation is able to reproduce results from equivalent simulations run with the existing functionality, but at 2.5 times greater computational efficiency. This new implementation is also able to run softcore simulations at the λ end states making direct calculation of free energies more accurate, compared to the extrapolation required in the existing implementation. The updated alchemical transformation functionality will be included in the next major release of Amber (scheduled for release in Q1 2014) and will be available at http://ambermd.org, under the Amber license.

Heme prevents amyloid beta peptide aggregation through hydrophobic interaction based on molecular dynamics simulation

Heme, which is abundant in hemoglobin and many other hemoproteins, is known to play an important role in electron transfer, oxygen transport, regulation of gene expression, and many other biological functions. With the belief that the aggregation of Aβ peptides forming higher order oligomers is one of the central pathological pathways in Alzheimer's disease, the formation of the Aβ-heme complex is essential as it inhibits Aβ aggregation and protects the neurons from degradation. In our studies, conventional molecular dynamics simulations were performed on the 1 Aβ + 1 heme and 2 Aβ + 4 hemes system, respectively, with the identification of several dominant binding motifs. We found that hydrophobic residues of the Aβ peptide have a high affinity to interact with heme instead of the histidine residue. We conclude that hydrophobic interaction plays a dominant role in the Aβ-heme complex formation which indirectly serves to physically prevent Aβ aggregation.

Study of the interaction of Huperzia saururus Lycopodium alkaloids with the acetylcholinesterase enzyme

In the present study, we describe and compare the binding modes of three Lycopodium alkaloids (sauroine, 6-hydroxylycopodine and sauroxine; isolated from Huperzia saururus) and huperzine A with the enzyme acetylcholinesterase. Refinement and rescoring of the docking poses (obtained with different programs) with an all atom force field helped to improve the quality of the protein-ligand complexes. Molecular dynamics simulations were performed to investigate the complexes and the alkaloid's binding modes. The combination of the latter two methodologies indicated that binding in the active site is favored for the active compounds. On the other hand, similar binding energies in both the active and the peripheral sites were obtained for sauroine, thus explaining its experimentally determined lack of activity. MM-GBSA predicted the order of binding energies in agreement with the experimental IC50 values.

An arm-swapped dimer of the Streptococcus pyogenes pilin specific assembly factor SipA

Streptococcus pyogenes (group A streptococcus [GAS]) is a major human pathogen. Attachment of GAS to host cells depends in large part on pili. These assemblies are built from multiple covalently linked subunits of a backbone protein (FctA), which forms the shaft of the pilus, and two minor pilin proteins, FctB anchoring the pilus to the cell wall and Cpa functioning as the adhesin at the tip. Polymerisation of the pilin subunits is mediated by a specific sortase, which catalyzes the formation of peptide bonds linking successive subunits. An additional gene, SipA, is also essential for GAS pilus polymerisation, but its function remains undefined. Here we report the crystal structure of a truncated SipA protein from GAS, determined at 1.67Å resolution. The structure reveals that SipA has the same core fold as the Escherichia coli type-I signal peptidase (SPase-I), but has a much smaller non-catalytic domain. The truncated protein, which lacks 9 N-terminal residues, forms an arm-swapped dimer in which the C-terminal β-strand of each monomer crosses over to interact with an N-terminal strand from the other monomer. In addition, there is no peptide binding cleft and significant differences in the putative membrane association region.

Protonation and pK changes in protein-ligand binding

Formation of protein-ligand complexes causes various changes in both the receptor and the ligand. This review focuses on changes in pK and protonation states of ionizable groups that accompany protein-ligand binding. Physical origins of these effects are outlined, followed by a brief overview of the computational methods to predict them and the associated corrections to receptor-ligand binding affinities. Statistical prevalence, magnitude and spatial distribution of the pK and protonation state changes in protein-ligand binding are discussed in detail, based on both experimental and theoretical studies. While there is no doubt that these changes occur, they do not occur all the time; the estimated prevalence varies, both between individual complexes and by method. The changes occur not only in the immediate vicinity of the interface but also sometimes far away. When receptor-ligand binding is associated with protonation state change at particular pH, the binding becomes pH dependent: we review the interplay between sub-cellular characteristic pH and optimum pH of receptor-ligand binding. It is pointed out that there is a tendency for protonation state changes upon binding to be minimal at physiologically relevant pH for each complex (no net proton uptake/release), suggesting that native receptor-ligand interactions have evolved to reduce the energy cost associated with ionization changes. As a result, previously reported statistical prevalence of these changes - typically computed at the same pH for all complexes - may be higher than what may be expected at optimum pH specific to each complex. We also discuss whether proper account of protonation state changes appears to improve practical docking and scoring outcomes relevant to structure-based drug design. An overview of some of the existing challenges in the field is provided in conclusion.

Effects of histidine protonation and rotameric states on virtual screening of M. tuberculosis RmlC

While it is well established that protonation and tautomeric states of ligands can significantly affect the results of virtual screening, such effects of ionizable residues of protein receptors are less well understood. In this study, we focus on histidine protonation and rotameric states and their impact on virtual screening of Mycobacterium tuberculosis enzyme RmlC. Depending on the net charge and the location of proton(s), a histidine can adopt three states: HIP (+1 charged, both δ- and ε-nitrogens protonated), HID (neutral, δ-nitrogen protonated), and HIE (neutral, ε-nitrogen protonated). Due to common ambiguities in X-ray crystal structures, a histidine may also be resolved as three additional states with its imidazole ring flipped. Here, we systematically investigate the predictive power of 36 receptor models with different protonation and rotameric states of two histidines in the RmlC active site by using results from a previous high-throughput screening. By measuring enrichment factors and area under the receiver operating characteristic curves, we show that virtual screening results vary depending on hydrogen bonding networks provided by the histidines, even in the cases where the ligand does not obviously interact with the side chain. Our results also suggest that, even with the help of widely used pKa prediction software, assigning histidine protonation and rotameric states for virtual screening can still be challenging and requires further examination and systematic characterization of the receptor-ligand complex.

Substrate placement influences reactivity in non-heme Fe(II) halogenases and hydroxylases

We employ error-corrected density functional theory methods to map out the dependence of reactivity on substrate position for SyrB2, a member of a family of non-heme iron halogenases and hydroxylases that are only reactive toward amino acid substrates delivered via prosthetic phosphopantetheine arms. For the initial hydrogen abstraction step, the inherent flexibility of the phosphopantetheine molecule weakens the position dependence for both the native substrate (threonine for SyrB2) and alternative substrates. Over a 5 Å window of substrate positions, the tethered hydrogen abstraction step proceeds with nearly identical activation energies and donor-acceptor distances in the transition state. The propensity of a particular substrate toward halogenation or hydroxylation is found to depend strongly on the substrate placement following hydrogen abstraction, with deeper substrate delivery into the active (for native substrates) site favoring halogenation and shallower substrate delivery favoring hydroxylation.

Disease-associated mutations in the coil 2B domain of human lamin A/C affect structural properties that mediate dimerization and intermediate filament formation

The lamin proteins are essential components of the nuclear lamina of eukaryotic cells, that are involved in a complex association mechanism to attain a functional supermolecular structure. Mutations of the lamin A/C gene are associated with several different neuromuscular diseases, and the detailed effect of disease-associated amino acid substitutions on the structure and stability of human lamin dimers is yet unknown. Here we present a structural and thermodynamic characterization by means of molecular dynamics simulations of the effect of pathological mutations (S326T, R331P, R331Q, E347K, E358K, M371K, and R377H) on the association of the coil 2B domains that mediate lamin A/C oligomerization. The structures attained during the simulations, along with the quantification of the contribution of each residue to the dimerization energies, support a lamin association mechanism mediated by homophilic intermolecular interactions promoted by dissociative conformational changes at distinct positions in the coiled coil. The pathogenic mutations can both increase or decrease the stability of lamin A/C dimers, and a possible correlation between the effect of the amino acid substitutions and disease onset and severity is presented.

Receptor-based virtual screening and biological characterization of human apurinic/apyrimidinic endonuclease (Ape1) inhibitors

The endonucleolytic activity of human apurinic/apyrimidinic endonuclease (AP endo, Ape1) is a major factor in maintaining the integrity of the genome. Conversely, as an undesired effect, Ape1 overexpression has been linked to resistance to radio- and chemotherapeutic treatments in several human tumors. Inhibition of Ape1 using siRNA or the expression of a dominant negative form of the protein has been shown to sensitize cells to DNA-damaging agents, including various chemotherapeutic agents. Therefore, inhibition of the enzymatic activity of Ape1 might result in a potent antitumor therapy. A number of small molecules have been described as Ape1 inhibitors; however, those compounds are in the early stages of development. Herein we report the identification of new compounds as potential Ape1 inhibitors through a docking-based virtual screening technique. Some of the compounds identified have in vitro activities in the low-to-medium micromolar range. Interaction of these compounds with the Ape1 protein was observed by mass spectrometry. These molecules also potentiate the cytotoxicity of the chemotherapeutic agent methyl methanesulfonate in fibrosarcoma cells. This study demonstrates the power of docking and virtual screening techniques as initial steps in the design of new drugs, and opens the door to the development of a new generation of Ape1 inhibitors.

Molecular modeling of the reaction pathway and hydride transfer reactions of HMG-CoA reductase

HMG-CoA reductase catalyzes the four-electron reduction of HMG-CoA to mevalonate and is an enzyme of considerable biomedical relevance because of the impact of its statin inhibitors on public health. Although the reaction has been studied extensively using X-ray crystallography, there are surprisingly no computational studies that test the mechanistic hypotheses suggested for this complex reaction. Theozyme and quantum mechanical (QM)/molecular mechanical (MM) calculations up to the B3LYP/6-31g(d,p)//B3LYP/6-311++g(2d,2p) level of theory were employed to generate an atomistic description of the enzymatic reaction process and its energy profile. The models generated here predict that the catalytically important Glu83 is protonated prior to hydride transfer and that it acts as the general acid or base in the reaction. With Glu83 protonated, the activation energies calculated for the sequential hydride transfer reactions, 21.8 and 19.3 kcal/mol, are in qualitative agreement with the experimentally determined rate constant for the entire reaction (1 s(-1) to 1 min(-1)). When Glu83 is not protonated, the first hydride transfer reaction is predicted to be disfavored by >20 kcal/mol, and the activation energy is predicted to be higher by >10 kcal/mol. While not involved in the reaction as an acid or base, Lys267 is critical for stabilization of the transition state in forming an oxyanion hole with the protonated Glu83. Molecular dynamics simulations and MM/Poisson-Boltzmann surface area free energy calculations predict that the enzyme active site stabilizes the hemithioacetal intermediate better than the aldehyde intermediate. This suggests a mechanism in which cofactor exchange occurs before the breakdown of the hemithioacetal. Slowing the conversion to aldehyde would provide the enzyme with a mechanism to protect it from solvent and explain why the free aldehyde is not observed experimentally. Our results support the hypothesis that the pK(a) of an active site acidic group is modulated by the redox state of the cofactor. The oxidized cofactor and deprotonated Glu83 are closer in space after hydride transfer, indicating that indeed the cofactor may influence the pK(a) of Glu83 through an electrostatic interaction. The enzyme is able to catalyze the transfer of a hydride to the structurally and electronically distinct substrates by maintaining the general shape of the active site and adjusting the electrostatic environment through acid-base chemistry. Our results are in good agreement with the well-studied hydride transfer reactions catalyzed by liver alcohol dehydrogenase in calculated energy profile and reaction geometries despite different mechanistic functionalities.

Enhancing Conformation and Protonation State Sampling of Hen Egg White Lysozyme Using pH Replica Exchange Molecular Dynamics

We evaluate the efficiency of the pH replica exchange molecular dynamics (pH-REMD) method proposed by Itoh et al. (Proteins2011, 79, 3420-3436) by using it to predict the pKa values of the titratable residues in hen egg white lysozyme (HEWL). pKa values predicted using pH-REMD converge significantly faster than those calculated using constant pH molecular dynamics (CpHMD). Furthermore, increasing the frequency between exchange attempts in pH-REMD simulations improves protonation and conformational state sampling. By enabling the simulation to sample both conformational and protonation states more rapidly, pH-REMD simulations provide valuable insight into the pH-dependence of HEWL that the CpHMD simulations failed to capture. We present an efficient and highly scalable implementation of pH-REMD as an attractive enhancement to traditional CpHMD methods.

3D-RISM-Dock: A New Fragment-Based Drug Design Protocol

We explore a new approach in the rational design of specificity in molecular recognition of small molecules based on statistical-mechanical integral equation theory of molecular liquids in the form of the three-dimensional reference interaction site model with the Kovalenko-Hirata closure (3D-RISM-KH). The numerically stable iterative solution of conventional 3D-RISM equations includes the fragmental decomposition of flexible ligands, which are treated as distinct species in solvent mixtures of arbitrary complexity. The computed density functions for solution (including ligand) molecules are obtained as a set of discrete spatial grids that uniquely describe the continuous solvent-site distribution around the protein solute. Potentials of mean force derived from these distributions define the scoring function interfaced with the AutoDock program for an automated ranking of docked conformations. As a case study in terms of solvent composition, we analyze cooperative interactions encountered in the binding of a flexible thiamine molecule to the prion protein at near-physiological conditions. The predicted location and residency times of computed binding modes are in excellent agreement with the available experimental data.

Selective monocationic inhibitors of neuronal nitric oxide synthase. Binding mode insights from molecular dynamics simulations

The reduction of pathophysiologic levels of nitric oxide through inhibition of neuronal nitric oxide synthase (nNOS) has the potential to be therapeutically beneficial in various neurodegenerative diseases. We have developed a series of pyrrolidine-based nNOS inhibitors that exhibit excellent potencies and isoform selectivities (J. Am. Chem. Soc. 2010, 132, 5437). However, there are still important challenges, such as how to decrease the multiple positive charges derived from basic amino groups, which contribute to poor bioavailability, without losing potency and/or selectivity. Here we present an interdisciplinary study combining molecular docking, crystallography, molecular dynamics simulations, synthesis, and enzymology to explore potential pharmacophoric features of nNOS inhibitors and to design potent and selective monocationic nNOS inhibitors. The simulation results indicate that different hydrogen bond patterns, electrostatic interactions, hydrophobic interactions, and a water molecule bridge are key factors for stabilizing ligands and controlling ligand orientation. We find that a heteroatom in the aromatic head or linker chain of the ligand provides additional stability and blocks the substrate binding pocket. Finally, the computational insights are experimentally validated with double-headed pyridine analogues. The compounds reported here are among the most potent and selective monocationic pyrrolidine-based nNOS inhibitors reported to date, and 10 shows improved membrane permeability.

Calmodulin readily switches conformation upon protonating high pKa acidic residues

We investigate protonation as a possible route for triggering conformational change in proteins by focusing on the calmodulin (CaM) example. Two hundred nanosecond molecular dynamics (MD) simulations are performed on both the extended and compact forms of calcium loaded CaM. The stability of both structures is confirmed under prevailing conditions. Protonation of nine acidic residues with upshifted pK(a) values leads to a large conformational change in less than 100 ns. The structure attained is consistent with fluorescence resonance energy transfer experimental results as well as structures from an ensemble compatible with NMR data. Analysis of the MD trajectories summing up to one microsecond implies that the key events leading to the completion of the conformational change begins with an initial formation of a salt bridge between the N-lobe and the linker, followed by the bending of the C-lobe and the organization of a stabilizing hydrophobic patch between the lobes. We find that CaM utilizes its Ca(2+) ions to harden/soften different regions so as to achieve various conformations. Thus, barrier crossing between extended and compact forms of CaM which is normally a rare event due to the repulsive electrostatic interactions between the two lobes is facilitated by protonation of high pK(a) residues. The results delineate how pH changes might be utilized in the cell to achieve different conformation-related functions.

NMR structure of a heterodimeric SAM:SAM complex: characterization and manipulation of EphA2 binding reveal new cellular functions of SHIP2

The sterile alpha motif (SAM) for protein-protein interactions is encountered in over 200 proteins, but the structural basis for its interactions is just becoming clear. Here we solved the structure of the EphA2-SHIP2 SAM:SAM heterodimeric complex by use of NMR restraints from chemical shift perturbations, NOE and RDC experiments. Specific contacts between the protein surfaces differ significantly from a previous model and other SAM:SAM complexes. Molecular dynamics and docking simulations indicate fluctuations in the complex toward alternate, higher energy conformations. The interface suggests that EphA family members bind to SHIP2 SAM, whereas EphB members may not; correspondingly, we demonstrate binding of EphA1, but not of EphB2, to SHIP2. A variant of EphB2 SAM was designed that binds SHIP2. Functional characterization of a mutant EphA2 compromised in SHIP2 binding reveals two previously unrecognized functions of SHIP2 in suppressing ligand-induced activation of EphA2 and in promoting receptor coordinated chemotactic cell migration.

H++ 3.0: automating pK prediction and the preparation of biomolecular structures for atomistic molecular modeling and simulations

The accuracy of atomistic biomolecular modeling and simulation studies depend on the accuracy of the input structures. Preparing these structures for an atomistic modeling task, such as molecular dynamics (MD) simulation, can involve the use of a variety of different tools for: correcting errors, adding missing atoms, filling valences with hydrogens, predicting pK values for titratable amino acids, assigning predefined partial charges and radii to all atoms, and generating force field parameter/topology files for MD. Identifying, installing and effectively using the appropriate tools for each of these tasks can be difficult for novice and time-consuming for experienced users. H++ (http://biophysics.cs.vt.edu/) is a free open-source web server that automates the above key steps in the preparation of biomolecular structures for molecular modeling and simulations. H++ also performs extensive error and consistency checking, providing error/warning messages together with the suggested corrections. In addition to numerous minor improvements, the latest version of H++ includes several new capabilities and options: fix erroneous (flipped) side chain conformations for HIS, GLN and ASN, include a ligand in the input structure, process nucleic acid structures and generate a solvent box with specified number of common ions for explicit solvent MD.

Bluues: a program for the analysis of the electrostatic properties of proteins based on generalized Born radii

BACKGROUND

The Poisson-Boltzmann (PB) equation and its linear approximation have been widely used to describe biomolecular electrostatics. Generalized Born (GB) models offer a convenient computational approximation for the more fundamental approach based on the Poisson-Boltzmann equation, and allows estimation of pairwise contributions to electrostatic effects in the molecular context.

RESULTS

We have implemented in a single program most common analyses of the electrostatic properties of proteins. The program first computes generalized Born radii, via a surface integral and then it uses generalized Born radii (using a finite radius test particle) to perform electrostaic analyses. In particular the ouput of the program entails, depending on user's requirement: 1) the generalized Born radius of each atom; 2) the electrostatic solvation free energy; 3) the electrostatic forces on each atom (currently in a developmental stage); 4) the pH-dependent properties (total charge and pH-dependent free energy of folding in the pH range -2 to 18; 5) the pKa of all ionizable groups; 6) the electrostatic potential at the surface of the molecule; 7) the electrostatic potential in a volume surrounding the molecule;

CONCLUSIONS

Although at the expense of limited flexibility the program provides most common analyses with requirement of a single input file in PQR format. The results obtained are comparable to those obtained using state-of-the-art Poisson-Boltzmann solvers. A Linux executable with example input and output files is provided as supplementary material.

Small angle neutron scattering reveals pH-dependent conformational changes in Trichoderma reesei cellobiohydrolase I: implications for enzymatic activity

Cellobiohydrolase I (Cel7A) of the fungus Trichoderma reesei (now classified as an anamorph of Hypocrea jecorina) hydrolyzes crystalline cellulose to soluble sugars, making it of key interest for producing fermentable sugars from biomass for biofuel production. The activity of the enzyme is pH-dependent, with its highest activity occurring at pH 4-5. To probe the response of the solution structure of Cel7A to changes in pH, we measured small angle neutron scattering of it in a series of solutions having pH values of 7.0, 6.0, 5.3, and 4.2. As the pH decreases from 7.0 to 5.3, the enzyme structure remains well defined, possessing a spatial differentiation between the cellulose binding domain and the catalytic core that only changes subtly. At pH 4.2, the solution conformation of the enzyme changes to a structure that is intermediate between a properly folded enzyme and a denatured, unfolded state, yet the secondary structure of the enzyme is essentially unaltered. The results indicate that at the pH of optimal activity, the catalytic core of the enzyme adopts a structure in which the compact packing typical of a fully folded polypeptide chain is disrupted and suggest that the increased range of structures afforded by this disordered state plays an important role in the increased activity of Cel7A through conformational selection.

Charge state of the globular histone core controls stability of the nucleosome

Presented here is a quantitative model of the wrapping and unwrapping of the DNA around the histone core of the nucleosome that suggests a mechanism by which this transition can be controlled: alteration of the charge state of the globular histone core. The mechanism is relevant to several classes of posttranslational modifications such as histone acetylation and phosphorylation; several specific scenarios consistent with recent in vivo experiments are considered. The model integrates a description based on an idealized geometry with one based on the atomistic structure of the nucleosome, and the model consistently accounts for both the electrostatic and nonelectrostatic contributions to the nucleosome free energy. Under physiological conditions, isolated nucleosomes are predicted to be very stable (38 +/- 7 kcal/mol). However, a decrease in the charge of the globular histone core by one unit charge, for example due to acetylation of a single lysine residue, can lead to a significant decrease in the strength of association with its DNA. In contrast to the globular histone core, comparable changes in the charge state of the histone tail regions have relatively little effect on the nucleosome's stability. The combination of high stability and sensitivity explains how the nucleosome is able to satisfy the seemingly contradictory requirements for thermodynamic stability while allowing quick access to its DNA informational content when needed by specific cellular processes such as transcription.

Identification and structural characterization of FYVE domain-containing proteins of Arabidopsis thaliana

BACKGROUND

FYVE domains have emerged as membrane-targeting domains highly specific for phosphatidylinositol 3-phosphate (PtdIns(3)P). They are predominantly found in proteins involved in various trafficking pathways. Although FYVE domains may function as individual modules, dimers or in partnership with other proteins, structurally, all FYVE domains share a fold comprising two small characteristic double-stranded beta-sheets, and a C-terminal alpha-helix, which houses eight conserved Zn2+ ion-binding cysteines. To date, the structural, biochemical, and biophysical mechanisms for subcellular targeting of FYVE domains for proteins from various model organisms have been worked out but plant FYVE domains remain noticeably under-investigated.

RESULTS

We carried out an extensive examination of all Arabidopsis FYVE domains, including their identification, classification, molecular modeling and biophysical characterization using computational approaches. Our classification of fifteen Arabidopsis FYVE proteins at the outset reveals unique domain architectures for FYVE containing proteins, which are not paralleled in other organisms. Detailed sequence analysis and biophysical characterization of the structural models are used to predict membrane interaction mechanisms previously described for other FYVE domains and their subtle variations as well as novel mechanisms that seem to be specific to plants.

CONCLUSIONS

Our study contributes to the understanding of the molecular basis of FYVE-based membrane targeting in plants on a genomic scale. The results show that FYVE domain containing proteins in plants have evolved to incorporate significant differences from those in other organisms implying that they play a unique role in plant signaling pathways and/or play similar/parallel roles in signaling to other organisms but use different protein players/signaling mechanisms.

Statistics and physical origins of pK and ionization state changes upon protein-ligand binding

This work investigates statistical prevalence and overall physical origins of changes in charge states of receptor proteins upon ligand binding. These changes are explored as a function of the ligand type (small molecule, protein, and nucleic acid), and distance from the binding region. Standard continuum solvent methodology is used to compute, on an equal footing, pK changes upon ligand binding for a total of 5899 ionizable residues in 20 protein-protein, 20 protein-small molecule, and 20 protein-nucleic acid high-resolution complexes. The size of the data set combined with an extensive error and sensitivity analysis allows us to make statistically justified and conservative conclusions: in 60% of all protein-small molecule, 90% of all protein-protein, and 85% of all protein-nucleic acid complexes there exists at least one ionizable residue that changes its charge state upon ligand binding at physiological conditions (pH = 6.5). Considering the most biologically relevant pH range of 4-8, the number of ionizable residues that experience substantial pK changes (DeltapK > 1.0) due to ligand binding is appreciable: on average, 6% of all ionizable residues in protein-small molecule complexes, 9% in protein-protein, and 12% in protein-nucleic acid complexes experience a substantial pK change upon ligand binding. These changes are safely above the statistical false-positive noise level. Most of the changes occur in the immediate binding interface region, where approximately one out of five ionizable residues experiences substantial pK change regardless of the ligand type. However, the physical origins of the change differ between the types: in protein-nucleic acid complexes, the pK values of interface residues are predominantly affected by electrostatic effects, whereas in protein-protein and protein-small molecule complexes, structural changes due to the induced-fit effect play an equally important role. In protein-protein and protein-nucleic acid complexes, there is a statistically significant number of substantial pK perturbations, mostly due to the induced-fit structural changes, in regions far from the binding interface.

Mechanism of Cdc25B phosphatase with the small molecule substrate p-nitrophenyl phosphate from QM/MM-MFEP calculations

Cdc25B is a dual-specificity phosphatase that catalyzes the dephosphorylation of the Cdk2/CycA protein complex. This enzyme is an important regulator of the human cell cycle and has been identified as a potential anticancer target. In general, protein tyrosine phosphatases are thought to bind the dianionic form of the phosphate and employ general acid catalysis via the Asp residue in the highly conserved WPD-loop. However, the Cdc25 phosphatases form a special subfamily based on their distinct differences from other protein tyrosine phosphatases. Although Cdc25B contains the (H/V)CX(5)R catalytic motif present in all other protein tyrosine phosphatases, it lacks an analogous catalytic acid residue. No crystallographic data currently exist for the complex of Cdc25B with Cdk2/CycA, so in addition to its natural protein substrate, experimental and theoretical studies are often carried out with small molecule substrates. In an effort to gain understanding of the dephosphorylation mechanism of Cdc25B with a commonly used small molecule substrate, we have performed simulations of the rate-limiting step of the reaction catalyzed by Cdc25B with the substrate p-nitrophenyl phosphate using the recently developed QM/MM Minimum Free Energy Path method (Hu et al. J. Chem. Phys. 2008, 034105). We have simulated the first step of the reaction with both the monoanionic and the dianionic forms of the substrate, and our calculations favor a mechanism involving the monoanionic form. Thus, Cdc25 may employ a unique dephosphorylation mechanism among protein tyrosine phosphatases, at least in the case of the small molecule substrate p-nitrophenyl phosphate.

Molecular dynamics simulation of the acidic compact state of apomyoglobin from yellowfin tuna

A molecular model of the acidic compact state of apomyoglobin (A-state) from yellowfin tuna was obtained using molecular dynamics simulations (MD) by calculating multiple trajectories. To cause partial unfolding within a reasonable amount of CPU time, both an acidic environment (pH 3 and 0.15M NaCl) and a temperature jump to 500 K were needed. Twenty-five acidic structures of apomyoglobin were generated by MD, 10 of them can be clustered by RMSD in an average structure having a common hydrophobic core as was reported for acidic sperm whale apomyoglobin, with shortened helices A,G,E, and H (the helix A appears to be translated along the sequence). Prolonging the MD runs at 500 K did not cause further substantial unfolding, suggesting that the ensemble of generated structures is indicative of a region of the conformational space accessible to the apoprotein at acidic pH corresponding to a local energy minimum. The comparison of experimentally determined values of specific spectroscopic properties of the apomyoglobin in acidic salt conditions with the expected ones on the basis of the MD generated structures shows a reasonable agreement considering the characteristic uncertainties of both experimental and simulation techniques. We used frequency domain fluorometry, acrylamide fluorescence quenching, and fluorescence correlation spectroscopy together with far UV circular dichroism to estimate the helical content, the Stern-Volmer quenching constant and the radius of gyration of the protein. Tuna apomyoglobin is a single tryptophan protein and thus, interpretation of its intrinsic fluorescence is simpler than for other proteins. The high sensitivity of the applied fluorescence techniques enabled experiments to be performed under very dilute conditions, that is, at concentrations of subnanomolar for the FCS measurements and 6 muM for the other fluorescence measurements. As high concentrations of proteins can strongly affect the association equilibrium among partially unfolded states, fluorescence techniques can provide complementary information with respect to other techniques requiring higher sample concentrations, such as NMR. The analysis of exposed hydrophobic regions in each of the MD-generated acidic structures reveals potential candidates involved in the aggregation processes of apomyoglobin in the acidic compact state. Our investigation represents an effective model system for studying amyloid fibril formation found in important diseases that are believed to proceed via aggregation of protein in the molten globule state.