Systematic Investigation of Tether Length and Phosphorus Configuration in Backbone Constrained Macrocyclic Nucleic Acids to Modulate Binding Kinetics for RNA.

We recently described a chemical strategy to pre-organize a trinucleotide subunit in a conformation suitable for Watson-Crick base pairing for modulating the binding kinetics of single-stranded oligonucleotides (ONs) using bis-phosphonate esters bridging hydrocarbon tethers to provide 11- and 15-membered macrocyclic analogues. In this manuscript, we describe the synthesis of all eight P-stereoisomers of macrocyclic 12-, 13-, 14-, and 16-membered hydrocarbon-bridged nucleotide trimers, their incorporation into ONs, and biophysical characterization of the modified ONs. The size of the macrocyclic tether and configuration at phosphorus had profound effects on hybridization kinetics. ONs containing 12- and 13-membered rings exhibited faster on-rates (up to 5-fold) and off-rates (up to 161-fold). In contrast, ONs using the larger ring size macrocycles generally exhibited smaller changes in binding kinetics relative to unmodified DNA. Interestingly, several of the analogues retained significant binding affinity for RNA based on their dissociation constants, despite being modestly destabilizing in the thermal denaturation experiments, highlighting the potential utility of measuring dissociation constants versus duplex thermal stability when evaluating novel nucleic acid analogues. Overall, our results provide additional insights into the ability of backbone-constrained macrocyclic nucleic acid analogues to modulate hybridization kinetics of modified ONs with RNA.

Ethidium bromide interactions with DNA: an exploration of a classic DNA-ligand complex with unbiased molecular dynamics simulations.

Visualization of double stranded DNA in gels with the binding of the fluorescent dye ethidium bromide has been a basic experimental technique in any molecular biology laboratory for >40 years. The interaction between ethidium and double stranded DNA has been observed to be an intercalation between base pairs with strong experimental evidence. This presents a unique opportunity for computational chemistry and biomolecular simulation techniques to benchmark and assess their models in order to see if the theory can reproduce experiments and ultimately provide new insights. We present molecular dynamics simulations of the interaction of ethidium with two different double stranded DNA models. The first model system is the classic sequence d(CGCGAATTCGCG)2 also known as the Drew-Dickerson dodecamer. We found that the ethidium ligand binds mainly stacked on, or intercalated between, the terminal base pairs of the DNA with little to no interaction with the inner base pairs. As the intercalation at the terminal CpG steps is relatively rapid, the resultant DNA unwinding, rigidification, and increased stability of the internal base pair steps inhibits further intercalation. In order to reduce these interactions and to provide a larger groove space, a second 18-mer DNA duplex system with the sequence d(GCATGAACGAACGAACGC) was tested. We computed molecular dynamics simulations for 20 independent replicas with this sequence, each with ∼27 µs of sampling time. Results show several spontaneous intercalation and base-pair eversion events that are consistent with experimental observations. The present work suggests that extended MD simulations with modern DNA force fields and optimized simulation codes are allowing the ability to reproduce unbiased intercalation events that we were not able to previously reach due to limits in computing power and the lack of extensively tested force fields and analysis tools.

The Amoebicidal Activity of Diferrocenyl Derivatives: A Significant Dependence on the Electronic Environment.

Amoebiasis is the second leading cause of death worldwide associated with parasitic disease and is becoming a critical health problem in low-income countries, urging new treatment alternatives. One of the most promising strategies is enhancing the redox imbalance within these susceptible parasites related to their limited antioxidant defense system. Metal-based drugs represent a perfect option due to their extraordinary capacity to stabilize different oxidation states and adopt diverse geometries, allowing their interaction with several molecular targets. This work describes the amoebicidal activity of five 2-(Z-2,3-diferrocenylvinyl)-4X-4,5-dihydrooxazole derivatives (X = H (3a), Me (3b), iPr (3c), Ph (3d), and benzyl (3e)) on Entamoeba histolytica trophozoites and the physicochemical, experimental, and theoretical properties that can be used to describe the antiproliferative activity. The growth inhibition capacity of these organometallic compounds is strongly related to a fine balance between the compounds' redox potential and hydrophilic character. The antiproliferative activity of diferrocenyl derivatives studied herein could be described either with the redox potential, the energy of electronic transitions, logP, or the calculated HOMO-LUMO values. Compound 3d presents the highest antiproliferative activity of the series with an IC50 of 23 µM. However, the results of this work provide a pipeline to improve the amoebicidal activity of these compounds through the directed modification of their electronic environment.

Riboflavin Stabilizes Abasic, Oxidized G-Quadruplex Structures.

The G-quadruplex is a noncanonical fold of DNA commonly found at telomeres and within gene promoter regions of the genome. These guanine-rich sequences are highly susceptible to damages such as base oxidation and depurination, leading to abasic sites. In the present work, we address whether a vacancy, such as an abasic site, in a G-quadruplex serves as a specific ligand recognition site. When the G-tetrad is all guanines, the vacant (abasic) site is recognized and bound by free guanine nucleobase. However, we aim to understand whether the preference for a specific ligand recognition changes with the presence of a guanine oxidation product 8-oxo-7,8-dihydroguanine (OG) adjacent to the vacancy in the tetrad. Using molecular dynamics simulation, circular dichroism, and nuclear magnetic resonance, we examined the ability for riboflavin to stabilize abasic site-containing G-quadruplex structures. Through structural and free energy binding analysis, we observe riboflavin's ability to stabilize an abasic site-containing G-quadruplex only in the presence of an adjacent OG-modified base. Further, when compared to simulation with the vacancy filled by free guanine, we observe that the free guanine nucleobase is pushed outside of the tetrad by OG to interact with other parts of the structure, including loop residues. These results support the preference of riboflavin over free guanine to fill an OG-adjacent G-quadruplex abasic vacancy.

Backbone Hydrocarbon-Constrained Nucleic Acids Modulate Hybridization Kinetics for RNA.

The binding affinity of therapeutic oligonucleotides (ONs) for their cognate RNA is determined by the rates of association (ka) and dissociation (kd). Single-stranded ONs are highly flexible and can adopt multiple conformations in solution, some of which may not be conducive for hybridization. We investigated if restricting rotation around the sugar-phosphate backbone, by tethering two adjacent backbone phosphonate esters using hydrocarbon bridges, can modulate hybridization kinetics of the modified ONs for complementary RNA. Given the large number of possible analogues with different tether lengths and configurations at the phosphorus atoms, we employed molecular dynamic simulations to optimize the size of the hydrocarbon bridge to guide the synthetic efforts. The backbone-constrained nucleotide trimers with stereodefined configurations at the contiguous backbone phosphorus atoms were assembled using a ring-closing metathesis reaction, then incorporated into oligonucleotides by an in situ synthesis of the phosphoramidites followed by coupling to solid supports. Evaluation of the modified oligonucleotides revealed that 15-membered macrocyclic-constrained analogues displayed similar or slightly improved on-rates but significantly increased off-rates compared to unmodified DNA ONs, resulting in reduced duplex stability. In contrast, LNA ONs with conformationally preorganized furanose rings showed similar on-rates to DNA ONs but very slow off-rates, resulting in net improvement in duplex stability. Furthermore, the experimental data generally supported the molecular dynamics simulation results, suggesting that this strategy can be used as a predictive tool for designing the next generation of constrained backbone ON analogues with improved hybridization properties.

Effect of 4-HNE Modification on ZU5-ANK Domain and the Formation of Their Complex with ß-Spectrin: A Molecular Dynamics Simulation Study.

4-HNE-modified ankyrins have been described in diseases such as diabetes, renal failure, G6PD deficient, sickle cell trait, and P. falciparum infected erythrocytes with different AB0 blood groups. However, effects at the atomic level of this carbonylation on structure and function of modified protein are not yet fully understood. We present a study based on molecular dynamics simulations of nine 4-HNE modified residues of the ZU5-ANK ankyrin domain with ß-spectrin and their binding energy profiles. Results show that 4-HNE induced local conformational changes over all protein systems evaluated, increased mobility in the modification sites, and localized structural changes between the positively charged patch of the ZU5-ANK domain. Carbonylation with 4-HNE on lysine residues decreased the affinity between ZU5-ANK and the 14-ß-spectrin repeat by reducing electrostatic and van der Waals interactions. The presented work provides further insight into understanding the loss of human erythrocyte deformation capacity under conditions of oxidative stress in different diseases.

Probing the influence of hypermodified residues within the tRNA3(Lys) anticodon stem loop interacting with the A-loop primer sequence from HIV-1.

Replication of the HIV-1 virus requires reverse transcription of the viral RNA genome, a process that is specifically initiated by human tRNA3(Lys) packaged within the infectious virion. The primary binding site for the tRNA involves the 3' 18 nucleotides with an additional interaction between an adenine rich loop (A-loop) in the template and the anticodon stem-loop region of the tRNA3(Lys). The loop of the tRNA primer contains two hypermodified base residues and a pseudouridine that are required for a proper binding and activity. Here, we investigate the influence on the structure, dynamics and binding stability of the three modified residues (mnm(5)s(2)U34, t(6)A37 and Ψ39) using extensive molecular dynamics and Quantum Theory of Atoms in Molecules (QTAIM) analysis. Consistent with experiment, the results suggest that the three modified residues are required for faithful binding. Residues mnm(5)s(2)U34 and Ψ39 have a major influence in stabilizing the anticodon loop whereas mnm(5)s(2)U34 and t(6)A37 appear to stabilize the formation of the complex of tRNA3(Lys) with the HIV-1 A-loop.

The Role of the DNA Backbone in Minor-Groove Ligand Binding.

Molecular recognition between ligands and nucleic acids plays a key role in therapeutic activity. Some molecules interact with DNA in a nonbonded manner through intercalation or through the DNA grooves. The recognition of minor-groove binders is attributed to a set of hydrogen-bonding interactions between the binders and the hydrogen-bond-acceptor groups on the groove floor and walls. It is commonly considered that interactions with the sugar groups of the DNA backbone are insignificant and do not contribute to the binding affinity or the specificity. However, our group has found that the deoxyribose rings have a central function in the recognition and the intercalation of metal complexes into DNA. Herein, we determined the specific interactions between the binder CGP 40215A and the minor-groove atoms, based on the local properties of electron density. We found that specific interactions between the deoxyribose moiety within the backbone of DNA and the binder are essential for molecular recognition, and they are responsible for one third of the interaction energy.

Structural and dynamical instability of DNA caused by high occurrence of d5SICS and dNaM unnatural nucleotides.

In vivo inclusion of unnatural base pairs (UBPs) into functional DNA was recently reported for compounds 2,6-dimethyl-2H-isoquiniline-1-thione (d5SICS, X) and 2-methoxy-3-methylnaphthalene (dNaM, Y) in a modified E. coli strand, for which high-fidelity replication was observed. Yet, little is known about possible genetic consequences they have in largely substituted DNA. Using a converged microsecond long molecular dynamics (MD) simulation of the sequences 5'-GCGCAAXTTGCGC-3' and 5'-GCGCXAXTXGCGC-3', where X represents the UBP, we show that in the system with only a single XY UBP pair present, the global RMS deviation from canonical B-DNA in our control simulations is ∼3 Å and a fully converged ensemble is achieved within 2 µs. With three UBPs, deviation increases to ∼5 Å and convergence is not achieved within 10 µs of sampling time. With five UBPs, no convergence is observed and the double helix collapses into a globular structure. A fully optimized structure of the trimer d(GXC) was obtained using the density functional theory method B97D/cc-pVTZ and resulted in an RMSD value of ∼2 Å when compared to the most populated structure obtained from the MD simulations. Their calculated interaction energy is -3.7 kcal mol-1. It is thus unlikely that d5SICS and dNaM could be useful as tools in DNA manipulation. This theoretical methodology can be used in the de novo design of UBPs.

Intercalation processes of copper complexes in DNA.

The family of anticancer complexes that include the transition metal copper known as Casiopeínas® shows promising results. Two of these complexes are currently in clinical trials. The interaction of these compounds with DNA has been observed experimentally and several hypotheses regarding the mechanism of action have been developed, and these include the generation of reactive oxygen species, phosphate hydrolysis and/or base-pair intercalation. To advance in the understanding on how these ligands interact with DNA, we present a molecular dynamics study of 21 Casiopeínas with a DNA dodecamer using 10 µs of simulation time for each compound. All the complexes were manually inserted into the minor groove as the starting point of the simulations. The binding energy of each complex and the observed representative type of interaction between the ligand and the DNA is reported. With this extended sampling time, we found that four of the compounds spontaneously flipped open a base pair and moved inside the resulting cavity and four compounds formed stacking interactions with the terminal base pairs. The complexes that formed the intercalation pocket led to more stable interactions.

Convergence and reproducibility in molecular dynamics simulations of the DNA duplex d(GCACGAACGAACGAACGC).

BACKGROUND: The structure and dynamics of DNA are critically related to its function. Molecular dynamics simulations augment experiment by providing detailed information about the atomic motions. However, to date the simulations have not been long enough for convergence of the dynamics and structural properties of DNA. METHODS: Molecular dynamics simulations performed with AMBER using the ff99SB force field with the parmbsc0 modifications, including ensembles of independent simulations, were compared to long timescale molecular dynamics performed with the specialized Anton MD engine on the B-DNA structure d(GCACGAACGAACGAACGC). To assess convergence, the decay of the average RMSD values over longer and longer time intervals was evaluated in addition to assessing convergence of the dynamics via the Kullback-Leibler divergence of principal component projection histograms. RESULTS: These molecular dynamics simulations-including one of the longest simulations of DNA published to date at ~44µs-surprisingly suggest that the structure and dynamics of the DNA helix, neglecting the terminal base pairs, are essentially fully converged on the ~1-5µs timescale. CONCLUSIONS: We can now reproducibly converge the structure and dynamics of B-DNA helices, omitting the terminal base pairs, on the µs time scale with both the AMBER and CHARMM C36 nucleic acid force fields. Results from independent ensembles of simulations starting from different initial conditions, when aggregated, match the results from long timescale simulations on the specialized Anton MD engine. GENERAL SIGNIFICANCE: With access to large-scale GPU resources or the specialized MD engine "Anton" it is possible for a variety of molecular systems to reproducibly and reliably converge the conformational ensemble of sampled structures. This article is part of a Special Issue entitled: Recent developments of molecular dynamics.

A mixed DFT-MD methodology for the in silico development of drug releasing macrocycles. Calix and thia-calix[N]arenes as carriers for Bosutinib and Sorafenib.

Interaction energies between a family of 36 calix[n]arenes, their corresponding thia- analogues, and two commercially available second generation tyrosine kinase III inhibitors-Bosutinib and Sorafenib-were calculated through DFT methods at the B97D/6-31G(d,p) level of theory, based on Natural Population Analysis, for the in silico development of suitable drug carriers based on the aforementioned macrocycles which can increase their bioavailability and in turn their pharmaceutical efficiency. Molecular Dynamics simulations (production runs: +500 ns) using the General Amber Force Field were also carried out in order to assess the releasing process of these drugs in an explicit aqueous environment. In total, 144 host-guest complexes are examined. According to our results, five-membered -SO3H and i-Pr functionalized-calixarenes are the best candidates for Sorafenib-carriers while six-membered ones -SO3H and C2H4NH2 functionalized- are the lead candidates for Bosutinib-carriers.

Using Wavelet Analysis To Assist in Identification of Significant Events in Molecular Dynamics Simulations.

Long time scale molecular dynamics (MD) simulations of biological systems are becoming increasingly commonplace due to the availability of both large-scale computational resources and significant advances in the underlying simulation methodologies. Therefore, it is useful to investigate and develop data mining and analysis techniques to quickly and efficiently extract the biologically relevant information from the incredible amount of generated data. Wavelet analysis (WA) is a technique that can quickly reveal significant motions during an MD simulation. Here, the application of WA on well-converged long time scale (tens of µs) simulations of a DNA helix is described. We show how WA combined with a simple clustering method can be used to identify both the physical and temporal locations of events with significant motion in MD trajectories. We also show that WA can not only distinguish and quantify the locations and time scales of significant motions, but by changing the maximum time scale of WA a more complete characterization of these motions can be obtained. This allows motions of different time scales to be identified or ignored as desired.

Potential amoebicidal activity of hydrazone derivatives: synthesis, characterization, electrochemical behavior, theoretical study and evaluation of the biological activity.

Four new hydrazones were synthesized by the condensation of the selected hydrazine and the appropriate nitrobenzaldehyde. A complete characterization was done employing 1H- and 13C-NMR, electrochemical techniques and theoretical studies. After the characterization and electrochemical analysis of each compound, amoebicidal activity was tested in vitro against the HM1:IMSS strain of Entamoeba histolytica. The results showed the influence of the nitrobenzene group and the hydrazone linkage on the amoebicidal activity. meta-Nitro substituted compound 2 presents a promising amoebicidal activity with an IC50 = 0.84 µM, which represents a 7-fold increase in cell growth inhibition potency with respect to metronidazole (IC50 = 6.3 µM). Compounds 1, 3, and 4 show decreased amoebicidal activity, with IC50 values of 7, 75 and 23 µM, respectively, as a function of the nitro group position on the aromatic ring. The observed differences in the biological activity could be explained not only by the redox potential of the molecules, but also by their capacity to participate in the formation of intra- and intermolecular hydrogen bonds. Redox potentials as well as the amoebicidal activity can be described with parameters obtained from the DFT analysis.

The π-back-bonding modulation and its impact in the electronic properties of Cu(II) antineoplastic compounds: an experimental and theoretical study.

A complete study of the electronic density distribution in antineoplastic mixed chelate complexes of the type [Cu(N-N)(glycinate)H2 O]NO3 (N-N=2,2'-bipyridine (bpy) (1), 4,4'-dimethyl-bpy (2), 5,5'-dimethyl-2,2'-bipyridine (3), 1,10-phenanthroline (phen) (4), 4-methyl-phen (5); 5-methyl-phen (6); 4,7-dimethyl-phen (7), 5,6-dimethyl-phen (8), and 3,4,7,8-tetramethyl-phen (9)), a family known as Casiopeínas, was carried out with cyclic voltammetry, EPR, and computational methods. Crystal structures of 1⋅H2 O, 2⋅H2 O, 3⋅H2 O, 6⋅H2 O, and 8⋅H2 O show the variability in the geometries adopted by the copper compounds in the solid state. Experimental properties are described employing electronic descriptors obtained from computational methods. The main descriptors found were: The total electronic population of the metal ion (N(Cu)), delocalization of the metal ion electrons over the donor atoms (Δ(Cu)), atomic dipolar moment (µ(Cu)), and the atomic quadrupole moment (Q(Cu)). It was found that π-back-bonding is the principal factor that modulates the distribution of the electron density around the metal ion. The electronic descriptors obtained from the computational approach can be used as electronic descriptors of inorganic compounds that have shown antiproliferative activities instead the experimental data, aiding the rational design of good candidates of metal-based drugs.

Molecular dynamics of structural effects of reactive carbonyl species derivate of lipid peroxidation on bovine serum albumin.

BACKGROUND: Serum albumin is the most abundant protein in the Mammalia blood plasma at where plays a decisive role in the transport wide variety of hydrophobic ligands. BSA undergoes oxidative modifications like the carbonylation by the reactive carbonyl species (RCSs) 4-hydroxy-2-nonenal (HNE), 4 hydroxy-2-hexenal (HHE), malondialdehyde (MDA) and 4-oxo-2-nonenal (ONE), among others. The structural and functional changes induced by protein carbonylation have been associated with the advancement of neurodegenerative, cardiovascular, metabolic and cancer diseases. METHODS: To elucidate structural effects of protein carbonylation with RCSs on BSA, parameters for six new non-standard amino acids were designated and molecular dynamics simulations of its mono­carbonylated-BSA systems were conducted in the AMBER force field. Trajectories were evaluated by RMSD, RMSF, PCA, RoG and SASA analysis. RESULTS: An increase in the conformational instability for all proteins modified with local changes were observed, without significant changes on the BSA global three-dimensional folding. A more relaxed compaction level and major solvent accessible surface area for modified systems was found. Four regions of high molecular fluctuation were identified in all modified systems, being the subdomains IA and IIIB those with the most remarkable local conformational changes. Regarding essential modes of domain movements, it was evidenced that the most representatives were those related to IA subdomain, while IIIB subdomain presented discrete changes. CONCLUSIONS: RCSs induces local structural changes on mono­carbonylated BSA. Also, this study extends our knowledge on how carbonylation by RCSs induce structural effects on proteins.

Structures and Dynamics of DNA Mini-Dumbbells Are Force Field Dependent.

Flexible nucleic acid structures can be challenging to accurately resolve with currently available experimental structural determination techniques. As an alternative, molecular dynamics (MD) simulations can provide a window into understanding the unique dynamics and population distributions of these biomolecules. Previously, molecular dynamics simulations of noncanonical (non-duplex) nucleic acids have proven difficult to accurately model. With a new influx of improved nucleic acid force fields, achieving an in-depth understanding of the dynamics of flexible nucleic acid structures may be achievable. In this project, currently available nucleic acid force fields are evaluated using a flexible yet stable model system: the DNA mini-dumbbell. Prior to MD simulations, nuclear magnetic resonance (NMR) re-refinement was accomplished using improved refinement techniques in explicit solvent to yield DNA mini-dumbbell structures with better agreement between the newly determined PDB snapshots, with the NMR data itself, as well as the unrestrained simulation data. Starting from newly determined structures, a total aggregate of over 800 µs of production data between 2 DNA mini-dumbbell sequences and 8 force fields was collected to compare to these newly refined structures. The force fields tested spanned from traditional Amber force fields: bsc0, bsc1, OL15, and OL21 to Charmm force fields: Charmm36 and the Drude polarizable force field, as well as force fields from independent developers: Tumuc1 and CuFix/NBFix. The results indicated slight variations not only between the different force fields but also between the sequences as well. Given our previous experiences with high populations of potentially anomalous structures in RNA UUCG tetraloops and in various tetranucleotides, we expected the mini-dumbbell system to be challenging to accurately model. Surprisingly, many of the recently developed force fields generated structures in good agreement with experiments. Yet, each of the force fields provided a different distribution of potentially anomalous structures.

Assessment of A- to B- DNA Transitions Utilizing the Drude Polarizable Force Field.

In addition to the well-characterized B-form of DNA, duplex DNA can adopt various conformations, such as A or Z-DNA. Though less common, these structures can be induced biologically through protein or ligand interactions or experimentally with niche environmental conditions, such as high salt concentrations or in mixed water-ethanol. Reproducing these alternate structures through molecular dynamics simulations in recent years has been quite challenging with the currently available force fields, simulation techniques, and time scales. In this study, the Drude polarizable force field is tested for its ability to facilitate transitions between A-DNA and B-DNA or maintain A-DNA. Though transitions away from B-DNA were observed in high concentrations of ethanol, the resulting structures had hybrid properties taken from both B-DNA and A-DNA structures. This was also true for A-DNA in ethanol, which lost some of the A-DNA properties that it was expected to maintain. When B-DNA was tested in high salt environments, the resulting B-DNA structures showed no distinguishable differences with the increasing salt concentrations tested. These results with the Drude FF and recent results with additive force fields suggest that at present the current additive and polarizable force fields do not facilitate a complete transition between B- to A-DNA conformations under the conditions simulated. At present, the Drude FF favors A-B DNA hybrid structures when simulated in nonphysiological conditions.

Assessing the Current State of Amber Force Field Modifications for DNA─2023 Edition.

Advances in molecular dynamics (MD) software alongside enhanced computational power and hardware have allowed for MD simulations to significantly expand our knowledge of biomolecular structure, dynamics, and interactions. Furthermore, it has allowed for the extension of conformational sampling times from nanoseconds to the microsecond level and beyond. This has not only made convergence of conformational ensembles through comprehensive sampling possible but consequently exposed deficiencies and allowed the community to overcome limitations in the available force fields. The reproducibility and accuracy of the force fields are imperative in order to produce biologically relevant data. The Amber nucleic acid force fields have been used widely since the mid-1980s, and improvement of these force fields has been a community effort with several artifacts revealed, corrected, and reevaluated by various research groups. Here, we focus on the Amber force fields for use with double-stranded DNA and present the assessment of two recently developed force field parameter sets (OL21 and Tumuc1). Extensive MD simulations were performed with six test systems and two different water models. We observe the improvement of OL21 and Tumuc1 compared to previous generations of the Amber DNA force. We did not detect any significant improvement in the performance of Tumuc1 compared to OL21 despite the reparameterization of bonded force field terms in the former; however, we did note discrepancies in Tumuc1 when modeling Z-DNA sequences.

Cysteine carbonylation with reactive carbonyl species from lipid peroxidation induce local structural changes on thioredoxin active site.

Protein oxidative modifications with reactive carbonyl species (RCS) is directly linked to metabolic processes in premature aging, cancer, neurodegenerative and infectious diseases. RCS as 4-Hydroxy-2-nonal (HNE), 4-Hydroxy-2-hexenal (HHE), 4-Oxo-2-nonenal (ONE) and Malondialdehyde (MDA) attack nucleophilic amino acids residues forming irreversible adducts with proteins as Thioredoxins (Trx). This is a class of small thiol oxide-reductases playing a key role in redox signaling and oxidative stress responses in mammals. Although proteomic studies have identified to Cys-32 residue as a target of HNE attack that inhibit its enzymatic activity, how this carbonylation affects its structure and dynamic behavior at the atomic level is unknown. Even more, the molecular bases for the atomistic behavior of these modified proteins have not been completely understood. We present molecular dynamics simulations of Trx-modified with four different RCS to analyze its global and local structural effects. For this, parameters supported in the AMBER force fields were built and validated for three non-natural cysteine residues modified with HHE, ONE and MDA. Results obtained showed a slight change in the global conformational stability of Trx modified with HNE and MDA, establishing that all modified proteins presented local regions of high mobility in the modified catalytic site and some regions far from the modification area. In addition, essential domain movement modes evidences that proteins modified with the RCS assayed induce changes in conformational flexibility. Finally, these data showed that the given conformational changes did not caused global changes in proteins but rather localized changes in particular regions.