Charge Engineering of the Nucleic Acid Binding Cleft of a Thermostable HIV-1 Reverse Transcriptase Reveals Key Interactions and a Novel Mechanism of RNase H Inactivation.

Coupled with PCR, reverse transcriptases (RTs) have been widely used for RNA detection and gene expression analysis. Increased thermostability and nucleic acid binding affinity are desirable RT properties to improve yields and sensitivity of these applications. The effects of amino acid substitutions in the RT RNase H domain were tested in an engineered HIV-1 group O RT, containing mutations K358R/A359G/S360A and devoid of RNase H activity due to the presence of E478Q (O3MQ RT). Twenty mutant RTs with Lys or Arg at positions interacting with the template-primer (i.e., at positions 473-477, 499-502 and 505) were obtained and characterized. Most of them produced significant amounts of cDNA at 37, 50 and 65 °C, as determined in RT-PCR reactions. However, a big loss of activity was observed with mutants A477K/R, S499K/R, V502K/R and Y505K/R, particularly at 65 °C. Binding affinity experiments confirmed that residues 477, 502 and 505 were less tolerant to mutations. Amino acid substitutions Q500K and Q500R produced a slight increase of cDNA synthesis efficiency at 50 and 65 °C, without altering the KD for model DNA/DNA and RNA/DNA heteroduplexes. Interestingly, molecular dynamics simulations predicted that those mutations inactivate the RNase H activity by altering the geometry of the catalytic site. Proof of this unexpected effect was obtained after introducing Q500K or Q500R in the wild-type HIV-1BH10 RT and mutant K358R/A359G/S360A RT. Our results reveal a novel mechanism of RNase H inactivation that preserves RT DNA binding and polymerization efficiency without substituting RNase H active site residues.

Proton Transfer in Guanine-Cytosine Base Pairs in B-DNA.

A double proton transfer reaction in a guanine-cytosine (GC) base pair has been proposed as a possible mechanism for rare tautomer (G*C*) formation and thus a source of spontaneous mutations. We analyze this system with free energy calculations based on extensive Quantum Mechanics/Molecular Mechanics simulations to properly consider the influence of the DNA biomolecular environment. We find that, although the G*C* rare tautomer is metastable in the gas phase, it is completely unstable in the conditions found in cells. Thus, our calculations show that a double proton reaction cannot be the source of spontaneous point mutations. We have also analyzed the intrabase H transfer reactions in guanine. Our results show that the DNA environment gives rise to a large free energy difference between the rare and canonical tautomers. These results show the key role of the DNA biological environment for the stability of the genetic code.

Heterochiral recognition among functionalized heptahelicenes on noble metal surfaces.

Chiral recognition among three differently functionalized heptahelicene derivatives on Ag(111) and Au(111) surfaces has been studied with scanning tunnelling microscopy. All three species were found to self-assemble into racemic zigzag structures, with alternation of (M)- and (P)-enantiomers.

Amino acid residues in HIV-2 reverse transcriptase that restrict the development of nucleoside analogue resistance through the excision pathway.

Nucleoside reverse transcriptase (RT) inhibitors (NRTIs) are the backbone of current antiretroviral treatments. However, the emergence of viral resistance against NRTIs is a major threat to their therapeutic effectiveness. In HIV-1, NRTI resistance-associated mutations either reduce RT-mediated incorporation of NRTI triphosphates (discrimination mechanism) or confer an ATP-mediated nucleotide excision activity that removes the inhibitor from the 3' terminus of DNA primers, enabling further primer elongation (excision mechanism). In HIV-2, resistance to zidovudine (3'-azido-3'-deoxythymidine (AZT)) and other NRTIs is conferred by mutations affecting nucleotide discrimination. Mutations of the excision pathway such as M41L, D67N, K70R, or S215Y (known as thymidine-analogue resistance mutations (TAMs)) are rare in the virus from HIV-2-infected individuals. Here, we demonstrate that mutant M41L/D67N/K70R/S215Y HIV-2 RT lacks ATP-dependent excision activity, and recombinant virus containing this RT remains susceptible to AZT inhibition. Mutant HIV-2 RTs were tested for their ability to unblock and extend DNA primers terminated with AZT and other NRTIs, when complexed with RNA or DNA templates. Our results show that Met73 and, to a lesser extent, Ile75 suppress excision activity when TAMs are present in the HIV-2 RT. Interestingly, recombinant HIV-2 carrying a mutant D67N/K70R/M73K RT showed 10-fold decreased AZT susceptibility and increased rescue efficiency on AZT- or tenofovir-terminated primers, as compared with the double-mutant D67N/K70R. Molecular dynamics simulations reveal that Met73influences ß3-ß4 hairpin loop conformation, whereas its substitution affects hydrogen bond interactions at position 70, required for NRTI excision. Our work highlights critical HIV-2 RT residues impeding the development of excision-mediated NRTI resistance.

Two-step ATP-driven opening of cohesin head.

The cohesin ring is a protein complex composed of four core subunits: Smc1A, Smc3, Rad21 and Stag1/2. It is involved in chromosome segregation, DNA repair, chromatin organization and transcription regulation. Opening of the ring occurs at the "head" structure, formed of the ATPase domains of Smc1A and Smc3 and Rad21. We investigate the mechanisms of the cohesin ring opening using techniques of free molecular dynamics (MD), steered MD and quantum mechanics/molecular mechanics MD (QM/MM MD). The study allows the thorough analysis of the opening events at the atomic scale: i) ATP hydrolysis at the Smc1A site, evaluating the role of the carboxy-terminal domain of Rad21 in the process; ii) the activation of the Smc3 site potentially mediated by the movement of specific amino acids; and iii) opening of the head domains after the two ATP hydrolysis events. Our study suggests that the cohesin ring opening is triggered by a sequential activation of the ATP sites in which ATP hydrolysis at the Smc1A site induces ATPase activity at the Smc3 site. Our analysis also provides an explanation for the effect of pathogenic variants related to cohesinopathies and cancer.

Effect of water-DNA interactions on elastic properties of DNA self-assembled monolayers.

DNA-water interactions have revealed as very important actor in DNA mechanics, from the molecular to the macroscopic scale. Given the particularly useful properties of DNA molecules to engineer novel materials through self-assembly and by bridging organic and inorganic materials, the interest in understanding DNA elasticity has crossed the boundaries of life science to reach also materials science and engineering. Here we show that thin films of DNA constructed through the self-assembly of sulfur tethered ssDNA strands demonstrate a Young's modulus tuning range of about 10 GPa by simply varying the environment relative humidity from 0% up to 70%. We observe that the highest tuning range occurs for ssDNA grafting densities of about 3.5 × 1013 molecules/cm 2, where the distance between the molecules maximizes the water mediated interactions between the strands. Upon hybridization with the complementary strand, the DNA self-assembled monolayers significantly soften by one order of magnitude and their Young's modulus dependency on the hydration state drastically decreases. The experimental observations are in agreement with molecular dynamics simulations.

Quantum Mechanics/Molecular Mechanics Free Energy Maps and Nonadiabatic Simulations for a Photochemical Reaction in DNA: Cyclobutane Thymine Dimer.

The absorption of ultraviolet radiation by DNA may result in harmful genetic lesions that affect DNA replication and transcription, ultimately causing mutations, cancer, and/or cell death. We analyze the most abundant photochemical reaction in DNA, the cyclobutane thymine dimer, using hybrid quantum mechanics/molecular mechanics (QM/MM) techniques and QM/MM nonadiabatic molecular dynamics. We find that, due to its double helix structure, DNA presents a free energy barrier between nonreactive and reactive conformations leading to the photolesion. Moreover, our nonadiabatic simulations show that most of the photoexcited reactive conformations return to standard B-DNA conformations after an ultrafast nonradiative decay to the ground state. This work highlights the importance of dynamical effects (free energy, excited-state dynamics) for the study of photochemical reactions in biological systems.

A Practical Quantum Mechanics Molecular Mechanics Method for the Dynamical Study of Reactions in Biomolecules.

Quantum mechanics/molecular mechanics (QM/MM) methods are excellent tools for the modeling of biomolecular reactions. Recently, we have implemented a new QM/MM method (Fireball/Amber), which combines an efficient density functional theory method (Fireball) and a well-recognized molecular dynamics package (Amber), offering an excellent balance between accuracy and sampling capabilities. Here, we present a detailed explanation of the Fireball method and Fireball/Amber implementation. We also discuss how this tool can be used to analyze reactions in biomolecules using steered molecular dynamics simulations. The potential of this approach is shown by the analysis of a reaction catalyzed by the enzyme triose-phosphate isomerase (TIM). The conformational space and energetic landscape for this reaction are analyzed without a priori assumptions about the protonation states of the different residues during the reaction. The results offer a detailed description of the reaction and reveal some new features of the catalytic mechanism. In particular, we find a new reaction mechanism that is characterized by the intramolecular proton transfer from O1 to O2 and the simultaneous proton transfer from Glu 165 to C2.

MEPSA: minimum energy pathway analysis for energy landscapes.

UNLABELLED: From conformational studies to atomistic descriptions of enzymatic reactions, potential and free energy landscapes can be used to describe biomolecular systems in detail. However, extracting the relevant data of complex 3D energy surfaces can sometimes be laborious. In this article, we present MEPSA (Minimum Energy Path Surface Analysis), a cross-platform user friendly tool for the analysis of energy landscapes from a transition state theory perspective. Some of its most relevant features are: identification of all the barriers and minima of the landscape at once, description of maxima edge profiles, detection of the lowest energy path connecting two minima and generation of transition state theory diagrams along these paths. In addition to a built-in plotting system, MEPSA can save most of the generated data into easily parseable text files, allowing more versatile uses of MEPSA's output such as the generation of molecular dynamics restraints from a calculated path. AVAILABILITY AND IMPLEMENTATION: MEPSA is freely available (under GPLv3 license) at: http://bioweb.cbm.uam.es/software/MEPSA/ CONTACT: pagomez@cbm.csic.es. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Torsion and curvature of FtsZ filaments.

FtsZ filaments participate in bacterial cell division, but it is still not clear how their dynamic polymerization and shape exert force on the underlying membrane. We present a theoretical description of individual filaments that incorporates information from molecular dynamic simulations. The structure of the crystallized Methanococcus jannaschii FtsZ dimer was used to model a FtsZ pentamer that showed a curvature and a twist. The estimated bending and torsion angles between monomers and their fluctuations were included in the theoretical description. The MD data also permitted positioning the curvature with respect to the protein coordinates and allowed us to explore the effect of the relative orientation of the preferred curvature with respect to the surface plane. We find that maximum tension is attained when filaments are firmly attached and oriented with their curvature perpendicular to the surface and that the twist serves as a valve to release or to tighten the tension exerted by the curved filaments on the membrane. The theoretical model also shows that the presence of torsion can explain the shape distribution of short filaments observed by Atomic Force Microscopy in previously published experiments. New experiments with FtsZ covalently attached to lipid membranes show that the filament on-plane curvature depends on lipid head charge, confirming the predicted monomer orientation effects. This new model underlines the fact that the combination of the three elements, filament curvature, twist and the strength and orientation of its surface attachment, can modulate the force exerted on the membrane during cell division.

Molecular basis of the association of H208Y and thymidine analogue resistance mutations M41L, L210W and T215Y in the HIV-1 reverse transcriptase of treated patients.

Thymidine analogue resistance mutations (TAMs) in HIV-1 reverse transcriptase (RT) associate in two clusters: (i) TAM1 (M41L, L210W and T215Y) and TAM2 (D67N, K70R, K219E/Q, and sometimes T215F). The amino acid substitution H208Y shows increased prevalence in patients treated with nucleoside analogues and is frequently associated with TAM1 mutations. We studied the molecular mechanism favoring the selection of H208Y in the presence of zidovudine, tenofovir and other nucleoside RT inhibitors (NRTIs). NRTI susceptibility was not affected by the addition of H208Y in phenotypic assays carried out in MT-4 cells using recombinant HIV-1 containing wild-type (subtype B, BH10), H208Y, M41L/L210W/T215Y or M41L/H208Y/L210W/T215Y RTs. However, enzymatic studies carried out with purified RTs revealed that in the presence of M41L/L210W/T215Y, H208Y increases the RT's ability to unblock and extend primers terminated with zidovudine, tenofovir and in a lesser extent, stavudine. These effects were observed with DNA/DNA complexes (but not with RNA/DNA) and resulted from the higher ATP-dependent excision activity of the M41L/H208Y/L210W/T215Y RT compared with the M41L/L210W/T215Y mutant. The increased rescue efficiency of the M41L/H208Y/L210W/T215Y RT was observed in the presence of ATP but not with GTP or ITP. Molecular dynamics studies predict an alteration of the stacking interactions between Tyr(215) and the adenine ring of ATP due to long-distance effects caused by tighter packaging of Tyr(208) and Trp(212). These studies provide a mechanistic explanation for the association of TAM-1 and H208Y mutations in viral isolates from patients treated with NRTIs.

fireball/amber: An Efficient Local-Orbital DFT QM/MM Method for Biomolecular Systems.

In recent years, quantum mechanics/molecular mechanics (QM/MM) methods have become an important computational tool for the study of chemical reactions and other processes in biomolecular systems. In the QM/MM technique, the active region is described by means of QM calculations, while the remainder of the system is described using a MM approach. Because of the complexity of biomolecules and the desire to achieve converged sampling, it is important that the QM method presents a good balance between accuracy and computational efficiency. Here, we report on the implementation of a QM/MM technique that combines a DFT approach specially designed for the study of complex systems using first-principles molecular dynamics simulations (fireball) with the amber force fields and simulation programs. We also present examples of the application of this QM/MM approach to three representative biomolecular systems: the analysis of the effect of electrostatic embedding in the behavior of a salt bridge between an aspartic acid and a lysine residue, a study of the intermediate states for the triosephosphate isomerase catalyzed conversion of dihydroxyacetone phosphate into glyceraldehyde 3-phosphate, and the detailed description, using DFT QM/MM molecular dynamics, of the cleavage of a phosphodiester bond in RNA catalyzed by the enzyme RNase A.

Simulation of catalytic water activation in mitochondrial F1-ATPase using a hybrid quantum mechanics/molecular mechanics approach: an alternative role for ß-Glu 188.

The use of quantum mechanics/molecular mechanics simulations to study the free energy landscape of the water activation at the catalytic site of mitochondrial F(1)-ATPase affords us insight into the generation of the nucleophile OH(-) prior to ATP hydrolysis. As a result, the ATP molecule was found to be the final proton acceptor. In the simulated pathway, the transfer of a proton to the nucleotide was not direct but occurred via a second water molecule in a manner similar to the Grotthuss mechanism proposed for proton diffusion. Residue ß-Glu 188, previously described as the putative catalytic base, was found to be involved in the stabilization of a transient hydronium ion during water activation. Simulations in the absence of the carboxylate moiety of ß-Glu 188 support this role.

Functional specificity of a protein-DNA complex mediated by two arginines bound to the minor groove.

A bacteriophage Ø29 transcriptional regulator, protein p4, interacts with its DNA target by employing two mechanisms: by direct readout of the chemical signatures of only one DNA base and by inducing local modification on the topology of short A tracts (indirect readout). p4 binds as a dimer to targets consisting of imperfect inverted repeats. Here we used molecular dynamic simulation to define interactions of a cluster of 12 positively charged amino acids of p4 with DNA and biochemical assays with modified DNA targets and mutated proteins to quantify the contribution of residues in the nucleoprotein complex. Our results show the implication of Arg54, with non-base-specific interaction in the central A tract, in p4 binding affinity. Despite being chemically equivalent and in identical protein monomers, the two Arg54 residues differed in their interactions with DNA. We discuss an indirect-readout mechanism for p4-DNA recognition mediated by dissimilar interaction of arginines penetrating the minor groove and the inherent properties of the A tract. Our findings extend the current understanding of protein-DNA recognition and contribute to the relevance of the sequence-dependent conformational malleability of the DNA, shedding light on the role of arginines in binding affinity. Characterization of mutant p4R54A shows that the residue is required for the activity of the protein as a transcriptional regulator.

Molecular dynamics simulation of GTPase activity in polymers of the cell division protein FtsZ.

FtsZ, the prokaryotic ortholog of tubulin, assembles into polymers in the bacterial division ring. The interfaces between monomers contain a GTP molecule, but the relationship between polymerization and GTPase activity is not unequivocally proven. A set of short FtsZ polymers were modelled and the formation of active GTPase structures was monitored using molecular dynamics. Only the interfaces nearest the polymer ends exhibited an adequate geometry for GTP hydrolysis. Simulated conversion of interfaces from close-to-end to internal position and vice versa resulted in their spontaneous rearrangement between active and inactive conformations. This predicted behavior of FtsZ polymer ends was supported by in vitro experiments.

Molecular dynamics analysis of conformational change of paramyxovirus F protein during the initial steps of membrane fusion.

The fusion of paramyxovirus to the cell membrane is mediated by fusion protein (F protein) present in the virus envelope, which undergoes a dramatic conformational change during the process. Unlike hemagglutinin in orthomyxovirus, this change is not mediated by an alteration of environmental pH, and its cause remains unknown. Steered molecular dynamics analysis leads us to suggest that the conformational modification is mediated only by stretching mechanical forces once the transmembrane fusion peptide of the protein is anchored to the cell membrane. Such elongating forces will generate major secondary structure rearrangement in the heptad repeat A region of the F protein; from ß-sheet conformation to an elongated coil and then spontaneously to an α-helix. In addition, it is proposed that the heptad repeat A region adopts a final three-helix coiled coil and that this structure appears after the formation of individual helices in each monomer.

The Role of Gln61 in HRas GTP hydrolysis: a quantum mechanics/molecular mechanics study.

Activation of the water molecule involved in GTP hydrolysis within the HRas·RasGAP system is analyzed using a tailored approach based on hybrid quantum mechanics/molecular mechanics (QM/MM) simulation. A new path emerges: transfer of a proton from the attacking water molecule to a second water molecule, then a different proton is transferred from this second water molecule to the GTP. Gln(61) will stabilize the transient OH(-) and H(3)O(+) molecules thus generated. This newly proposed mechanism was generated by using, for the first time to our knowledge, the entire HRas-RasGAP protein complex in a QM/MM simulation context. It also offers a rational explanation for previous experimental results regarding the decrease of GTPase rate found in the HRas Q61A mutant and the increase exhibited by the HRas Q61E mutant.

Technical phosphoproteomic and bioinformatic tools useful in cancer research.

Reversible protein phosphorylation is one of the most important forms of cellular regulation. Thus, phosphoproteomic analysis of protein phosphorylation in cells is a powerful tool to evaluate cell functional status. The importance of protein kinase-regulated signal transduction pathways in human cancer has led to the development of drugs that inhibit protein kinases at the apex or intermediary levels of these pathways. Phosphoproteomic analysis of these signalling pathways will provide important insights for operation and connectivity of these pathways to facilitate identification of the best targets for cancer therapies. Enrichment of phosphorylated proteins or peptides from tissue or bodily fluid samples is required. The application of technologies such as phosphoenrichments, mass spectrometry (MS) coupled to bioinformatics tools is crucial for the identification and quantification of protein phosphorylation sites for advancing in such relevant clinical research. A combination of different phosphopeptide enrichments, quantitative techniques and bioinformatic tools is necessary to achieve good phospho-regulation data and good structural analysis of protein studies. The current and most useful proteomics and bioinformatics techniques will be explained with research examples. Our aim in this article is to be helpful for cancer research via detailing proteomics and bioinformatic tools.

Thymidine analogue excision and discrimination modulated by mutational complexes including single amino acid deletions of Asp-67 or Thr-69 in HIV-1 reverse transcriptase.

Single amino acid deletions in the ß3-ß4 hairpin loop of human immunodeficiency virus type 1 (HIV-1) reverse transcriptase (RT) have been identified in heavily treated patients. The deletion of Asp-67 together with mutations T69G and K70R (Δ67 complex) are usually associated with thymidine analog resistance mutations (TAMs) (e.g. M41L, T215Y, etc.) while the deletion of Thr-69 (Δ69) is rarely found in isolates containing TAMs. Here, we show that the complex Δ67/T69G/K70R enhances ATP-dependent phosphorolytic activity on primers terminated with 3'-azido-3'-deoxythymidine (AZT) or 2',3'-didehydro-2',3'-dideoxythymidine (d4T) both in the presence or absence of TAMs (i.e. M41L/T215Y), while Δ69 (or the complex S68G/Δ69/K70G) antagonize the effects of TAMs in ATP-mediated excision. These effects are consistent with AZT susceptibility data obtained with recombinant HIV-1 bearing the relevant RTs. Molecular dynamics studies based on models of wild-type HIV-1 RT and mutant Δ69, Δ67/T69G/K70R, and D67N/K70R RTs support a relevant role for Lys/Arg-70 in the excision reaction. In Δ69 RT, the side chain of Lys-70 locates away from the putative pyrophosphate binding site. Therefore, its participation in interactions required for the excision reaction is unlikely. Our theoretical studies also suggest a role for Lys-219 in thymidine analog excision/discrimination. However, pre-steady-state kinetics revealed only minor differences in selectivity of AZT-triphosphate versus dTTP between deletion-containing RTs and their homologous enzymes having the K219E mutation. K219E reduced both ATP- and pyrophosphate-mediated excision of primers terminated with thymidine analogues, only when introduced in RTs bearing Δ69 or S68G/Δ69/K70G, providing further biochemical evidence that explains the lack of association of Δ69 and TAMs in HIV-1 isolates.

Nanomechanics of the cadherin ectodomain: "canalization" by Ca2+ binding results in a new mechanical element.

Cadherins form a large family of calcium-dependent cell-cell adhesion receptors involved in development, morphogenesis, synaptogenesis, differentiation, and carcinogenesis through signal mechanotransduction using an adaptor complex that connects them to the cytoskeleton. However, the molecular mechanisms underlying mechanotransduction through cadherins remain unknown, although their extracellular region (ectodomain) is thought to be critical in this process. By single molecule force spectroscopy, molecular dynamics simulations, and protein engineering, here we have directly examined the nanomechanics of the C-cadherin ectodomain and found it to be strongly dependent on the calcium concentration. In the presence of calcium, the ectodomain extends through a defined ("canalized") pathway that involves two mechanical resistance elements: a mechanical clamp from the cadherin domains and a novel mechanostable component from the interdomain calcium-binding regions ("calcium rivet") that is abolished by magnesium replacement and in a mutant intended to impede calcium coordination. By contrast, in the absence of calcium, the mechanical response of the ectodomain becomes largely "decanalized" and destabilized. The cadherin ectodomain may therefore behave as a calcium-switched "mechanical antenna" with very different mechanical responses depending on calcium concentration (which would affect its mechanical integrity and force transmission capability). The versatile mechanical design of the cadherin ectodomain and its dependence on extracellular calcium facilitate a variety of mechanical responses that, we hypothesize, could influence the various adhesive properties mediated by cadherins in tissue morphogenesis, synaptic plasticity, and disease. Our work represents the first step toward the mechanical characterization of the cadherin system, opening the door to understanding the mechanical bases of its mechanotransduction.