Conformational Energy Research Articles

Structures and conformational dynamics of trans- and cis-crotonaldehydes in excited electronic states: A quantum-chemical study

For trans- and cis-crotonaldehydes CH3CHCHCHO in the triplet and singlet lowest excited electronic states of n,π* and π,π* types, the geometries of conformers, conformer energy differences, and conformational dynamics have been addressed using ab initio multi-reference techniques (CASSCF and CASPT2). The conformers of investigated molecules under study have been shown to possess a symmetry plane in 1,3(n,π*) states and to be significantly non-planar in 1,3(π,π*) states. Couplings of asymmetric internal rotations around CC and CC bonds with others large amplitude molecular motions such as methyl torsion, non-planar distortions of CCHO, CCHC, and CCH(CH3) fragments are discussed. The type of coupling between large-amplitude motions has been shown to be connected with the localization of electronic excitation.

Differences in the Conformational Energy Landscape of CDK1 and CDK2 Suggest a Mechanism for Achieving Selective CDK Inhibition

SummaryDysregulation of the cell cycle characterizes many cancer subtypes, providing a rationale for developing cyclin-dependent kinase (CDK) inhibitors. Potent CDK2 inhibitors might target certain cancers in which CCNE1 is amplified. However, current CDK2 inhibitors also inhibit CDK1, generating a toxicity liability. We have used biophysical measurements and X-ray crystallography to investigate the ATP-competitive inhibitor binding properties of cyclin-free and cyclin-bound CDK1 and CDK2. We show that these kinases can readily be distinguished by such inhibitors when cyclin-free, but not when cyclin-bound. The basis for this discrimination is unclear from either inspection or molecular dynamics simulation of ligand-bound CDKs, but is reflected in the contacts made between the kinase N- and C-lobes. We conclude that there is a subtle but profound difference between the conformational energy landscapes of cyclin-free CDK1 and CDK2. The unusual properties of CDK1 might be exploited to differentiate CDK1 from other CDKs in future cancer therapeutic design.

Using iSpartan To Support a Student-Centered Activity on Alkane Conformations

Strategies employed in the organic chemistry course must have the shared benefit of promoting concept mastery and visual literacy. With this in mind, an iSpartan-enabled visualization and computational activity was developed to complement learning in the first-semester, introductory organic chemistry course. The activity provides students an opportunity to draw molecules, collect structure-based data, and engage in group discussions related to acyclic alkane conformations. Students do not receive a lecture or complete required reading on conformations. Instead, the activity along with group- and instructor-facilitated discussions serves as the method of content delivery for the topic. Pretest and post-test scores suggest improvements in each student’s ability to interpret Newman Projections and determine relative conformation energies and stabilities. A survey was administered to determine the degree to which each student valued the activity and the associated technology, the response to which suggests th...

Conjugated Polymers Based on Thiazole Flanked Naphthalene Diimide for Unipolar n-Type Organic Field-Effect Transistors

This paper reports the rational design and synthesis of a novel electron-withdrawing building block, thiazole flanked naphthalene diimide (TzNDI), which offers a coplanar conformation and deep-lying highest-occupied molecular orbitals energy level in resulting conjugated polymers. A series of conjugated polymers (PTzNDI-2FT, PTzNDI-T, PTzNDI-Se, and PTzNDI-2T) consisting of TzNDI and different donor units were synthesized and characterized. The polymers all possess a high molecular weight and excellent thermal property. Their intense light absorption in low energy bands suggests an enhanced intramolecular charge transfer. The organic field-effect transistors (OFETs) based on these polymers exhibit unipolar n-type transport characteristics with low off current and high on–off current ratio. More importantly, all the devices exhibit near ideal transfer curves with kink-free transfer characteristics. Among these polymers, PTzNDI-2FT exhibits the highest electron mobility (μe) of 0.57 cm2 V–1 s–1, outperformi...

Archaic and alternative chaperones preserve pilin folding energy by providing incomplete structural information

Adhesive pili are external component of fibrous adhesive organelles and help bacteria attach to biotic or abiotic surfaces. The biogenesis of adhesive pili via the chaperone-usher pathway (CUP) is independent of external energy sources. In the classical CUP, chaperones transport assembly-competent pilins in a folded but expanded conformation. During donor-strand exchange, pilins subsequently collapse, producing a tightly packed hydrophobic core and releasing the necessary free energy to drive fiber formation. Here, we show that pilus biogenesis in non-classical, archaic, and alternative CUPs uses a different source of conformational energy. High-resolution structures of the archaic Csu-pili system from Acinetobacter baumannii revealed that non-classical chaperones employ a short donor strand motif that is insufficient to fully complement the pilin fold. This results in chaperone-bound pilins being trapped in a substantially unfolded intermediate. The exchange of this short motif with the longer donor strand from adjacent pilin provides the full steric information essential for folding, and thereby induces a large unfolded-to-folded conformational transition to drive assembly. Our findings may inform the development of anti-adhesion drugs (pilicides) to combat bacterial infections.

Raman, infrared and NMR spectra, vibrational assignments and quantum mechanical calculations of centrosymmetric 3,6-Dicholoro-1,2,4,5-tetrazine

The Raman spectrum (2000−100 cm−1) of 3,6-Dicholoro-1,2,4,5-tetrazine (DCTZ, C2N4Cl2) in solid phase has been recorded. Moreover, the mid-IR (2000–400 cm−1) and ATR-IR (2000-480 cm−1), 13C (0–200 ppm) and 15N (0–440 ppm) NMR spectra were recorded. Initially, optimization followed by frequency calculations were carried out for DCTZ restricted to Ci and D2h point groups to work out the least energy conformer with real frequencies. DFT(B3LYP) and MP2(full) computations utilizing 6-31G(d) and 6–311++G(d,p) basis sets favor Ci symmetry conformer by 781/598 cm−1 (9.3/7.2 kjmol−1) and 745/698 cm−1 (8.9/8.3 kjmol−1). The simulated spectra are well compiled to the experimentally recorded IR, Raman, 13C and 15N regardless of few coincident vibrational frequencies that belong to D2h symmetry conformer. Thus, we have carried out additional CASTEP solid state optimization for one molecule per unit cell (Ci symmetry) employing Materials Studio using DFT-PBE method. Similar calculations were completed for D2h molecular symmetry for one molecule and repeated molecules per unit cell in which Ci symmetry is favored by 341 cm−1 (4.1 kjmol−1). It is worth mentioning that computations for both Ci and D2h conformers produce real vibrational frequencies. Assisted by normal coordinate analysis, potential energy distributions (PEDs), complete and confident vibrational assignments are provided herein for Ci conformer. The results are conferred herein are compared with comparable molecules whenever applicable.

Efficiently Capturing Weak Interactions in ab Initio Molecular Dynamics with on-the-Fly Basis Set Extrapolation.

Weak interactions have a critical role in accurately portraying conformational change. However, the computational study of these often requires large basis electronic structure calculations that are generally cost-prohibitive within ab initio molecular dynamics. Here, we present a new approach to efficiently obtain AIMD trajectories in agreement with large, triple-ζ, polarized valence basis functions, at much reduced computational cost. For example, it follows from our studies that AIMD trajectories can indeed be constructed in agreement with basis sets such as 6-311++G(2df,2pd) with computational effort commensurate with those from much smaller basis sets such as 6-31+G(d), for polypeptide systems with 100+ atoms. The method is based on molecular fragmentation and allows a range-specified repartitioning of intramolecular (and potentially intermolecular) interactions where noncovalent interactions are selectively assembled using a piece-wise reconstruction based on a set-theoretic inclusion-exclusion principle generalization of ONIOM. Through a simplex decomposition of molecular systems the approach efficiently provides the necessary many-body interactions to faithfully represent noncovalent interactions at the large basis limit. Conformational stabilization energies are provided at close to the complete-basis limit at much reduced cost, and similarly AIMD trajectories (both Born-Oppenheimer and Car-Parrinello-type) are obtained in agreement with very large basis set sizes, in an extremely efficient and accurate manner. The method is demonstrated through simulations on polypeptide fragments of a variety of sizes.

Low Energy Conformations for S100 Binding Peptide from the Negative Regulatory Domain of p53.

We have computed the low energy conformations of the negative regulatory domain of p53, residues 374-388 using Empirical Conformational Energies of Peptides Program including solvation and computed the statistical weights of distinct conformational states. We find that there are two high probability conformations, one an α-helix from Lys 374-Lys 381, followed by another helical structure involving Lys 382-Glu 388 (statistical weight of 0.48) and an all-α-helix for the entire sequence (statistical weight of 0.23). Both structure are superimposable on the NMR structure of this sequence bound to the S100 protein. The global minimum structure (statistical weight of 0.014) is a beta structure from Gly 374-Arg 379 followed by an α-helix from His 380-Glu 388. Based on these results, we propose a possible strategy for enhancement of p53 anti-tumor activity in cancer cells. Since the structure of this sequence bound to the sirtuin protein, Sir2, is a β-sheet, we further propose that the global minimum may be an intermediate on the α-β structure transition.

Investigation of structural stability and functionality of homodimeric gramicidin towards peptide-based drug: a molecular simulation approach.

Increasing death rates due to antibiotic resistance deteriorate the existing treatment measures. Antimicrobial peptides have turned into the emerging cure for multidrug resistance. However, the stability and functionality determine an antimicrobial peptide as a drug. Analyses of the homodimeric β-helical peptide, gramicidin have suggested the significant role of gramicidin-A, gramicidin-B, and gramicidin-C as antimicrobial compounds, but the structural basis for understanding the stability and functionality is insufficient to resolve multidrug resistance. To identify the best template among gramicidin types as a therapeutic product, we combined a detailed comparative static analysis and dynamic analysis along with conformational free energy and secondary structure prediction. We observed that the high intramolecular interactions and the geometrical features favored gramicidin-A among other types of gramicidin. Our analyses further revealed that the secondary structure of gramicidin-A showed β sheets with coils along the conformations without any disruption, thereby enhanced its membrane interactions in terms of binding free energy. In conclusion, gramicidin-A has definitely showed enhanced structural stability and functionality; this could be considered the best template for a potential therapeutic product.

A prelude to building mathematical models for polypeptide folding: analysis on the conformational potential energy hypersurface cross-sections of N-acetyl-glycyl-glycine-N′-methylamide

Finding a relationship on how a three-dimensional protein folds from its linear amino acid chain gets more complex with increasing chain length, so working on a smaller peptide conformational problem can provide initial ideas on what are the main molecular forces and how these influence the folding process. Following the study of conformations of amino acid units entering the proteins to understand the secondary structure of small peptides, this paper proposes mathematical models for the several two-rotor cross-sections of the five-dimensional N-acetyl-glycyl-glycine-N′-methylamide potential energy hypersurface (PEHS). These cross-sections are extracted along the first glycine subunit, with its coordinates fixed at the five energy minima of the glycine diamide. The resulting mathematical models yield an average RMSE of 1.36 kJ mol−1 and an average R2 of 0.9923 with respect to energy values obtained from DFT calculations. The minima geometries obtained from these models are also in good agreement with DFT-optimized energy minima conformers. An important aspect of this study also tackles the relationship between the PEHS of the glycyl-glycine diamide and its glycine subunits. It has been observed that there are deviations up to 28.35 kJ mol−1 and 29.52 kJ mol−1 between the PEHS cross-sections along γL and γD conformations, respectively, in the first glycine subunit. This may suggest that there are significant backbone–backbone intermolecular forces acting on the dipeptide. The abovementioned findings can help in developing more complex mathematical models for polypeptide folding from amino acid subunits.

Conformational aspects of polymorphs and phases of 2-propyl-1H-benzimidazole.

This paper reports on the polymorphism of 2-propyl-1H-benzimidazole (2PrBzIm) induced by temperature change. Upon heating, an irreversible reconstructive-type phase transition at T = 384 K from the ordered form I (P212121) to a new polymorph, form II HT (Pcam), was observed. The structural transformation between forms I and II involves significant changes in the crystal packing, as well as a key conformational variation around the propyl chain of the molecule. After the first irreversible phase transition, the II HT form undergoes two further (reversible) phase transitions upon cooling at 361 K (II RT) and 181 K (II LT). All three phases (forms II HT, II RT and II LT) have almost identical crystal packing and, given the reversibility of the conversions as a function of temperature, they are referred to as form II temperature phases. They differ, however, with respect to conformational variations around the propyl chain of 2PrBzIm. Energy calculations of the gas-phase conformational energy landscape of this compound about its flexible bonds allowed us to classify the observed conformational variations of all forms into changes and adjustments of conformers. This reveals that forms I and II are related by conformational change, and that two of the form II phases (HT and RT) are related by conformational adjustment, whilst the other two (RT and LT) are related by conformational change. We introduce the term 'conformational phases' for different crystal phases with almost identical packing but showing changes in conformation.

Structural basis of urea‐induced unfolding of Fasciola gigantica glutathione S‐transferase

Glutathione S-transferases (GSTs) are enzymes that are involved in the detoxification of harmful electrophilic endogenous and exogenous compounds by conjugating with glutathione (GSH). The liver fluke GSTs have multifunctional roles in the host-parasite interaction, such as general detoxification and bile acid sequestration to synthase activity. The GSTs have been highlighted as vaccine candidates towards parasitic flukes. In this study, we have thoroughly examined the urea-induced unfolding of a mu-class Fasciola gigantica GST1 (FgGST1) using spectroscopic techniques and molecular dynamic simulations. FgGST1 is a highly cooperative molecule, because during urea-induced equilibrium unfolding, a concurrent unfolding of the protein without stabilization of any folded intermediate was observed. The protein was stabilized with conformational free energy of about ~12.36 kcal/mol. The protein loses its activity with increasing urea concentration, as the GSH molecule is not able to bind to the protein. We also studied the fluorescence quenching of Trp residues and the obtained K SV data that provided additional information on the unfolding of FgGST1. Molecular dynamic trajectories simulated in different urea concentrations and temperatures indicated that urea destabilizes FgGST1 structure by weakening hydrophobic interactions and the hydrogen bond network. We observed a precise correlation between the in vitro and in silico studies.

Electronic Energies Are Not Enough: An Ion Mobility-Aided, Quantum Chemical Benchmark Analysis of H+GPGG Conformers.

Peptides and proteins exist not as single, lowest energy structures but as ensembles of states separated by small barriers. In order to study these species we must be able to correctly identify their gas phase conformational distributions, and ion mobility spectrometry (IMS) has arisen as an experimental method for assessing the gas phase energetics of flexible peptides. Here, we present a thorough exploration and benchmarking of the low energy conformers of the small, hairpin peptide H+GPGG with the aid of ion mobility spectrometry against a wide swath of density functionals (35 dispersion-corrected and uncorrected functionals represented by 21 unique exchange-correlation functionals) and wave function theory methods (15 total levels of theory). The three experimentally resolved IMS peaks were found to correspond to three distinct pairs of conformers, each pair composed of species differing only by the chair or boat configuration of the proline. Two of the H+GPGG conformer pairs possess a cis configuration about the Pro-Gly1 peptide bond while the other adopts a trans configuration. While the experimental spectrum reports a higher intensity for the cis-1 conformer than the trans conformer, 13 WFT levels of theory, including a complete basis set CCSD(T) extrapolation, obtain trans to be favored in terms of the electronic energy. This same effect is seen in 14 of the 18 dispersion-corrected density functionals studied, whereas the remaining 17 functionals show more variety. Only when Gibbs free energies are considered do the WFT methods and dispersion-corrected functionals reflect the experimental distribution. CAM-B3LYP-D3BJ emerges as the best-performing density functional, matching the experimental distribution and the CCSD(T)/CBS relative energies within 6%. Further analysis reveals the trans conformer to be favored electronically, but entropically disfavored, leading to the experimental preference of the cis-1 conformer. These results highlight the danger in considering only the electronic energies, which is common practice in electronic structure theory predictions of conformational energy distributions. Additionally, the effects of temperature and scaling of the frequencies used in obtaining the Gibbs free energies are explored.

Identification of Effective Dimeric Gramicidin-D Peptide as Antimicrobial Therapeutics over Drug Resistance: In-Silico Approach.

Discovering and developing the antimicrobial peptides are recently focused on pharmaceutical firm, since they serve as complementary to antibiotics in prevailing over drug resistance by eliciting the disruption of microbial membrane. Still, there are lots of challenges to bring up the structurally stable and functionally efficient antimicrobial peptides. It is well known that gramicidin D is the prominent antimicrobial peptide that exists as g-AB, g-BC, and g-AC. This study analyzes the structural stability and the functional activity of hetero-dimeric double-stranded gramicidin-D peptides, thereby demonstrating its potent antimicrobial activity against antibiotic-resistant micro-organisms. To investigate the structural stability and functionality of gramicidin D, we performed static and dynamic analysis. Initially, we observed a maximum number of intermolecular interactions and membrane penetration in g-AB as compared to g-BC and g-AC. To substantiate further, the geometrical and thermodynamic parameters revealed the retention of maximum stability in g-AB than g-AC and g-BC. Thus, the conformational free energy and the binding free energy showed the variation among gramicidin-D peptides for the prediction of increased stability and functionality. In conclusion, g-AB peptide has definitely demonstrated adequate structural stability and functionality and this work will need to be considered in peptide-based drug discovery.

The Electron-Conformational Model of Ryanodine Receptors of the Heart Cell

The main ideas of an electron-conformational (EC) model of ligand-activated ryanodine channels (RyR) that reduces many degrees of freedom of this gigantic nanoscopic molecular complex to two so-called electron and conformational degrees of freedom are presented. In the toy model, the conformational potential or energy profile of the RyR channel represents two branches that describe the dependence of energy on the conformational coordinate (reaction coordinate) in the initial and ligand-activated states, respectively. A model extension that takes into account both the tetrameric structure of the RyR channel and additional (orthogonal) rotational conformational mode has been considered. The EC model has been demonstrated to give a biophysical basis to the traditional phenomenological model of Markov chains. To demonstrate the possibilities of the EC model, we refer to examples of model description of the dynamics of isolated RyR channels and clusters of RyR channels in release units of the heart cell.

Hybrid parallelization of a multi-tree path search algorithm: Application to highly-flexible biomolecules

The study of the conformational energy landscape of a molecule is essential for the understanding of its physicochemical properties. This requires the exploration of a continuous, high-dimensional space to identify the most probable conformations and the transition paths between them. The problem is computationally difficult, in particular for highly-flexible biomolecules such as Intrinsically Disordered Proteins (IDPs). In recent years, a robotics-inspired algorithm called Transition-based Rapidly-exploring Random Tree (TRRT) has been proposed to solve this problem, and has been shown to provide good results with small and middle-sized biomolecules. Aiming to treat larger systems, we propose a hybrid strategy for the efficient parallelization of a multi-tree variant of TRRT, called Multi-TRRT, enabling an efficient execution in (possibly large) computer clusters. The parallel algorithm uses OpenMP multi-threading for computation inside each multi-core processor and MPI to perform the communication between processors. Results show a near-linear speedup for a wide range of cluster configurations. Although the paper mainly deals with the application of the proposed parallel algorithm to the investigation of biomolecules, the explanations concerning the methods are general, aiming to inspire future work on the parallelization of related algorithms.

Revisiting the Case of an Intramolecular Hydrogen Bond Network Forming Four- and Five-Membered Rings in d-Glucose.

The conformational behavior of cyclic monosaccharides has been widely studied over the past years, but there is no general agreement about which effects are in fact responsible for the observed conformational preferences. A recent microwave spectroscopy study determined the conformational equilibrium of d-glucose in the gas phase with a preference for a counterclockwise arrangement of the hydroxyl groups. Nevertheless, the effects that control this orientation are still uncertain because the role of intramolecular hydrogen bonds (IHBs), electrostatic and steric repulsions is not clear. This work reports a density functional theory approach based on the conformational energies of d-glucose and of some derivatives in which the anomeric hydroxyl is replaced with hydrogen (H, small and not prone to participate in proton transfer), fluorine (F, small, electronegative, and as capable as OH of forming hydrogen bonds as a proton acceptor), and chlorine (Cl, big and not anticipated to be involved in effective hydrogen bond formation) to obtain insights into the effects of the substituent at the anomeric carbon on the arrangement of the hydroxyl groups in d-glucose. The nature of the substituents at this position is crucial to determine the orientation of the remaining hydroxyl groups. Natural bond orbital (NBO) and quantum theory of atoms in molecules (QTAIM) analyses, in addition to NMR chemical shift calculations, have been provided to support the conformational energy outcomes. Overall, the results agree with the lack of IHBs forming four- and five-membered rings in d-glucose and emphasize that steric and electrostatic repulsions involving the hydroxyl groups in the clockwise orientation are driving forces of the conformational behavior.

Entropy Production Associated with Aggregation into Granules in a Subdiffusive Environment

We study the entropy production that is associated with the growing or shrinking of a small granule in, for instance, a colloidal suspension or in an aggregating polymer chain. A granule will fluctuate in size when the energy of binding is comparable to , which is the “quantum” of Brownian energy. Especially for polymers, the conformational energy landscape is often rough and has been commonly modeled as being self-similar in its structure. The subdiffusion that emerges in such a high-dimensional, fractal environment leads to a Fokker–Planck Equation with a fractional time derivative. We set up such a so-called fractional Fokker–Planck Equation for the aggregation into granules. From that Fokker–Planck Equation, we derive an expression for the entropy production of a growing granule.

Deuterium Rotating Frame NMR Relaxation Measurements in the Solid State under Static Conditions for Quantification of Dynamics.

The feasibility of static deuterium rotating frame NMR relaxation measurements for characterization of slow timescale motions in powder systems is demonstrated. Using a model compound dimethyl sulfone-d6 , we show that these measurements yield conformational exchange rates and activation energy values in accordance with results obtained with other techniques. Furthermore, we demonstrate that the full Liouvillian approach as opposed to the Redfield approximation is necessary to analyze the experimental data.

Unraveling the mechanism of the cadherin-catenin-actin catch bond.

The adherens junctions between epithelial cells involve a protein complex formed by E-cadherin, β-catenin, α-catenin and F-actin. The stability of this complex was a puzzle for many years, since in vitro studies could reconstitute various stable subsets of the individual proteins, but never the entirety. The missing ingredient turned out to be mechanical tension: a recent experiment that applied physiological forces to the complex with an optical tweezer dramatically increased its lifetime, a phenomenon known as catch bonding. However, in the absence of a crystal structure for the full complex, the microscopic details of the catch bond mechanism remain mysterious. Building on structural clues that point to α-catenin as the force transducer, we present a quantitative theoretical model for how the catch bond arises, fully accounting for the experimental lifetime distributions. The underlying hypothesis is that force induces a rotational transition between two conformations of α-catenin, overcoming a significant energy barrier due to a network of salt bridges. This transition allosterically regulates the energies at the interface between α-catenin and F-actin. The model allows us to predict these energetic changes, as well as highlighting the importance of the salt bridge rotational barrier. By stabilizing one of the α-catenin states, this barrier could play a role in how the complex responds to additional in vivo binding partners like vinculin. Since significant conformational energy barriers are a common feature of other adhesion systems that exhibit catch bonds, our model can be adapted into a general theoretical framework for integrating structure and function in a variety of force-regulated protein complexes.