Structural Rearrangements of Rhodopsin Subunits in a Dimer Complex: a Molecular Dynamics Simulation Study

Investigating the Role of Rhodopsin F45L Mutation in Mouse Rod Photoreceptor Signaling and Survival

Rhodopsin is the critical receptor molecule which enables vertebrate rod photoreceptor cells to detect a single photon of light and initiate a cascade of molecular events leading to visual perception. Recently, it has been suggested that the F45L mutation in the transmembrane helix of rhodopsin disrupts its dimerization in vitro To determine whether this mutation of rhodopsin affects its signaling properties in vivo, we generated knock-in mice expressing the rhodopsin F45L mutant. We then examined the function of rods in the mutant mice versus wild-type controls, using in vivo electroretinography and transretinal and single cell suction recordings, combined with morphologic analysis and spectrophotometry. Although we did not evaluate the effect of the F45L mutation on the state of dimerization of the rhodopsin in vivo, our results revealed that F45L-mutant mice exhibit normal retinal morphology, normal rod responses as measured both in vivo and ex vivo, and normal rod dark adaptation. We conclude that the F45L mutation does not affect the signaling properties of rhodopsin in its natural setting.

Prediction and targeting of GPCR oligomer interfaces

GPCR oligomerization has emerged as a hot topic in the GPCR field in the last years. Receptors that are part of these oligomers can influence each other's function, although it is not yet entirely understood how these interactions work. The existence of such a highly complex network of interactions between GPCRs generates the possibility of alternative targets for new therapeutic approaches. However, challenges still exist in the characterization of these complexes, especially at the interface level. Different experimental approaches, such as FRET or BRET, are usually combined to study GPCR oligomer interactions. Computational methods have been applied as a useful tool for retrieving information from GPCR sequences and the few X-ray-resolved oligomeric structures that are accessible, as well as for predicting new and trustworthy GPCR oligomeric interfaces. Machine-learning (ML) approaches have recently helped with some hindrances of other methods. By joining and evaluating multiple structure-, sequence- and co-evolution-based features on the same algorithm, it is possible to dilute the issues of particular structures and residues that arise from the experimental methodology into all-encompassing algorithms capable of accurately predict GPCR-GPCR interfaces. All these methods used as a single or a combined approach provide useful information about GPCR oligomerization and its role in GPCR function and dynamics. Altogether, we present experimental, computational and machine-learning methods used to study oligomers interfaces, as well as strategies that have been used to target these dynamic complexes.

Oligomerization and cooperativity in GPCRs from the perspective of the angiotensin AT1 and dopamine D2 receptors

G Protein-Coupled Receptors (GPCRs) can form homo- and heterodimers or constitute higher oligomeric clusters with other heptahelical GPCRs. In this article, multiscale molecular modeling approaches as well as experimental techniques which are used to study oligomerization of GPCRs are reviewed. In particular, the effect of dimerization/oligomerization to the ligand binding affinity of individual protomers and also on the efficacy of the oligomer are discussed by including diverse examples from the literature. In addition, possible allosteric effects that may emerge upon interaction of GPCRs with membrane components, like cholesterol, is also discussed. Investigation of these above-mentioned interactions may greatly contribute to the candidate molecule screening studies and development of novel therapeutics with fewer adverse effects.

Characterization of Ofloxacin Interaction with Mutated (A91V) Quinolone Resistance Determining Region of DNA Gyrase in Mycobacterium Leprae through Computational Simulation

Mycobacterium leprae, the causal agent of leprosy is non-cultivable in vitro. Thus, the assessment of antibiotic activity against Mycobacterium leprae depends primarily upon the time-consuming mouse footpad system. The GyrA protein of Mycobacterium leprae is the target of the antimycobacterial drug, Ofloxacin. In recent times, the GyrA mutation (A91V) has been found to be resistant to Ofloxacin. This phenomenon has necessitated the development of new, long-acting antimycobacterial compounds. The underlying mechanism of drug resistance is not completely known. Currently, experimentally crystallized GyrA-DNA-OFLX models are not available for highlighting the binding and mechanism of Ofloxacin resistance. Hence, we employed computational approaches to characterize the Ofloxacin interaction with both the native and mutant forms of GyrA complexed with DNA. Binding energy measurements obtained from molecular docking studies highlights hydrogen bond-mediated efficient binding of Ofloxacin to Asp47 in the native GyrA-DNA complex in comparison with that of the mutant GyrA-DNA complex. Further, molecular dynamics studies highlighted the stable binding of Ofloxacin with native GyrA-DNA complex than with the mutant GyrA-DNA complex. This mechanism provided a plausible reason for the reported, reduced effect of Ofloxacin to control leprosy in individuals with the A91V mutation. Our report is the first of its kind wherein the basis for the Ofloxacin drug resistance mechanism has been explored with the help of ternary Mycobacterium leprae complex, GyrA-DNA-OFLX. These structural insights will provide useful information for designing new drugs to target the Ofloxacin-resistant DNA gyrase.

An Ensemble-Based Protocol for the Computational Prediction of Helix-Helix Interactions in G Protein-Coupled Receptors using Coarse-Grained Molecular Dynamics

The accurate identification of the specific points of interaction between G protein-coupled receptor (GPCR) oligomers is essential for the design of receptor ligands targeting oligomeric receptor targets. A coarse-grained molecular dynamics computer simulation approach would provide a compelling means of identifying these specific protein–protein interactions and could be applied both for known oligomers of interest and as a high-throughput screen to identify novel oligomeric targets. However, to be effective, this in silico modeling must provide accurate, precise, and reproducible information. This has been achieved recently in numerous biological systems using an ensemble-based all-atom molecular dynamics approach. In this study, we describe an equivalent methodology for ensemble-based coarse-grained simulations. We report the performance of this method when applied to four different GPCRs known to oligomerize using error analysis to determine the ensemble size and individual replica simulation time required. Our measurements of distance between residues shown to be involved in oligomerization of the fifth transmembrane domain from the adenosine A_2A receptor are in very good agreement with the existing biophysical data and provide information about the nature of the contact interface that cannot be determined experimentally. Calculations of distance between rhodopsin, CXCR4, and β₁AR transmembrane domains reported to form contact points in homodimers correlate well with the corresponding measurements obtained from experimental structural data, providing an ability to predict contact interfaces computationally. Interestingly, error analysis enables identification of noninteracting regions. Our results confirm that GPCR interactions can be reliably predicted using this novel methodology.

Methods used to study the oligomeric structure of G-protein-coupled receptors

G-protein-coupled receptors (GPCRs), which constitute the largest family of cell surface receptors, were originally thought to function as monomers, but are now recognized as being able to act in a wide range of oligomeric states and indeed, it is known that the oligomerization state of a GPCR can modulate its pharmacology and function. A number of experimental techniques have been devised to study GPCR oligomerization including those based upon traditional biochemistry such as blue-native PAGE (BN-PAGE), co-immunoprecipitation (Co-IP) and protein-fragment complementation assays (PCAs), those based upon resonance energy transfer, FRET, time-resolved FRET (TR-FRET), FRET spectrometry and bioluminescence resonance energy transfer (BRET). Those based upon microscopy such as FRAP, total internal reflection fluorescence microscopy (TIRFM), spatial intensity distribution analysis (SpIDA) and various single molecule imaging techniques. Finally with the solution of a growing number of crystal structures, X-ray crystallography must be acknowledged as an important source of discovery in this field. A different, but in many ways complementary approach to the use of more traditional experimental techniques, are those involving computational methods that possess obvious merit in the study of the dynamics of oligomer formation and function. Here, we summarize the latest developments that have been made in the methods used to study GPCR oligomerization and give an overview of their application.

Basic Pharmacological and Structural Evidence for Class A G-Protein-Coupled Receptor Heteromerization

Cell membrane receptors rarely work on isolation, often they form oligomeric complexes with other receptor molecules and they may directly interact with different proteins of the signal transduction machinery. For a variety of reasons, rhodopsin-like class A G-protein-coupled receptors (GPCRs) seem an exception to the general rule of receptor-receptor direct interaction. In fact, controversy surrounds their potential to form homo- hetero-dimers/oligomers with other class A GPCRs; in a sense, the field is going backward instead of forward. This review focuses on the convergent, complementary and telling evidence showing that homo- and heteromers of class A GPCRs exist in transfected cells and, more importantly, in natural sources. It is time to decide between questioning the occurrence of heteromers or, alternatively, facing the vast scientific and technical challenges that class A receptor-dimer/oligomer existence pose to Pharmacology and to Drug Discovery.

Using computer-aided drug design and medicinal chemistry strategies in the fight against diabetes

The aim of this work is to present a simple, practical and efficient protocol for drug design, in particular Diabetes, which includes selection of the illness, good choice of a target as well as a bioactive ligand and then usage of various computer aided drug design and medicinal chemistry tools to design novel potential drug candidates in different diseases. We have selected the validated target dipeptidyl peptidase IV (DPP-IV), whose inhibition contributes to reduce glucose levels in type 2 diabetes patients. The most active inhibitor with complex X-ray structure reported was initially extracted from the BindingDB database. By using molecular modification strategies widely used in medicinal chemistry, besides current state-of-the-art tools in drug design (including flexible docking, virtual screening, molecular interaction fields, molecular dynamics, ADME and toxicity predictions), we have proposed 4 novel potential DPP-IV inhibitors with drug properties for Diabetes control, which have been supported and validated by all the computational tools used herewith.

Thermal stability and unfolding pathways of Sso7d and its mutant F31A: insight from molecular dynamics simulation

The thermo-stability and unfolding behaviors of a small hyperthermophilic protein Sso7d as well as its single-point mutation F31A are studied by molecular dynamics simulation at temperatures of 300 K, 371 K and 500 K. Simulations at 300 K show that the F31A mutant displays a much larger flexibility than the wild type, which implies that the mutation obviously decreases the protein's stability. In the simulations at 371 K, although larger fluctuations were observed, both of these two maintain their stable conformations. High temperature simulations at 500 K suggest that the unfolding of these two proteins evolves along different pathways. For the wild-type protein, the C-terminal alpha-helix is melted at the early unfolding stage, whereas it is destroyed much later in the unfolding process of the F31A mutant. The results also show that the mutant unfolds much faster than its parent protein. The deeply buried aromatic cluster in the F31A mutant dissociates quickly relative to the wild-type protein at high temperature. Besides, it is found that the triple-stranded antiparallel β-sheet in the wild-type protein plays an important role in maintaining the stability of the entire structure.

Microsecond scale replica exchange molecular dynamic simulation of villin headpiece: an insight into the folding landscape

Reaching the experimental time scale of millisecond is a grand challenge for protein folding simulations. The development of advanced Molecular Dynamics techniques like Replica Exchange Molecular Dynamics (REMD) makes it possible to reach these experimental timescales. In this study, an attempt has been made to reach the multi microsecond simulation time scale by carrying out folding simulations on a three helix bundle protein, Villin, by combining REMD and Amber United Atom model. Twenty replicas having different temperatures ranging from 295 K to 390 K were simulated for 1.5 µs each. The lowest Root Mean Square Deviation (RMSD) structure of 2.5 Å was obtained with respect to native structure (PDB code 1VII), with all the helices formed. The folding population landscapes were built using segment-wise RMSD and Principal Components as reaction coordinates. These analyses suggest the two-stage folding for Villin. The combination of REMD and Amber United Atom model may be useful to understand the folding mechanism of various fast folding proteins.

Protein-protein interactions: principles, techniques, and their potential role in new drug development

A vast network of genes is inter-linked through protein-protein interactions and is critical component of almost every biological process under physiological conditions. Any disruption of the biologically essential network leads to pathological conditions resulting into related diseases. Therefore, proper understanding of biological functions warrants a comprehensive knowledge of protein-protein interactions and the molecular mechanisms that govern such processes. The importance of protein-protein interaction process is highlighted by the fact that a number of powerful techniques/methods have been developed to understand how such interactions take place under various physiological and pathological conditions. Many of the key protein-protein interactions are known to participate in disease-associated signaling pathways, and represent novel targets for therapeutic intervention. Thus, controlling protein-protein interactions offers a rich dividend for the discovery of new drug targets. Availability of various tools to study and the knowledge of human genome have put us in a unique position to understand highly complex biological network, and the mechanisms involved therein. In this review article, we have summarized protein-protein interaction networks, techniques/methods of their binding/kinetic parameters, and the role of these interactions in the development of potential tools for drug designing.

Potent inhibitor design against H1N1 swine influenza: structure-based and molecular dynamics analysis for M2 inhibitors from traditional Chinese medicine database

The rapid spread of influenza virus subtype H1N1 poses a great threat to million lives worldwide. To search for new anti-influenza compounds, we performed molecular docking and molecular dynamics simulation to identify potential traditional Chinese medicine (TCM) constituents that could block influenza M2 channel activity. Quinic acid, genipin, syringic acid, cucurbitine, fagarine, and methyl isoferulate all have extremely well docking results as compared to control amantadine. Further de novo drug design suggests that derivatives of genipin and methyl isoferulate could have enhanced binding affinity towards M2 channel. Selected molecular dynamics simulations of M2-derivative complexes show stable hydrogen bond interactions between the derivatives and M2 residues, Ser10 and Ala9. To our best knowledge, this is the first study on the anti-viral activity of the above listed TCM compounds.

Conformational analysis of genotoxic benzo[a]pyrene-7,8-dione-duplex DNA adducts using a molecular dynamics method (II)

The conformations of the benzo[a]pyrene-7,8-quinone (BPQ) modified oligonucleotide were investigated using molecular dynamic simulation. In the initial structures, the central guanine base was modified with BPQ resulting in the formation of four structurally distinguishable 10-(N2-deoxyguanosyl)-9,10-dihydro-9-hydroxy benzo[a]pyrene-7,8-dione adducts (BPQ-G3,4). Each of the oligonucleotide adduct consisted of two conformers, namely syn and anti conformations, depending on the rotation around the glycosidic bond between BPQ and the guanine base. The results revealed that the BPQ moiety was located in the major groove for all four syn conformers. The relative energies of these conformers were high, and the backbone largely deviated from the B-form. On the other hand, BPQ was located in the minor groove with relatively low energies, and backbone was retained in all of the anti conformer cases. The most conceivable BPQ-modified double stranded oligonucleotide structure was proposed from the energy calculation and the structural analysis.

RNA stability under different combinations of amber force fields and solvation models

The proper matching of force field and solvent is critical to obtain correct result in molecular dynamics simulation of bio-molecules. This problem has been intensively investigated for protein but not for RNA yet. In this paper, we use standard molecular dynamics and replica exchange molecular dynamics to take a series of tests on the RNA stability under different combinations of Amber force field parameters (ff98, ff99 and ff99bsc0) and the general Born implicit solvent models (igb1, igb2 and igb5). It is found that only ff98 and ff99bsc0 with igb1 can keep the native conformations of RNA hairpin and duplex. Our results suggest that ff98 plus igb1 may be reasonable choice for molecular dynamics simulation of RNA dynamics.

Influence of the acetylcholinesterase active site protonation on omega loop and active site dynamics

Existence of alternative entrances in acetylcholinesterase (AChE) could explain the contrast between the very high AChE catalytic efficiency and the narrow and long access path to the active site revealed by X-ray crystallography. Alternative entrances could facilitate diffusion of the reaction products or at least water and ions from the active site. Previous molecular dynamics simulations identified side door and back door as the most probable alternative entrances. The simulations of non-inhibited AChE suggested that the back door opening events occur only rarely (0.8% of the time in the 10ns trajectory). Here we present a molecular dynamics simulation of non-inhibited AChE, where the back door opening appears much more often (14% of the time in the 12ns trajectory) and where the side door opening was observed quite frequently (78% of trajectory time). We also present molecular dynamics, where the back door does not open at all, or where large conformational changes of the AChE omega loop occur together with alternative passage opening events. All these differences in AChE dynamical behavior are caused by different protonation states of two glutamate residues located on bottom of the active site gorge (Glu202 and G450 in Mus musculus AChE). Our results confirm the results of previous molecular dynamics simulations, expand the view and suggest the probable reasons for the overall conformational behavior of AChE omega loop.

A novel strategy for designing the selective PPAR agonist by the "sum of activity" model

Peroxisome proliferator-activated receptors alpha, delta and gamma are a collection of ligand-activated transcription factors crucial in lipid and glucose homeostasis. The involvement of these receptors in lipid metabolism makes them perfect therapeutic target for treating obesity and stroke. In this study, 'sum of activity' model was employed to design multi-target agonists. We used a new strategy to design agonists that fit both alpha and delta but not gamma, to avoid side effect. The CoMFA and CoMSIA models were used to explore the pharmacophore features by constructing three individual models: (a) alpha-model, (b) delta-model and (c) gamma-model, and two sum models: (d) alpha, delta- model, and (e) alpha, delta and gamma-model. The CoMFA model yielded a significant cross validation value, q(2), of 0.729 and non-cross validation value, r(2), of 0.933 in the alpha, delta-model. The CoMSIA studies yielded the best predictive models with q(2) of 0.622 in A+S and with r(2) of 0.911 in the alpha, delta-model. Finally, we proposed that distinct features shown in models (a), (b), (d) but not (c) and (e) should be accounted in designing weight-controlling drugs.

A theoretical elucidation of glucose interaction with HSA's domains

The interaction of different domains belonging to Human Serum Albumin (HSA) with open form of glucose have been investigated using molecular dynamics simulation methods. Applying docking, primary structures involving interaction of some residues with glucose have been obtained. Subsequently, equilibrium geometries at 300 K and minimum geometries have been determined for each of aforementioned structures by employing MD simulation and simulated annealing. The stability of species has been evaluated using a SAWSA v2.0 model. Ultimately, NBO analysis has been carried out to specify possible hydrogen bonding regarding the HSA interaction with glucose. Results obtained show that glucose can interact with Lys195, Lys199, and Glu153. In these interactions, each lysine forms an H-bonding with glucose. The H-bonding is obtained by stretching of N-H bond belonging to NH(3)(+) group of lysine along an oxygen atom of glucose. In addition, the above mentioned lysines are protonated, and there is an electrostatic interaction between glucose with Lys195 or Lys199. In addition, an H-bonding is formed between O atom of -COO group belonging to Glu153 and H atom of OH group belonging to glucose. Because, the N-H group of Lys195 interacts with the O atom of latter OH group, reaction of Lys195 is more desirable than that of Lys199. In fact, glucose is placed in the vicinity of Lys195 along with electrostatic interaction and H-bonding to Lys195 and Lys199 as well as H-bonding with Glu153, which subsequently reacts with Lys195. Thus, Lys195 is the primary site in reaction of glucose with HSA.

Insight derived from molecular dynamics simulation into substrate-induced changes in protein motions of proteinase K

Because of the significant industrial, agricultural and biotechnological importance of serine protease proteinase K, it has been extensively investigated using experimental approaches such as X-ray crystallography, site-directed mutagenesis and kinetic measurement. However, detailed aspects of enzymatic mechanism such as substrate binding, release and relevant regulation remain unstudied. Molecular dynamics (MD) simulations of the proteinase K alone and in complex with the peptide substrate AAPA were performed to investigate the effect of substrate binding on the dynamics/molecular motions of proteinase K. The results indicate that during simulations the substrate-complexed proteinase K adopt a more compact and stable conformation than the substrate-free form. Further essential dynamics (ED) analysis reveals that the major internal motions are confined within a subspace of very small dimension. Upon substrate binding, the overall flexibility of the protease is reduced; and the noticeable displacements are observed not only in substrate-binding regions but also in regions opposite the substrate-binding groove/pockets. The dynamic pockets caused by the large concerted motions are proposed to be linked to the substrate recognition, binding, orientation and product release; and the significant displacements in regions opposite the binding groove/pockets are considered to play a role in modulating the dynamics of enzyme-substrate interaction. Our simulation results complement the biochemical and structural studies, highlighting the dynamic mechanism of the functional properties of proteinase K.

Molecular dynamics studies on troponin (TnI-TnT-TnC) complexes: insight into the regulation of muscle contraction

Mutations of any subunit of the troponin complex may lead to serious disorders. Rational approaches to managing these disorders require knowledge of the complex interactions among the three subunits that are required for proper function. Molecular dynamics (MD) simulations were performed for both skeletal (sTn) and cardiac (cTn) troponin. The interactions and correlated motions among the three components of the troponin complex were analyzed using both Molecular Mechanics-Generalized Born Surface Area (MMGBSA) and cross-correlation techniques. The TnTH2 helix was strongly positively correlated with the two long helices of TnI. The C domain of TnC was positively correlated with TnI and TnT. The N domain of TnC was negatively correlated with TnI and TnT in cTn, but not in sTn. The two C-domain calcium-binding sites of TnC were dynamically correlated. The two regulatory N-domain calcium-binding sites of TnC were dynamically correlated, even though the calcium-binding site I is dysfunctional. The strong interaction residue pairs and the strong dynamically correlated residues pairs among the three components of troponin complexes were identified. These correlated motions are consistent with the idea that there is a high degree of cooperativity among the components of the regulatory complex in response to Ca(2+) and other effectors. This approach may give insight into the mechanism by which mutations of troponin cause disease. It is interesting that some observed disease causing mutations fall within regions of troponin that are strongly correlated or interacted.

The therapeutically anti-prion active antibody-fragment scFv-W226: paramagnetic relaxation-enhanced NMR spectroscopy aided structure elucidation of the paratope-epitope interface

Antibodies have become indispensable reagents with numerous applications in biological and biotechnical analysis, in diagnostics as well as in therapy. In all cases, selective interaction with an epitope is crucial and depends on the conformation of the paratope. While epitopes are routinely mapped at high throughput, methods revealing structural insights on a rather short timescale are rare. We here demonstrate paramagnetic relaxation-enhanced (PRE) NMR spectroscopy to be a powerful tool unraveling structural information about epitope-orientation in a groove spanned by the complementary determining regions. In particular, we utilize the spin label TOAC, which is fused to the peptidic epitope using standard solid-phase chemistry and which is characterized by a reduced mobility compared to, e.g., spin labels attached to the side-chain functionalities of cysteine or lysine residues. We apply the method to determine the orientation of helix 1 of the prion protein, which is the epitope for the therapeutically anti-prion active scF(v) fragment W226.

Conformational and oligomeric effects on the cysteine pK(a) of tryparedoxin peroxidase

Typical 2-Cys peroxiredoxins (Prxs) are peroxidases which regulate cell signaling pathways, apoptosis, and differentiation. These enzymes are obligate homodimers, and can form decamers in solution. During catalysis, Prxs exhibit cysteine-dependent reactivity which requires the deprotonation of the peroxidatic cysteine (C(p)) supported by a lowered pK(a) in the initial step. We present the results of molecular dynamics simulations combined with pKa calculations on the monomeric, dimeric and decameric forms of one typical 2-Cys Prx, the tryparedoxin peroxidase from Trypanosoma cruzi (PDB id, 1uul). The calculations indicate that C(p) (C52) pK(a) values are highly affected by oligomeric state; an unshifted C(p) pK(a) (approximately 8.3, comparable to the pK(a) of isolated cysteine) is calculated for the monomer. In the dimers, starting with essentially identical structures, the C(p)s evolve dynamically asymmetric pK(a)s during the simulations; one subunit's C(p) pK(a) is shifted downward at a time, while the other is unshifted. However, when averaged over time, or multiple simulations, the two subunits within a dimer exhibit the same C(p), showing no preference for a lowered pK(a) in either subunit. Two conserved pathways that communicate the asymmetric pK(a)s between C(p)s of different subunits can be identified. In the decamer, all the C(p) pK(a)s are shifted downward, with slight asymmetry in the dimers which form the decamers. Structural analyses implicate oligomerization effects as responsible for these oligomeric state-dependent C(p) pK(a) shifts. The intra-dimer and the inter-dimer subunit contacts in the decamer restrict the conformations of the side chains of several residues (T49, T54 and E55) calculated to be key in shifting the C(p) pK(a). In addition, the backbone fluctuations of a few residues (M46, D47 and F48) result in a different electrostatic environment for the C(p) in dimers relative to the monomers. These side chain and backbone interactions which contribute to pK(a) modulation indicate the importance of oligomerization to the function of the typical 2-Cys Prxs.

Molecular dynamics simulations of cyclohexyl modified peptide nucleic acids (PNA)

Peptide Nucleic Acids (PNA) that bind sequence specifically to DNA/RNA are of major interest in the field of molecular biology and could form the basis for gene-targeted drugs. Molecular dynamics simulations are aimed to characterize the structural and dynamical features to understand the effect of backbone modification on the structure and dynamics along with the stability of the resulting 10mer complexes of PNA with DNA/RNA. Twelve Molecular Dynamics (MD) simulations of duplexes and triplexes with and without cyclohexyl modification were carried out for 10ns each. The simulations indicate that the cyclohexyl modification with different stereoisomers has influenced all the PNA-DNA/RNA complexes. Modification has added rigidity to backbone by restricting beta to +60 in case of (1R,2S) cyclohexyl PNA and to -60 in case of (1S,2R) cyclohexyl PNA. The results of MD simulations were able to show the backbone rigidification and preference for RNA complexes over DNA due to presence of cyclohexyl ring in the PNA backbone.