Adolescents' support for an outdoor smoke-free policy at sports clubs in the Netherlands.

INTRODUCTION: Outdoor smoke-free policies (SFPs) at sports clubs may contribute to the prevention of smoking among adolescents. Adolescents' support for such policy is important to its success. The aim of this study is to explore adolescents' perceptions with regard to an outdoor SFP at sports clubs in the Netherlands. METHODS: Focus group discussions (n=27) were held with 180 participants (aged 13-18 years) at 16 sports clubs. Thematic analysis was used to analyze the data. RESULTS: Participants generally supported an outdoor SFP at sports clubs. Five reasons for this support were reported: 1) children should not be exposed to smoking, 2) smoking and sports (clubs) do not fit together, 3) secondhand smoke is undesirable, 4) an outdoor SFP may enhance a sports club's image, and 5) an outdoor SFP contributes to the prevention of smoking. Some participants voiced considerations against an outdoor SFP, arguing amongst others, that smokers need to be taken into account, and that problems may occur with compliance and enforcement. Support for an outdoor SFP was stronger among participants at clubs with an outdoor SFP than among those without such policy. CONCLUSIONS: This study shows that adolescents generally support an outdoor smoke-free policy at sports clubs. After implementation, the outdoor SFP was generally experienced as a normal practice. These results could encourage sports clubs without an outdoor SFP to become smoke-free as well.

The screening of SLC6A8 deficiency among Estonian families with X-linked mental retardation.

The urinary creatine:creatinine (Cr:Crn) ratio was measured in males from 49 families with a family history compatible with X-linked mental retardation (XLMR) in order to estimate the prevalence of SLC6A8 deficiency in Estonia. We identified 11 boys from 9 families with an increased urinary Cr:Crn ratio (18%). In three related boys, a hemizygous missense mutation (c.1271G>A; p.Gly424Asp) was identified. Their mother was heterozygous for the same mutation. Although many missense mutations have been described, the p.Gly424Asp mutation has not been previously reported. The clinical expression varied widely among affected males of this family. Patients 1 and 3 had relatively mild clinical expression (mild mental retardation (MR) and attention deficit disorder), but patient 2 had all typical clinical signs of SLC6A8 defect such as moderate MR, autistic features, expressive dysphasia and epilepsy. Among our patients, we saw significant problems in speech and language development combined with attention and behavioural difficulties. The number of false-positive biochemical results with increased urinary Cr:Crn ratio was higher (18%) in our study than in previous reports (1.8­10%). We therefore suggest that repeated biochemical testing should be performed before DNA sequencing analysis. Our study suggests that 2% (95% confidence limits: 0.05­11.1%) of this Estonian XLMR panel are due to mutations in the SLC6A8, which is similar to the prevalence reported in other populations. We therefore conclude that creatine transporter deficiency is a relatively common genetic disorder in males with sporadic or familiar MR and diagnostic screening of them should always include screening for SLC6A8 deficiency.

Sudden death of a girl with Prader-Willi syndrome.

We report on the sudden death of a 3.5-year-old girl with Prader-Willi syndrome (PWS) and 15q11-q13 deletion. She suffered from severe chronic breathing disturbances and recurrent bronchitis. During an episode of acute bronchitis she had a cardiac arrest and died two months later of the sequelae. Brain CT imaging three weeks after the arrest showed bilateral symmetrical haemorrhages in the basal ganglia region. The spatial distribution of the haemorrhages can possibly suggest that the basal ganglia in PWS may be especially susceptible to hypoxemia.

[Effect of various suture strength factors on behavior of meniscus sutures in cyclic loading conditions]. / Der Einfluss verschiedener Nahtstärken auf das Verhalten von Meniskusnähten unter zyklischer Zugbelastung.

The aim of this study was to analyze meniscal sutures under cyclic loading conditions for different suture types (vertical and horizontal mattress sutures) and suture materials (absorbable monofilament sutures: PDS 2-0; PDS-0, and PDS-1 USP). Testing was performed on medial porcine menisci, using a well-established biomechanical testing model with a complete longitudinal tear 3 mm from the periphery of the meniscus. Sixty specimens were used. One suture was tested at a time. During cyclic testing 100 load cycles were applied with a crosshead speed of 50 mm/min. Three different maximum loads (10 N, 20 N, and 40 N) were used. The preload was set at 5 N. After cyclic loading, the specimens were loaded until failure. During cyclic loading, a gap appeared between the two parts of the meniscus, and partial tissue failures were observed at the surface of the meniscus. Gapping was more marked with higher loads and with the weaker suture material (p < 0.001). Using PDS 0 and PDS 1 sutures, less partial tissue failures were observed compared to PDS 2-0 (p < 0.001). The ultimate failure loads after cyclic loading were higher with PDS-0 and PDS-1 sutures. With these suture materials vertical sutures were stronger than horizontal sutures (p < 0.05). Using PDS 2-0 this difference could not be found. These results show that the primary strength of meniscal sutures depends on the suture material. The frequency and the amount of gapping and partial tissue failures, which can be observed under cyclic loading, are less distinct with PDS-0 and PDS-1 compared to PDS 2-0. From a biomechanical point of view, PDS 0 and PDS 1 sutures are recommended for meniscal sutures to guarantee a high primary stability, a small amount of gapping, and few partial tissue failures.

Familial Williams-Beuren syndrome.

Williams-Beuren syndrome (WBS) occurs sporadically; however, at least four familial cases of WBS have been described previously. We describe a mother and her son with typical WBS. The diagnosis of WBS in the son was confirmed by molecular cytogenetic analysis fluorescence in situ hybridization. He had a deletion of 7q11.23 at the ELN locus. The mother was diagnosed after the identification of WBS in her affected son. She is deceased and was thus not studied by FISH. However, her combined symptoms make it very clear that she had WBS. Two traits uncommon in WBS were observed, unilateral renal hypoplasia in the mother and a hemivertebra at L5 in the son.

Conformational requirements of suramin to target angiogenic growth factors.

The conformational requirements of suramin for its binding to basic fibroblast growth factor (bFGF) and platelet-derived growth factor (PDGF) were examined by molecular modeling and docking simulations using the conformational features of suramin determined by the present proton nuclear magnetic resonance (1H-NMR) studies and the crystal structures of growth factors reported previously. The assignment of resonances of suramin to individual protons was accomplished by the combined analysis of the coupling constants, two-dimensional correlated spectroscopy (COSY) and nuclear Overhauser effect spectroscopy (NOESY). The NOESY data obtained for suramin were utilized in a conformational search algorithm with constraints to generate a family of conformers which were further refined by restrained energy minimization. A family of nine conformers generated by restrained modeling falls primarily into one of two categories, either the conformer's two naphthyl rings are far apart, approximately 28-30 A, from one another or the conformer's two naphthyl rings are relatively close, approximately 16-18 A. The NMR data provide evidence for the presence of more than one conformer in solution. The modeling and docking simulation studies suggest that suramin binds efficiently to bFGF and PDGF by an induced-fit mechanism, wherein suramin complements bFGF or PDGF by adjusting its conformational freedom around the two pairs of single bonds that link the middle phenyl rings to the secondary amide backbone. The interaction of suramin with bFGF or PDGF primarily involves ion-pair, hydrophobic and hydrogen bonding interactions, in addition to van der Waals' contacts. The results indicate that suramin not only sterically blocks the receptor binding loop of the growth factors, but also competes for the binding sites of agonists such as heparin. The results suggest that suramin's propensity to bind to several polypeptides of varying size and structure is due to its conformational flexibility. Collectively, the data emphasize that conformationally constrained suramin analogs that selectively and competitively target angiogenic growth factors could be designed to reduce non-specific binding and, accordingly, toxicity.

Structure of an anti-HIV-1 hammerhead ribozyme complex with a 17-mer DNA substrate analog of HIV-1 gag RNA and a mechanism for the cleavage reaction: 750 MHz NMR and computer experiments.

The structure of an anti-HIV-1 ribozyme-DNA abortive substrate complex was investigated by 750 MHz NMR and computer modeling experiments. The ribozyme was a chimeric molecule with 30 residues-18 DNA nucleotides, and 12 RNA residues in the conserved core. The DNA substrate analog had 17 residues. The chimeric ribozyme and the DNA substrate formed a shortened ribozyme-abortive substrate complex of 47 nucleotides with two DNA stems (stems I and III) and a loop consisting of the conserved core residues. Circular dichroism spectra showed that the DNA stems assume A-family conformation at the NMR concentration and a temperature of 15 degrees C, contrary to the conventional wisdom that DNA duplexes in aqueous solution populate entirely in the B-form. It is proposed that the A-family RNA residues at the core expand the A-family initiated at the core into the DNA stems because of the large free energy requirement for the formation of A/B junctions. Assignments of the base H8/H6 protons and H1' of the 47 residues were made by a NOESY walk. In addition to the methyl groups of all T's, the imino resonances of stems I and III and AH2's were assigned from appropriate NOESY walks. The extracted NMR data along with available crystallographic data, were used to derive a structural model of the complex. Stems I and III of the final model displayed a remarkable similarity to the A form of DNA; in stem III, a GC base pair was found to be moving into the floor of the minor groove defined by flanking AT pairs; data suggest the formation of a buckled rhombic structure with the adjacent pair; in addition, the base pair at the interface of stem III and the loop region displayed deformed geometry. The loop with the catalytic core, and the immediate region of the stems displayed conformational multiplicity within the NMR time scale. A catalytic mechanism for ribozyme action based on the derived structure, and consistent with biochemical data in the literature, is proposed. The complex between the anti HIV-1 gag ribozyme and its abortive DNA substrate manifests in the detection of a continuous track of A.T base pairs; this suggests that the interaction between the ribozyme and its DNA substrate is stronger than the one observed in the case of the free ribozyme where the bases in stem I and stem III regions interact strongly with the ribozyme core region (Sarma, R. H., et al. FEBS Letters 375, 317-23, 1995). The complex formation provides certain guidelines in the design of suitable therapeutic ribozymes. If the residues in the ribozyme stem regions interact with the conserved core, it may either prevent or interfere with the formation of a catalytically active tertiary structure.

[Echinococcosis of the heart]. / Echinokokkose des Herzens.

Human disease caused by echinococcus granulosus is rare in Northern Europe, but is seen increasingly due to migration of labour and of tourism to endemic areas e. g. Southern Europe and the Middle East. Echinococcosis of the heart occurs in only 0.5 to 2% of cases of echinococcosis. Diagnosis is difficult because of the long latency between infection and manifestation of disease, and also because symptoms are unspecific. It is a diagnosis "to think of". Apart from a thorough history serological tests and medical imaging (ultrasound, X-rays, computed tomography) are used in the diagnosis of echinococcosis. Surgery is the preferred therapy, if not feasible medical treatment with benzimidazoles (Albendazole and Mebendazole) is effective. We describe a case of echinococcosis of the heart as the only manifestation of the disease: a 34-year-old male Turkish patient came in severe cardiac shock due to cardiac tamponade to the emergency department. Immediate heart surgery revealed a pericardium filled with cysts of echinococcus granulosus, infiltrating the left ventricle of the heart (Figure 1). The cysts were removed and the patient was put on a 2-year course of Mebendazole. The patient is free of relapse for more than 6 years now and he is considered to be cured.

Modeling DNA hydration: comparison of calculated and experimental hydration properties of nuclic acid bases.

Hydration properties of individual nucleic acid bases were calculated and compared with the available experimental data. Three sets of classical potential functions (PF) used in simulations of nucleic acid hydration were juxtaposed: (i) the PF developed by Poltev and Malenkov (PM), (ii) the PF of Weiner and Kollman (WK), which together with Jorgensen's TIP3P water model are widely used in the AMBER program, and (iii) OPLS (optimized potentials for liquid simulations) developed by Jorgensen (J). The global minima of interaction energy of single water molecules with all the natural nucleic acid bases correspond to the formation of two water-base hydrogen bonds (water bridging of two hydrophilic atoms of the base). The energy values of these minima calculated via PM potentials are in somewhat better conformity with mass-spectrometric data than the values calculated via WK PF. OPLS gave much weaker water-base interactions for all compounds considered, thus these PF were not used in further computations. Monte Carlo simulations of the hydration of 9-methyladenine, 1-methyluracil and 1-methylthymine were performed in systems with 400 water molecules and periodic boundary conditions. Results of simulations with PM potentials give better agreement with experimental data on hydration energies than WK PF. Computations with PM PF of the hydration energy of keto and enol tautomers of 9-methylguanine can account for the shift in the tautomeric equilibrium of guanine in aqueous media to a dominance of the keto form in spite of nearly equal intrinsic stability of keto and enol tautomers. The results of guanine hydration computations are discussed in relation to mechanisms of base mispairing errors in nucleic acid biosynthesis. The data presented in this paper along with previous results on simulation of hydration shell structures in DNA duplex grooves provide ample evidence for the advantages of PM PF in studies of nucleic-acid hydration.

Modeling of a possible conformational change associated with the catalytic mechanism in the hammerhead ribozyme.

Here we describe a possible model of the cleavage mechanism in the hammerhead ribozyme. In this model, the 2' hydroxyl of C17 is moved into an appropriate orientation for an in-line attack on the G1.1 phosphate through a change in its sugar pucker from C3' endo to C2' endo. This conformational change in the active site is caused by a change in the uridine turn placing the N2 and N3 atoms of G5 of the conserved core in hydrogen bonding geometry with the N3 and N2 atoms on the conserved G16.2 residue. The observed conformational change in the uridine turn suggests an explanation for the conservation of G5. In the crystal structure of H.M. Pley et al., Nature 372, 68-74 (1994), G5 is situated 5.3A away from G16.2. However, the uridine turn is sufficiently flexible to allow this conformational change with relatively modest changes in the backbone torsion angles (average change of 14.2 degrees). Two magnesium ions were modeled into the active site with positions analogous to those described in the functionally similar Klenow fragment 3'-5' exonuclease (L.S. Beese and T.A. Steitz, EMBO J. 10, 25-33 (1991)), the Group I intron (T.A. Steitz and J.A. Steitz, P.N.A.S. U.S.A. 90, 6498-6502 (1993); R.F. Setlik et al., J. Biomol. Str. Dyn. 10, 945-972 (1993)) and other phosphotransferases. Comparison of this model with one in which the uridine turn conformation was not changed showed that although the changes in the C17 sugar pucker could be modeled, insufficient space existed for the magnesium ions in the active site.

The digoxigenin (DIG) system for non-radioactive labelling and detection of nucleic acids--an overview.

The use of non-radioactive nucleic acid probes has increased dramatically over the last years. The convenience of not having to deal with radioactive isotopes combined with the stability and economy of the non-radioactive systems has led to applications in various techniques. One of the most successful labelling and detection systems is based on the hapten digoxigenin. Here, the different methods for labelling nucleic acids with digoxigenin as well as the various possibilities for detection are described. Some typical applications illustrate the utility of the DIG system.

Finding the global minimum: a fuzzy end elimination implementation.

The 'fuzzy end elimination theorem' (FEE) is a mathematically proven theorem that identifies rotameric states in proteins which are incompatible with the global minimum energy conformation. While implementing the FEE we noticed two different aspects that directly affected the final results at convergence. First, the identification of a single dead-ending rotameric state can trigger a 'domino effect' that initiates the identification of additional rotameric states which become dead-ending. A recursive check for dead-ending rotameric states is therefore necessary every time a dead-ending rotameric state is identified. It is shown that, if the recursive check is omitted, it is possible to miss the identification of some dead-ending rotameric states causing a premature termination of the elimination process. Second, we examined the effects of removing dead-ending rotameric states from further considerations at different moments of time. Two different methods of rotameric state removal were examined for an order dependence. In one case, each rotamer found to be incompatible with the global minimum energy conformation was removed immediately following its identification. In the other, dead-ending rotamers were marked for deletion but retained during the search, so that they influenced the evaluation of other rotameric states. When the search was completed, all marked rotamers were removed simultaneously. In addition, to expand further the usefulness of the FEE, a novel method is presented that allows for further reduction in the remaining set of conformations at the FEE convergence. In this method, called a tree-based search, each dead-ending pair of rotamers which does not lead to the direct removal of either rotameric state is used to reduce significantly the number of remaining conformations. In the future this method can also be expanded to triplet and quadruplet sets of rotameric states. We tested our implementation of the FEE by exhaustively searching ten protein segments and found that the FEE identified the global minimum every time. For each segment, the global minimum was exhaustively searched in two different environments: (i) the segments were extracted from the protein and exhaustively searched in the absence of the surrounding residues; (ii) the segments were exhaustively searched in the presence of the remaining residues fixed at crystal structure conformations. We also evaluated the performance of the method for accurately predicting side chain conformations. We examined the influence of factors such as type and accuracy of backbone template used, and the restrictions imposed by the choice of potential function, parameterization and rotamer database. Conclusions are drawn on these results and future prospects are given.

Secondary structure in solution of two anti-HIV-1 hammerhead ribozymes as investigated by two-dimensional 1H 500 MHz NMR spectroscopy in water.

Two hammerhead chimeric RNA/DNA ribozymes (HRz) were synthesized in pure form. Both were 30 nucleotides long, and the sequences were such that they could be targeted to cleave the HIV-1 gag RNA. Named HRz-W and HRz-M, the former had its invariable core region conserved, the latter had a uridine in the invariable region replaced by a guanine. Their secodary structures were determined by 2D NOESY 1H 500 MHz NMR spectroscopy in 90% water and 10% D2(0), following the imino protons. The data show that both HRz-M and HRz-W form identical secondary structures with stem regions consisting of continuous stacks of AT and GT pairs. An energy minimized computer model of this stem region is provided. The results suggest that the loss of catalytic activity that is known to result when an invariant core residue is replaced is not related to the secondary structure of the ribozymes in the absence of substrate.

A full-coordinate model of the polymerase domain of HIV-1 reverse transcriptase and its interaction with a nucleic acid substrate.

We present a full-coordinate model of residues 1-319 of the polymerase domain of HIV-I reverse transcriptase. This model was constructed from the x-ray crystallographic structure of Jacobo-Molina et al. (Jacobo-Molina et al., P.N.A.S. USA 90, 6320-6324 (1993)) which is currently available to the degree of C-coordinates. The backbone and side-chain atoms were constructed using the MAXSPROUT suite of programs (L. Holm and C. Sander, J. Mol. Biol. 218, 183-194 (1991)) and refined through molecular modeling. A seven base pair A-form dsDNA was positioned in the nucleic acid binding cleft to represent the template-primer complex. The orientation of the template-primer complex in the nucleic acid binding cleft was guided by the positions of phosphorus atoms in the crystal structure.

Applicability of PM3 to transphosphorylation reaction path: toward designing a minimal ribozyme.

A growing body of evidence shows that RNA can catalyze many of the reactions necessary both for replication of genetic material and the possible transition into the modern protein-based world. However, contemporary ribozymes are too large to have self-assembled from a prebiotic oligonucleotide pool. Still, it is likely that the major features of the earliest ribozymes have been preserved as molecular fossils in the catalytic RNA of today. Therefore, the search for a minimal ribozyme has been aimed at finding the necessary structural features of a modern ribozyme (Beaudry and Joyce, 1990). Both a three-dimensional model and quantum chemical calculations are required to quantitatively determine the effects of structural features of the ribozyme on the reaction it catalyzes. Using this model, quantum chemical calculations must be performed to determine quantitatively the effects of structural features on catalysis. Previous studies of the reaction path have been conducted at the ab initio level, but these methods are limited to small models due to enormous computational requirements. Semiempirical methods have been applied to large systems in the past; however, the accuracy of these methods depends largely on the system under investigation. In the present study we assess the validity of the MNDO/PM3 method on a simple model of the ribozyme-catalyzed reaction, or hydrolysis of phosphoric acid. We find that the results are qualitatively similar to ab initio results using large basis sets. Therefore, PM3 is suitable for studying the reaction path of the ribozyme-catalyzed reaction.

[Genotoxic modification of nucleic acid bases and biological consequences of it. Review and prospects of experimental and computational investigations]. / Genotoksicheskaia modifikatsiia osnovanii nukleinovykh kislot i ee biologicheskie posledstviia. Obzor i perspektivy éksperimental'nykh i raschetno-teoreticheskikh issledovanii.

The review is presented of experimental and computational data on the influence of genotoxic modification of bases (deamination, alkylation, oxidation) on the structure and biological functioning of nucleic acids. Pathways are discussed for the influence of modification on coding properties of bases, on possible errors of nucleic acid biosynthesis, and on configurations of nucleotide mispairs. The atomic structure of nucleic acid fragments with modified bases and the role of base damages in mutagenesis and carcinogenesis are considered.

Modeling study on the cleavage step of the self-splicing reaction in group I introns.

A three-dimensional model of the Tetrahymena thermophila group I intron is used to further explore the catalytic mechanism of the transphosphorylation reaction of the cleavage step. Based on the coordinates of the catalytic core model proposed by Michel and Westhof (Michel, F., Westhof, E. J. Mol. Biol. 216, 585-610 (1990)), we first converted their ligation step model into a model of the cleavage step by the substitution of several bases and the removal of helix P9. Next, an attempt to place a trigonal bipyramidal transition state model in the active site revealed that this modified model for the cleavage step could not accommodate the transition state due to insufficient space. A lowering of P1 helix relative to surrounding helices provided the additional space required. Simultaneously, it provided a better starting geometry to model the molecular contacts proposed by Pyle et al. (Pyle, A. M., Murphy, F. L., Cech, T. R. Nature 358, 123-128. (1992)), based on mutational studies involving the J8/7 segment. Two hydrated Mg2+ complexes were placed in the active site of the ribozyme model, using the crystal structure of the functionally similar Klenow fragment (Beese, L.S., Steitz, T.A. EMBO J. 10, 25-33 (1991)) as a guide. The presence of two metal ions in the active site of the intron differs from previous models, which incorporate one metal ion in the catalytic site to fulfill the postulated roles of Mg2+ in catalysis. The reaction profile is simulated based on a trigonal bipyramidal transition state, and the role of the hydrated Mg2+ complexes in catalysis is further explored using molecular orbital calculations.

Cell-associated haemolytic activity of Helicobacter pylori.

Helicobacter pylori cells cultured on solid medium were quantitatively tested for haemolytic activity against erythrocytes of man, sheep, the guinea pig and rabbit. Using 4-day and 8-day cultures of two standard strains (ATCC 43504, IMMi 676), human erythrocytes were not lysed by 10% bacterial suspensions. Rabbit erythrocytes were the most sensitive to 8-day cultures. Hot-cold incubation yielded the highest haemolysis titres. The extent of haemolysis strongly correlated with the number of bacterial cells. Supplementation of the test medium (PBS, pH 7.4) with L-cysteine, dithiothreitol, MgCl2, EDTA, cholesterol, lecithin or sphingomyelin did not influence the haemolysis titres. They were significantly reduced in the presence of pronase E, human serum, bovine serum albumin or CaCl2, and by heat treatment of the bacteria. Supplementation of the test medium with cardiolipin strongly increased the haemolysis titres. Comparing the cell-associated haemolytic activity of 18 strains, the titres ranged from < 2 to 64, with a median titre of 16. No correlation was found between the haemolytic activity and phospholipase C activity of the cell suspensions. It was concluded that the formation of lysophosphatides and non-enzymatic factors rather than a sulphydryl-activated cytolysin or phospholipase C are responsible for the cell-associated haemolytic activity. This property may be involved in the pathogenicity and virulence of Helicobacter pylori.

Multipole correction of atomic monopole models of molecular charge distribution. I. Peptides.

The defects in atomic monopole models of molecular charge distribution have been analyzed for several model-blocked peptides and compared with accurate quantum chemical values. The results indicate that the angular characteristics of the molecular electrostatic potential around functional groups capable of forming hydrogen bonds can be considerably distorted within various models relying upon isotropic atomic charges only. It is shown that these defects can be corrected by augmenting the atomic point charge models by cumulative atomic multipole moments (CAMMs). Alternatively, sets of off-center atomic point charges could be automatically derived from respective multipoles, providing approximately equivalent corrections. For the first time, correlated atomic multipoles have been calculated for N-acetyl, N'-methylamide-blocked derivatives of glycine, alanine, cysteine, threonine, leucine, lysine, and serine using the MP2 method. The role of the correlation effects in the peptide molecular charge distribution are discussed.

Point charge representation of multicenter multipole moments in calculation of electrostatic properties.

Distributed Point Charge Models (PCM) for CO, (H2O)2, and HS-SH molecules have been computed from analytical expressions using multi-center multipole moments. The point charges (set of charges including both atomic and non-atomic positions) exactly reproduce both molecular and segmental multipole moments, thus constituting an accurate representation of the local anisotropy of electrostatic properties. In contrast to other known point charge models, PCM can be used to calculate not only intermolecular, but also intramolecular interactions. Comparison of these results with more accurate calculations demonstrated that PCM can correctly represent both weak and strong (intramolecular) interactions, thus indicating the merit of extending PCM to obtain improved potentials for molecular mechanics and molecular dynamics computational methods.