Conformational Energy Research Articles

DFT computational correlations on conformational barriers of Zn2+ and Ni2+ chiral meso-(α,β-unsaturated)- porphyrins.

Correlations between DFT and experimental measurements on Zn2+ and Ni2+ chiral meso-(α,β-unsaturated)- porphyrins were performed using Kohn-Sham methodology. The exchange-correlation Becke88-Perdew86 functional was used in conjunction with double-zeta Slater basis sets. An accurate description of the electronic processes depending on the metal ion (Zn, Ni) or ligand (perilaldehyde and myrtenal) was made, confirming experimental results in terms of structural and electronic modifications. Moreover, this theoretical study provides a stronger knowledge and interpretation of the dynamical conformational features of the free base, Zn and Ni structures. Fundamental links between the central metallic atom and distortions of the porphyrinic core and ligands were demonstrated, in agreement with experimental data. We observed that the core in ZnPeriP and ZnMyrtP species is almost flat, in comparison with the Ni porphyrinic core, which appeared much more distorted. The type of distortion differs between PeriP and MyrtP ligands, with a combined saddled-ruffled characteristic with the former and a pronounced ruffled twisting for the latter. Finally, conformational energy barriers were extracted by spinning one of the arms in steps of 20° in a 360° dihedral angle. The resulted conformational barriers for NiPeriP or NiMyrtP are lower in energy than for ZnPeriP or ZnMyrtP, in agreement with experimental data.

Carba-cyclophellitols Are Neutral Retaining-Glucosidase Inhibitors.

The conformational analysis of glycosidases affords a route to their specific inhibition through transition-state mimicry. Inspired by the rapid reaction rates of cyclophellitol and cyclophellitol aziridine—both covalent retaining β-glucosidase inhibitors—we postulated that the corresponding carba “cyclopropyl” analogue would be a potent retaining β-glucosidase inhibitor for those enzymes reacting through the 4H3 transition-state conformation. Ab initio metadynamics simulations of the conformational free energy landscape for the cyclopropyl inhibitors show a strong bias for the 4H3 conformation, and carba-cyclophellitol, with an N-(4-azidobutyl)carboxamide moiety, proved to be a potent inhibitor (Ki = 8.2 nM) of the Thermotoga maritimaTmGH1 β-glucosidase. 3-D structural analysis and comparison with unreacted epoxides show that this compound indeed binds in the 4H3 conformation, suggesting that conformational strain induced through a cyclopropyl unit may add to the armory of tight-binding inhibitor designs.

Conformational analysis of N-methylacetamide molecule in the ground and excited electronic states

Equilibrium geometry parameters, harmonic vibrational frequencies and barriers to conformational transitions and conformer energy differences for N-methylacetamide (CH3CONHCH3) molecule in S0, S1 and T1 electronic states were calculated in MP2 (S0) and CASPT2 (S1, T1) levels of theory. It is known that there are two planar (or effectively planar) conformers of N-methylacetamide in the ground electronic state. It is proved that the electronic excitation to aforementioned electronic states resulted in both CCON and HNCC fragments becoming pyramidal and rotated around «central» CN bond. One-, two- and three-dimensional potential energy surface (PES) sections by different large amplitude motions (LAM) coordinates were obtained by MP2 (S0) and CASPT2 (S1, T1) methods with cc-pVTZ basis set. On 2D PES sections by «central» torsion (CN) and one of the CCON or HNCC fragments’ inversion coordinates there are six different minima. Using PES sections as the potentials in model Hamiltonians, different anharmonic vibrational problems were solved and the frequencies of large amplitude vibrations were estimated variationally.

Simple Design of an Enzyme-Inspired Supported Catalyst Based on a Catalytic Triad

Summary Enzyme active sites afford an intricate interplay of functional groups to mediate complex organic and inorganic reactions. Many hydrolytic enzymes use a catalytic triad comprising three different functional residues—(Ser(-OH), Hist(-imidazole), Asp(-CO 2 H))—that catalyze the hydrolysis of numerous unique substrates. Inspired by this design, we have developed a simple one-step synthesis for preparing a new supported catalytic system in which the three reactive groups of the catalytic triad (alcohol, imidazole, and carboxylate) are incorporated into a single functional unit. These artificial active sites can be coupled to a solid-phase support (Merrifield resin) by copper(I)-catalyzed azide-alkyne cycloaddition "click chemistry," and their effectiveness as esterolysis catalysts was demonstrated. Furthermore, tuning the local hydrophobicity of the resin particles with an approach analogous to the native enzyme hydrophobic pocket increased the catalytic efficiency. Quantum mechanics and molecular dynamics computational modeling were used to probe the catalytic effect and suggested a concerted two-step mechanism and hydrophobic nanoenvironment similar to that of hydrolytic enzymes.

Fast de novo discovery of low-energy protein loop conformations.

In the prediction of protein structure from amino acid sequence, loops are challenging regions for computational methods. Since loops are often located on the protein surface, they can have significant roles in determining protein functions and binding properties. Loop prediction without the aid of a structural template requires extensive conformational sampling and energy minimization, which are computationally difficult. In this article we present a new de novo loop sampling method, the Parallely filtered Energy Targeted All-atom Loop Sampler (PETALS) to rapidly locate low energy conformations. PETALS explores both backbone and side-chain positions of the loop region simultaneously according to the energy function selected by the user, and constructs a nonredundant ensemble of low energy loop conformations using filtering criteria. The method is illustrated with the DFIRE potential and DiSGro energy function for loops, and shown to be highly effective at discovering conformations with near-native (or better) energy. Using the same energy function as the DiSGro algorithm, PETALS samples conformations with both lower RMSDs and lower energies. PETALS is also useful for assessing the accuracy of different energy functions. PETALS runs rapidly, requiring an average time cost of 10 minutes for a length 12 loop on a single 3.2 GHz processor core, comparable to the fastest existing de novo methods for generating an ensemble of conformations. Proteins 2017; 85:1402-1412. © 2017 Wiley Periodicals, Inc.

Polarizable Force Field for DNA Based on the Classical Drude Oscillator: I. Refinement Using Quantum Mechanical Base Stacking and Conformational Energetics.

Empirical force fields seek to relate the configuration of a set of atoms to its energy, thus yielding the forces governing its dynamics, using classical physics rather than more expensive quantum mechanical calculations that are computationally intractable for large systems. Most force fields used to simulate biomolecular systems use fixed atomic partial charges, neglecting the influence of electronic polarization, instead making use of a mean-field approximation that may not be transferable across environments. Recent hardware and software developments make polarizable simulations feasible, and to this end, polarizable force fields represent the next generation of molecular dynamics simulation technology. In this work, we describe the refinement of a polarizable force field for DNA based on the classical Drude oscillator model by targeting quantum mechanical interaction energies and conformational energy profiles of model compounds necessary to build a complete DNA force field. The parametrization strategy employed in the present work seeks to correct weak base stacking in A- and B-DNA and the unwinding of Z-DNA observed in the previous version of the force field, called Drude-2013. Refinement of base nonbonded terms and reparametrization of dihedral terms in the glycosidic linkage, deoxyribofuranose rings, and important backbone torsions resulted in improved agreement with quantum mechanical potential energy surfaces. Notably, we expand on previous efforts by explicitly including Z-DNA conformational energetics in the refinement.

The effect of mechanical loads on the degradation of aliphatic biodegradable polyesters.

Aliphatic biodegradable polyesters have been the most widely used synthetic polymers for developing biodegradable devices as alternatives for the currently used permanent medical devices. The performances during biodegradation process play crucial roles for final realization of their functions. Because physiological and biochemical environment in vivo significantly affects biodegradation process, large numbers of studies on effects of mechanical loads on the degradation of aliphatic biodegradable polyesters have been launched during last decades. In this review article, we discussed the mechanism of biodegradation and several different mechanical loads that have been reported to affect the biodegradation process. Other physiological and biochemical factors related to mechanical loads were also discussed. The mechanical load could change the conformational strain energy and morphology to weaken the stability of the polymer. Besides, the load and pattern could accelerate the loss of intrinsic mechanical properties of polymers. This indicated that investigations into effects of mechanical loads on the degradation should be indispensable. More combination condition of mechanical loads and multiple factors should be considered in order to keep the degradation rate controllable and evaluate the degradation process in vivo accurately. Only then can the degradable devise achieve the desired effects and further expand the special applications of aliphatic biodegradable polyesters.

Mapping the conformational free energy of aspartic acid in the gas phase and in aqueous solution

The conformational free energy landscape of aspartic acid, a proteogenic amino acid involved in a wide variety of biological functions, was investigated as an example of the complexity that multiple rotatable bonds produce even in relatively simple molecules. To efficiently explore such a landscape, this molecule was studied in the neutral and zwitterionic forms, in the gas phase and in water solution, by means of molecular dynamics and the enhanced sampling method metadynamics with classical force-fields. Multi-dimensional free energy landscapes were reduced to bi-dimensional maps through the non-linear dimensionality reduction algorithm sketch-map to identify the energetically stable conformers and their interconnection paths. Quantum chemical calculations were then performed on the minimum free energy structures. Our procedure returned the low energy conformations observed experimentally in the gas phase with rotational spectroscopy [M. E. Sanz et al., Phys. Chem. Chem. Phys. 12, 3573 (2010)]. Moreover, it provided information on higher energy conformers not accessible to experiments and on the conformers in water. The comparison between different force-fields and quantum chemical data highlighted the importance of the underlying potential energy surface to accurately capture energy rankings. The combination of force-field based metadynamics, sketch-map analysis, and quantum chemical calculations was able to produce an exhaustive conformational exploration in a range of significant free energies that complements the experimental data. Similar protocols can be applied to larger peptides with complex conformational landscapes and would greatly benefit from the next generation of accurate force-fields.

A Study about Regioisomeric Hydroquinones with Multiple Intramolecular Hydrogen Bonding.

A theoretical exploration about hydrogen bonding in a series of synthetic regioisomeric antitumor tricyclic hydroquinones is presented. The stabilization energy for the intramolecular hydrogen bond (IHB) formation in four structurally different situations were evaluated: (a) IHB between the proton of a phenolic hydroxyl group and an ortho-carbonyl group (forming a six-membered ring); (b) between the oxygen atom of a phenolic hydroxyl group and the proton of an hydroxyalkyl group (seven membered ring); (c) between the proton of a phenolic hydroxyl group with the oxygen atom of the hydroxyl group of a hydroxyalkyl moiety (seven-membered ring); and (d) between the proton of a phenolic hydroxyl group and an oxygen atom directly bonded to the aromatic ring in ortho position (five-membered ring). A conformational analysis for the rotation around the hydroxyalkyl substituent is also performed. It is observed that there is a correspondence between the conformational energies and the IHB. The strongest intramolecular hydrogen bonds are those involving a phenolic proton and a carbonyl oxygen atom, forming a six-membered ring, and the weakest are those involving a phenolic proton with the oxygen atom of the chromenone, forming five-membered rings. Additionally, the synthesis and structural assignment of two pairs of regioisomeric hydroquinones, by 2D-NMR experiments, are reported. These results can be useful in the design of biologically-active molecules.

Nucleotide Dependent Switching in Rho GTPase: Conformational Heterogeneity and Competing Molecular Interactions

Ras superfamily of GTPases regulate myriad cellular processes through a conserved nucleotide (GTP/GDP) dependent switching mechanism. Unlike Ras family of GTPases, for the Rho GTPases, there is no clear evidence for the existence of “sub-states” such as state 1 & state 2 in the GTP bound form. To explore the nucleotide dependent conformational space of the Switch I loop and also to look for existence of state 1 like conformations in Rho GTPases, atomistic molecular dynamics and metadynamics simulations on RhoA were performed. These studies demonstrate that both the nucleotide-free state and the GDP bound “OFF” state have very similar conformations, whereas the GTP bound “ON” state has unique conformations with signatures of two intermediate states. The conformational free energy landscape for these systems suggests the presence of multiple intermediate states. Interestingly, the energetic penalty of exposing the non-polar residues in the GTP bound form is counter balanced by the favourable hydrogen bonded interactions between the γ-phosphate group of GTP with the highly conserved Tyr34 and Thr37 residues. These competing molecular interactions lead to a tuneable energy landscape of the Switch I conformation, which can undergo significant changes based on the local environment including changes upon binding to effectors.

Computational Approaches to Understanding the Biological Behaviour of Intrinsically Disordered Proteins

Intrinsically disordered proteins (IDPs) represent a class of proteins that lack a persistent folded conformation and exist as dynamic ensembles in their native state. Inherent lack of a well-defined structure and remarkable structural plasticity have facilitated their functioning in a wide range of crucial cellular processes such as signalling transduction and cell cycle regulation as well as responsible for their aberrant toxic amyloidogenic conformations implicated in a wide range of neurodegenerative diseases, cancer, etc. Their ubiquitous presence in nature, role in biological function and diseases have spurred interest in the biophysical and conformational characterization of IDPs. Conventional methods of structure determination are less feasible owing to structural and spatiotemporal heterogeneity of IDPs, which demand the development of novel biophysical methods as well as rigorous computational techniques for their characterization. In this review, we provide a brief overview of the widely used computational techniques probing the rugged conformational energy landscape of IDPs, their kinetics of structural transitions and molecular interactions key to their functions. Advances in the development of calibrated computational approaches for statistical representation of highly dynamic structural ensemble of IDPs are provided with examples. Challenges in modelling this unique class of proteins as well as the existing and futuristic avenues are also discussed.

Temperature-induced first-order displacive phase transition of isonicotinamide-4-methoxybenzoic acid co-crystal

Isonicotinamide-4-methoxybenzoic acid co-crystal (1), C6H6N2O·C8H8O3, is formed through slow evaporation from methanol solution and it undergoes a first-order isosymmetry (monoclinic I2/a ↔ monoclinic I2/a) structural phase transition at Tc = 142.5 (5) K, which has been confirmed by an abrupt jump of crystallographic interaxial angle β from variable-temperature single-crystal XRD and small heat hysteresis (6.25 K) in differential scanning calorimetry measurement. The three-dimensional X-ray crystal structures of (1) at the low-temperature phase (LTP) (100, 140 and 142 K) and the high-temperature phase (HTP) (143, 150, 200, 250 and 300 K) were solved and refined as a simple non-disordered model with final R[F2 > 2σ(F2)] ≃ 0.05. The asymmetric unit of (1) consists of crystallographically independent 4-methoxybenzoic acid (A) and isonicotinamide (B) molecules in both enantiotropic phases. Molecule A adopts a `near-hydroxyl' conformation in which the hydroxyl and methoxy groups are positioned on the same side. Both `near-hydroxyl' and `near-carbonyl' molecular conformations possess minimum conformational energies with an energy difference of < 0.15 kJ mol-1 from a potential energy surface scan. In the crystal, molecules are joined into linear ABBA arrays by intermolecular N-H...O and O-H...N hydrogen bonds which were preserved in both phases. However, these ABBA arrays are displaced from planarity upon LTP-to-HTP transition and the changes in inter-array interactions are observed in two-dimensional fingerprint plots of their Hirshfeld surfaces. The PIXEL energies of each molecular pair in both phases were calculated to investigate the difference in intermolecular interaction energies before and after the displacement of ABBA arrays from planarity, which directly leads to the single-crystal-to-single-crystal phase transition of (1).

Dynamic Descriptions of Highly Flexible Molecules from NMR Dipolar Couplings: Physical Basis and Limitations.

Biomolecules that control physiological function by changing their conformation play key roles in biology and remain poorly characterized. NMR dipolar couplings (DCs) depend intrinsically on both molecular shape and structural fluctuations, thereby providing the enticing prospect of tracking these conformational changes at atomic detail. Although this dual dependence has until now severely complicated analysis of DCs from highly dynamic systems, general approaches have recently been proposed that simplify interpretation of experimental DCs, by entirely eliminating molecular alignment from the analysis. Using simple and intuitive simulation of target ensembles, we investigate the impact of such approaches on the resulting descriptions of the conformational energy landscape. We find that ensemble descriptions of highly flexible systems derived from DCs without explicit consideration of the alignment properties of the constituent conformations can be compromised and inaccurate, despite exhibiting high correlation with experimental measurement.

The Michaelis Complex of Arginine Kinase Samples the Transition State at a Frequency That Matches the Catalytic Rate.

Arginine kinase (AK), which is a member of the phosphagen kinase family, serves as a model system for studying the structural and dynamic determinants of biomolecular enzyme catalysis of all major states involved of the enzymatic cycle. These states are the apo state (substrate free), the Michaelis complex analogue AK:Arg:Mg·AMPPNP (MCA), a product complex analogue AK:pAIE:Mg·ADP (PCA), and the transition state analogue AK:Arg:Mg·ADP:NO3- (TSA). The conformational dynamics of these states have been studied by NMR relaxation dispersion measurements of the methyl groups of the Ile, Leu, and Val residues at two static magnetic fields. Although all states undergo significant amounts of μs-ms time scale dynamics, only the MCA samples a dominant excited state that resembles the TSA, as evidenced by the strong correlation between the relaxation dispersion derived chemical shift differences Δω and the equilibrium chemical shift differences Δδ of these states. The average lifetime of the MCA is 36 ms and the free energy difference to the TSA-like form is 8.5 kJ/mol. It is shown that the conformational energy landscape of the Michaelis complex analogue is shaped in a way that at room temperature it channels passage to the transition state, thereby determining the rate-limiting step of the phosphorylation reaction of arginine. Conversely, relaxation dispersion experiments of the TSA reveal that it samples the structures of the Michaelis complex analogue or the apo state as its dominant excited state. This reciprocal behavior shows that the free energy of the TSA, with all ligands bound, is lower by only about 8.9 kJ/mol than that of the Michaelis or apo complex conformations with the TSA ligands present.

Solid-State Characterization of 2-[(2,6-Dichlorophenyl)amino]-Benzaldehyde: An Experimental and Theoretical Investigation

The solid-state properties of 2-[(2,6-dichlorophenyl)amino]-benzaldehyde (DCABA) were investigated. Unlike its precursor diclofenac acid, for which three polymorphs are currently known, only one crystalline form of DCABA was found. It was further characterized by other spectroscopic and spectrometric methods including IR, Raman, and UV–vis spectroscopy and powder X-ray diffraction (PXRD). The thermal behavior of the crystalline form was studied by differential scanning calorimetry (DSC). Theoretical studies, including Hirshfeld surface analysis and conformational energy searches, were performed to provide insight into the factors contributing to the stability of the crystal and to assess the possibility of additional polymorphs. The full characterization of this compound can help fast and accurate identification of DCABA both in dosage forms and in the environment.

Thermodynamics of globular proteins

The analysis of temperature-induced unfolding of proteins in aqueous solutions was performed. Based on the data of thermodynamic parameters of protein unfolding and using the method of semi-empirical calculations of hydration parameters at reference temperature 298 K, we obtained numerical values of enthalpy, free energy, and entropy which characterize the unfolding of proteins in the ‘gas phase’. It was shown that specific values of the energy of weak intramolecular bonds (∆Hint), conformational free energy (∆Gconf) and entropy (∆Sconf) are the same for proteins with molecular weight 7–25 kDa. Using the energy value (∆Hint) and the proposed approach for estimation of the conformational entropy of native protein (SNC), numerical values of the absolute free energy (GNC) were obtained.

Experimental conformational energy maps of proteins and peptides.

We have presented an extensive analysis of the peptide backbone dihedral angles in the PDB structures and computed experimental Ramachandran plots for their distributions seen under a various constraints on X-ray resolution, representativeness at different sequence identity percentages, and hydrogen bonding distances. These experimental distributions have been converted into isoenergy contour plots using the approach employed previously by F. M. Pohl. This has led to the identification of energetically favored minima in the Ramachandran (ϕ, ψ) plots in which global minima are predominantly observed either in the right-handed α-helical or the polyproline II regions. Further, we have identified low energy pathways for transitions between various minima in the (ϕ,ψ) plots. We have compared and presented the experimental plots with published theoretical plots obtained from both molecular mechanics and quantum mechanical approaches. In addition, we have developed and employed a root mean square deviation (RMSD) metric for isoenergy contours in various ranges, as a measure (in kcal.mol-1 ) to compare any two plots and determine the extent of correlation and similarity between their isoenergy contours. In general, we observe a greater degree of compatibility with experimental plots for energy maps obtained from molecular mechanics methods compared to most quantum mechanical methods. The experimental energy plots we have investigated could be helpful in refining protein structures obtained from X-ray, NMR, and electron microscopy and in refining force field parameters to enable simulations of peptide and protein structures that have higher degree of consistency with experiments. Proteins 2017; 85:979-1001. © 2017 Wiley Periodicals, Inc.

Electrostatic Interactions as Mediators in the Allosteric Activation of Protein Kinase A RIα.

Close-range electrostatic interactions that form salt bridges are key components of protein stability. Here we investigate the role of these charged interactions in modulating the allosteric activation of protein kinase A (PKA) via computational and experimental mutational studies of a conserved basic patch located in the regulatory subunit's B/C helix. Molecular dynamics simulations evidenced the presence of an extended network of fluctuating salt bridges spanning the helix and connecting the two cAMP binding domains in its extremities. Distinct changes in the flexibility and conformational free energy landscape induced by the separate mutations of Arg239 and Arg241 suggested alteration of cAMP-induced allosteric activation and were verified through in vitro fluorescence polarization assays. These observations suggest a mechanical aspect to the allosteric transition of PKA, with Arg239 and Arg241 acting in competition to promote the transition between the two protein functional states. The simulations also provide a molecular explanation for the essential role of Arg241 in allowing cooperative activation, by evidencing the existence of a stable interdomain salt bridge with Asp267. Our integrated approach points to the role of salt bridges not only in protein stability but also in promoting conformational transition and function.

Conformational energy range of ligands in protein crystal structures: The difficult quest for accurate understanding.

In this review, we address a fundamental question: What is the range of conformational energies seen in ligands in protein-ligand crystal structures? This value is important biophysically, for better understanding the protein-ligand binding process; and practically, for providing a parameter to be used in many computational drug design methods such as docking and pharmacophore searches. We synthesize a selection of previously reported conflicting results from computational studies of this issue and conclude that high ligand conformational energies really are present in some crystal structures. The main source of disagreement between different analyses appears to be due to divergent treatments of electrostatics and solvation. At the same time, however, for many ligands, a high conformational energy is in error, due to either crystal structure inaccuracies or incorrect determination of the reference state. Aside from simple chemistry mistakes, we argue that crystal structure error may mainly be because of the heuristic weighting of ligand stereochemical restraints relative to the fit of the structure to the electron density. This problem cannot be fixed with improvements to electron density fitting or with simple ligand geometry checks, though better metrics are needed for evaluating ligand and binding site chemistry in addition to geometry during structure refinement. The ultimate solution for accurately determining ligand conformational energies lies in ultrahigh-resolution crystal structures that can be refined without restraints.

Augmenting the anisotropic network model with torsional potentials improves PATH performance, enabling detailed comparison with experimental rate data.

PATH algorithms for identifying conformational transition states provide computational parameters—time to the transition state, conformational free energy differences, and transition state activation energies—for comparison to experimental data and can be carried out sufficiently rapidly to use in the “high throughput” mode. These advantages are especially useful for interpreting results from combinatorial mutagenesis experiments. This report updates the previously published algorithm with enhancements that improve correlations between PATH convergence parameters derived from virtual variant structures generated by RosettaBackrub and previously published kinetic data for a complete, four-way combinatorial mutagenesis of a conformational switch in Tryptophanyl-tRNA synthetase.