Biomedical subjects

Leandro Martínez

Publications and source records attributed to Leandro Martínez.

4 recordsLinked to original sources

PubMed2006

Molecular dynamics simulations of ligand dissociation from thyroid hormone receptors: evidence of the likeliest escape pathway and its implications for the design of novel ligands.

Steered molecular dynamics simulations of ligand dissociation from Thyroid hormone receptors indicate that dissociation is favored via rearrangements in a mobile part of the LBD comprising H3, the loop between H1 and H2, and nearby beta-sheets, contrary to current models in which the H12 is mostly involved. Dissociation is facilitated in this path by the interaction of the hydrophilic part of the ligand with external water molecules, suggesting strategies to enhance ligand binding affinity.

Computer Simulation↗

PubMed2005

Molecular dynamics simulations reveal multiple pathways of ligand dissociation from thyroid hormone receptors.

Nuclear receptor (NR) ligands occupy a pocket that lies within the core of the NR ligand-binding domain (LBD), and most NR LBDs lack obvious entry/exit routes upon the protein surface. Thus, significant NR conformational rearrangements must accompany ligand binding and release. The precise nature of these processes, however, remains poorly understood. Here, we utilize locally enhanced sampling (LES) molecular dynamics computer simulations to predict molecular motions of x-ray structures of thyroid hormone receptor (TR) LBDs and determine events that permit ligand escape. We find that the natural ligand 3,5,3'-triiodo-L-thyronine (T(3)) dissociates from the TRalpha1 LBD along three competing pathways generated through i), opening of helix (H) 12; ii), separation of H8 and H11 and the Omega-loop between H2 and H3; and iii), opening of H2 and H3, and the intervening beta-strand. Similar pathways are involved in dissociation of T(3) and the TRbeta-selective ligand GC24 from TRbeta; the TR agonist IH5 from the alpha- and beta-TR forms; and Triac from two natural human TRbeta mutants, A317T and A234T, but are detected with different frequencies in simulations performed with the different structures. Path I was previously suggested to represent a major pathway for NR ligand dissociation. We propose here that Paths II and III are also likely ligand escape routes for TRs and other NRs. We also propose that different escape paths are preferred in different situations, implying that it will be possible to design NR ligands that only associate stably with their cognate receptors in specific cellular contexts.

Algorithms↗

PubMed2004

Globally convergent trust-region methods for self-consistent field electronic structure calculations.

As far as more complex systems are being accessible for quantum chemical calculations, the reliability of the algorithms used becomes increasingly important. Trust-region strategies comprise a large family of optimization algorithms that incorporates both robustness and applicability for a great variety of problems. The objective of this work is to provide a basic algorithm and an adequate theoretical framework for the application of globally convergent trust-region methods to electronic structure calculations. Closed shell restricted Hartree-Fock calculations are addressed as finite-dimensional nonlinear programming problems with weighted orthogonality constraints. A Levenberg-Marquardt-like modification of a trust-region algorithm for constrained optimization is developed for solving this problem. It is proved that this algorithm is globally convergent. The subproblems that ensure global convergence are easy-to-compute projections and are dependent only on the structure of the constraints, thus being extendable to other problems. Numerical experiments are presented, which confirm the theoretical predictions. The structure of the algorithm is such that accelerations can be easily associated without affecting the convergence properties.

Journal Article↗

PubMed2003

Packing optimization for automated generation of complex system's initial configurations for molecular dynamics and docking.

Molecular Dynamics is a powerful methodology for the comprehension at molecular level of many chemical and biochemical systems. The theories and techniques developed for structural and thermodynamic analyses are well established, and many software packages are available. However, designing starting configurations for dynamics can be cumbersome. Easily generated regular lattices can be used when simple liquids or mixtures are studied. However, for complex mixtures, polymer solutions or solid adsorbed liquids (for example) this approach is inefficient, and it turns out to be very hard to obtain an adequate coordinate file. In this article, the problem of obtaining an adequate initial configuration is treated as a "packing" problem and solved by an optimization procedure. The initial configuration is chosen in such a way that the minimum distance between atoms of different molecules is greater than a fixed tolerance. The optimization uses a well-known algorithm for box-constrained minimization. Applications are given for biomolecule solvation, many-component mixtures, and interfaces. This approach can reduce the work of designing starting configurations from days or weeks to few minutes or hours, in an automated fashion. Packing optimization is also shown to be a powerful methodology for space search in docking of small ligands to proteins. This is demonstrated by docking of the thyroid hormone to its nuclear receptor.

Journal Article↗