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.

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.

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.

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.

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.

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.

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.

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.

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.

DNA sequence-specific recognition by a transcriptional regulator requires indirect readout of A-tracts.

The bacteriophage Ø29 transcriptional regulator p4 binds to promoters of different intrinsic activities. The p4-DNA complex contains two identical protomers that make similar interactions with the target sequence 5'-AACTTTTT-15 bp-AAAATGTT-3'. To define how the various elements in the target sequence contribute to p4's affinity, we studied p4 binding to a series of mutated binding sites. The binding specificity depends critically on base pairs of the target sequence through both direct as well as indirect readout. There is only one specific contact between a base and an amino acid residue; other contacts take place with the phosphate backbone. Alteration of direct amino acid-base contacts, or mutation of non-contacted A.T base pairs at A-tracts abolished binding. We generated three 5 ns molecular dynamics (MD) simulations to investigate the basis for the p4-DNA complex specificity. Recognition is controlled by the protein and depends on DNA dynamic properties. MD results on protein-DNA contacts and the divergence of p4 affinity to modified binding sites reveal an inherent asymmetry, which is required for p4-specific binding and may be crucial for transcription regulation.

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.

A Mg2+-induced conformational switch rendering a competent DNA polymerase catalytic complex.

The structural and dynamical changes occurring before nucleotide addition were studied using molecular dynamics (MD) simulations of human immunodeficiency virus type 1 (HIV-1) reverse transcriptase (RT) complexes containing one or two Mg2+ ions in the presence of dNTP. Our models revealed that the formation of a catalytically competent DNA polymerase complex required subtle rearrangements at the catalytic site A, which occurred only when an Mg2+ ion was bound. This model has been validated using pre-steady-state kinetics to show that free Mg2+ is necessary to obtain a catalytically competent polymerase. Kinetic studies carried out with Be2+ as a cofactor permitted the functional discrimination between metal sites A and B. At low concentrations, Be2+ increased the catalytic efficiency of the polymerase, while at higher concentrations, it competed with Mg2+ for binding to site A, and inhibited DNA polymerization. In agreement with experimental data, MD simulations revealed that the catalytic attack distance between the 3-OH of the primer and the phosphorus in complexes containing Be2+ instead of Mg2+ at site A was above 4.5 A. Our findings provide a detailed description of the mechanism of DNA polymerization and should be helpful to understand the molecular basis of DNA replication fidelity.

Residues in human respiratory syncytial virus P protein that are essential for its activity on RNA viral synthesis.

Human respiratory syncytial virus (HRSV) P protein, 241 amino acid long, is a structural homotetrameric phosphoprotein. Viral transcription and replication processes are dependent on functional P protein interactions inside viral ribonucleoprotein complexes (RNPs). Binding capacity to RNPs proteins and transcription and replication complementation analyses, using inactive P protein variants, have identified residues essential for functional interactions with itself, L, N and M2-1 proteins. P protein may establish some of these interactions as monomer, but efficient viral transcription and replication requires P protein oligomerization through the central region of the molecule. A structurally stable three-dimensional model has been generated in silico for this region (residues 98-158). Our analysis has indicated that P protein residues L135, D139, E140 and L142 are involved in homotetramerization. Additionally, the residues D136, S156, T160 and E179 appear to be essential for P protein activity on viral RNA synthesis and very high turnover phosphorylation at S143, T160 and T210 could regulate it. Thus, compounds targeted to those of these residues, located in the modeled three-dimensional structure, could have specific anti-HRSV effect.

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.

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.

Phosphorylation modulates the alpha-helical structure and polymerization of a peptide from the third tau microtubule-binding repeat.

Paired helical filaments (PHFs) isolated from patients with Alzheimer's disease (AD) mainly consist of the microtubule-associated protein tau in a hyperphosphorylated form. It has been found that PHFs are the first example of pathological protein aggregation associated with formation of alpha-helices [Biochemistry (2002) 41, 7150-5]. In an effort to investigate the interplay between phosphorylation and the putative role of short regions of alpha-helix in the polymerization of tau, we have focused on the region of tau encompassing residues 317 to 335. This region is able to form protein fibrils in vitro and has two serines that are often found phosphorylated in PHFs. Using trifluoroethanol as an indicator of the alpha-helix, we find that the stability of the alpha-helix conformation is enhanced by phosphorylation. Circular dichroism data show that the phosphorylated peptide in water presents a content in alpha-helix similar to the unphosphorylated peptide at 40% of trifluoroethanol. Phosphorylation also stimulates the effect of juglone in promoting the in vitro polymerization. Furthermore, Fourier transformed infrared spectroscopy of samples of phosphorylated peptide polymerized with juglone renders a spectrum with maxima at approximately 1665 and approximately 1675 cm(-1), which are suggestive of a mixture of turns and alpha-helix conformations. Our results provide a direct mechanistic connection between phosphorylation and polymerization in tau. The connection between phosphorylation and polymerization appears to involve formation of alpha-helix structure.

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.