Protein effects in non-heme iron enzyme catalysis: insights from multiscale models

Many non-heme iron enzymes have similar sets of ligands but still catalyze widely different reactions. A key question is, therefore, the role of the protein in controlling reactivity and selectivity. Examples from multiscale simulations, primarily QM/MM, of both mono- and binuclear non-heme iron enzymes are used to analyze the stability of these models and what they reveal about the protein effects. Consistent results from QM/MM modeling are the importance of the hydrogen bond network to control reactivity and electrostatic stabilization of electron transfer from second-sphere residues. The long-range electrostatic effects on reaction barriers are small for many systems. In the systems where large electrostatic effects have been reported, these lead to higher barriers. There is thus no evidence of any significant long-range electrostatic effects contributing to the catalytic efficiency of non-heme iron enzymes. However, the correct evaluation of electrostatic contributions is challenging, and the correlation between calculated residue contributions and the effects of mutation experiments is not very strong. The largest benefits of QM/MM models are thus the improved active-site geometries, rather than the calculation of accurate energies. Reported differences in mechanistic predictions between QM and QM/MM models can be explained by differences in hydrogen bonding patterns in and around the active site. Correctly constructed cluster models can give results with similar accuracy as those from multiscale models, but the latter reduces the risk of drawing the wrong mechanistic conclusions based on incorrect geometries and are preferable for all types of modeling, even when using very large QM parts.     

Anesthetic interaction with ketosteroid isomerase: insights from molecular dynamics simulations

The nature and the sites of interactions between anesthetic halothane and homodimeric Delta5-3-ketosteroid isomerase (KSI) are characterized by flexible ligand docking and confirmed by 1H-15N NMR. The dynamics consequence of halothane interaction and the implication of the dynamic changes to KSI function are studied by multiple 5-ns molecular dynamics simulations in the presence and absence of halothane. Both docking and MD simulations show that halothane prefer the amphiphilic dimeric interface to the hydrophobic active site of KSI. Halothane occupancy at the dimer interface disrupted the intersubunit hydrogen bonding formed either directly through side chains of polar residues or indirectly through the mediation of the interfacial water molecules. Moreover, in the presence of halothane, the exchange rate of the bound waters with bulk water was increased. Halothane perturbation to the dimer interface affected the overall flexibility of the active site. This action is likely to contribute to the halothane-induced reduction of the KSI activity. The allosteric halothane modulation of the dynamics-function relationship of KSI without direct competition at the enzymatic active sites may be generalized to offer a unifying explanation of anesthetic action on a diverse range of multidomain neuronal proteins that are potentially relevant to clinical general anesthesia.     

Molecular dynamics insights into protein-glycosaminoglycan systems from microsecond-scale simulations

Heparin is a key player in cell signaling via its physical interactions with protein targets in the extracellular matrix. However, basic molecular level understanding of these highly biologically relevant intermolecular interactions is still incomplete. In this study, for the first time, microsecond-scale MD simulations are reported for a complex between fibroblast growth factor 1 and heparin. We rigorously analyze this molecular system in terms of the conformational space, structural, energetic, and dynamic characteristics. We reveal that the conformational selection mechanism of binding denotes a recognition specificity determinant. We conclude that the length of the simulation could be crucial for evaluation of some of the analyzed parameters. Our data provide novel significant insights into the interactions in the fibroblast growth factor 1 complex with heparin, in particular, and into the physical-chemical nature of protein-glycosaminoglycan systems in general, which have potential applicability for biomaterials development in the area of regenerative medicine.     

Sequence dependence of amyloid fibril formation: insights from molecular dynamics simulations

The clarification of the physico-chemical determinants underlying amyloid deposition is critical for our understanding of misfolding diseases. With this purpose we have performed a systematic all-atom molecular dynamics (MD) study of a series of single point mutants of the de novo designed amyloidogenic peptide STVIIE. Sixteen different 50ns long simulations using explicit solvent have been carried out starting from four different conformations of a polymeric six-stranded beta-sheet. The simulations have provided evidence for the influence of a small number of site-specific hydrophobic interactions on the packing and stabilization of nascent aggregates, as well as the interplay between side-chain interactions and the net charge of the molecule on the strand arrangement of polymeric beta-sheets. This MD analysis has also shed light into the origin of the position dependence on mutation of beta-sheet polymerization that was found experimentally for this model system. Our results suggest that MD can be applied to detect critical positions for beta-sheet aggregation within a given amyloidogenic stretch. Studies similar to the one presented here can guide site-directed mutations or the design of drugs that specifically disrupt the key stabilizing interactions of beta-sheet aggregates.     

Protein structure and dynamics in nonaqueous solvents: insights from molecular dynamics simulation studies

Protein structure and dynamics in nonaqueous solvents are here investigated using molecular dynamics simulation studies, by considering two model proteins (ubiquitin and cutinase) in hexane, under varying hydration conditions. Ionization of the protein groups is treated assuming "pH memory," i.e., using the ionization states characteristic of aqueous solution. Neutralization of charged groups by counterions is done by considering a counterion for each charged group that cannot be made neutral by establishing a salt bridge with another charged group; this treatment is more physically reasonable for the nonaqueous situation, contrasting with the usual procedures. Our studies show that hydration has a profound effect on protein stability and flexibility in nonaqueous solvents. The structure becomes more nativelike with increasing values of hydration, up to a certain point, when further increases render it unstable and unfolding starts to occur. There is an optimal amount of water, approximately 10% (w/w), where the protein structure and flexibility are closer to the ones found in aqueous solution. This behavior can explain the experimentally known bell-shaped dependence of enzyme catalysis on hydration, and the molecular reasons for it are examined here. Water and counterions play a fundamental and dynamic role on protein stabilization, but they also seem to be important for protein unfolding at high percentages of bound water.     

Regulatory phosphorylation of cyclin-dependent kinase 2: insights from molecular dynamics simulations

The structures of fully active cyclin-dependent kinase-2 (CDK2) complexed with ATP and peptide substrate, CDK2 after the catalytic reaction, and CDK2 inhibited by phosphorylation at Thr14/Tyr15 were studied using molecular dynamics (MD) simulations. The structural details of the CDK2 catalytic site and CDK2 substrate binding box were described. Comparison of MD simulations of inhibited complexes of CDK2 was used to help understand the role of inhibitory phosphorylation at Thr14/Tyr15. Phosphorylation at Thr14/Tyr15 causes ATP misalignment for the phosphate-group transfer, changes in the Mg²⁺ coordination sphere, and changes in the H-bond network formed by CDK2 catalytic residues (Asp127, Lys129, Asn132). The inhibitory phosphorylation causes the G-loop to shift from the ATP binding site, which leads to opening of the CDK2 substrate binding box, thus probably weakening substrate binding. All these effects explain the decrease in kinase activity observed after inhibitory phosphorylation at Thr14/Tyr15 in the G-loop. Interaction of the peptide substrate, and the phosphorylated peptide product, with CDK2 was also studied and compared. These results broaden hypotheses drawn from our previous MD studies as to why a basic residue (Arg/Lys) is preferred at the P₊₂ substrate position. Figure View of the substrate binding site of the fully active cyclin-dependent kinase-2 (CDK2) (pT160-CDK2/cyclin A/ATP). The pThr160 activation site is located in the T-loop (yellow secondary structure). The G-loop, which partly forms the ATP binding site, is shown in blue. The Thr14 and Tyr15 inhibitory phosphorylation sites located in the G-loop are shown in licorice representation     

Protein flexibility using constraints from molecular dynamics simulations

Proteins are held together in the native state by hydrophobic interactions, hydrogen bonds and interactions with the surrounding water, whose strength as well as spatial and temporal distribution affects protein flexibility and hence function. We study these effects using 10 ns molecular dynamics simulations of pure water and of two proteins, the glutamate receptor ligand binding domain and barnase. We find that most of the noncovalent interactions flicker on and off over typically nanoseconds, and so we can obtain good statistics from the molecular dynamics simulations. Based on this information, a topological network of rigid bonds corresponding to a protein structure with covalent and noncovalent bonds is constructed, with account being taken of the influence of the flickering hydrogen bonds. We define the duty cycle for the noncovalent interactions as the percentage of time a given interaction is present, which we use as an input to investigate flexibility/rigidity patterns, in the algorithm FIRST which constructs and analyses topological networks.     

Protein dynamics simulations from nanoseconds to microseconds.

There have been a number of advances in atomic resolution simulations of biomolecules during the past few years. These have arisen partly from improvements to computer power and partly from algorithmic improvements. There have also been advances in measuring time-dependent fluctuations in proteins using NMR spectroscopy, revealing the importance of fluctuations in the microsecond to millisecond time range. Progress has also been made in measuring how far the simulations are able to represent the accessible phase space that is available to the protein in its native state, in solution, at room temperature. Another area of development is the simulation of protein unfolding at atomic resolution.     

Structure and dynamics of Mycobacterium tuberculosis truncated hemoglobin N: insights from NMR spectroscopy and molecular dynamics simulations

The potent nitric oxide dioxygenase (NOD) activity (trHbN-Fe2?-O? + (?)NO → trHbN-Fe3?-OH? + NO??) of Mycobacterium tuberculosis truncated hemoglobin N (trHbN) protects aerobic respiration from inhibition by (?)NO. The high activity of trHbN has been attributed in part to the presence of numerous short-lived hydrophobic cavities that allow partition and diffusion of the gaseous substrates (?)NO and O? to the active site. We investigated the relation between these cavities and the dynamics of the protein using solution NMR spectroscopy and molecular dynamics (MD). Results from both approaches indicate that the protein is mainly rigid with very limited motions of the backbone N-H bond vectors on the picoseconds-nanoseconds time scale, indicating that substrate diffusion and partition within trHbN may be controlled by side-chains movements. Model-free analysis also revealed the presence of slow motions (microseconds-milliseconds), not observed in MD simulations, for many residues located in helices B and G including the distal heme pocket Tyr33(B10). All currently known crystal structures and molecular dynamics data of truncated hemoglobins with the so-called pre-A N-terminal extension suggest a stable α-helical conformation that extends in solution. Moreover, a recent study attributed a crucial role to the pre-A helix for NOD activity. However, solution NMR data clearly show that in near-physiological conditions these residues do not adopt an α-helical conformation and are significantly disordered and that the helical conformation seen in crystal structures is likely induced by crystal contacts. Although this lack of order for the pre-A does not disagree with an important functional role for these residues, our data show that one should not assume an helical conformation for these residues in any functional interpretation. Moreover, future molecular dynamics simulations should not use an initial α-helical conformation for these residues in order to avoid a bias based on an erroneous initial structure for the N-termini residues. This work constitutes the first study of a truncated hemoglobin dynamics performed by solution heteronuclear relaxation NMR spectroscopy.     

Protein dynamics and function: insights from the energy landscape and solvent slaving

Protein motions are complex and a good way to describe them is in terms of a very high-dimensional conformation space. We give here a simple explanation of the conformation space and the energy landscape, the conformational motions and protein reactions, based on an analogy to a traffic problem. The analogy provides insight into the slaving of protein processes to bulk solvent fluctuations, in both the native and unfolded states.     

Folding and stability of the three-stranded beta-sheet peptide Betanova: insights from molecular dynamics simulations

The dynamics of the three-stranded beta-sheet peptide Betanova has been studied at four different temperatures (280, 300, 350, and 450 K by molecular dynamics simulation techniques, in explicit water. Two 20-ns simulations at 280 K indicate that the peptide remains very flexible under "folding" conditions sampling a range of conformations that together satisfy the nuclear magnetic resonance (NMR)-derived experimental constraints. Two simulations at 300 K (above the experimental folding temperature) of 20 ns each show partial formation of "native"-like structure, which also satisfies most of the NOE constraints at 280 K. At higher temperature, the presence of compact states, in which a series of hydrophobic contacts remain present, are observed. This is consistent with experimental observations regarding the role of hydrophobic contacts in determining the peptide's stability and in initiating the formation of turns and loops. A set of different structures is shown to satisfy NMR-derived distance restraints and a possible mechanism for the folding of the peptide into the NMR-determined structure is proposed.     

The pH dependence of flavivirus envelope protein structure: insights from molecular dynamics simulations

The flavivirus membrane fusion is triggered by the acid pH of the endosomes after virus endocytosis. The proposed mechanism involves changes in the protonation state of conserved histidine residues of the E protein present in the viral surface that undergoes a series of structural rearrangements that result in the fusion between the endosome and viral bilayers. We studied the pH dependence of E protein rearrangements of dengue virus type 2, used as a model, in the pH range experimented by the virus along the fusion process. We employed a low computational cost scheme to explore the behavior of the E protein by molecular dynamics (MD) simulations of complete systems that include the protein, the solvent, and ions. The procedure alternates cyclically the update of the ionization states of the protein residues with common MD steps applied to the new ionization configuration. Important pH-dependent protein structure rearrangements consistent with the changes of the protonation states of conserved histidine residues were observed. The involvement of other conserved residues in the flavivirus in the rearrangements was also identified. The results show interesting correlations with a proposed model for the fusion mechanism, as well as the experimentally identified key residues, contributing to a better understanding of the structural changes in protein E that lead to the fusion process.     

Investigating the mechanism of peptide aggregation: insights from mixed monte carlo-molecular dynamics simulations

The early stages of peptide aggregation are currently not accessible by experimental techniques at atomic resolution. In this article, we address this problem through the application of a mixed simulation scheme in which a preliminary coarse-grained Monte Carlo analysis of the free-energy landscape is used to identify representative conformations of the aggregates and subsequent all-atom molecular dynamics simulations are used to analyze in detail possible pathways for the stabilization of oligomers. This protocol was applied to systems consisting of multiple copies of the model peptide GNNQQNY, whose detailed structures in the aggregated state have been recently solved in another study. The analysis of the various trajectories provides dynamical and structural insight into the details of aggregation. In particular, the simulations suggest a hierarchical mechanism characterized by the initial formation of stable parallel β-sheet dimers and identify the formation of the polar zipper motif as a fundamental feature for the stabilization of initial oligomers. Simulation results are consistent with experimentally derived observations and provide an atomically detailed view of the putative initial stages of fibril formation.     

New insights into the structure of abasic DNA from molecular dynamics simulations

Abasic (AP) sites constitute a common form of DNA damage, arising from the spontaneous or enzymatic breakage of the N-glycosyl bond and the loss of a nucleotide base. To examine the effects of such damage on DNA structure, especially in the vicinity of the abasic sugar, four 1.5 ns molecular dynamics simulations of double-helical DNA dodecamers with and without a single abasic (tetrahydrofuran, X) lesion in a 5′-d(CXT) context have been performed and analyzed. The results indicate that the abasic site does not maintain a hole or gap in the DNA, but instead perturbs the canonical structure and induces additional flexibility close to the abasic site. In the apurinic simulations (i.e., when a pyrimidine is opposite the AP site), the abasic sugar flipped in and out of the minor groove, and the gap was water filled, except during the occurrence of a novel non-Watson–Crick C-T base pair across the abasic site. The apyrimidinic gap was not penetrated by water until the abasic sugar flipped out and remained extrahelical. Both AP helices showed kinks of 20–30° at the abasic site. The Watson–Crick hydrogen bonds are more transient throughout the DNA double helices containing an abasic site. The abasic sugar displayed an unusually broad range of sugar puckers centered around the northern pucker. The increased motion of the bases and backbone near the abasic site appear to correlate with sequence-dependent helical stability. The data indicate that abasic DNA contorts more easily and in specific ways relative to unmodified DNA, an aspect likely to be important in abasic site recognition and hydrolysis.     

Self-assembly of cationic surfactants on the carbon nanotube surface: insights from molecular dynamics simulations

The insolubility of carbon nanotubes (CNTs) in aqueous media has been a limitation for the practical application of this unique material. Recent studies have demonstrated that the suspend ability of CNT can be substantially improved by employing appropriate surfactants. Although various surfactants have been tested, the exact mechanism by which carbon nanotubes and the different surfactants interact is not fully understood. To deepen the understanding of molecular interaction between CNT and surfactants, as well as to investigate the influence of the surfactant tail length on the adsorption process, we report here the first detailed large-scale all-atomistic molecular dynamics simulation study of the adsorption and morphology of aggregates of the cationic surfactants containing trimethylammonium headgroups (C₁₂TAB and C₁₆TAB) on single-walled carbon nanotube (SWNT) surfaces. We find that the aggregation morphology of both C₁₂TAB and C₁₆TAB on the SWNT is dependent upon the number of the surfactants in the simulation box. As the number of the surfactants increases the random monolayer structure gradually changes to the cylinder-like monolayer structure. Moreover, we make a comparison between the C₁₂TAB and C₁₆TAB adsorption onto SWNTs to clarify the role of the surfactant tail length on the adsorption process. This comparison indicates that by increasing the number of surfactant molecules, the larger number of the C₁₆TAB molecules tend to adsorb onto SWNTs. Further, our results show that a longer chain yields the higher packed aggregates in which the surfactant heads are extended far into the aqueous phase, which in turn may increase the SWNTs stabilization in aqueous suspensions.     

Molecular insights of SAH enzyme catalysis and implication for inhibitor design

Biological transmethylation reaction is a key step in the duplication of virus life cycle, in which S-adenosylmethionine plays as the methyl donor. The product of this reactions, S-adenosylhomocysteine (AdoHcy) inhibits the transmethylation process. AdoHcy is hydrolysed to adenosine and L-homocysteine by the action of S-adenosylhomocysteine hydrolase (SAH). Thus the virus life cycle should be cut off once the action of SAH is inhibited. Our study was focussed on the discovery of potential inhibitor against SAH. We performed a similarity search in Traditional Chinese Medicine Database and retrieved 17 hits with high similarity. After that we virtually docked the 17 compounds as well as the natural substrates to the hydrolase using Autodock 3.0.1 software. Then we discussed about the mechanism of the inhibition reaction, followed by proposing the potential inhibitors by comparing best docked solutions and possible modification for the best inhibitors.     

PIK3CA somatic mutations in breast cancer: Mechanistic insights from Langevin dynamics simulations

Somatic mutations in PIK3CA (phosphatidylinositol-3 kinase, catalytic subunit, alpha isoform) are reported in breast and other human cancers to concentrate at hotspots within its kinase and helical domains. Most of these mutations cause kinase gain of function in vitro and are associated with oncogenicity in vivo. However, little is known about the mechanisms driving tumor development. We have performed computational structural studies on a homology model of wildtype PIK3CA plus recurrent H1047R, H1047L, and P539R mutations, located in the kinase and helical domains, respectively. The time evolution of the structures show that H1047R/L mutants exhibit a larger area of the catalytic cleft between the kinase N- and C-lobes compared with the wildtype that could facilitate the entrance of substrates. This larger area might yield enhanced substrate-to-product turnover associated with oncogenicity. In addition, the H1047R/L mutants display increased kinase activation loop mobility, compared with the wildtype. The P539R mutant forms more hydrogen bonds and salt-bridge interactions than the wildtype, properties that are associated with enhanced thermostability. Mutant-specific differences in the catalytic cleft and activation loop behavior suggest that structure-based mutant-specific inhibitors can be designed for PIK3CA-positive breast cancers.     

Structure and folding of disulfide-rich miniproteins: insights from molecular dynamics simulations and MM-PBSA free energy calculations

The fold of small disulfide-rich proteins largely relies on two or more disulfide bridges that are main components of the hydrophobic core. Because of the small size of these proteins and their high cystine content, the cysteine connectivity has been difficult to ascertain in some cases, leading to uncertainties and debates in the literature. Here, we use molecular dynamics simulations and MM-PBSA free energy calculations to compare similar folds with different disulfide pairings in two disulfide-rich miniprotein families, namely the knottins and the short-chain scorpion toxins, for which the connectivity has been discussed. We first show that the MM-PBSA approach is able to discriminate the correct knotted topology of knottins from the laddered one. Interestingly, a comparison of the free energy components for kalata B1 and MCoTI-II suggests that cyclotides and squash inhibitors, although sharing the same scaffold, are stabilized through different interactions. Application to short-chain scorpion toxins suggests that the conventional cysteine pairing found in many homologous toxins is significantly more stable than the unconventional pairing reported for maurotoxin and for spinoxin. This would mean that native maurotoxin and spinoxin are not at the lowest free energy minimum and might result from kinetically rather than thermodynamically driven oxidative folding processes. For both knottins and toxins, the correct or conventional disulfide connectivities provide lower flexibilities and smaller deviations from the initial conformations. Overall, our work suggests that molecular dynamics simulations and the MM-PBSA approach to estimate free energies are useful tools to analyze and compare disulfide bridge connectivities in miniproteins.     

OmpT: molecular dynamics simulations of an outer membrane enzyme

Five molecular dynamics simulations (total duration >25 ns) have been performed on the Escherichia coli outer membrane protease OmpT embedded in a dimyristoylphosphatidylcholine lipid bilayer. Globally the protein is conformationally stable. Some degree of tilt of the beta-barrel is observed relative to the bilayer plane. The greatest degree of conformational flexibility is seen in the extracellular loops. A complex network of fluctuating H-bonds is formed between the active site residues, such that the Asp210-His212 interaction is maintained throughout, whereas His212 and Asp83 are often bridged by a water molecule. This supports a catalytic mechanism whereby Asp83 and His212 bind a water molecule that attacks the peptide carbonyl. A configuration yielded by docking calculations of OmpT simulation snapshots and a model substrate peptide Ala-Arg-Arg-Ala was used as the starting point for an extended Huckel calculation on the docked peptide. These placed the lowest unoccupied molecular orbital mainly on the carbon atom of the central C=O in the scissile peptide bond, thus favoring attack on the central peptide by the water held by residues Asp83 and His212. The trajectories of water molecules reveal exchange of waters between the intracellular face of the membrane and the interior of the barrel but no exchange at the extracellular mouth. This suggests that the pore-like region in the center of OmpT may enable access of water to the active site from below. The simulations appear to reveal the presence of specific lipid interaction sites on the surface of the OmpT barrel. This reveals the ability of extended MD simulations to provide meaningful information on protein-lipid interactions.