Serine proteases: an ab initio molecular dynamics study

In serine proteases (SPs), the H-bond between His57 and Asp102 and that between Gly193 and the transition state intermediate play a crucial role in enzymatic function. To shed light on the nature of these interactions, we have carried out ab initio molecular dynamics simulations on complexes representing adducts between the reaction intermediate and elastase (one protein belonging to the SP family). Our calculations indicate the presence of a low-barrier H-bond between His57 and Asp102, in complete agreement with NMR experiments on enzyme-transition state analogue complexes. Comparison with an ab initio molecular dynamics simulation on a model of the substrate-enzyme adduct indicates that the Gly193-induced strong stabilization of the intermediate is accomplished by charge/dipole interactions and not by H-bonding as previously suggested. Inclusion of the protein electric field in the calculations does not affect significantly the charge distribution.     

Structural phase transition of CdTe: an ab initio study

A constant pressure ab initio MD technique and density functional theory with a generalized gradient approximation (GGA) was used to study the pressure-induced phase transition in zinc-blende CdTe. We found that CdTe undergoes a structural first-order phase transition to $ {\text{I}}\overline 4 {\text{m2}} $ (binary β-tin) tetragonal structure in the constant pressure molecular dynamics simulation at 20 GPa. When the pressure was increased to 50 GPa, the phase of tetragonal structure converted to a new Imm2 orthorhombic structure. These phase transformations were also calculated by using the enthalpy calculations. Transition phases, lattice parameters and bulk properties we attained are comparable with experimental and theoretical data.     

Pressure-induced phase transition in wurtzite ZnTe: an ab initio study

A constant pressure ab initio MD technique and density functional theory with a generalized gradient approximation (GGA) was used to study the pressure-induced phase transition in wurtzite ZnTe. A first-order phase transition from the wurtzite structure to a Cmcm structure was successfully observed in a constant-pressure molecular dynamics simulation. This phase transformation was also analyzed using enthalpy calculations. We also investigated the stability of wurtzite (WZ) and zinc-blende (ZB) phases from energy–volume calculations, and found that both structures show quite similar equations of state and transform into a Cmcm structure at 16 GPa using enthalpy calculations, in agreement with experimental observations. The transition phase, lattice parameters and bulk properties we obtained are comparable with experimental and theoretical data.     

Contribution of hydrogen bonds to protein stability

Our goal was to gain a better understanding of the contribution of the burial of polar groups and their hydrogen bonds to the conformational stability of proteins. We measured the change in stability, Δ(ΔG), for a series of hydrogen bonding mutants in four proteins: villin headpiece subdomain (VHP) containing 36 residues, a surface protein from Borrelia burgdorferi (VlsE) containing 341 residues, and two proteins previously studied in our laboratory, ribonucleases Sa (RNase Sa) and T1 (RNase T1). Crystal structures were determined for three of the hydrogen bonding mutants of RNase Sa: S24A, Y51F, and T95A. The structures are very similar to wild type RNase Sa and the hydrogen bonding partners form intermolecular hydrogen bonds to water in all three mutants. We compare our results with previous studies of similar mutants in other proteins and reach the following conclusions. (1) Hydrogen bonds contribute favorably to protein stability. (2) The contribution of hydrogen bonds to protein stability is strongly context dependent. (3) Hydrogen bonds by side chains and peptide groups make similar contributions to protein stability. (4) Polar group burial can make a favorable contribution to protein stability even if the polar groups are not hydrogen bonded. (5) The contribution of hydrogen bonds to protein stability is similar for VHP, a small protein, and VlsE, a large protein.     

Structural phase transition of BeTe: an ab initio molecular dynamics study

Beryllium telluride (BeTe) with cubic zinc-blende (ZB) structure was studied using ab initio constant pressure method under high pressure. The ab initio molecular dynamics (MD) approach for constant pressure was studied and it was found that the first order phase transition occurs from the ZB structure to the nickel arsenide (NiAs) structure. It has been shown that the MD simulation predicts the transition pressure P _T more than the value obtained by the static enthalpy and experimental data. The structural pathway reveals MD simulation such as cubic → tetragonal → orthorhombic → monoclinic → orthorhombic → hexagonal, leading the ZB to NiAs phase. The phase transformation is accompanied by a 10% volume drop and at 80 GPa is likely to be around 35 GPa in the experiment. In the present study, our obtained values can be compared with the experimental and theoretical results.

Open image in new window

Graphical abstract The energy-volume relation and ZB phase for the BeTe

     

A partially buried site in homologous HPr proteins is not optimized for stability

The energetic consequences of site-specific replacement of a residue at a partially buried site in the two homologous HPr proteins from Escherichia coli and Bacillus subtilis is described. We determined previously that the replacement of a partially buried Lys residue with Glu at position 49 in E.coli HPr increased the conformational stability of the protein substantially because the side-chain of the latter residue could act as a hydrogen-bond acceptor. Here, we extend this analysis to other side-chains with different chemical properties and abilities to form hydrogen bonds to compare the properties of this position in the backgrounds of two different homologous HPr proteins. We find that the variants with polar residues that can form a tertiary hydrogen bond with a nearby site in the protein are more stable than either hydrophobic residues or polar residues that become buried yet are incapable of forming a new hydrogen bond. Furthermore, the protein with the wild-type residue in each HPr variant is not among the most stable of the proteins studied. These results suggest a general strategy for designing variants in which the overall stability of a protein can be modulated in a defined fashion.     

Effect of microsolvation on zwitterionic glycine: an ab initio and density functional theory study

The effect of microsolvation on zwitterionic glycine, considering both (-NH3(+)) as proton donor and (-COO(-)) as proton acceptor at correlated ab initio (MP2) level and density functional methods (B3LYP, PW91, MPW1PW91 and PBE) using 6-311++G** basis set has been reported. DFT methods have been employed so as to compare the performance/quality of different gradient-corrected correlation functionals (PW91, PBE), hybrid functionals (B3LYP, MPW1PW91) and to predict the near quantitative structural and vibrational properties, at reduced computational cost. B3LYP method outperforms among the different DFT methods for the computed hydrogen bond distances and found closer to the value obtained by correlated MP2 level, whereas MPW1PW91 and PBE methods shows very similar values but approximately 0.03 A less, compared to B3LYP method. MP2 calculation and single point CCSD(T)//MP2 calculation have been considered to decompose the interaction energy, including corrections for basis set superposition error (BSSE). Moreover, charge distribution analysis has also been carried out to understand the long raised questions, how and why the two body energies have significant contribution to the total binding energy.     

Structural phase transition of SnSe under uniaxial stress and hydrostatic pressure: an ab initio study

We study the structural behavior of SnSe under the hydrostatic pressure using a constant pressure ab initio technique. We find SnSe undergoes a structural second order phase transition from the orthorhombic (Pnma) structure to orthorhombic (Cmcm) structure in the constant pressure simulation at 7 GPa which is in good agreement with the recent experimental study. The Cmcm structure is fivefold coordinated. This phase transition is also analyzed from the total energy calculations. Besides, we study the behavior of SnSe under uniaxial stress.     

Theoretical ab initio study of the effects of methylation on structure and stability of G:C Watson-Crick base pair

Methylation of DNA occurs most readily at N(3), N(7), and O(6) of purine bases and N(3) and O(2) of pyrimidines. Methylated bases are continuously formed through endogenous and exogenous mechanisms. The results of a theoretical ab initio study on the methylation of G:C base pair components are reported. The geometries of the local minima were optimized without symmetry restrictions by the gradient procedure at DFT level of theory and were verified by energy second derivative calculations. The standard 6-31G(d) basis set was used. The single-point calculations have been performed at the MP2/6-31G(d,p), MP2/6-31++G(d,p), and MP2/6-311++G(2d,2p) levels of theory. The geometrical parameters, relative stability and counterpoise corrected interaction energies are reported. Also, using a variation-perturbation energy decomposition scheme we have found the vital contributions to the total interaction energy.     

NADH interactions with WT- and S94A-acyl carrier protein reductase from Mycobacterium tuberculosis: an ab initio study

We present an ab initio molecular dynamics study of the complex between acyl carrier protein reductase InhA from M. tuberculosis and isonicotinic acid hydrazide-NADH. We focus on wild-type (WT) InhA and a mutant causing drug resistance (S94A) for which structural information is available (Rozwarski et al., 1998;279:98--102; Dessen et al., 1995;267:1638--1641). Our calculations suggest that the water-mediated H-bond interactions between Ser94 side chain and NADH, present in WT InhA X-ray structure, can be lost during the dynamics. This conformational change is accompanied by a structural rearrangement of Gly14. The calculated structure of WT is rather similar to the X-ray structure of the S94A mutant in terms of geometrical parameters and chemical bonding. Further evidence for the mobility of Ser94 is provided by a 1-ns-long classical molecular dynamics on the entire protein. The previously unrecognized high mobility of Ser94 can provide a rationale of the small change in free binding energies on passing from WT to S94A InhA.     

3- and 5-isoxazolol zwitterions: an ab initio molecular orbital study relating to GABA agonism and antagonism

The Hartree-Fock ab initio molecular orbital method has been applied to eight compounds: GABA (gamma-amino butyric acid) (1), its partially rigidified analog, TACA (trans-4-aminocrotonic acid) (2), six isoxazolol analogs; muscimol (5-aminomethylisoxazol-3-ol (3), THIP (4,5,6,7-tetrahydroisoxazolo[5,4-c]pyridin-3-ol) (4), THAZ (5,6,7,8-tetrahydro-4H-isoxazolo[4,5-d]azepin-3-ol) (5), isomuscimol (3-aminomethylisoxazol-5-ol) (6), iso-THIP (4,5,6,7-tetrahydroisoxazolo[3,4-c] pyridin-5-ol) (7), and iso-THAZ (5,6,7,8-tetrahydro-4H-isoxazolo[3,4-d]azepin-5-ol) (8). GABA is an endogenous inhibitory transmitter. The four following molecules (2), (3), (4) and (5) are agonist: they bind themselves to the GABA receptors and induce approximately the same effect as GABA. (6) is lightly agonist, presenting a lower affinity. Compounds (7) and (8) are antagonists, giving rise to convulsion. Optimized molecular conformations of GABA (1), muscimol (3) and isomuscimol (6) are discussed. Geometric and electronic parameters showing the presence of intramolecular hydrogen bonds are presented. The permutation of the heteroatoms in the isoxazole ring has no effect on the side-chain orientation explaining maybe the agonist character of isomuscimol, being able to adopt easily and exactly the active conformation. Atomic charge distributions and electronic overlap populations for all compounds have been computed in order to try to understand why their GABAergic activities can be so different. The computed values show that the 3-isoxazolol ring mimics in a good way the carboxylic function of GABA. They also illustrate the larger electronic delocalization within the 5-isoxazolol ring and therefore the resulting antagonist character, except for isomuscimol.     

A Birch-like mechanism in enzymatic benzoyl-CoA reduction: a kinetic study of substrate analogues combined with an ab initio model

Benzoyl-CoA reductase from the anaerobic bacterium Thauera aromatica catalyzes the ATP-driven two-electron reduction of the aromatic moiety of benzoyl-CoA. A Birch mechanism involving alternate one-electron and one-proton transfer steps to the aromatic ring was previously proposed for benzoyl-CoA reductase. Due to the high redox barrier, the first electron transfer step yielding a radical anion is considered the rate-limiting step in this reaction. Focusing on the mechanism of substrate reduction, this work combines the kinetic analysis of a number of substrate analogues with a model based on the ab initio calculated electron density of the radical anion of benzoyl-CoA, a transition state model of the proposed Birch mechanism. Both K(m) and k(cat) of ortho-substituted benzoyl-CoA increased in parallel with the substituent's acceptor strength (F > Cl = H > OH > NH(2)). Among the isomers of monofluorobenzoyl-CoA, reduction rates decreased in the following order: ortho > meta > para; the K(m) values increased in the following order: meta > ortho > para. Five-ring heteroaromatic acid thiol esters were reduced in the following order: thiophene > furan > pyrrole; the 2-isomers are reduced much faster than the 3-isomers. Most of these results could be rationalized by the model. A Hammett plot indicated that the reaction mechanism is only slightly polar, suggesting the involvement of a partial protonation of the carbonyl oxygen of benzoyl-CoA and/or a simultaneous transfer of the first electron and proton. Surprisingly, benzoyl-CoA reductase exhibited a hydrogen kinetic isotope effect on k(cat) for pyridine-2-carbonyl-CoA (2.1) but only a negligible one for benzoyl-CoA (1.2), indicating that pyridine-2-carbonyl-CoA reduction proceeds according to a varied mechanism.     

Structure,stability and water permeation of ([D-Leu-L-Lys-(D-Gln-L-Ala)3]) cyclic peptide nanotube: a molecular dynamics study

The structural stability of 8 × ([D-Leu-L-Lys-(D-Gln-L-Ala)₃]) cyclic peptide nanotube (CPN) in water and different phospholipid bilayers were explored by 100 ns independent molecular dynamics (MD) simulations. The role of non-bonded interaction energy between the side and main chains of cyclic peptide rings in different membrane environments assessed, wherein the repulsive electrostatic interaction energy between neighbouring cyclic peptide rings was found adequate to break hydrogen bond energy thereby to crumple CPN. Further, the water permeation across the CPN channel was studied in four types of phospholipid bilayers- DMPG (1,2-Dimyristoyl-sn-glycero-3-phosphorylglycerol), DMPS (1,2-Dimyristoyl-sn-glycero-3-phosphoserine), POPC (1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) and POPE (1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine) from MD simulations. DMPS membrane shows higher non-bonded interaction energies (?1913.06 kJ/mol of electrostatic interaction energy and ?994.13 kJ/mol of van der Waals interaction energy) with CPN due to the presence of polar molecules in lipid structure. Thusly, the non-bonded interaction energies were essential towards the stability of CPN than hydrogen bonds between the nearby cyclic peptides. The result also reveals the role of side chains, hydrogen bonds and non-bonded interaction energies in an aqueous environment. The diffusion coefficient of water obtained from means square deviation calculation shows similar coefficients irrespective of the lipid surroundings. However, the permeation coefficients demonstrate water flow in the channel relies upon the environment.     

Asp79 makes a large, unfavorable contribution to the stability of RNase Sa

The two most buried carboxyl groups in ribonuclease Sa (RNase Sa) are Asp33 (99% buried; pK 2.4) and Asp79 (85% buried; pK 7.4). Above these pK values, the stability of the D33A variant is 6kcal/mol less than wild-type RNase Sa, and the stability of the D79A variant is 3.3kcal/mol greater than wild-type RNase Sa. The key structural difference between the carboxyl groups is that Asp33 forms three intramolecular hydrogen bonds, and Asp79 forms no intramolecular hydrogen bond. Here, we focus on Asp79 and describe studies of 11 Asp79 variants. Most of the variants were at least 2kcal/mol more stable than wild-type RNase Sa, and the most interesting was D79F. At pH 3, below the pK of Asp79, RNase Sa is 0.3kcal/mol more stable than the D79F variant. At pH 8.5, above the pK of Asp79, RNase Sa is 3.7kcal/mol less stable than the D79F variant. The unfavorable contribution of Asp79 to the stability appears to result from the Born self-energy of burying the charge and, more importantly, from unfavorable charge-charge interactions. To counteract the effect of the negative charge on Asp79, we prepared the Q94K variant and the crystal structure showed that the amino group of the Lys formed a hydrogen-bonded ion pair (distance, 2.71A; angle, 100 degrees ) with the carboxyl group of Asp79. The stability of the Q94K variant was about the same as the wild-type at pH 3, where Asp79 is uncharged, but 1kcal/mol greater than that of wild-type RNase Sa at pH 8.5, where Asp79 is charged. Differences in hydrophobicity, steric strain, Born self-energy, and electrostatic interactions all appear to contribute to the range of stabilities observed in the variants. When it is possible, replacing buried, non-hydrogen bonded, ionizable side-chains with non-polar side-chains is an excellent means of increasing protein stability.     

基于密度泛函理论研究Watson–Crick碱基对中的氢键特征(英文)

基于密度泛函理论(DFT)研究腺嘌呤、胸腺嘧啶、鸟嘌呤、胞嘧啶以及腺嘌呤胸腺嘧啶碱基对、鸟嘌呤胞嘧啶碱基对。在DFT-B3LYP/6-31G**水平上利用自然键轨道理论分析研究结果显示,互补碱基对的结构和电子特征有利于氢键的形成。本文中讨论几何结构、电子结构、分子轨道和能量对于氢键形成的影响。此研究结果将有助于更好的理解AT和GC碱基对中氢键与它们的结构特性之间的关系。     

The contribution of cross-links to protein stability: A normal mode analysis of the configurational entropy of the native state

The vibrational entropy of native BPTI, with three disulfide bonds, was determined by use of normal mode calculations and compared with that of folded variants having either one less disulfide bond or lacking a peptide bond at the trypsin-reactive site. Favorable contributions to the free energy of 2.5–5.1 kcal/mol at 300 K were calculated for the reduction of disulfide bonds in the folded state, whereas no favorable contribution was found for the hydrolysis of the peptide bond cleaved by trypsin. This is on the order of the effect of disulfides in the unfolded state. The implications of these results for the stabilization of a folded protein by the introduction of crosslinks are discussed. © 1993 Wiley-Liss, Inc.     

Structural and electronic properties of diazonium functionalized (4, 4) single walled carbon nanotube: an ab initio study

In this work, ab initio density functional theory (DFT) calculations are performed to study the structural and electronic properties of diazonium reagent functionalized (4, 4) single-walled carbon nanotube (SWCNT). We find the aryl group covalently bonds with SWCNT and prefers to be perpendicular to the side wall of nanotube. It has a rotational barrier of 0.35 eV around the formed aryl-tube bond axis and should be thermodynamically stable at room temperature. Additionally, new peaks appeared around the Fermi energy in the density of state (DOS) due to the weak band dispersion. Increasing of the coverage of the functional group will result in significant upshift of the Fermi level.