Biomedical subjects

N Foloppe

Publications and source records attributed to N Foloppe.

8 recordsLinked to original sources

PubMed2006

Towards predictive ligand design with free-energy based computational methods?

The accurate prediction of ligand-biopolymer binding affinities is of general interest to medicinal chemistry, as well as to the broader field of molecular recognition. The ability to predict computationally the thermodynamics of these molecular recognition processes has been relatively weak until recently, however, continued developments on several fronts are extending the scope of applicability of these methods. The rapid growth in the number of protein-ligand structures has initially led to the development of a range of empirical scoring functions based on relatively simple descriptions of intermolecular interactions. These methods have had some success in ranking binding affinities when tuned to particular protein systems or in rather qualitative estimates of molecular fit in fast docking calculations. However, they are too unreliable for more detailed, quantitative, assessment and comparison of binding affinities. Physics-based free energy calculations are in principle more general and have the potential to be significantly more accurate. These approaches have seen steady development over many years and rely on carefully calibrated molecular energy functions (force-fields), simulations of the systems with explicit solvent, and the coming-of-age of continuum solvation models. In addition to the initially developped Free Energy Perturbation (FEP) and Thermodynamic Integration (TI) methods, new approaches include the Molecular Mechanics-Poisson-Boltzmann Surface Area (MM-PBSA) and the Linear Interaction Energy (LIE) approaches. This review concentrates on MM-PBSA and LIE, and their variants. The routine application of these calculations is becoming possible because of enhanced computational hardware and the development of a range of computational chemistry tools. This review addresses: i) the basic principles behind free energy calculations ii) recent methodological advances iii) comparisons of predicted and experimentally determined affinities iv) the uncertainties and limitations of both the computational and experimental data v) areas where progress can be made vi) the practicality of applying the methods at the different stages of the drug discovery and optimization process.

Animals↗

PubMed2001

Structure, dynamics and electrostatics of the active site of glutaredoxin 3 from Escherichia coli: comparison with functionally related proteins.

The chemistry of active-site cysteine residues is central to the activity of thiol-disulfide oxidoreductases of the thioredoxin superfamily. In these reactions, a nucleophilic thiolate is required, but the associated pK(a) values differ vastly in the superfamily, from less than 4 in DsbA to greater than 7 in Trx. The factors that stabilize this thiolate are, however, not clearly established. The glutaredoxins (Grxs), which are members of this superfamily, contain a Cys-Pro-Tyr-Cys motif in their active site. In reduced Grxs, the pK(a) of the N-terminal active-site nucleophilic cysteine residue is lowered significantly, and the stabilization of the corresponding thiolate is expected to influence the redox potential of these enzymes. Here, we use a combination of long molecular dynamics (MD) simulations, pK(a) calculations, and experimental investigations to derive the structure and dynamics of the reduced active site from Escherichia coli Grx3, and investigate the factors that stabilize the thiolate. Several different MD simulations converged toward a consensus conformation for the active-site cysteine residues (Cys11 and Cys14), after a number of local conformational changes. Key features of the model were tested experimentally by measurement of NMR scalar coupling constants, and determination of pK(a) values of selected residues. The pK(a) values of the Grx3 active-site residues were calculated during the MD simulations, and support the underlying structural model. The structure of Grx3, in combination with the pK(a) calculations, indicate that the pK(a) of the N-terminal active-site cysteine residue in Grx3 is intermediate between that of its counterpart in DsbA and Trx. The pK(a) values in best agreement with experiment are obtained with a low (<4) protein dielectric constant. The calculated pK(a) values fluctuate significantly in response to protein dynamics, which underscores the importance of the details of the underlying structures when calculating pK(a) values. The thiolate of Cys11 is stabilized primarily by direct hydrogen bonding with the amide protons of Tyr13 and Cys14 and the thiol proton of Cys14, rather than by long-range interactions from charged groups or from a helix macrodipole. From the comparison of reduced Grx3 with other members of the thioredoxin superfamily, a unifying theme for the structural basis of thiol pK(a) differences in this superfamily begins to emerge.

Bacterial Proteins↗

PubMed2000

New insights into the structure of abasic DNA from molecular dynamics simulations.

Abasic (AP) sites constitute a common form of DNA damage, arising from the spontaneous or enzymatic breakage of the N-glycosyl bond and the loss of a nucleotide base. To examine the effects of such damage on DNA structure, especially in the vicinity of the abasic sugar, four 1.5 ns molecular dynamics simulations of double-helical DNA dodecamers with and without a single abasic (tetrahydrofuran, X) lesion in a 5'-d(CXT) context have been performed and analyzed. The results indicate that the abasic site does not maintain a hole or gap in the DNA, but instead perturbs the canonical structure and induces additional flexibility close to the abasic site. In the apurinic simulations (i.e., when a pyrimidine is opposite the AP site), the abasic sugar flipped in and out of the minor groove, and the gap was water filled, except during the occurrence of a novel non-Watson-Crick C-T base pair across the abasic site. The apyrimidinic gap was not penetrated by water until the abasic sugar flipped out and remained extrahelical. Both AP helices showed kinks of 20-30 degrees at the abasic site. The Watson-Crick hydrogen bonds are more transient throughout the DNA double helices containing an abasic site. The abasic sugar displayed an unusually broad range of sugar puckers centered around the northern pucker. The increased motion of the bases and backbone near the abasic site appear to correlate with sequence-dependent helical stability. The data indicate that abasic DNA contorts more easily and in specific ways relative to unmodified DNA, an aspect likely to be important in abasic site recognition and hydrolysis.

Apurinic Acid↗

PubMed1999

Intrinsic conformational properties of deoxyribonucleosides: implicated role for cytosine in the equilibrium among the A, B, and Z forms of DNA.

Structural properties of biomolecules are dictated by their intrinsic conformational energetics in combination with environmental contributions. Calculations using high-level ab initio methods on the deoxyribonucleosides have been performed to investigate the influence of base on the intrinsic conformational energetics of nucleosides. Energy minima in the north and south ranges of the deoxyribose pseudorotation surfaces have been located, allowing characterization of the influence of base on the structures and energy differences between those minima. With all bases, chi values associated with the south energy minimum are lower than in canonical B-DNA, while chi values associated with the north energy minimum are close to those in canonical A-DNA. In deoxycytidine, chi adopts an A-DNA conformation in both the north and south energy minima. Energy differences between the A and B conformations of the nucleosides are <0.5 kcal/mol in the present calculations, except with deoxycytidine, where the A form is favored by 2.3 kcal/mol, leading the intrinsic conformational energetics of GC basepairs to favor the A form of DNA by 1.5 kcal/mol as compared with AT pairs. This indicates that the intrinsic conformational properties of cytosine at the nucleoside level contribute to the A form of DNA containing predominately GC-rich sequences. In the context of a B versus Z DNA equilibrium, deoxycytidine favors the Z form over the B form by 1.6 kcal/mol as compared with deoxythymidine, suggesting that the intrinsic conformational properties of cytosine also contribute to GC-rich sequences occurring in Z DNA with a higher frequency than AT-rich sequences. Results show that the east pseudorotation energy barrier involves a decrease in the furanose amplitude and is systematically lower than the inversion barrier, with the energy differences influenced by the base. Energy barriers going from the south (B form) sugar pucker to the east pseudorotation barrier are lower in pyrimidines as compared with purines, indicating that the intrinsic conformational properties associated with base may also influence the sugar pseudorotational population distribution seen in DNA crystal structures and the kinetics of B to A transitions. The present work provides evidence that base composition, in addition to base sequence, can influence DNA conformation.

Base Composition↗

PubMed1995

Structural model of the photosynthetic reaction center of Rhodobacter capsulatus.

The reaction center (RC) from the photosynthetic bacterium Rhodobacter (Rb.) capsulatus has been the subject of a considerable amount of molecular biological and spectroscopic work aimed at improving our understanding of the primary steps of photosynthesis. However, no three-dimensional structure is available for this protein. We present here a model obtained by combining information from the structure of the highly homologous RC from Rhodopseudomonas (Rps.) viridis with molecular mechanics and simulated annealing calculations. In the Rb. capsulatus model the orientations of the bacteriochlorophyll monomer and the bacteriopheophytin on the branch inactive in electron transfer differ significantly from those in the RCs of Rps. viridis and Rb. sphaeroides. The bacteriopheophytin orientational difference is in good accord with previous linear dichroism measurements. A comparison is made of interactions between the pigments and the protein environment that may be of functional significance in Rps. viridis, Rb. sphaeroides, and Rb. capsulatus.

Amino Acids↗

PubMedDate not supplied

Development and current status of the CHARMM force field for nucleic acids.

The CHARMM27 all-atom force field for nucleic acids represents a highly optimized model for investigations of nucleic acids via empirical force field calculations. The force field satisfactorily treats the A, B, and Z forms of DNA as well as RNA, and it also useful for nucleosides and nucleotides. In addition, it is compatible with the CHARMM force fields for proteins and lipids, allowing for simulation studies of heterogeneous systems.

Biophysics↗

PubMedDate not supplied

Ab initio conformational analysis of nucleic acid components: intrinsic energetic contributions to nucleic acid structure and dynamics.

In recent years, the use of high-level ab initio calculations has allowed for the intrinsic conformational properties of nucleic acid building blocks to be revisited. This has provided new insights into the intrinsic conformational energetics of these compounds and its relationship to nucleic acids structure and dynamics. In this article we review recent developments and present new results. New data include comparison of various levels of theory on conformational properties of nucleic acid building blocks, calculations on the abasic sugar, known to occur in vivo in DNA, on the TA conformation of DNA observed in the complex with the TATA box binding protein, and on inosine. Tests of the Hartree-Fock (HF), second-order Møller-Plesset (MP2), and Density Functional Theory/Becke3, Lee, Yang and Par (DFT/B3LYP) levels of theory show the overall shape of backbone torsional energy profiles (for gamma, epsilon, and chi) to be similar for the different levels, though some systematic differences are identified between the MP2 and DFT/B3LYP profiles. The east pseudorotation energy barrier in deoxyribonucleosides is also sensitive to the level of theory, with the HF and DFT/B3LYP east barriers being significantly lower (approximately 2.5 kcal/mol) than the MP2 counterpart (approximately 4.0 kcal/mol). Additional calculations at various levels of theory suggest that the east barrier in deoxyribonucleosides is between 3.0 and 4.0 kcal/mol. In the abasic sugar, the west pseudorotation energy barrier is found to be slightly lower than the east barrier and the south pucker is favored more than in standard nucleosides. Results on the TA conformation suggest that, at the nucleoside level, this conformation is significantly destabilized relative to the global energy minimum, or relative to the A- and B-DNA conformations. Deoxyribocytosine would destabilize the TA conformation more than other bases relative to the A-DNA conformation, but not relative to the B-DNA conformation.

DNA↗

PubMedDate not supplied

Drug-like annotation and duplicate analysis of a 23-supplier chemical database totalling 2.7 million compounds.

We have implemented five drug-like filters, based on 1D and 2D molecular descriptors, and applied them to characterize the drug-like properties of commercially available chemical compounds. In addition to previously published filters (Lipinski and Veber), we implemented a filter for medicinal chemistry tractability based on lists of chemical features drawn up by a panel of medicinal chemists. A filter based on the modeling of aqueous solubility (>1 microM) was derived in-house, as well as another based on the modeling of Caco-2 passive membrane permeability (>10 nm/s). A library of 2.7 million compounds was collated from the 23 compound suppliers and analyzed with these filters, highlighting a tendency toward highly lipophilic compounds. The library contains 1.6 M unique structures, of which 37% (607,223) passed all five drug-like filters. None of the 23 suppliers provides all the members of the drug-like subset, emphasizing the benefit of considering compounds from various compound suppliers as a source of diversity for drug discovery.

Journal Article↗