Jesús I. Mendieta-Moreno

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.

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.

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 MEPSAs 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.

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.

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 × 1013molecules/cm2, 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.

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.

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.

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.

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.

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.